default to constant PML in har3d

simulation/HarocheHelper.cpp

@@ -41,24 +41,24 @@ Complex PML::dampingX(fullVector<double>& xyz){
   double     f = Xmax + SizeX - fabs(xyz(0));
   double overF = 1 / (f) - 1 / (SizeX);
 
-  return Complex(1, -overF / kHaroche);
-  //return a;
+  //return Complex(1, -overF / kHaroche);
+  return a;
 }
 
 Complex PML::dampingY(fullVector<double>& xyz){
   double     f = Ymax + SizeY - fabs(xyz(1));
   double overF = 1 / (f) - 1 / (SizeY);
 
-  return Complex(1, -overF / kHaroche);
-  //return a;
+  //return Complex(1, -overF / kHaroche);
+  return a;
 }
 
 Complex PML::dampingZ(fullVector<double>& xyz){
   double     f = Zmax + SizeZ - fabs(xyz(2));
   double overF = 1 / (f) - 1 / (SizeZ);
 
-  return Complex(1, -overF / kHaroche);
-  //return a;
+  //return Complex(1, -overF / kHaroche);
+  return a;
 }
 
 // Air