public void TemporalFrequencyPointSourceGreensFunction_Test()
{
double[] realGreensFunctionValues = new double[] { 0.195264, 0.0411688, 0.00145147 };
double[] imaginaryGreensFunctionValues = new double[] { -0.0212481, -0.0138797, -0.00274176 };
for (int iR = 0; iR < rTestValues.Length; iR++)
{
var tfpsGF = DiffusionGreensFunctions.TemporalFrequencyPointSourceGreensFunction(dp, rTestValues[iR], k);
var relDiffReal = Math.Abs(tfpsGF.Real - realGreensFunctionValues[iR]) / realGreensFunctionValues[iR];
Assert.IsTrue(relDiffReal < thresholdValue, "Test failed for r =" + rTestValues[iR] +
"mm, with relative difference " + relDiffReal + " for the real compoment");
var relDiffImag = Math.Abs(tfpsGF.Imaginary - imaginaryGreensFunctionValues[iR]) / imaginaryGreensFunctionValues[iR];
Assert.IsTrue(relDiffImag < thresholdValue, "Test failed for r =" + rTestValues[iR] +
"mm, with relative difference " + relDiffReal + " for the imaginary compoment");
}
}
