sym mesh

NonLinearEVP/NonLinearEVP.geo

@@ -17,11 +17,11 @@
 
 Include "NonLinearEVP_data.geo";
 lc_cell = a_lat/paramaille;
-lc_sq   = lc_cell/4;
+lc_sq   = lc_cell/6;
 lc_pml  = lc_cell*1.2;
 lc_sqa  = lc_sq;
 lc_sq_inside = lc_sq*1.4;
-epsc         = lc_sq*0.8;
+epsc         = lc_sq*0.9;
 
 
 If (flag_Tmesh==0)
@@ -171,10 +171,10 @@ If (flag_Tmesh==1)
   Point(3)  = { a_lat/2. , d_sq/2.+space2pml, 0. , lc_cell};
   Point(4)  = {-a_lat/2. , d_sq/2.+space2pml, 0. , lc_cell};
 
-  Point(5)  = {-d_sq/2. ,-d_sq/2., 0. , lc_sq*15};
-  Point(6)  = { d_sq/2. ,-d_sq/2., 0. , lc_sq*15};
-  Point(7)  = { d_sq/2. , d_sq/2., 0. , lc_sq*15};
-  Point(8)  = {-d_sq/2. , d_sq/2., 0. , lc_sq*15};
+  Point(5)  = {-d_sq/2. ,-d_sq/2., 0. , lc_sq};
+  Point(6)  = { d_sq/2. ,-d_sq/2., 0. , lc_sq};
+  Point(7)  = { d_sq/2. , d_sq/2., 0. , lc_sq};
+  Point(8)  = {-d_sq/2. , d_sq/2., 0. , lc_sq};
 
   Point(9)  = {-a_lat/2. ,-d_sq/2.-space2pml-pmlsize, 0. , lc_pml};
   Point(10) = { a_lat/2. ,-d_sq/2.-space2pml-pmlsize, 0. , lc_pml};
@@ -419,17 +419,17 @@ If (flag_Tmesh==1)
   Periodic Line {95} = {38}  Translate {0,epsc,0} ;
 
   // Transfinite Surface { expression-list } | "*" < = { expression-list } > < Left | Right | Alternate | AlternateRight | AlternateLeft > ;
-  Transfinite Surface { 9 } Left;
-  Transfinite Surface { 10 } Right;
+  Transfinite Surface { 9 }  AlternateLeft;
+  Transfinite Surface { 10 } AlternateRight;
 
-  Transfinite Surface { 19 } Left;
-  Transfinite Surface { 20 } Left;
+  Transfinite Surface { 19 } AlternateLeft;
+  Transfinite Surface { 20 } AlternateLeft;
 
-  Transfinite Surface { 29 } Left;
-  Transfinite Surface { 30 } Left;
+  Transfinite Surface { 29 } AlternateLeft;
+  Transfinite Surface { 30 } AlternateLeft;
 
-  Transfinite Surface { 39 } Left;
-  Transfinite Surface { 40 } Right;
+  Transfinite Surface { 39 } AlternateLeft;
+  Transfinite Surface { 40 } AlternateRight;
 
 
   Physical Surface(100) = {7,8,9,17,18,19,23,24,29,37,38,40,46};   // 1 dom in

NonLinearEVP/NonLinearEVP_data.geo

@@ -53,8 +53,8 @@ DefineConstant[
   eig_min_im     = {0.0005, Name StrCat[pp3 , "5EV imag min [2Pic\a]"] , Highlight Str[colorpp3]  , Closed close_menu} , 
   eig_max_im     = {4.    , Name StrCat[pp3 , "6EV imag max [2Pic\a]"] , Highlight Str[colorpp3]  , Closed close_menu} , 
 
-  paramaille     = {4     , Name StrCat[pp4 , "0number of mesh elements per period []"]  , Highlight Str[colorpp4], Min 2, Max 10} , 
-  flag_Tmesh     = {0     , Name StrCat[pp4 , "2locally structured mesh?"] , Choices {0="unstruct",1="struct"} },
+  paramaille     = {5     , Name StrCat[pp4 , "0number of mesh elements per period []"]  , Highlight Str[colorpp4], Min 2, Max 10} , 
+  flag_Tmesh     = {1     , Name StrCat[pp4 , "2locally structured mesh?"] , Choices {0="unstruct",1="struct"} },
   flag_o2g       = {1     , Name StrCat[pp4 , "3Geometrical order"] , Choices {0="order 1 (linear)",1="order 2 (curved)"} },
   flag_o2i       = {1     , Name StrCat[pp4 , "4Interpolation order"] , Choices {0="order 1",1="full order 2"}, ServerAction "ResetDatabase"},
   flag_rounding  = {1     , Name StrCat[pp4 , "5Corner rounding/0Enable rounding of corners!"], Choices{0,1}, ServerAction "ResetDatabase"},