|
From: | Yves Renard |
Subject: | Re: [Getfem-users] divergence of strain tensor |
Date: | Fri, 13 Dec 2013 15:28:19 +0100 |
User-agent: | Mozilla/5.0 (X11; Linux x86_64; rv:24.0) Gecko/20100101 Thunderbird/24.1.1 |
Dear Wen, If you really need to compute the divergence operator (it is rare to really have to compute it, except for intance for residual a posteriori estimators) you can use compute_hessian in getfem_derivatives.h, yes. The derivatives are stored in Fortran order (\partial^2 u / \partial_1 \partial_1, \partial^2 u / \partial_1 \partial_2, ... \partial^2 u / \partial_1 \partial_N, \partial^2 u / \partial_2 \partial_1) on each finite element node on the finite element method it is interpolated and gathered on the whole vector (i.e. derivative on the first dof, then derivatives on the second dof, and so on). Yves. Le 12/12/2013 16:46, Wen Jiang 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] |