Hi Ben, Yes, my mesh has mixed elements, and some of the elements need two different kinds of fe.
I do not create new fe objects, just reinitialize them with new elements of the same kind. Thanks, Manav On Oct 18, 2008, at 3:11 PM, Kirk, Benjamin (JSC-EG) wrote: > 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
