42a inaccurate fdpm amplitude and phase derivative values with respect to mu a and mu s

src/Vts/Factories/ComputationFactory.cs

@@ -146,8 +146,15 @@ public static double[] ComputeReflectance(
             //    new[] { op.Mua, op.Musp, op.G, op.N }
             //double[] optimizationParameters = GetOptimizationParameters(forwardSolver, solutionDomainType, independentAxisTypes); // will need this for inverse solver
 
-            // create a list of inputs (besides optical properties) that corresponds to the behavior of the function above
+            if (IsComplexSolver(solutionDomainType) && forwardAnalysisType != ForwardAnalysisType.R) // complex derivative
+            {
+                // create a list of inputs (besides optical properties) that corresponds to the behavior of the function above
+                // create Complex derivative function with re/im and dre/dim in it
+                return func(independentValues)
+                    .Concat(func.GetDerivativeFunc(forwardAnalysisType)(independentValues)).ToArray();
+            }
 
+            // prior code
             return forwardAnalysisType == ForwardAnalysisType.R
                 ? func(independentValues)
                 : func.GetDerivativeFunc(forwardAnalysisType)(independentValues);