kscience/kmath
9
0
Fork 1
Code Issues Pull Requests 6 Packages Projects Releases 16 Wiki Activity

fixed error

Browse Source
This commit is contained in:
margarita0303 2022-07-26 21:56:50 +03:00
parent 9a46dd8966
commit f0e445587d
1 changed files with 2 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
kmath-tensors/src/commonMain/kotlin/space/kscience/kmath/tensors/core/internal/linUtils.kt
View File

@ -549,12 +549,13 @@ internal fun MutableStructure2D<Double>.svdHelper(u: MutableStructure2D<Double>,
if (flag != 0) {
c = 0.0
s = 1.0
for (i in l until k) {
for (i in l until k + 1) {
f = s * rv1[i]
rv1[i] = c * rv1[i]
if (abs(f) + anorm == anorm) {
break
}
g = w.mutableBuffer.array()[w.bufferStart + i]
h = pythag(f, g)
w.mutableBuffer.array()[w.bufferStart + i] = h
h = 1.0 / h
@ -579,11 +580,6 @@ internal fun MutableStructure2D<Double>.svdHelper(u: MutableStructure2D<Double>,
break
}
// надо придумать, что сделать - выкинуть ошибку?
// if (its == 30) {
// return
// }
x = w.mutableBuffer.array()[w.bufferStart + l]
nm = k - 1
y = w.mutableBuffer.array()[w.bufferStart + nm]