The approach should work, but I have not tried it. Are your meshes mixed element types? If so you likely also need to store quadrature rules as well or things will get very confused. If all elements are the same type you should be able to attach all the fe's to one qrule.
Have you tried just moving your fe object definition outside of the assemble loop so they are not created and destroyed at each iteration? This will not reduce the computational overhead but will drastically reduce memory allocation insode reinit. -Ben ----- Original Message ----- From: Manav Bhatia <[EMAIL PROTECTED]> To: [email protected] <[email protected]> Sent: Sat Oct 18 13:04:06 2008 Subject: [Libmesh-users] fe reinit Hi, For my application requiring the solution of a nonlinear transient system, I am doing a reinit of the fe object per elem per nonlinear iteration. For larger systems this has started to be a major CPU time expense. I am now considering saving one fe per elem in memory so that I do not have to do these reinits. Ofcourse, I will be committing a considerable amount of memory as well. I am writing to ask if anyone has tried this, and could share his/ her experiences or comment on this. Thanks, Manav ------------------------------------------------------------------------- This SF.Net email is sponsored by the Moblin Your Move Developer's challenge Build the coolest Linux based applications with Moblin SDK & win great prizes Grand prize is a trip for two to an Open Source event anywhere in the world http://moblin-contest.org/redirect.php?banner_id=100&url=/ _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users ------------------------------------------------------------------------- This SF.Net email is sponsored by the Moblin Your Move Developer's challenge Build the coolest Linux based applications with Moblin SDK & win great prizes Grand prize is a trip for two to an Open Source event anywhere in the world http://moblin-contest.org/redirect.php?banner_id=100&url=/ _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
