|
From: | Yves Renard |
Subject: | Re: [Getfem-users] Curvilinear structures in Getfem |
Date: | Tue, 21 Nov 2017 10:12:44 +0100 |
User-agent: | Mozilla/5.0 (X11; Linux x86_64; rv:52.0) Gecko/20100101 Thunderbird/52.4.0 |
Dear Jean-François, The problem with the _expression_ "0.5*sqr(Normalized(element_K).(Grad_u*Normalized(element_K)))" was that element_K is in fact declared as a matrix ( a 2x1 matrix here) and the dot product is restricted to vectors. I commited a fix to transform element_K as a vector for 1D elements. You can download this fix or use ":" instead of "." for the moment. In fact, in that case, I think taht your _expression_ is equivalent to "Grad_u:Grad_Test_u" since the gradient is along the tangent. Just a remark. It is of course allowed to use potential (it will be derived twice to obtain the tangent matrix), but I recommend to gives the weak formulation _expression_ not to have some surprise, especially for coupled problems (for instance if a coefficient is depending on another variable). Best regards, Yves Le 21/11/2017 à 00:41, Jean-François Barthélémy a écrit :
-- Yves Renard (address@hidden) tel : (33) 04.72.43.87.08 Pole de Mathematiques, INSA-Lyon fax : (33) 04.72.43.85.29 20, rue Albert Einstein 69621 Villeurbanne Cedex, FRANCE http://math.univ-lyon1.fr/~renard --------- |
[Prev in Thread] | Current Thread | [Next in Thread] |