public double Interp(double mu0, double lambda0)
{
PolynomInterpolator pi;
double deltaLambda = this.lambda[1] - this.lambda[0];
int i0 = (int)((lambda0 - this.lambda[0]) / deltaLambda);
int i1 = i0 + 1;
double[] intenSet0 = new double[3];
double[] intenSet1 = new double[3];
intenSet0[0] = this.specGrid[0][i0];
intenSet0[1] = this.specGrid[1][i0];
intenSet0[2] = this.specGrid[2][i0];
intenSet1[0] = this.specGrid[0][i1];
intenSet1[1] = this.specGrid[1][i1];
intenSet1[2] = this.specGrid[2][i1];
pi = new PolynomInterpolator(this.mu, intenSet0);
double inten0 = pi.Interp(mu0);
double inten1 = pi.Interp(mu0);
return(inten0 + (inten1 - inten0) * (lambda0 - this.lambda[i0]) / deltaLambda);
}
public double Interp(double mu0, double lambda0)
{
PolynomInterpolator pi0, pi1;
double deltaLambda = this.lambda[1] - this.lambda[0];
int i0 = (int)((lambda0 - this.lambda[0]) / deltaLambda);
int i1 = i0 + 1;
double[] intenSet0 = new double[3];
double[] intenSet1 = new double[3];
if (i0 < 0)
{
intenSet0[0] = this.specGrid[0][0];
intenSet0[1] = this.specGrid[1][0];
intenSet0[2] = this.specGrid[2][0];
pi0 = new PolynomInterpolator(this.mu, intenSet0);
return(pi0.Interp(mu0));
}
if (i1 > specGrid[0].Length - 1)
{
intenSet0[0] = this.specGrid[0][specGrid[0].Length - 1];
intenSet0[1] = this.specGrid[1][specGrid[1].Length - 1];
intenSet0[2] = this.specGrid[2][specGrid[2].Length - 1];
pi0 = new PolynomInterpolator(this.mu, intenSet0);
return(pi0.Interp(mu0));
}
intenSet0[0] = this.specGrid[0][i0];
intenSet0[1] = this.specGrid[1][i0];
intenSet0[2] = this.specGrid[2][i0];
intenSet1[0] = this.specGrid[0][i1];
intenSet1[1] = this.specGrid[1][i1];
intenSet1[2] = this.specGrid[2][i1];
pi0 = new PolynomInterpolator(this.mu, intenSet0);
pi1 = new PolynomInterpolator(this.mu, intenSet1);
double inten0 = pi0.Interp(mu0);
double inten1 = pi1.Interp(mu0);
return(inten0 + (inten1 - inten0) * (lambda0 - this.lambda[i0]) / deltaLambda);
}