I'm working some optimizations with how libMesh stores & allocates degree of freedom indices when there are multiple variables associated with a given System. These generically fall into two categories -
(1) homogeneous types, where there are multiple variables but they all have the same finite element type, and (2) mixed types, where there are different variables with different finite element types as well (like Taylor-Hood Q2/Q1 incompressible NS, for example) for both cases, we could do some optimization in how we store degree of freedom indices & do things like compute the sparse matrix graph. for case (1) in particular, we could also take advantage of the identical block-structure of the linear system as well and use e.g. PETSc's block matrix and vector support. So my question is, for your multivariable systems, are they typically of type (1) or (2)? How many variables and what types? My reason in asking is because some of the optimizations are easier for case (1), which is where I plan to start, but don't want to do anything that would preclude additional optimization for (2), especially if there are people with ~5 or more variables of different finite element types in a system. Thanks, -Ben ------------------------------------------------------------------------------ LogMeIn Rescue: Anywhere, Anytime Remote support for IT. Free Trial Remotely access PCs and mobile devices and provide instant support Improve your efficiency, and focus on delivering more value-add services Discover what IT Professionals Know. Rescue delivers http://p.sf.net/sfu/logmein_12329d2d _______________________________________________ Libmesh-devel mailing list Libmesh-devel@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-devel
