Blowup with OpenMP
Dear CROCO users,
I'm running the test case of RIP. Using the default cpp setup (undef OPENMP and undef MPI), the codes runs well. But when I define OPENMP, the codes blow up after a few steps with the identical croco.in.Rip file. For example, set NPP=2, NSUB_X=1, NSUB_E=NPP, and export OMP_NUM_THREADS=2. Then I checked the output, the solutions have extremes around the boundary between two tiles. So I believe the blowup is caused by some parallel bugs. Could someone please take a look?
Thank you, Peng