PyTrilinos: Initializing an Epetra.CrsMatrix using Epetra.Vector.Map() does not work.
Created by: hakeemo
If I run
from PyTrilinos import Epetra
comm = Epetra.PyComm()
my_id = comm.MyPID()
num_proc = comm.NumProc()
num_cells = 1000
row_map = Epetra.Map(num_cells, 0 , comm)
V = Epetra.Vector(row_map)
A = Epetra.CrsMatrix(Epetra.Copy, V.Map(), 1) #This line raises an error
I get a NotImplementedError. This is strange since V.Map() and row_map are both BlockMaps, but using the later to initialize a CrsMatrix works fine.