// release memory before reallocation
matrix.reset ();
- temp_vector.clear ();
matrix.reset (new Epetra_FECrsMatrix(*m.matrix));
}
const bool exchange_data)
{
// release memory before reallocation
- temp_vector.clear();
matrix.reset();
// if we want to exchange data, build a usual Trilinos sparsity pattern
void
SparseMatrix::reinit (const SparsityPattern &sparsity_pattern)
{
- temp_vector.clear ();
matrix.reset ();
// reinit with a (parallel) Trilinos
SparseMatrix::reinit (const SparseMatrix &sparse_matrix)
{
column_space_map.reset (new Epetra_Map (sparse_matrix.domain_partitioner()));
- temp_vector.clear ();
matrix.reset ();
matrix.reset (new Epetra_FECrsMatrix
(Copy, sparse_matrix.trilinos_sparsity_pattern(), false));
const Epetra_CrsGraph *graph = &input_matrix.Graph();
- temp_vector.clear ();
matrix.reset ();
matrix.reset (new Epetra_FECrsMatrix(Copy, *graph, false));
// empty matrix.
column_space_map.reset (new Epetra_Map (0, 0,
Utilities::Trilinos::comm_self()));
- temp_vector.clear();
matrix.reset (new Epetra_FECrsMatrix(View, *column_space_map, 0));
matrix->FillComplete();