aboutsummaryrefslogtreecommitdiffhomepage
path: root/Eigen/src/Core/StableNorm.h
diff options
context:
space:
mode:
authorGravatar Gael Guennebaud <g.gael@free.fr>2009-07-20 10:57:31 +0200
committerGravatar Gael Guennebaud <g.gael@free.fr>2009-07-20 10:57:31 +0200
commit4c85fa8c734c8d103343c163fe95e9ffdeb9f5aa (patch)
tree39e21c0afa5de0a26be5a5747e97e78926ad29ac /Eigen/src/Core/StableNorm.h
parentc10b919edb9a362388a98f620f8cef7292ff496a (diff)
compilation fix (sun CC)
Diffstat (limited to 'Eigen/src/Core/StableNorm.h')
-rw-r--r--Eigen/src/Core/StableNorm.h10
1 files changed, 5 insertions, 5 deletions
diff --git a/Eigen/src/Core/StableNorm.h b/Eigen/src/Core/StableNorm.h
index 383c64f01..22809633d 100644
--- a/Eigen/src/Core/StableNorm.h
+++ b/Eigen/src/Core/StableNorm.h
@@ -115,17 +115,17 @@ MatrixBase<Derived>::blueNorm() const
ei_assert(false && "the algorithm cannot be guaranteed on this computer");
}
iexp = -((1-iemin)/2);
- b1 = std::pow(ibeta, iexp); // lower boundary of midrange
+ b1 = std::pow(double(ibeta),iexp); // lower boundary of midrange
iexp = (iemax + 1 - it)/2;
- b2 = std::pow(ibeta,iexp); // upper boundary of midrange
+ b2 = std::pow(double(ibeta),iexp); // upper boundary of midrange
iexp = (2-iemin)/2;
- s1m = std::pow(ibeta,iexp); // scaling factor for lower range
+ s1m = std::pow(double(ibeta),iexp); // scaling factor for lower range
iexp = - ((iemax+it)/2);
- s2m = std::pow(ibeta,iexp); // scaling factor for upper range
+ s2m = std::pow(double(ibeta),iexp); // scaling factor for upper range
overfl = rbig*s2m; // overfow boundary for abig
- eps = std::pow(ibeta, 1-it);
+ eps = std::pow(double(ibeta), 1-it);
relerr = ei_sqrt(eps); // tolerance for neglecting asml
abig = 1.0/eps - 1.0;
if (RealScalar(nbig)>abig) nmax = abig; // largest safe n