[Trilinos-Users] Having Different Row maps and A Transposition : What To Do ?
Davood Ansari
david.ansari at gmail.com
Tue Jun 17 01:33:12 MDT 2008
Hi all
Besides assembling my *A* and *B* matrices (as in *A *x* = *lambda* B *y
for an FE eigen-problem)
a third (constraints) matrix is needed to be assembled (call this one *C*).
For A and B, I have opted for the *usual* linearly sliced row map.
For C, the assembling process is much different (kind of recursive) from
ordinary FE matrices. Yet, there is a natural map that can partition *C*'s
rows into independent
slices. In the Solver stage, it is needed to multiply the transpose of
*C*many times
(like A and B) with the vectors that are constructed based on the
*usual*row map.
Now, which one is possibly the right way to go:
Option I.
1.Assemble *C*, using its natural map, thus requiring no off process
communications (requires Epetra_Crs_Matrix).
2.Then, transform *C* into a new matrix which is the transpose of *C* and
complies to the *usual* map(don't even know how)
Option II.
1.Assemble *C*, using the *usual* map, thus requiring off process
communications (definitely Epetra_FECrs_Matrix)
2.Then transpose it.
Kindly comment on this
Davood
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://software.sandia.gov/mailman/private/trilinos-users/attachments/20080617/2e7afa79/attachment.html
More information about the Trilinos-Users
mailing list