Ask Your Question
0

Blowup with OpenMP

asked 2019-02-27 23:34:21 +0200

Peng gravatar image

updated 2019-03-05 07:28:01 +0200

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

edit retag flag offensive close merge delete

5 answers

Sort by ยป oldest newest most voted
0

answered 2019-03-20 09:54:20 +0200

Are you using the released version available on CROCO web site or the Git Master version?

edit flag offensive delete link more

Comments

It was downloaded from CROCO web site https://www.croco-ocean.org/download/croco-project/

Peng gravatar imagePeng ( 2019-03-21 05:36:21 +0200 )edit
0

answered 2019-03-15 02:27:49 +0200

Hi Peng and Pmarches (and everyone),

I encountered a similar problem with OPENMP and gfortran in UBUNTU. Step2d velocity blowed-up between two tiles. When I turned off OPENMP, this blowing-up was gone.

Do you have more clue about this issue? Thank you. wf

edit flag offensive delete link more
0

answered 2019-03-13 15:51:22 +0200

Hi Peng, I don't see the problem. The RIP case is routinely tested for parallelization (with automatic procedure). I would advise to make sure that your jobcomp, cppdefs.h and param.h are taken from the OCEAN directory. What linux distribution are you on? ... Patrick

edit flag offensive delete link more
0

answered 2019-03-06 19:57:01 +0200

Peng gravatar image

Hi Pmarches, thank you for your reply. I'm using Linux and ifort. I did not modify jobcomp except changing the compiler to ifort. (gfortran also has this blowup issue.) I did a few tests and found that the blowup issue depends on how to divide the subdomains and number of threads. For example, set NPP=2, NSUB_X=1, NSUB_E=NPP, and export OMP_NUM_THREADS=2, the codes blow up. But set NPP=4, NSUB_X=1, NSUB_E=NPP, and export OMP_NUM_THREADS=2, the codes had no problem. I am not sure why this happens.

edit flag offensive delete link more
0

answered 2019-03-06 12:08:29 +0200

Hi Peng, I just did a test of RIP case with OpenMP (NPP=2 or 4) and had no problem. What computer and compiler are you using? Did you modify jobcomp?

edit flag offensive delete link more

Comments

Hi Pmarches, I also have similar problem. I use ubuntu 18.04 . thank you. Apri

Apriansyah gravatar imageApriansyah ( 2019-10-02 16:51:05 +0200 )edit
Login/Signup to Answer

Question Tools

2 followers

Stats

Asked: 2019-02-27 23:34:21 +0200

Seen: 10,461 times

Last updated: Mar 20 '19