diff --git a/tbox/src/wDss.h b/tbox/src/wDss.h
index da42e4ac37db46fd320c6bef2e4a8814db7278d2..f5a27ae2e7d86307356654500b09b4d053b748c9 100644
--- a/tbox/src/wDss.h
+++ b/tbox/src/wDss.h
@@ -23,7 +23,6 @@
 #include "wLinearSolver.h"
 #include <mkl_dss.h>
 #include <Eigen/Sparse>
-#include <limits>
 
 namespace tbox
 {
@@ -48,7 +47,7 @@ public:
     virtual void solve(Eigen::Map<Eigen::VectorXd> const &b, Eigen::Map<Eigen::VectorXd> &x) override;
     virtual void compute(Eigen::SparseMatrix<double> const &A, Eigen::Map<Eigen::VectorXd> const &b, Eigen::Map<Eigen::VectorXd> &x) override;
 
-    virtual double getError() override { return std::numeric_limits<double>::epsilon(); }
+    virtual double getError() override { return -1.; }
     virtual int getIterations() override { return 1; }
 
     virtual void write(std::ostream &out) const override;
diff --git a/tbox/src/wMumps.h b/tbox/src/wMumps.h
index dae62bd54973e6bec0fa4b1f7ebaf7172b308377..b3f14923750f05ddc4a4a68ae5dc7eb59c2b3df9 100644
--- a/tbox/src/wMumps.h
+++ b/tbox/src/wMumps.h
@@ -23,7 +23,6 @@
 #include "wLinearSolver.h"
 #include <dmumps_c.h>
 #include <Eigen/Sparse>
-#include <limits>
 
 namespace tbox
 {
@@ -61,7 +60,7 @@ public:
     virtual void solve(Eigen::Map<Eigen::VectorXd> const &b, Eigen::Map<Eigen::VectorXd> &x) override;
     virtual void compute(Eigen::SparseMatrix<double> const &A, Eigen::Map<Eigen::VectorXd> const &b, Eigen::Map<Eigen::VectorXd> &x) override;
 
-    virtual double getError() override { return std::numeric_limits<double>::epsilon(); }
+    virtual double getError() override { return -1.; }
     virtual int getIterations() override { return 1; }
 
     virtual void write(std::ostream &out) const override;
diff --git a/tbox/src/wPardiso.h b/tbox/src/wPardiso.h
index 3e6c57657b4cded842547860c4c336036443ba66..8cd7be4c74660660a59339e2e78cd82a2255b108 100644
--- a/tbox/src/wPardiso.h
+++ b/tbox/src/wPardiso.h
@@ -22,7 +22,6 @@
 
 #include "wLinearSolver.h"
 #include <Eigen/PardisoSupport>
-#include <limits>
 
 namespace tbox
 {
@@ -46,7 +45,7 @@ public:
     virtual void solve(Eigen::Map<Eigen::VectorXd> const &b, Eigen::Map<Eigen::VectorXd> &x) override;
     virtual void compute(Eigen::SparseMatrix<double> const &A, Eigen::Map<Eigen::VectorXd> const &b, Eigen::Map<Eigen::VectorXd> &x) override;
 
-    virtual double getError() override { return std::numeric_limits<double>::epsilon(); }
+    virtual double getError() override { return -1.; }
     virtual int getIterations() override { return 1; }
 
     virtual void write(std::ostream &out) const override;
diff --git a/tbox/src/wSparseLu.h b/tbox/src/wSparseLu.h
index 6055f3815bd616844bafa843ec794c42e62bf45e..402dffe920ff8c1945f91047e79835f59df9ca98 100644
--- a/tbox/src/wSparseLu.h
+++ b/tbox/src/wSparseLu.h
@@ -20,7 +20,6 @@
 #include "tbox.h"
 #include "wLinearSolver.h"
 #include <Eigen/Sparse>
-#include <limits>
 
 namespace tbox
 {
@@ -46,7 +45,7 @@ public:
     virtual void solve(Eigen::Map<Eigen::VectorXd> const &b, Eigen::Map<Eigen::VectorXd> &x) override;
     virtual void compute(Eigen::SparseMatrix<double> const &A, Eigen::Map<Eigen::VectorXd> const &b, Eigen::Map<Eigen::VectorXd> &x) override;
 
-    virtual double getError() override { return std::numeric_limits<double>::epsilon(); }
+    virtual double getError() override { return -1.; }
     virtual int getIterations() override { return 1; }
 
     virtual void write(std::ostream &out) const override;