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

[Tetsuo Koyama]

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