Bug fix

termoshtt · termoshtt · commit c5db4b8a424b · 2019-05-05T18:28:45.000+09:00
diff --git a/src/krylov/householder.rs b/src/krylov/householder.rs
@@ -21,7 +21,7 @@ impl<A: Scalar> Householder<A> {
         let w = self.v[k].slice(s![k..]);
         let c = A::from(2.0).unwrap() * w.inner(&a.slice(s![k..]));
         for l in k..self.dim {
-            a[l] -= c * w[l];
+            a[l] -= c * w[l - k];
         }
     }
 }

Original file line number	Diff line number	Diff line change
`@@ -21,7 +21,7 @@ impl<A: Scalar> Householder<A> {`
`21`	`21`	`let w = self.v[k].slice(s![k..]);`
`22`	`22`	`let c = A::from(2.0).unwrap() * w.inner(&a.slice(s![k..]));`
`23`	`23`	`for l in k..self.dim {`
`24`		`- a[l] -= c * w[l];`
	`24`	`+ a[l] -= c * w[l - k];`
`25`	`25`	`}`
`26`	`26`	`}`
`27`	`27`	`}`