// is initialized with both row and
// column map. Maybe find something
// more out about this...
- //reinit (input_map, input_map, sparsity_pattern);
+ reinit (input_map, input_map, sparsity_pattern);
- matrix.reset();
+ /* matrix.reset();
unsigned int n_rows = sparsity_pattern.n_rows();
(new Epetra_FECrsMatrix(Copy, *row_map, &n_entries_per_row[0],
false));
- reinit (sparsity_pattern);
+ reinit (sparsity_pattern);*/
}
{
// flush buffers
int ierr;
- if (row_map->SameAs(*col_map))
- ierr = matrix->GlobalAssemble ();
- else
+ //if (row_map->SameAs(*col_map))
+ //ierr = matrix->GlobalAssemble ();
+ //else
ierr = matrix->GlobalAssemble (*col_map, *row_map);
AssertThrow (ierr == 0, ExcTrilinosError(ierr));
AssertThrow ((static_cast<signed int>(index) >= vector.map->MinMyGID()) &&
(static_cast<signed int>(index) <= vector.map->MaxMyGID()),
ExcAccessToNonLocalElement (index, vector.map->MinMyGID(),
- vector.map->MaxMyGID()-1));
+ vector.map->MaxMyGID()));
return (*(vector.vector))[0][index];
}
{
// Extract local indices in
// the vector.
- int trilinos_i = map->LID(index);
+ int trilinos_i = vector->Map().LID(index);
TrilinosScalar value = 0.;
if (trilinos_i == -1 )
{