[Trilinos-Users] errors build MPI tpetra CrsMatrix

[David Hysom]

I'm getting an error building a Tpetra::CrsMatrix;
The attached is a slightly edited version of the code
used to build the matrix.

When I invoke our driver with a single processor,
everything works fine. When I use two or more
processors I get a traceback (using valgrind; see attached)
That shows an error in the call to A->fillComplete();
this subsequently shows a problem when makeColMap()
is invoked.

My debug output shows that the correct rows/columns are being
inserted on the appropriate processors.

Can you spot an error, or suggest ways of debugging
to find out what's going on?

thanks, David


-------------- next part --------------
Conditional jump or move depends on uninitialised value(s)           std::_Rb_tree<int, int, std::_Identity<int>, std::less<int>, std::allocator<int> >::_M_insert_unique(int const&) (stl_tree.h:1167)
  at 0x41FF7E: std::_Rb_tree<int, int, std::_Identity<int>, std::less<int>, std::allocator<int> >::_M_insert_unique(int const&) (stl_tree.h:1167)
  by 0x4204B5: std::set<int, std::less<int>, std::allocator<int> >::insert(int const&) (stl_set.h:411)
  by 0x58F0E5: Tpetra::CrsGraph<int, int, KokkosClassic::SerialNode, KokkosClassic::AltSparseOps<void, int, KokkosClassic::SerialNode, KokkosClassic::details::AltSparseOpsDefaultAllocator<int, KokkosClassic::SerialNode> > >::makeColMap() (Tpetra_CrsGraph_def.hpp:3253)
  by 0x54792B: Tpetra::CrsMatrix<double, int, int, KokkosClassic::SerialNode, KokkosClassic::AltSparseOps<void, int, KokkosClassic::SerialNode, KokkosClassic::details::AltSparseOpsDefaultAllocator<int, KokkosClassic::SerialNode> > >::fillComplete(Teuchos::RCP<Tpetra::Map<int, int, KokkosClassic::SerialNode> const> const&, Teuchos::RCP<Tpetra::Map<int, int, KokkosClassic::SerialNode> const> const&, Teuchos::RCP<Teuchos::ParameterList> const&) (Tpetra_CrsMatrix_def.hpp:3061)
  by 0x548F4B: Tpetra::CrsMatrix<double, int, int, KokkosClassic::SerialNode, KokkosClassic::AltSparseOps<void, int, KokkosClassic::SerialNode, KokkosClassic::details::AltSparseOpsDefaultAllocator<int, KokkosClassic::SerialNode> > >::fillComplete(Teuchos::RCP<Teuchos::ParameterList> const&) (Tpetra_CrsMatrix_def.hpp:2937)
  by 0x460444: MatrixUtil::buildLaplacian_stupid(Teuchos::RCP<std::map<int, std::set<int, std::less<int>, std::allocator<int> >, std::less<int>, std::allocator<std::pair<int const, std::set<int, std::less<int>, std::allocator<int> > > > > >&, int, int) (MatrixUtil.cpp:66)
  by 0x474132: main (example_2.cpp:41)
Uninitialised value was created by a stack allocation
  at 0x58E6F8: Tpetra::CrsGraph<int, int, KokkosClassic::SerialNode, KokkosClassic::AltSparseOps<void, int, KokkosClassic::SerialNode, KokkosClassic::details::AltSparseOpsDefaultAllocator<int, KokkosClassic::SerialNode> > >::makeColMap() (Tpetra_CrsGraph_def.hpp:3173)
U
-------------- next part --------------
A non-text attachment was scrubbed...
Name: buildMat.cpp
Type: text/x-c++src
Size: 1615 bytes
Desc: not available
URL: <http://software.sandia.gov/pipermail/trilinos-users/attachments/20140903/23168bab/attachment.bin>