[Getfem-commits] [getfem-commits] branch master updated: Making work the directional derivative with interpolate transformation in a very special case

[Yves Renard]

This is an automated email from the git hooks/post-receive script.

renard pushed a commit to branch master
in repository getfem.

The following commit(s) were added to refs/heads/master by this push:
     new 2f17bad  Making work the directional derivative with interpolate 
transformation in a very special case
2f17bad is described below

commit 2f17bad1068571c33123aadbdbba5dd2ffa532a0
Author: Yves Renard <Yves.Renard@insa-lyon.fr>
AuthorDate: Sun Jan 10 20:11:44 2021 +0100

    Making work the directional derivative with interpolate transformation in a 
very special case
---
 .gitignore                              |  1 +
 src/getfem_generic_assembly_semantic.cc | 41 +++++++++++++++++++++++----------
 2 files changed, 30 insertions(+), 12 deletions(-)

diff --git a/.gitignore b/.gitignore
index 3dc5d59..ad82f14 100644
--- a/.gitignore
+++ b/.gitignore
@@ -299,6 +299,7 @@ Makefile
 /contrib/xfem_contact/xfem_dirichlet
 /contrib/xfem_contact/xfem_stokes
 /contrib/xfem_stab_unilat_contact/xfem_stab_unilat_contact
+/contrib/continuum_mechanics/results/
 /interface/tests/*/*.vtk
 /interface/tests/*/*.pos
 /interface/tests/*/*.dx
diff --git a/src/getfem_generic_assembly_semantic.cc 
b/src/getfem_generic_assembly_semantic.cc
index b2df289..dbb9029 100644
--- a/src/getfem_generic_assembly_semantic.cc
+++ b/src/getfem_generic_assembly_semantic.cc
@@ -257,28 +257,45 @@ namespace getfem {
         break;
       case GA_NODE_INTERPOLATE_VAL_TEST: case GA_NODE_INTERPOLATE_GRAD_TEST:
       case GA_NODE_INTERPOLATE_HESS_TEST: case GA_NODE_INTERPOLATE_DIVERG_TEST:
-        GMM_ASSERT1(false,
-                    "Sorry, directional derivative do not work for the moment "
-                    "with interpolate transformations. Future work.");
+        if (pexpr->node_type == GA_NODE_VAL_TEST) {
+          pnode->name = pexpr->name;
+        } else {
+          GMM_ASSERT1(false,
+                      "Sorry, directional derivative do not work for the "
+                      "moment with interpolate transformations. Future work.");
+        }
       case GA_NODE_ELEMENTARY_VAL_TEST: case GA_NODE_ELEMENTARY_GRAD_TEST:
       case GA_NODE_ELEMENTARY_HESS_TEST: case GA_NODE_ELEMENTARY_DIVERG_TEST:
-        GMM_ASSERT1(false,
-                    "Sorry, directional derivative do not work for the moment "
-                    "with elementary transformations. Future work.");
+        if (pexpr->node_type == GA_NODE_VAL_TEST) {
+          pnode->name = pexpr->name;
+        } else {
+          GMM_ASSERT1(false,
+                      "Sorry, directional derivative do not work for the "
+                      "moment with elementary transformations. Future work.");
+        }
       case GA_NODE_SECONDARY_DOMAIN_VAL_TEST:
       case GA_NODE_SECONDARY_DOMAIN_GRAD_TEST:
       case GA_NODE_SECONDARY_DOMAIN_HESS_TEST:
       case GA_NODE_SECONDARY_DOMAIN_DIVERG_TEST:
-        GMM_ASSERT1(false,
-                    "Sorry, directional derivative do not work for the moment "
-                    "with secondary domains. Future work.");
+        if (pexpr->node_type == GA_NODE_VAL_TEST) {
+          pnode->name = pexpr->name;
+        } else {
+          GMM_ASSERT1(false,
+                      "Sorry, directional derivative do not work for the "
+                      "moment with secondary domains. Future work.");
+        }
       case GA_NODE_XFEM_PLUS_VAL_TEST: case GA_NODE_XFEM_PLUS_GRAD_TEST:
       case GA_NODE_XFEM_PLUS_HESS_TEST: case GA_NODE_XFEM_PLUS_DIVERG_TEST:
       case GA_NODE_XFEM_MINUS_VAL_TEST: case GA_NODE_XFEM_MINUS_GRAD_TEST:
       case GA_NODE_XFEM_MINUS_HESS_TEST: case GA_NODE_XFEM_MINUS_DIVERG_TEST:
-        GMM_ASSERT1(false,
-                    "Sorry, directional derivative do not work for the moment "
-                    "with Xfem_plus and Xfem_minus operations. Future work.");
+        if (pexpr->node_type == GA_NODE_VAL_TEST) {
+          pnode->name = pexpr->name;
+        } else {
+          GMM_ASSERT1(false,
+                      "Sorry, directional derivative do not work for the "
+                      "moment with Xfem_plus and Xfem_minus operations. "
+                      "Future work.");
+        }
       default: break;
       }
     }