Re: [Libmesh-users] fe reinit

2008-10-18 Thread Roy Stogner
On Sat, 18 Oct 2008, Derek Gaston wrote: > What I'm getting at here is that there are lots of ways to spend > computational time... and it is _very_ application dependent. This is definitely true. > I'm personally not worried about how long FE::reinit() is taking > because I know for a fact tha

Re: [Libmesh-users] fe reinit

2008-10-18 Thread Manav Bhatia
hmmm. actually, that might be it. My FEBase::build is inside the assembly routine. I will move this out of the routine and update you on how things change. Thanks, Manav On Oct 18, 2008, at 5:24 PM, Benjamin Kirk wrote: >> I reinit a fe object with only the same elem kind. Hence, that

Re: [Libmesh-users] fe reinit

2008-10-18 Thread Benjamin Kirk
> I reinit a fe object with only the same elem kind. Hence, that aspect > of reallocating memory space due to changing elem types does not seem > to be a problem. What I was specifically asking is where code like AutoPtr fe (FEBase::build(dim, fe_type)); Sits. The issue is if it is inside your

Re: [Libmesh-users] fe reinit

2008-10-18 Thread Derek Gaston
On Oct 18, 2008, at 2:13 PM, Manav Bhatia wrote: > I am curious about this: doesn't any nonlinear transient computation > require this information per element per iteration? Am I the first one > to consider caching this information? How do the CPU/memory overheads > for your problems work out? Th

Re: [Libmesh-users] fe reinit

2008-10-18 Thread Manav Bhatia
I reinit a fe object with only the same elem kind. Hence, that aspect of reallocating memory space due to changing elem types does not seem to be a problem. Why would non-affine maps affect the memory footprint regarding caching of information? Since my mesh geometry does not change during

Re: [Libmesh-users] fe reinit

2008-10-18 Thread John Peterson
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 > ex

Re: [Libmesh-users] fe reinit

2008-10-18 Thread Manav Bhatia
- Original Message - > From: Manav Bhatia <[EMAIL PROTECTED]> > To: libmesh-users@lists.sourceforge.net > > Sent: Sat Oct 18 13:04:06 2008 > Subject: [Libmesh-users] fe reinit > > Hi, > > For my application requiring the solution of a nonlinear transien

Re: [Libmesh-users] fe reinit

2008-10-18 Thread Kirk, Benjamin (JSC-EG)
t;[EMAIL PROTECTED]> To: libmesh-users@lists.sourceforge.net 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 lar

[Libmesh-users] fe reinit

2008-10-18 Thread Manav Bhatia
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 d