试图解决负的特征值

MechDancer · Mar 13, 2020 · 798126f · 798126f
1 parent 71ff758
commit 798126f
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 2 deletions.
diff --git a/build.gradle.kts b/build.gradle.kts
@@ -8,7 +8,7 @@ plugins {
 }
 
 group = "org.mechdancer"
-version = "0.2.8-dev-2"
+version = "0.2.8-snapshot-3"
 
 repositories {
     mavenCentral()

diff --git a/src/main/kotlin/org/mechdancer/algebra/function/matrix/Jocabi.kt b/src/main/kotlin/org/mechdancer/algebra/function/matrix/Jocabi.kt
@@ -89,7 +89,7 @@ fun Matrix.svd(epsilon: Double = 1e-8): Triple<Matrix, Matrix, Matrix> {
 
     fun partEvd(a: Matrix, at: Matrix): Triple<Matrix, Matrix, Matrix> {
         val (m, square) = (a * at).evd(epsilon)!!
-        val s = square.diagonal.map { if (it < DOUBLE_PRECISION) .0 else sqrt(it) }
+        val s = square.diagonal.map { if (it < DOUBLE_PRECISION) DOUBLE_PRECISION else sqrt(it) }
         val w = listMatrixOf(row, column) { r, c -> if (r == c) s[r] else .0 }
         val others = (at * m * DiagonalMatrix(s.map { 1 / it }))
         return Triple(m, w, others)