From 0e163e0f41146477c7c1dbeeedae6f79b7a22cdb Mon Sep 17 00:00:00 2001
From: Guillaume Demesy <guillaume.demesy@fresnel.fr>
Date: Thu, 19 Sep 2019 11:26:55 +0200
Subject: [PATCH] bug: is it due to BF_PerpendicularEdge* - adding Form0

---
 NonLinearEVP/NonLinearEVP.pro | 24 +++++++++++++++++++++++-
 1 file changed, 23 insertions(+), 1 deletion(-)

diff --git a/NonLinearEVP/NonLinearEVP.pro b/NonLinearEVP/NonLinearEVP.pro
index 5edd798..0e6956a 100644
--- a/NonLinearEVP/NonLinearEVP.pro
+++ b/NonLinearEVP/NonLinearEVP.pro
@@ -181,6 +181,27 @@ FunctionSpace {
     }
   }
 
+  { Name Hgrad; Type Form0;
+    BasisFunction {
+      { Name un;  NameOfCoef un;  Function BF_Node   ; Support Region[Om]; Entity NodesOf[All]; }
+      { Name un2; NameOfCoef un2; Function BF_Node_2E; Support Region[Om]; Entity EdgesOf[All]; }
+      If(flag_o2i==1)
+        { Name un3; NameOfCoef un3; Function BF_Node_2F; Support Region[Om]; Entity FacetsOf[All]; }
+        { Name un4; NameOfCoef un4; Function BF_Node_3E; Support Region[Om]; Entity EdgesOf[All]; }
+        { Name un5; NameOfCoef un5; Function BF_Node_3F; Support Region[Om]; Entity FacetsOf[All]; }
+      EndIf
+     }
+    Constraint {
+      { NameOfCoef un;  EntityType NodesOf ; NameOfConstraint BlochX; }
+      { NameOfCoef un2; EntityType EdgesOf ; NameOfConstraint BlochX; }
+      If(flag_o2i==1)
+        { NameOfCoef un3; EntityType FacetsOf ; NameOfConstraint BlochX; }
+        { NameOfCoef un4; EntityType EdgesOf  ; NameOfConstraint BlochX; }
+        { NameOfCoef un5; EntityType FacetsOf ; NameOfConstraint BlochX; }
+      EndIf
+    }
+  }
+
   { Name Hcurl; Type Form1;
     BasisFunction {
       { Name sn;  NameOfCoef un;  Function BF_Edge   ; Support Region[Om]; Entity EdgesOf[All]; }
@@ -333,7 +354,8 @@ Formulation {
     }
   }
   {Name modal_helmholtz_PEP_h; Type FemEquation;
-    Quantity {{ Name u   ; Type Local; NameOfSpace Hgrad_perp  ;}}
+    // Quantity {{ Name u   ; Type Local; NameOfSpace Hgrad_perp  ;}}
+    Quantity {{ Name u   ; Type Local; NameOfSpace Hgrad  ;}}
     Equation {
       Galerkin {    [-1/epsr_nod[]* om_d_1^2     * Dof{Curl u}, {Curl u} ];          In Om_2; Jacobian JSur; Integration Int_1;}
       Galerkin { Eig[ 1/epsr_nod[]*eps_oo_1*gam_1* Dof{Curl u}, {Curl u} ]; Order 1; In Om_2; Jacobian JSur; Integration Int_1;}
-- 
GitLab