Commit 1f60bd18 authored by Guillaume Demesy's avatar Guillaume Demesy

Revert "test"

This reverts commit dc040487.
parent dc040487
......@@ -3,6 +3,8 @@
//// _data.geo ////
///////////////////////////////////
deg2rad = Pi/180;
pp0 = "1Geometry/0";
pp1 = "2Polarization-Bloch/0";
pp3 = "3Eigenvalue problem parameters/0";
......@@ -30,9 +32,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 , "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} ,
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} ,
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} ,
......@@ -51,7 +53,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;
......@@ -60,9 +62,7 @@ epsilon0 = epsf*nm;
mu0 = muf*nm;
norm = a_lat/(2.*Pi*cel);
Printf("norm %f",norm);
// Normalization
// Normalize
d_sq = d_sq * a_lat;
space2pml = space2pml * a_lat;
pmlsize = pmlsize * a_lat;
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment