[Trilinos-Users] bug in EpetraExt_HDF5?

Heroux, Michael A maherou at sandia.gov
Fri Aug 8 07:31:56 MDT 2008


These changes will be part of the Trilinos 9.0 release in September.  If you need them sooner, just let me know.

Mike


On 8/8/08 3:08 AM, "Simone Deparis" <simone.deparis at epfl.ch> wrote:




Simone Deparis wrote:
> 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

Thanks to Mike Heroux, the problem is fixed:
changes in EpetraExt_HDF5.cpp

902c902
<    Epetra_Map LinearMap(X.GlobalLength(), X.Map().IndexBase(), Comm_);
---
>    Epetra_Map LinearMap(X.GlobalLength(), 0, Comm_);

930c930
<   hsize_t offset[] = {0, LinearX->Map().GID(0)-X.Map().IndexBase()};
---
>   hsize_t offset[] = {0, LinearX->Map().GID(0)};


As far as I understood, this checnges are already in the development
branch. When will they be part of the distribution?

Thank you again
best
Simone

_______________________________________________
Trilinos-Users mailing list
Trilinos-Users at software.sandia.gov
http://software.sandia.gov/mailman/listinfo/trilinos-users


-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://software.sandia.gov/mailman/private/trilinos-users/attachments/20080808/fbf9ee9e/attachment.html 


More information about the Trilinos-Users mailing list