Pytrilinos Epetra map: Cannot create overlapping maps
Created by: hakeemo
I cannot seem to create a nonunique Epetra map in PyTrilinos. When I run the following code
from PyTrilinos import Epetra
import numpy as np
len_global = 100
comm = Epetra.PyComm()
unique_map = Epetra.Map(len_global, 0 , comm)
print unique_map.MyGlobalElements()
overlap_length = 10
num_proc = comm.NumProc()
my_id = comm.MyPID()
my_range = np.arange(len_global//num_proc*my_id,
min(len_global//num_proc*(my_id+1) + overlap_length, len_global,len_global), dtype = np.int32)
my_range= my_range.tolist()
nonunique_map = Epetra.Map(len_global, my_range , 0, comm)
print nonunique_map.MyGlobalElements()
I get
File "/home/tigon/trilinos-install/lib/python2.7/site-packages/PyTrilinos/Epetra.py", line 3222, in init
this = _Epetra.new_Map(*args)
Epetra.Error: Error code = -4
Unless I set overlap_length=0, which works fine.