[Trilinos-Users] Scalar Epetra_Maps and Vector Epetra_Maps and linear solves

[Williams, Alan B]

John,
Which method to you use to obtain the domain-map from the vbr matrix?
If you use the method OperatorDomainMap(), you'll get the point-entry version of the map (with element sizes == 1). But if you use DomainMap(), you *should* get the correct block-map.
If not, I'll look into it deeper and see what's going wrong.
Alan


From: trilinos-users-bounces at software.sandia.gov [mailto:trilinos-users-bounces at software.sandia.gov] On Behalf Of Mitchell, John A
Sent: Tuesday, March 29, 2011 5:42 PM
To: 'trilinos-users at software.sandia.gov'
Subject: [Trilinos-Users] Scalar Epetra_Maps and Vector Epetra_Maps and linear solves

I need help / guidance on Epetra_BlockMaps, Epetra_FEVbrMatrix and linear solves.

I construct my matrix with an Epetra_BlockMap that has elements of length 3 - I am solving a problem that has 3 degrees of freedom per node.

I use 'BeginInsertGlobalValues' for each row and SubmitBlockEntry(k) where 'k' is a 3x3 block entry ( Epetra_SerialDenseMatrix)  for each column of the row;

I finish the row with EndSubmitEntries()

I finish all rows with FillComplete(rowMap,rowMap) where rowMap is the Epetra_BlockMap that has my owned rows and has elements of length 3.  I used this map to construct my matrix.  I was hoping that this version FillComplete would keep my rowMap as the domainMap and rangeMap but no such luck.

I construct my right hand side vector with 'rowMap'.  It's most logical to compute my rhs based upon the rowMap.  However, the domain and range of the matrix have Epetra_Maps that have element length=1.  These maps are incompatible with my rhs map.

How do I map data from my 'rowMap' to the 'domainMap' of the operator?

Thanks,
John

-------------- next part --------------
An HTML attachment was scrubbed...
URL: https://software.sandia.gov/pipermail/trilinos-users/attachments/20110330/665b9b42/attachment.html