[Trilinos-Users] Big Big Error on generating Tridiag matrix, am I wrong?

john ax.equal.b at gmail.com
Thu Oct 29 21:02:01 MDT 2009

Dear all,

I updated to Trililnos 10.0.2 yesterday.
But I found that the generated Tridiag matrix in Both triutils and galeri
packages are wrong.

In source Trilinos_Util_CrsMatrixGallery.cpp reads following lines
    } else {
      // off-diagonal for internal row
      Indices[0] = MyGlobalElements_[i]-1;
     * Values[1] = b_;// Should the index be **Values[0] = b_*
      Indices[1] = MyGlobalElements_[i]+1;
     * Values[0] = c_;// Should the index be **Values[1] = c_*
      NumEntries = 2;
Here b_ is the subdiagonal and c_ is the superdiagonal.

The same code appears in Galeri_Tridiag.h.

Does my observation is wrong?

-------------- next part --------------
An HTML attachment was scrubbed...
URL: https://software.sandia.gov/pipermail/trilinos-users/attachments/20091030/2d8dc802/attachment.html 

More information about the Trilinos-Users mailing list