Commit e78eba44 authored by Axel Modave's avatar Axel Modave

analytics: remove useless function that I added a while ago

parent 56d7baaf
Pipeline #1884 passed with stage
in 47 seconds
......@@ -1146,7 +1146,6 @@ struct StringXFunction2Nbr F_Function[] = { /* #Par #Arg */
{"AcousticFieldHardSphere", (CAST)F_AcousticFieldHardSphere, 2, 1 },
{"RCSHardSphere", (CAST)F_RCSHardSphere, 2, 1 },
{"AcousticFieldSoftCylinder", (CAST)F_AcousticFieldSoftCylinder, -1, 1 },
{"CylindricalHarmonic", (CAST)F_CylindricalHarmonic, 3, 1 },
{"AcousticFieldSoftCylinderABC", (CAST)F_AcousticFieldSoftCylinderABC, 5, 1 },
{"DrAcousticFieldSoftCylinder", (CAST)F_DrAcousticFieldSoftCylinder, 2, 1 },
{"RCSSoftCylinder", (CAST)F_RCSSoftCylinder, 2, 1 },
......
......@@ -42,7 +42,6 @@ void F_DrAcousticFieldSoftSphere(F_ARG) ;
void F_RCSSoftSphere(F_ARG) ;
void F_AcousticFieldHardSphere(F_ARG) ;
void F_RCSHardSphere(F_ARG) ;
void F_CylindricalHarmonic(F_ARG) ;
void F_AcousticFieldSoftCylinder(F_ARG) ;
void F_AcousticFieldSoftCylinderABC(F_ARG) ;
void F_DrAcousticFieldSoftCylinder(F_ARG) ;
......
......@@ -1635,26 +1635,6 @@ void F_JFIE_TransZPolCyl(F_ARG)
V->Type = SCALAR ;
}
/* Cylindrical harmonics normalized by the max. value at radius R.
Returns the mode 'm' */
void F_CylindricalHarmonic(F_ARG)
{
double theta = atan2(A->Val[1], A->Val[0]);
double r = sqrt(A->Val[0]*A->Val[0] + A->Val[1]*A->Val[1]);
double k = Fct->Para[0];
double R = Fct->Para[1];
int m = Fct->Para[2];
double kr = k*r;
double kR = k*R;
std::complex<double> HnkR( jn(m,kR), yn(m,kR) );
std::complex<double> Hnkr( jn(m,kr), yn(m,kr) );
std::complex<double> val = cos(m*theta) * Hnkr/HnkR;
V->Val[0] = std::real(val);
V->Val[MAX_DIM] = std::imag(val);
V->Type = SCALAR ;
}
/* Scattering by acoustically soft circular cylinder of radius R,
under plane wave incidence e^{ikx}. Returns scatterered field
outside */
......
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