[Trilinos-Users] Replacing maps in Epetra CRS Matrix

[Pavel Jiránek]

Hello,

I have a distributed CRS matrix (call it A) with a given row map. For some
reason, I would need to make a permutation of the global IDs so that the
row map is linear, if possible without a need to create another matrix and
fill it again using different (linear) row map.

What I do at the moment is:
1) Create a linear row map of the same local dimension as the row map of A.
2) Make an int vector view to its GIDs (on the row map of A).
3) Make an int vector on the col map of A and import the vector from the
previous step.
4) Create a col map deduced from the imported GIDs.

Then I call ReplaceRowMap(), ReplaceColMap(), and
ReplaceDomainMapAndImporter() on A. I suppose that this is correct (no
error happens) but still, the matrix still has an incorrect range map
equivalent to the old (nonlinear) row map. Is there any way to replace also
that one?

Thanks in advance for any feedback.

Regards,

Pavel
-------------- next part --------------
An HTML attachment was scrubbed...
URL: https://software.sandia.gov/pipermail/trilinos-users/attachments/20131001/d28006f5/attachment.html