FD uses sparse matrices