Epetra_CrsMatrix::Multiply with locally replicated range map question

Created by: bartgol

@trilinos/epetra

Current Behavior

Currently, if the range map is locally replicated, the result of A.Apply(X,Y) is not what I expected. In particular, Epetra does a final call Y.Reduce(), which ends up giving a Y vector that is numProc times larger than what I expected. Why is Epetra doing this? I.e., why are we reducing the output vector, in case of a locally replicated range map? Isn't the result (before the Reduce() call) already what one most likely wants? The user can always call Y.Reduce() if he is interested in that.

Possible Solution

Assuming my concern is correct, simply remove line 3212 (and the likes in other branches of if statements) in Epetra_CrsMatrix.cpp