Re: [Getfem-users] Where is the source of heat_equation

[Yves Renard]

Dear Tetsuo,

No, it's not the expected behavior. I do not see why a term would be avoided if it is added after the addition of the Newmark scheme.

May be if you can reduce your program as the minimal program that exhibit the error, I could have a look to it.

Best Regards,

Yves

Le 17/03/2019 à 07:07, Tetsuo Koyama a écrit :

Dear Yves and GetFEM++ Project

Thank you for your support. I could solve the problem.
The reason is that I added the source term brick before "add_Newmark_scheme".
After I added the source term brick before adding Newmark scheme, I
could assemble it.
Is this a specification or a bug? Can you give me your opinion?

Best Regard Tetsuo.

2019年3月16日(土) 12:09 Tetsuo Koyama <address@hidden>:

Dear Yves

Thank you for your advice.
But I'm sorry that rhs is still zero (I also checked the case which is
not use the region).
I will continue investigation. If you notice something, please contact me.

Best Regard Tetsuo.

2019年3月15日(金) 14:55 Tetsuo Koyama <address@hidden>:

Dear Yves

Thank you for always teaching me getfem.

"[1,0,0].Test_u"

I'll try it. Sorry for my misunderstood and email subject was no correct.

Best Regard,Tetsuo

2019年3月15日(金) 23:39 Yves Renard <address@hidden>:

Dear Tetsuo,

Your term

getfem::add_source_term_generic_assembly_brick(model, mim, "1.0*Test_u",
TOP);

is a priori correct and should add a source term on the region TOP
(Neumann condition), of course if the region is non-empty.

Apart from the fact that u is a vector field in your case, so that it
should complain on the fact that "1*Test_u" is not a scalar. A correct
_expression_ should be for instance "[1,0,0].Test_u"

Best Regards,

Yves


Le 15/03/2019 à 10:38, Tetsuo Koyama a écrit :

Dear GetFEM++ Project

Hello. Thank you for your support. I could assemble the matrix in model object.
Today's question is that I couldn't build the rhs of equation.

I coded rhs's source_term_generic_assembly_brick by reading the
following document.
http://getfem.org/userdoc/model_time_integration.html#small-example-heat-equation
I added the source_term_generic_assembly_brick at line 432 for test.
But, the debugging output is zero in line 550. And the message's seem
not build the source_term_generic_assembly_brick.

I think my use it not correct. Could you teach me why it not work?
Thank you for reading.

Best Regard Tetsuo.

--

   Yves Renard (address@hidden)       tel : (33) 04.72.43.87.08
   INSA-Lyon
   20, rue Albert Einstein
   69621 Villeurbanne Cedex, FRANCE
   http://math.univ-lyon1.fr/~renard

---------

-- 

  Yves Renard (address@hidden)       tel : (33) 04.72.43.87.08
  INSA-Lyon 
  20, rue Albert Einstein
  69621 Villeurbanne Cedex, FRANCE
  http://math.univ-lyon1.fr/~renard

---------