recon.h

Go to the documentation of this file.

#ifndef DVRLIB_RECON_H
#define DVRLIB_RECON_H
 
namespace dvrlib {
 
class vector;
class matrix;
 
template <class argtype, class restype>
class func {
public:
  virtual ~func() = default;
  virtual restype operator()(const argtype& arg) const = 0;
};
 
void lin_cov_update(const matrix& S_x,
                    const matrix& F,
                    matrix& S_v);
 
void lin_recon(const vector& r,
               const matrix& S_x,
               const matrix& F,
               vector& v);
 
void lin_recon_update(const vector& r,
                      const matrix& S_x_inv,
                      const matrix& F,
                      const vector& v,
                      vector& dv);
 
int recon(const vector& x,
          const matrix& S_x,
          const func<vector, vector>& f,
          const func<vector, matrix>& J,
          vector& v,
          matrix& S_v,
          double eps = 1e-6,
          int maxiter = 50);
 
void extract_confidence(const matrix& S_xnew,
                        vector& conf_results);
 
double confint2var(double confint);
 
double var2confint(double var);
 
}  // namespace dvrlib
 
#endif  // DVRLIB_RECON_H