Okay, so dxyz seems to be appropriate for storing the metric - that problem is solved then. Now, one last thing: All DOFs are at the vertices (and we have only one DOF per vertex). Ie the global vector and matrices have dimension M or MxM if M is the number of mesh vertices. However, each DOF is not only shared with all the triangles incident to that vertex, but also with all elements which have this vertex in their patch (ie in their 1-neighborhood). E.g., for the element e shown below, the element matrix contribution Ke(i,j) is a (N+6)x(N+6) matrix, where only the local indices 0,1 and N directly belong to element e, while the other DOFs belong to the elements around e. I don't think that dof_map can provide such a map out of the box, or am I wrong?
I was therefore thinking of the following: Create a custom subdiv_dof_map, which collects the dof indices of all vertices in the patch of e and then (in the correct order) combines them into one dof_indices vector. This one can then be used as before. I guess it would be minimally invasive and only need small changes in the application code. I just hope there's no hidden catch... -Norbert >> >> N+4 - N+1 - N+2 >> / \ / \ / \ >> / \ / \ / \ >> N+5 -- N --- 1 -- N+3 >> \ / \ e / \ / >> \ / \ / \ / >> N-1--- 0 --- 2 >> \ /|\ / >> \ / | \ / >> 5 -4- 3 >> ------------------------------------------------------------------------- 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
