Fixes #41 and corrects bug from #40

ChristopherRabotin · Apr 7, 2017 · 89ebb62 · 89ebb62
1 parent c25e79c
commit 89ebb62
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 3 deletions.
diff --git a/c2d.go b/c2d.go
@@ -14,7 +14,7 @@ func VanLoan(A, Γ, W *mat64.Dense, Δt float64) (*mat64.Dense, *mat64.SymDense,
 	var err error
 	// Check aliasing
 	λ := mat64.Eigen{}
-	λ.Factorize(A, false)
+	λ.Factorize(A, false, false)
 	λmaxImag := -math.MaxFloat64
 	var λmax complex128
 	for _, λ := range λ.Values(nil) {

diff --git a/truth.go b/truth.go
@@ -21,7 +21,7 @@ func (t *BatchGroundTruth) Error(k int, est Estimate) Estimate {
 func (t *BatchGroundTruth) ErrorWithOffset(k int, est Estimate, offset *mat64.Vector) Estimate {
 	esR, _ := est.State().Dims()
 	estState := mat64.NewVector(esR, nil)
-	if k > 0 {
+	if k >= 0 {
 		estState.CopyVec(est.State())
 		if offset != nil {
 			estState.AddVec(estState, offset)
@@ -39,7 +39,7 @@ func (t *BatchGroundTruth) ErrorWithOffset(k int, est Estimate, offset *mat64.Ve
 
 	emR, _ := est.Measurement().Dims()
 	estMeas := mat64.NewVector(emR, nil)
-	if k > 0 {
+	if k >= 0 {
 		estMeas.CopyVec(est.Measurement())
 		trueMeas := mat64.NewVector(esR, nil)
 		if t.states != nil {