sparse_matrix_ptrtype
newZeroMatrix( datamap_ptrtype const& d1, datamap_ptrtype const& d2 )
{
auto A = sparse_matrix_ptrtype( new eigen_sparse_matrix_type( d1->nGlobalElements(), d2->nGlobalElements() ) );
//A->setMatrixProperties( matrix_properties );
return A;
}
sparse_matrix_ptrtype
newMatrix( datamap_ptrtype const& d1, datamap_ptrtype const& d2, size_type matrix_properties = NON_HERMITIAN, bool init = true )
{
auto A = sparse_matrix_ptrtype( new eigen_sparse_matrix_type( d1->nGlobalElements(), d2->nGlobalElements() ) );
A->setMatrixProperties( matrix_properties );
return A;
}
vector_ptrtype
newVector( datamap_ptrtype const& d )
{
return vector_ptrtype( new eigen_vector_type( d->nGlobalElements() ) );
}
