public static IEnumerable <double> TemporalPointSourceCenterlinePHD(
this IForwardSolver myForwardSolver,
IEnumerable <OpticalProperties> ops, IEnumerable <double> rProbes,
IEnumerable <double> rhos, IEnumerable <double> zs, IEnumerable <double> ts)
{
foreach (var op in ops)
{
DiffusionParameters dp = DiffusionParameters.Create(op, ForwardModel.SDA);
foreach (var rProbe in rProbes)
{
foreach (var rho in rhos)
{
foreach (var z in zs)
{
var r11 = CalculatorToolbox.GetRadius(rProbe, z - dp.zp);
var r12 = CalculatorToolbox.GetRadius(rProbe, z + dp.zp + 2 * dp.zb);
var r21 = CalculatorToolbox.GetRadius(rho - rProbe, z);
var r22 = CalculatorToolbox.GetRadius(rho - rProbe, z + 2 * dp.zb);
foreach (var t in ts)
{
Func <double, double> integrandConvolve =
tau =>
{
return
(DiffusionGreensFunctions.TemporalPointSourceImageGreensFunction
(dp, r11, r12, tau) *
DiffusionGreensFunctions.TemporalPointSourceImageGreensFunction
(dp, r21, r22, t - tau));
};
yield return(FunctionMath.Integrate(
integrandConvolve, Meta.Numerics.Interval.FromEndpoints(0.0, t)));
}
}
}
}
}
}
