Dear Yves Thank you for your reply.
I will try to isolate the problem and investigate the cause. I will contact you when I can identify the problem. Best Regards, Tetsuo 2019年3月20日(水) 0:28 Yves Renard <[email protected]>: > > > 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 <[email protected]>: > > 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 > > --------- > > > -- > > 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 > > ---------
