GitHub-Proxy/eigen
1
0
Fork 0
mirror of https://gitlab.com/libeigen/eigen.git synced 2025-05-01 16:24:28 +08:00
Code Issues Packages Projects Releases Wiki Activity

compilation fix (sun CC)

Browse Source
This commit is contained in:
Gael Guennebaud 2009-07-20 10:57:31 +02:00
parent c10b919edb
commit 4c85fa8c73
1 changed files with 5 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
Eigen/src/Core/StableNorm.h
View File

@ -115,17 +115,17 @@ MatrixBase<Derived>::blueNorm() const
ei_assert(false && "the algorithm cannot be guaranteed on this computer"); ei_assert(false && "the algorithm cannot be guaranteed on this computer");
} }
iexp = -((1-iemin)/2); 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; 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; 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); 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 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 relerr = ei_sqrt(eps); // tolerance for neglecting asml
abig = 1.0/eps - 1.0; abig = 1.0/eps - 1.0;
if (RealScalar(nbig)>abig) nmax = abig; // largest safe n if (RealScalar(nbig)>abig) nmax = abig; // largest safe n