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 <[email protected]>: > > 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 <[email protected]>: > > > > > > 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 ([email protected]) tel : (33) 04.72.43.87.08 > > INSA-Lyon > > 20, rue Albert Einstein > > 69621 Villeurbanne Cedex, FRANCE > > http://math.univ-lyon1.fr/~renard > > > > --------- > >
