test

NonLinearEVP/NonLinearEVP_data.geo

@@ -3,8 +3,6 @@
 ////   _data.geo               ////
 ///////////////////////////////////
 
-deg2rad  = Pi/180;
-
 pp0        = "1Geometry/0";
 pp1        = "2Polarization-Bloch/0";
 pp3        = "3Eigenvalue problem parameters/0";
@@ -32,9 +30,9 @@ DefineConstant[
   flag_Hparallel = {1     , Name StrCat[pp1 , "1polarization case"] , Choices {0="E //",1="H //"} },
   kx             = {0.75  , Name StrCat[pp1 , "2kx [Pi\a]"] , Highlight Str[colorpp1]  , Closed close_menu} , 
 
-  eps_oo_1       = {1     , Name StrCat[pp2 , "0epsilon_inf_Drude [ - ]"] , Highlight Str[colorpp2]  , Closed close_menu} , 
-  om_d_1         = {1.1   , Name StrCat[pp2 , "1Omega_Drude [2Pic\a]"] , Highlight Str[colorpp2]  , Closed close_menu} , 
-  gam_1          = {0.05  , Name StrCat[pp2 , "2Gamma_Drude [2Pic\a]"] , Highlight Str[colorpp2]  , Closed close_menu} , 
+  eps_oo_1       = {1     , Name StrCat[pp2 , "0eps_oo_1 [ - ]"] , Highlight Str[colorpp2]  , Closed close_menu} , 
+  om_d_1         = {1.1   , Name StrCat[pp2 , "1om_d_1 [2Pic\a]"] , Highlight Str[colorpp2]  , Closed close_menu} , 
+  gam_1          = {0.05  , Name StrCat[pp2 , "2gam_1 [2Pic\a]"] , Highlight Str[colorpp2]  , Closed close_menu} , 
   
   neig          =  {1     , Name StrCat[pp3 , "0Number of eigenvalues [int]"] , Highlight Str[colorpp3]  , Closed close_menu} ,
   eig_target_re  = {0.0077, Name StrCat[pp3 , "1EV real part target [2Pic\a]"] , Highlight Str[colorpp3]  , Closed close_menu} , 
@@ -53,7 +51,7 @@ DefineConstant[
   flag_outEigvec = {1     , Name StrCat[pp4, "output eigenvector?"], Choices{0,1}}
 ];
 
-// normalized units so eigenvalues are close to unity
+
 cel      = a_lat/(2*Pi);
 epsf     = 8.854187817e-3;
 muf      = 400.*Pi;
@@ -62,7 +60,9 @@ epsilon0 = epsf*nm;
 mu0      = muf*nm;
 norm     = a_lat/(2.*Pi*cel);
 
-// Normalize
+Printf("norm %f",norm);    
+
+// Normalization
 d_sq           = d_sq          * a_lat;
 space2pml      = space2pml     * a_lat;
 pmlsize        = pmlsize       * a_lat;