Merge branch 'master' of https://github.com/VirtualPhotonics/vts

VirtualPhotonics · Mar 27, 2019 · 2481087 · 2481087
2 parents 8d48b01 + 33b28c7
commit 2481087
Showing 1 changed file with 8 additions and 1 deletion.
diff --git a/src/Vts/Modeling/ForwardSolvers/TwoLayerSDAForwardSolver.cs b/src/Vts/Modeling/ForwardSolvers/TwoLayerSDAForwardSolver.cs
@@ -281,11 +281,14 @@ public override IEnumerable<double> FluenceOfRhoAndZ(
             IEnumerable<double> rhos,
             IEnumerable<double> zs)
         {
-
             foreach (var region in regions)
             {
                 var dp = GetDiffusionParameters(region);
                 var layerThicknesses = GetLayerThicknesses(region);
+                if (dp[0].zp > layerThicknesses[0])
+                {
+                    throw new ArgumentException("Top layer thickness must be greater than l* = 1/(mua+musp)");
+                }
                 foreach (var rho in rhos)
                 {
                     foreach (var z in zs)
@@ -436,6 +439,10 @@ public override IEnumerable<Complex> FluenceOfRhoAndZAndFt(
             {
                 var dps  = GetDiffusionParameters(region);
                 var layerThicknesses = GetLayerThicknesses(region);
+                if (dps[0].zp > layerThicknesses[0])
+                {
+                    throw new ArgumentException("Top layer thickness must be greater than l* = 1/(mua+musp)");
+                }
                 foreach (var rho in rhos)
                 {
                     foreach (var z in zs)