[Trilinos-Users] Having Different Row maps and A Transposition : What To Do ?
david.ansari at gmail.com
Tue Jun 17 01:33:12 MDT 2008
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
(like A and B) with the vectors that are constructed based on the
Now, which one is possibly the right way to go:
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)
1.Assemble *C*, using the *usual* map, thus requiring off process
communications (definitely Epetra_FECrs_Matrix)
2.Then transpose it.
Kindly comment on this
-------------- next part --------------
An HTML attachment was scrubbed...
More information about the Trilinos-Users