[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
## Re: [Getfem-users] Where is the source of heat_equation

**From**: |
Yves Renard |

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

**Date**: |
Fri, 15 Mar 2019 15:39:20 +0100 |

**User-agent**: |
Mozilla/5.0 (X11; Linux x86_64; rv:60.0) Gecko/20100101 Thunderbird/60.5.1 |

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