Christophe Ortiz <[email protected]> writes: > For the moment, dof is small (dof=4). Still doing some tests with the > classes and methods. But should reach 1000-10000 in production.
This is a huge difference. It's a waste of time to profile cases that don't matter, so try to profile the cases that matter (even if the "physics" is mocked). Note that for 4x4, you can use PetscScalar emat[4][4]; > Ahhh...did not think about it ! This would allow to create the vector only > once in the main (after dof is determined) and pass it as argument to > FormIJacobian(). I will try. Thanks ! Yes, that's what the context is for.
pgpgTdXibLSlM.pgp
Description: PGP signature
