cargo run --release --features intel_mkl --bin brusselator_pde -- --npoint=151 --blas-nt=24 -g klu -o metis Second-book problem = true Number of points along x and y = 151 Tolerance (abs_tol = rel_tol) = 1.00E-04 Concurrent real and complex sys = true Problem dimension (ndim) = 45602 Number of non-zeros (jac_nnz) = 319214 Number of BLAS threads = 24 Linear solver = Klu Radau5: Radau method (Radau IIA) (implicit, order 5, embedded) Number of function evaluations = 247 Number of Jacobian evaluations = 22 Number of factorizations = 40 Number of lin sys solutions = 69 Number of performed steps = 40 Number of accepted steps = 33 Number of rejected steps = 7 Number of iterations (maximum) = 3 Number of iterations (last step) = 1 Last accepted/suggested stepsize = 0.03483119190513051 Max time spent on a step = 2.897575677s Max time spent on the Jacobian = 1.600501ms Max time spent on factorization = 3.298239682s Max time spent on lin solution = 25.142284ms Total time = 1m55.943576188s cargo run --release --features intel_mkl --bin brusselator_pde -- --npoint=151 --blas-nt=24 -g umfpack -o metis Second-book problem = true Number of points along x and y = 151 Tolerance (abs_tol = rel_tol) = 1.00E-04 Concurrent real and complex sys = true Problem dimension (ndim) = 45602 Number of non-zeros (jac_nnz) = 319214 Number of BLAS threads = 24 Linear solver = Umfpack Radau5: Radau method (Radau IIA) (implicit, order 5, embedded) Number of function evaluations = 247 Number of Jacobian evaluations = 22 Number of factorizations = 40 Number of lin sys solutions = 69 Number of performed steps = 40 Number of accepted steps = 33 Number of rejected steps = 7 Number of iterations (maximum) = 3 Number of iterations (last step) = 1 Last accepted/suggested stepsize = 0.034831191910329906 Max time spent on a step = 231.65154ms Max time spent on the Jacobian = 1.614761ms Max time spent on factorization = 334.260584ms Max time spent on lin solution = 36.129836ms Total time = 10.957887941s cargo run --release --features intel_mkl --bin brusselator_pde -- --npoint=151 --blas-nt=24 -g mumps -o metis Second-book problem = true Number of points along x and y = 151 Tolerance (abs_tol = rel_tol) = 1.00E-04 Concurrent real and complex sys = true Problem dimension (ndim) = 45602 Number of non-zeros (jac_nnz) = 319214 Number of BLAS threads = 24 Linear solver = Mumps Radau5: Radau method (Radau IIA) (implicit, order 5, embedded) Number of function evaluations = 247 Number of Jacobian evaluations = 22 Number of factorizations = 40 Number of lin sys solutions = 69 Number of performed steps = 40 Number of accepted steps = 33 Number of rejected steps = 7 Number of iterations (maximum) = 3 Number of iterations (last step) = 1 Last accepted/suggested stepsize = 0.0348311919052513 Max time spent on a step = 181.65336ms Max time spent on the Jacobian = 1.726563ms Max time spent on factorization = 406.78126ms Max time spent on lin solution = 16.975554ms Total time = 8.175700517s