gsl_integration_qawo_table(double omega, double L,
gsl_integration_qawo_enum sine,
int n)
{
if (n == 0)
{
throw new ArgumentOutOfRangeException("table length n must be positive integer");
}
chebmo = new double[25 * n];
this.n = n;
this.sine = sine;
this.omega = omega;
this.L = L;
par = 0.5 * omega * L;
/* precompute the moments */
{
int i;
double scale = 1.0;
for (i = 0; i < this.n; i++)
{
compute_moments(par * scale, chebmo, 25 * i);
scale *= 0.5;
}
}
}
public void set(double omega, double L, gsl_integration_qawo_enum sine)
{
this.omega = omega;
this.sine = sine;
this.L = L;
par = 0.5 * omega * L;
/* recompute the moments */
{
int i;
double scale = 1.0;
for (i = 0; i < n; i++)
{
compute_moments(par * scale, chebmo, 25 * i);
scale *= 0.5;
}
}
return; //GSL_SUCCESS;
}
public void set(double omega, double L, gsl_integration_qawo_enum sine)
{
this.omega = omega;
this.sine = sine;
this.L = L;
this.par = 0.5 * omega * L;
/* recompute the moments */
{
int i;
double scale = 1.0;
for (i = 0; i < this.n; i++)
{
compute_moments(this.par * scale, this.chebmo, 25 * i);
scale *= 0.5;
}
}
return; //GSL_SUCCESS;
}
gsl_integration_qawo_table(double omega, double L,
gsl_integration_qawo_enum sine,
int n)
{
if (n == 0)
{
throw new ArgumentOutOfRangeException("table length n must be positive integer");
}
this.chebmo = new double[25 * n];
this.n = n;
this.sine = sine;
this.omega = omega;
this.L = L;
this.par = 0.5 * omega * L;
/* precompute the moments */
{
int i;
double scale = 1.0;
for (i = 0; i < this.n; i++)
{
compute_moments(this.par * scale, this.chebmo, 25 * i);
scale *= 0.5;
}
}
}
