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
>>>
>>

Reply via email to