Re: Nonlinear term MPI running does not end

[Konstantinos Poulios]

Dear Tetsuo,

The relevant code is in the void model::assembly function in getfem_models.cc. The relevant code assembling the term you add with md.add_nonlinear_term(..) must be executed inside the if condition

if (generic_expressions.size()) {...}

You can have a look there and ask for further help if it looks too complex. You should also check if the test works when you run it with md.add_nonlinear_term but setting the number of MPI processes to one.

BR

Kostas

On Tue, May 11, 2021 at 10:44 AM Tetsuo Koyama <tkoyama010@gmail.com> wrote:

Dear Kostas

Thank you for your reply.

> Interesting. In order to isolate the issue, can you also check with

> md.add_linear_term(..)

> ?
It ends when using md.add_linear_term(..).

It seems that it is a problem of md.add_nonlinear_term(..).

Is there a point which I can check?

Best regards Tetsuo.

2021年5月11日(火) 17:19 Konstantinos Poulios <logari81@googlemail.com>:

Dear Tetsuo,

Interesting. In order to isolate the issue, can you also check with
md.add_linear_term(..)

?

Best regards

Kostas

On Tue, May 11, 2021 at 12:22 AM Tetsuo Koyama <tkoyama010@gmail.com> wrote:

Dear GetFEM community

I am running MPI Parallelization of GetFEM.The running command is

$ git clone https://git.savannah.nongnu.org/git/getfem.git
$ cd getfem
$ bash autogen.sh
$ ./configure --with-pic --enable-paralevel=2
$ make
$ make install
$ mpirun -n 2 python demo_parallel_laplacian.py

The python script ends correctly. But when I changed the following linear term to nonlinear term the script did not end.

-md.add_Laplacian_brick(mim, 'u')
+md.add_nonlinear_term(mim, "Grad_u.Grad_Test_u")

Do you know the reason?
Best regards Tetsuo