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

[Tetsuo Koyama]

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