On Wed, 10 Sep 2008, Kirk, Benjamin (JSC-EG) wrote:
In the projection part we loop over var/elem under the assumption that
the var types are different. That way the fe objects are not created
over and over...
Yes, but we could just keep a small cache of FE objects, like
FEMSystem does. They're not usually too big, they're often sharable
between variables, and there's not too many of them.
In my case (and I suspect most frequently) the vars are all the same, so
it really should just be a loop over the elements, get the interpolation
weights, and then apply them to all the vars. In my case that world
speed this up by 5x.
Exactly.
Also, if each system projects all of its N vectors at once that could be
a big win - potentially another factor of N. The only drawback in this
case is all the old and new vectors must exist concurrently instead of
sequentially, so it costs more in RAM.
That cost could bother people; more CPU time means more CPU time, but
more memory sometimes means you just can't solve the problem.
Projecting vectors simultaneously (within a System or within a whole
EquationSystems) should be a non-default option at best, I think.
---
Roy
-------------------------------------------------------------------------
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-devel mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/libmesh-devel
