On Sat, Oct 18, 2008 at 1:04 PM, Manav Bhatia <[EMAIL PROTECTED]> wrote: > 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.
I think it depends which part of reinit is taking up the most time. If it is the derivative calculations then I'm not sure how much info you can really cache. I'm assuming here you have elements with non-affine maps... If you have a hybrid mesh (a mixture of geometric element types) you will probably gain some performance if you loop over all elements of a given type instead of switching back and forth repeatedly between geometric element types. -- John ------------------------------------------------------------------------- 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
