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

[Yves Renard]

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

---------