[Trilinos-Users] Need Advice On A Special Kind of FE Matrix Assebmly

Davood Ansari david.ansari at gmail.com
Wed May 28 01:41:35 MDT 2008


I am assembling a FEM matrix. And of course I am using Epetra's
FE_CrsMatrix class.

Unlike ordinary FE matrix assembly which is an iteration over the entire
elements, I have my elements sorted in a binary tree. First I have to
the deepest level elements (branch ends) into the global matrix. Then, I
have to
assemble the elements that are located one level higher (higher means
towards the tree root).
At each level, of the assembly I need to have access to the data that has
been previously
(only form the previous level) assembled into the global matrix.

If I were going to read the required data directly from the global matrix,
then I'd have to call
FillComplete and GlobalAssemble methods repeatedly. I am not sure if such an
approach would even
work (probably it will not without sophisticated synchronizations between
different processes in parallel)
and predict that it will turn out to be very inefficient.
I understand that this can be done using a recursive process but I am
avoiding it simply because it is
not a good parallel solution.

I wonder if you guys have a better suggestion in mind.
Please Comment

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://software.sandia.gov/mailman/private/trilinos-users/attachments/20080528/a747a3eb/attachment.html 

More information about the Trilinos-Users mailing list