From 0abb148b7dcefc1d7ba0f39ca98b708df6152d8e Mon Sep 17 00:00:00 2001 From: Thomas Capricelli Date: Fri, 21 Aug 2009 00:27:11 +0200 Subject: [PATCH] use ei_sqrt instead of sqrt --- unsupported/Eigen/src/NonLinear/lmder.h | 2 +- unsupported/Eigen/src/NonLinear/lmdif.h | 2 +- unsupported/Eigen/src/NonLinear/lmstr.h | 2 +- 3 files changed, 3 insertions(+), 3 deletions(-) diff --git a/unsupported/Eigen/src/NonLinear/lmder.h b/unsupported/Eigen/src/NonLinear/lmder.h index 57c1038fa..f134927b1 100644 --- a/unsupported/Eigen/src/NonLinear/lmder.h +++ b/unsupported/Eigen/src/NonLinear/lmder.h @@ -261,7 +261,7 @@ L200: /* L230: */ } temp1 = ei_enorm(n, &wa3[1]) / fnorm; - temp2 = sqrt(par) * pnorm / fnorm; + temp2 = ei_sqrt(par) * pnorm / fnorm; /* Computing 2nd power */ prered = temp1 * temp1 + temp2 * temp2 / p5; dirder = -(temp1 * temp1 + temp2 * temp2); diff --git a/unsupported/Eigen/src/NonLinear/lmdif.h b/unsupported/Eigen/src/NonLinear/lmdif.h index 27b725171..3db36fcbf 100644 --- a/unsupported/Eigen/src/NonLinear/lmdif.h +++ b/unsupported/Eigen/src/NonLinear/lmdif.h @@ -263,7 +263,7 @@ L200: /* L230: */ } temp1 = ei_enorm(n, &wa3[1]) / fnorm; - temp2 = sqrt(par) * pnorm / fnorm; + temp2 = ei_sqrt(par) * pnorm / fnorm; /* Computing 2nd power */ prered = temp1 * temp1 + temp2 * temp2 / p5; dirder = -(temp1 * temp1 + temp2 * temp2); diff --git a/unsupported/Eigen/src/NonLinear/lmstr.h b/unsupported/Eigen/src/NonLinear/lmstr.h index d2810ab31..e18ee4cb0 100644 --- a/unsupported/Eigen/src/NonLinear/lmstr.h +++ b/unsupported/Eigen/src/NonLinear/lmstr.h @@ -285,7 +285,7 @@ L240: /* L270: */ } temp1 = ei_enorm(n, &wa3[1]) / fnorm; - temp2 = sqrt(par) * pnorm / fnorm; + temp2 = ei_sqrt(par) * pnorm / fnorm; /* Computing 2nd power */ prered = temp1 * temp1 + temp2 * temp2 / p5; dirder = -(temp1 * temp1 + temp2 * temp2);