BTW, without telling me so much that you have to kill me...

What are you solving with 9 variables?  More importantly, are they all the same 
FE type?

When I first started putting the DofMap together Taylor-Hood elements were 
typical, with varying orders for dfferent variables.  But in *all* the 
applications since equal order for all variables has worked fine.

The reason  I ask is 'cus if you dig around in the DofMap, especially the 
distribution part, there is a lot of code which could be simplified if you knew 
a priori all variable types were identical.  This is especially the case when 
we build the graph of the sparse matrix - it should instead be done for one 
prototypical variable and the pattern for the rest of the matrix can be 
inferred.  That is an n_vars^2 savings on something that happens at each 
adaptive step.

Similarly, constraints are calculated n_vars times too many...

I've been thinking about a derived specialization of the DofMap for this 
special (but common) case.

-Ben
-------------------------------------------------------------------------
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

Reply via email to