invalid read with fortran example, dense package 1.1.1 about strumpack HOT 3 OPEN

I can reproduce the problem when running with 4 mpi processes, but not with 1.
I will look into it.

We will try to add a Fortran interface to the newest version. We no longer further develop the 1.1.1 version of the code. However, since it might take a while before we have the Fortran interface to the newest code ready, I recommend you keep using the 1.1.1 version for now. So I will check this invalid read issue as soon as possible.

Pieter

from strumpack.

The problem is on line 298 in HSS_par_tools.hpp:
pgemr2d(m,n,A,child1->tree_ib,child2->tree_ib,descA,A12[my_i],IONE,IONE,desc,globctxt);
replace that with:
pgemr2d(m,n,A,child1->tree_ib,child2->tree_ib,descA,(desc[BLACSctxt] != -1) ? A12[my_i] : NULL,IONE,IONE,desc,globctxt);

Same for line 314 in HSS_par_tools.hpp:
pgemr2d(n,m,A,child2->tree_ib,child1->tree_ib,descA,A21[my_i],IONE,IONE,desc,globctxt);
replace that with:
pgemr2d(n,m,A,child2->tree_ib,child1->tree_ib,descA,(desc[BLACSctxt] != -1) ? A21[my_i] : NULL,IONE,IONE,desc,globctxt);

from strumpack.

Thanks Pieter, this fixed the problem indeed. I tried to debug it myself, but I am not a C man.

I will use version 1.1.1 for now, but am indeed interested to switching to 2 when a Fortran wrapper becomes available.

Thanks for the effort to make this library!

from strumpack.