Best performance for initialising a variable number of matrices/vectors within a Crank-Nicolson scheme

5 posts were split to a new topic: Why re-use factorization(A) rather than inv(A)