Changeset 42827 for trunk/psLib/src/math/psMatrix.c

Timestamp:

May 8, 2025, 4:49:24 PM (12 months ago)

Author:

eugene

Message:

allow larger dynamic range for F64 matrix inversion

File:

• trunk/psLib/src/math/psMatrix.c (modified) (3 diffs)

trunk/psLib/src/math/psMatrix.c

@@ -459 +459 @@
 // (based on version by William Kahan -- see Ohana/src/libohana/doc/kahan-gji.pdf)
 
-# define MAX_RANGE 1.0e7
+# define MAX_RANGE_F64 1.0e16
+# define MAX_RANGE_F32 1.0e8
 // MAX_RANGE is used to test for ill-conditioned input matrices.  For an ill-conditioned
 // matrix, one or more of the pivots trends towards zero, and growth goes to infinity.  Rather
@@ -474 +475 @@
     PS_ASSERT_INT_EQUAL(a->numCols, a->numRows, false);
     PS_ASSERT_VECTOR_SIZE(b, (long int)a->numCols, false);
+
+    psF32 MAX_RANGE = MAX_RANGE_F32;
 
     // Check for non-finite entries in matrix
@@ -499 +502 @@
               }
           }
-          break;
+          MAX_RANGE = MAX_RANGE_F64;
+          break;
       }
       default: