79 int_type NumGlobalElements = (int_type) OldDomainMap.NumGlobalElements64();
80 assert( orig.RowMap().NumMyElements() ==
NewRowMap_.NumMyElements() );
82 if (NumGlobalElements == 0 && orig.RowMap().NumGlobalElements64() == 0 )
95 Epetra_Map tmpColMap( NumGlobalElements, NumMyElements, 0, OldDomainMap.Comm() );
97 for(
int i = 0; i < NumMyElements; ++i )
98 Cols[i] = (int_type) tmpColMap.GID64(i);
100 NewCols.Import( Cols, Importer,
Insert );
102 std::vector<int_type*> NewColIndices(1);
103 NewCols.ExtractView( &NewColIndices[0] );
105 int NumMyColElements = OldColMap.NumMyElements();
106 int_type NumGlobalColElements = (int_type) OldColMap.NumGlobalElements64();
108 NewColMap_ =
new Epetra_Map( NumGlobalColElements, NumMyColElements, NewColIndices[0], (int_type) OldColMap.IndexBase64(), OldColMap.Comm() );
117 int numMyRows = NewMatrix->NumMyRows();
118 for(
int i = 0; i < numMyRows; ++i )
120 orig.ExtractMyRowView( i, indicesCnt, myValues, myIndices );
121 NewMatrix->InsertMyValues( i, indicesCnt, myValues, myIndices );
124 NewMatrix->FillComplete();