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 <tkoyama...@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 <logar...@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 <tkoyama...@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 >>> >>