[Trilinos-Users] bug in EpetraExt_HDF5?

Simone Deparis simone.deparis at epfl.ch
Thu Aug 7 08:59:49 MDT 2008


Hi,

(trilinos 8.0.7)
I suspect a bug in EpetraExt_HDF5.cpp, line 902:
    Epetra_Map LinearMap(X.GlobalLength(), 0, Comm_);

if the roiginal multivector has indexBase == 1, the temporary
multivector LinearX has one missing entry.

Can someone confirm, please? If I put
    Epetra_Map LinearMap(X.GlobalLength(), X.Map().IndexBase(), Comm_);

should work?

Thank you
best
Simone



More information about the Trilinos-Users mailing list