[deal.II] Re: Trilinos SparseMatrix mmult bug

2016-09-29 Thread Ben Shields
So this turned out to be a bug in my own code after I made the original change in trilinos_sparse_matrix.cc. When I was constructing the sparsity patterns for the matrices using DoFTools::make_sparsity_pattern, there is a flag I had set keep_constrained_dofs = false. Setting this flag to true ha

[deal.II] Re: Trilinos SparseMatrix mmult bug

2016-09-23 Thread Ben Shields
Hello Bruno, Sorry for the vagueness, I was not confident if this was an actual bug in deal.ii (I've never encountered one before!). Here I've attached a bare-bones version of my code which reproduces the bug I'm encountering. I basically look at 4 cases: 1) deal.ii's TrilinosWrappers::SparseM

[deal.II] Re: Trilinos SparseMatrix mmult bug

2016-09-23 Thread Bruno Turcksin
Ben, On Thursday, September 22, 2016 at 5:41:26 PM UTC-4, Ben Shields wrote: > > Hey there, thanks for the response. I made the recommended change to > trilinos_sparse_matrix.cc, and the corresponding test code returns the > appropriate output. However, I'm still getting an issue in my particula

[deal.II] Re: Trilinos SparseMatrix mmult bug

2016-09-22 Thread Ben Shields
I've just updated to 8.4.2 and repeated the change to trilinos_sparse_matrix.cc, and the matrices appear the same as my previous post. So I've at least ruled out that it's not an issue specific to my version number. -- The deal.II project is located at http://www.dealii.org/ For mailing list/f

[deal.II] Re: Trilinos SparseMatrix mmult bug

2016-09-22 Thread Ben Shields
Hey there, thanks for the response. I made the recommended change to trilinos_sparse_matrix.cc, and the corresponding test code returns the appropriate output. However, I'm still getting an issue in my particular case. It's improved in that the resulting matrix is now square and the dimensions

[deal.II] Re: Trilinos SparseMatrix mmult bug

2016-09-21 Thread Daniel Jodlbauer
Sounds like we ran into the same problem https://github.com/dealii/dealii/pull/2536, although is also experienced this issue in serial. -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You receiv