ConceptPythonIntermediate Finite differences: the heat equation as a matrix, and why the time step has a limit solve = factorized(sparse.eye_array(m, format="csc") - dt * D * A_sp)
ConceptPythonAdvanced The adjoint method: why 200 derivatives cost two solves, not 201 solve = factorized(K_of(k)) # one LU factorization for both solves