restore la/matrix_ops.v

spytheman · spytheman · commit 7b1688b4ceb1 · 2024-09-14T20:31:54.000+03:00
diff --git a/la/matrix_ops.v b/la/matrix_ops.v
@@ -90,7 +90,7 @@ pub fn matrix_svd(mut s []f64, mut u Matrix[f64], mut vt Matrix[f64], mut a Matr
 	if copy_a {
 		acpy = a.clone()
 	}
-	vlas.dgesvd(c'A', c'A', a.m, a.n, acpy.data, 1, s, u.data, a.m, vt.data, a.n, superb)
+	vlas.dgesvd(`A`, `A`, a.m, a.n, acpy.data, 1, s, u.data, a.m, vt.data, a.n, superb)
 }
 
 // matrix_inv computes the inverse of a general matrix (square or not). It also computes the

Original file line number	Diff line number	Diff line change
`@@ -90,7 +90,7 @@ pub fn matrix_svd(mut s []f64, mut u Matrix[f64], mut vt Matrix[f64], mut a Matr`
`90`	`90`	`if copy_a {`
`91`	`91`	`acpy = a.clone()`
`92`	`92`	`}`
`93`		`- vlas.dgesvd(c'A', c'A', a.m, a.n, acpy.data, 1, s, u.data, a.m, vt.data, a.n, superb)`
	`93`	+ vlas.dgesvd(`A`, `A`, a.m, a.n, acpy.data, 1, s, u.data, a.m, vt.data, a.n, superb)
`94`	`94`	`}`
`95`	`95`
`96`	`96`	`// matrix_inv computes the inverse of a general matrix (square or not). It also computes the`