Bug fix to convert from Pozar coordinates to Kass coordinates.

project8 · Jul 17, 2023 · 92d49cc · 92d49cc
1 parent b8c3700
commit 92d49cc
Show file tree

Hide file tree

Showing 2 changed files with 8 additions and 8 deletions.
diff --git a/Source/Fields/LMCPozarCylindricalCavity.cc b/Source/Fields/LMCPozarCylindricalCavity.cc
@@ -25,7 +25,7 @@ namespace locust
 
     	double z = zKass + L/2.;
 
-    	if ((r > R) || (fabs(z) > L/2.))  // outside the cavity
+    	if ((r > R) || (fabs(zKass) > L/2.))  // outside the cavity
     	{
     		return {0.,0.};
     	}
@@ -66,7 +66,7 @@ namespace locust
 
     	double z = zKass + L/2.;
 
-    	if ((r > R) || (fabs(z) > L/2.))  // outside the cavity
+    	if ((r > R) || (fabs(zKass) > L/2.))  // outside the cavity
     	{
     		return {0.,0.,0.};
     	}
@@ -105,7 +105,7 @@ namespace locust
     {
     	double z = zKass + L/2.;
 
-    	if ((r > R) || (fabs(z) > L/2.))  // outside the cavity
+    	if ((r > R) || (fabs(zKass) > L/2.))  // outside the cavity
     	{
     		return {0.,0.,0.};
     	}
@@ -145,7 +145,7 @@ namespace locust
     {
     	double z = zKass + L/2.;
 
-    	if ((r > R) || (fabs(z) > L/2.))  // outside the cavity
+    	if ((r > R) || (fabs(zKass) > L/2.))  // outside the cavity
     	{
     		return {0.,0.};
     	}

diff --git a/Source/Fields/LMCPozarRectangularWaveguide.cc b/Source/Fields/LMCPozarRectangularWaveguide.cc
@@ -23,7 +23,7 @@ namespace locust
     	double x = xKass + dimX/2.;
     	double y = yKass + dimY/2.;
 
-    	if ((fabs(x) > dimX/2.) || (fabs(y) > dimY/2.))
+    	if ((fabs(xKass) > dimX/2.) || (fabs(yKass) > dimY/2.))
     	{
     		return {0.,0.};
     	}
@@ -51,7 +51,7 @@ namespace locust
     	double x = xKass + dimX/2.;
     	double y = yKass + dimY/2.;
 
-    	if ((fabs(x) > dimX/2.) || (fabs(y) > dimY/2.))
+    	if ((fabs(xKass) > dimX/2.) || (fabs(yKass) > dimY/2.))
     	{
     		return {0.,0.};
     	}
@@ -81,7 +81,7 @@ namespace locust
     	double x = xKass + dimX/2.;
     	double y = yKass + dimY/2.;
 
-    	if ((fabs(x) > dimX/2.) || (fabs(y) > dimY/2.))
+    	if ((fabs(xKass) > dimX/2.) || (fabs(yKass) > dimY/2.))
     	{
     		return {0.,0.};
     	}
@@ -109,7 +109,7 @@ namespace locust
     	double x = xKass + dimX/2.;
     	double y = yKass + dimY/2.;
 
-    	if ((fabs(x) > dimX/2.) || (fabs(y) > dimY/2.))
+    	if ((fabs(xKass) > dimX/2.) || (fabs(yKass) > dimY/2.))
     	{
     		return {0.,0.};
     	}