LinearSolver::getError returns -1 for direct solver

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;

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;

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;

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;