TimeDomainLeastSquareDecon.h

#ifndef __TIME_DOMAIN_LEAST_SQUARE_DECON_H__
#define __TIME_DOMAIN_LEAST_SQUARE_DECON_H__
#include "mspass/algorithms/deconvolution/ScalarDecon.h"
#include "mspass/algorithms/deconvolution/ShapingWavelet.h"
#include "mspass/seismic/CoreTimeSeries.h"
#include "mspass/utility/Metadata.h"
#include <boost/archive/text_iarchive.hpp>
#include <boost/archive/text_oarchive.hpp>
#include <boost/serialization/base_object.hpp>
#include <vector>
 
namespace mspass::algorithms::deconvolution {
class TimeDomainLeastSquareDecon : public ScalarDecon {
public:
  TimeDomainLeastSquareDecon();
  TimeDomainLeastSquareDecon(const TimeDomainLeastSquareDecon &parent);
  TimeDomainLeastSquareDecon(const mspass::utility::Metadata &md);
  TimeDomainLeastSquareDecon(const mspass::utility::Metadata &md,
                             const std::vector<double> &wavelet,
                             const std::vector<double> &data);
  void changeparameter(const mspass::utility::Metadata &md);
  void process();
  mspass::seismic::CoreTimeSeries output_shaping_wavelet();
  mspass::seismic::CoreTimeSeries actual_output();
  mspass::seismic::CoreTimeSeries inverse_wavelet(const double t0parent);
  mspass::seismic::CoreTimeSeries inverse_wavelet();
  mspass::utility::Metadata QCMetrics();
 
private:
  void read_metadata(const mspass::utility::Metadata &md);
  int diagnostic_length() const;
  std::vector<double> solve_for(const std::vector<double> &rhs_data) const;
  std::vector<double> apply_wavelet(const std::vector<double> &model) const;
  std::vector<double> apply_shaping_wavelet(
      const std::vector<double> &model);
  int output_length;
  int sample_shift;
  double damp;
  double dt;
  mutable double residual_norm;
  mutable double data_norm;
  mutable double regularization_parameter;
 
  friend boost::serialization::access;
  template <class Archive>
  void serialize(Archive &ar, const unsigned int version) {
    ar &boost::serialization::base_object<ScalarDecon>(*this);
    ar & output_length;
    ar & sample_shift;
    ar & damp;
    ar & dt;
    ar & residual_norm;
    ar & data_norm;
    ar & regularization_parameter;
  }
};
} // namespace mspass::algorithms::deconvolution
#endif