[Getfem-commits] r5227 - in /trunk/getfem: doc/sphinx/source/userdoc/ src/

[Yves . Renard]

Author: renard
Date: Wed Jan 20 11:13:44 2016
New Revision: 5227

URL: http://svn.gna.org/viewcvs/getfem?rev=5227&view=rev
Log:
minor fixes

Modified:
    trunk/getfem/doc/sphinx/source/userdoc/gasm_high.rst
    trunk/getfem/doc/sphinx/source/userdoc/rmesh.rst
    trunk/getfem/src/getfem_nonlinear_elasticity.cc

Modified: trunk/getfem/doc/sphinx/source/userdoc/gasm_high.rst
URL: 
http://svn.gna.org/viewcvs/getfem/trunk/getfem/doc/sphinx/source/userdoc/gasm_high.rst?rev=5227&r1=5226&r2=5227&view=diff
==============================================================================
--- trunk/getfem/doc/sphinx/source/userdoc/gasm_high.rst        (original)
+++ trunk/getfem/doc/sphinx/source/userdoc/gasm_high.rst        Wed Jan 20 
11:13:44 2016
@@ -671,7 +671,7 @@
 
 **CAUTION**: You have to think that when some variables are used in the 
transformation, the computation of the tangent system takes into account these 
dependence. However, the second derivative of a transformation with respect to 
a variable used has not been implemented. Thus, such a transformation is not 
allowed in the definition of a potential since it cannot be derived twice.
 
-.. _ud-gasm-high-elem-trans:
+.. _ud-gasm-high-inter-elt-disc:
 
 Evaluating discontinuities across inter-element edges/faces
 -----------------------------------------------------------
@@ -709,23 +709,6 @@
 
 Compared to other interpolate transformations, this transformation is more 
optimized and benefits from finite element and geometric transformation 
pre-computations.
 
-Evaluating discontinuities across inter-element edges/faces
------------------------------------------------------------
-
-A specific interpolate transformation (see previous section), called 
``neighbour_elt`` is defined by default in all models. This transformation can 
only be used when a computation is made on an internal edge/face of a mesh, 
i.e. an element face shared at least by two elements. It aims to compute 
discontinuity jumps of a variable across inter-element faces. It is 
particularly suitable to implement Discontinuous Galerkin and interior penalty 
methods, Ghost penalty terms or a posteriori estimators. The expressions::
-
-  Interpolate(Normal, neighbour_elt)
-  Interpolate(X, neighbour_elt)
-  Interpolate(u, neighbour_elt)
-  Interpolate(Grad_u, neighbour_elt)
-  Interpolate(Div_u, neighbour_elt)
-  Interpolate(Hess_u, neighbour_elt)
-  Interpolate(Test_u, neighbour_elt)
-  Interpolate(Grad_Test_u, neighbour_elt)
-  Interpolate(Div_Test_u, neighbour_elt)
-  Interpolate(Hess_Test_u, neighbour_elt)
-
-
 .. _ud-gasm-high-elem-trans:
 
 Elementary transformations

Modified: trunk/getfem/doc/sphinx/source/userdoc/rmesh.rst
URL: 
http://svn.gna.org/viewcvs/getfem/trunk/getfem/doc/sphinx/source/userdoc/rmesh.rst?rev=5227&r1=5226&r2=5227&view=diff
==============================================================================
--- trunk/getfem/doc/sphinx/source/userdoc/rmesh.rst    (original)
+++ trunk/getfem/doc/sphinx/source/userdoc/rmesh.rst    Wed Jan 20 11:13:44 2016
@@ -51,4 +51,4 @@
    \int_e |\hspace{0.01em}[\hspace{-0.12em}[
    \partial_n u ]\hspace{-0.12em}]\hspace{0.01em}|^2 d \Gamma,
 
-where :math:`[\hspace{-0.12em}[\partial_n u]\hspace{-0.12em}]` is the jump of 
the normal derivative. Then, the error estimate for a given element is the sum 
of the computed quantities on each internal face multiplied by the element 
diameter. This basic error estimate can be taken as a model for more elaborated 
ones. It uses the high-level generic assembly and the ``neighbour_elt`` 
interpolate transformation (see :ref:`ud-gasm-high-elem-trans`).
+where :math:`[\hspace{-0.12em}[\partial_n u]\hspace{-0.12em}]` is the jump of 
the normal derivative. Then, the error estimate for a given element is the sum 
of the computed quantities on each internal face multiplied by the element 
diameter. This basic error estimate can be taken as a model for more elaborated 
ones. It uses the high-level generic assembly and the ``neighbour_elt`` 
interpolate transformation (see :ref:`ud-gasm-high-inter-elt-disc`).

Modified: trunk/getfem/src/getfem_nonlinear_elasticity.cc
URL: 
http://svn.gna.org/viewcvs/getfem/trunk/getfem/src/getfem_nonlinear_elasticity.cc?rev=5227&r1=5226&r2=5227&view=diff
==============================================================================
--- trunk/getfem/src/getfem_nonlinear_elasticity.cc     (original)
+++ trunk/getfem/src/getfem_nonlinear_elasticity.cc     Wed Jan 20 11:13:44 2016
@@ -603,7 +603,7 @@
   scalar_type Neo_Hookean_hyperelastic_law::strain_energy
   (const base_matrix &E, const base_vector &params,
    scalar_type det_trans) const {
-    if (compressible && det_trans <= scalar_type(0)) return 1e200;
+    if (det_trans <= scalar_type(0)) return 1e200;
     size_type N = gmm::mat_nrows(E);
     GMM_ASSERT1(N == 3, "Neo Hookean hyperelastic law only defined "
                 "on dimension 3, sorry");