Tpetra::CrsMatrix: Possible bug in transposed apply
Created by: mhoemmen
@trilinos/tpetra
@ikalash reported the following:
I am trying to use Tpetra::CrsMatrix apply method with a Teuchos::TRANS combine mode, and the method does not appear to be working correctly. I end up with a vector of zeros even though the operator is nonzero, nor is the input vector. If I set the combine mode to Teuchos::NO_TRANS, things work correctly. I assume Teuchos::TRANS has been tested, so perhaps I am doing something wrong, although I am not sure what it could be. Is there some caveat about the method’s usage with the TRANS combine mode? I printed the Boolean returned when calling hasTransposeApply() and it prints true.
@ikalash also sent data, which I'll post here.