Tiny bug Report : GMResMPI or BicgstabMPI running in subMPI domain failed with unsuccessful ' norm2() '.

[Zhouccc]

Dear sir,

STRUMPACK/src/dense/DistributedVector.hpp

Line 41 in e0fc9d5

return std::sqrt(MPIComm().all_reduce(vdotv*vdotv, MPI_SUM));

 return std::sqrt(MPIComm().all_reduce(vdotv*vdotv, MPI_SUM));

I have decomposited my MPI_COMMON_WORLD(size 16) into 4 domains with 4 MPIs every region.
Then here, run ' spss_ -> solve() ' in 4 MPI domians
and found others will be barried when the 1st of 4 solve() finish.

Result : Used by GMResMPI or BicgstabMPI, the norm2() contains ' MPIComm() ' , which default to MPI_COMMON_WORLD , rather than given ' comm' .

FIX code :

return std::sqrt(comm.all_reduce(vdotv*vdotv, MPI_SUM));

After a few test by solving FEM matrix , it works well :)