[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: [Getfem-users] Where is the source of heat_equation
From: |
Tetsuo Koyama |
Subject: |
Re: [Getfem-users] Where is the source of heat_equation |
Date: |
Wed, 20 Mar 2019 10:05:47 +0900 |
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 <address@hidden>:
>
>
> 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 <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
>
> ---------
>
>
> --
>
> 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
>
> ---------