Skip to content
Snippets Groups Projects
Commit 9aad97c5 authored by Nicolas Marsic's avatar Nicolas Marsic
Browse files

default to constant PML in har3d

parent a9e8c149
No related branches found
No related tags found
No related merge requests found
......@@ -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
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment