Changeset 23348

Timestamp:

09/28/18 07:59:34 (6 years ago)

Author:

Mathieu Morlighem

Message:

BUG: fixed solver in parallel with CoDiPack

File:

• issm/trunk-jpl/src/c/toolkits/mumps/MumpsSolve.cpp (modified) (1 diff)

issm/trunk-jpl/src/c/toolkits/mumps/MumpsSolve.cpp

@@ -525 +525 @@
 
   updateVectorAdjoint(tape, indexB, adjX, n);
+
+  // bcast dp_y (the solution of the forward system)
+  ISSM_MPI_Bcast(valueX,n,ISSM_MPI_PDOUBLE,0,IssmComm::GetComm());
+  // bcast the adjoint of the right-hand-side, i.e. this solution
+  ISSM_MPI_Bcast(adjX,n,ISSM_MPI_PDOUBLE,0,IssmComm::GetComm());
+
   for(int i=0; i<local_nnz; ++i) {
     // we access the transposed matrix here because we stored the indices in a transposed way