--- trunk/src/lib/common/eigen3x3.c 2011/12/18 13:17:09 1678
+++ trunk/src/lib/common/eigen3x3.c 2011/12/19 20:31:04 1679
@@ -200,7 +200,7 @@
D = QQQ - R*R;
if (D > epsilon) {
/* three distinct roots- this is the most common case */
- double mm, ss, cc;
+ Diderot_real_t mm, ss, cc;
theta = ATAN2(SQRT(D), R)/3.0;
mm = SQRT(Q);
ss = SIN(theta);
@@ -212,7 +212,7 @@
}
/* else D is near enough to zero */
else if (R < -epsilon || epsilon < R) {
- double U;
+ Diderot_real_t U;
/* one double root and one single root */
U = CBRT(R); /* cube root function */
if (U > 0) {
@@ -300,7 +300,7 @@
D = QQQ - R*R;
if (D > epsilon) {
/* three distinct roots- this is the most common case */
- double mm, ss, cc;
+ Diderot_real_t mm, ss, cc;
theta = ATAN2(SQRT(D), R)/3.0;
mm = SQRT(Q);
ss = SIN(theta);
@@ -312,7 +312,7 @@
}
/* else D is near enough to zero */
else if (R < -epsilon || epsilon < R) {
- double U;
+ Diderot_real_t U;
/* one double root and one single root */
U = CBRT(R); /* cube root function */
if (U > 0) {