Merge pull request #70 from vanroekel/evaluateCubicBug

Fixes a bug in evaluate cubic
CVMix · Apr 12, 2018 · 534fc38 · 534fc38
2 parents 04136cf + f6746fd
commit 534fc38
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/shared/cvmix_math.F90 b/src/shared/cvmix_math.F90
@@ -244,7 +244,7 @@ function cvmix_math_evaluate_cubic(coeffs, x_in, fprime)
       cvmix_math_evaluate_cubic = cvmix_math_evaluate_cubic +                 &
                                   coeffs(i)*(x_in**(i-1))
       if (present(fprime).and.(i.gt.2)) &
-        fprime = fprime + real(i-1,cvmix_r8)*(x_in**(i-2))
+        fprime = fprime + coeffs(i)*real(i-1,cvmix_r8)*(x_in**(i-2))
     end do
 
   end function cvmix_math_evaluate_cubic