LinearObsOperatorBase.h

Go to the documentation of this file.

/*
 * (C) Copyright 2017-2021 UCAR
 *
 * This software is licensed under the terms of the Apache Licence Version 2.0
 * which can be obtained at http://www.apache.org/licenses/LICENSE-2.0.
 */
 
#ifndef UFO_LINEAROBSOPERATORBASE_H_
#define UFO_LINEAROBSOPERATORBASE_H_
 
#include <map>
#include <memory>
#include <string>
#include <vector>
 
#include <boost/noncopyable.hpp>
 
#include "eckit/config/Configuration.h"
#include "ioda/ObsSpace.h"
#include "oops/base/Variables.h"
#include "oops/util/AssociativeContainers.h"
#include "oops/util/parameters/ConfigurationParameter.h"
#include "oops/util/parameters/HasParameters_.h"
#include "oops/util/parameters/Parameters.h"
#include "oops/util/parameters/ParametersOrConfiguration.h"
#include "oops/util/parameters/RequiredPolymorphicParameter.h"
#include "oops/util/Printable.h"
#include "ufo/ObsOperatorParametersBase.h"
 
namespace ioda {
class ObsVector;
}
 
namespace ufo {
class GeoVaLs;
class ObsDiagnostics;
 
// -----------------------------------------------------------------------------
/// Base class for linear observation operators
///
/// Note: subclasses can opt to extract their settings either from
/// a Configuration object or from a subclass of ObsOperatorParametersBase.
///
/// In the former case, they should provide a constructor with the following signature:
///
///    SubclassName(const ioda::ObsSpace &, const eckit::Configuration &);
///
/// In the latter case, the implementer should first define a subclass of ObsOperatorParametersBase
/// holding the settings of the operator in question. The LinearObsOperatorBase subclass should
/// then typedef `Parameters_` to the name of the ObsOperatorParametersBase subclass and provide a
/// constructor with the following signature:
///
///    SubclassName(const ioda::ObsSpace &, const Parameters_ &);
class LinearObsOperatorBase : public util::Printable,
                              private boost::noncopyable {
 public:
  explicit LinearObsOperatorBase(const ioda::ObsSpace & odb)
    : odb_(odb) {}
  virtual ~LinearObsOperatorBase() {}
 
/// Obs Operator
  virtual void setTrajectory(const GeoVaLs &, ObsDiagnostics &) = 0;
  virtual void simulateObsTL(const GeoVaLs &, ioda::ObsVector &) const = 0;
  virtual void simulateObsAD(GeoVaLs &, const ioda::ObsVector &) const = 0;
 
/// Operator input required from Model
  virtual const oops::Variables & requiredVars() const = 0;
 
/// \brief List of variables simulated by this operator.
///
/// The default implementation returns the list of all simulated variables in the ObsSpace.
  virtual oops::Variables simulatedVars() const;
 
/// \brief The space containing the observations to be simulated by this operator.
  const ioda::ObsSpace &obsspace() const { return odb_; }
 
 private:
  virtual void print(std::ostream &) const = 0;
 
 private:
  const ioda::ObsSpace & odb_;
};
 
// -----------------------------------------------------------------------------
 
class LinearObsOperatorFactory;
 
// -----------------------------------------------------------------------------
 
/// \brief Contains a polymorphic parameter holding an instance of a subclass of
/// ObsOperatorParametersBase.
class LinearObsOperatorParametersWrapper : public oops::Parameters {
  OOPS_CONCRETE_PARAMETERS(LinearObsOperatorParametersWrapper, Parameters)
 public:
  /// After deserialization, holds an instance of a subclass of ObsOperatorParametersBase
  /// controlling the behavior of a linear observation operator. The type of the subclass is
  /// determined by the value of the "name" key in the Configuration object from which this object
  /// is deserialized.
  oops::RequiredPolymorphicParameter<ObsOperatorParametersBase, LinearObsOperatorFactory>
    operatorParameters{"name", this};
};
 
// -----------------------------------------------------------------------------
 
/// Linear obs operator factory
class LinearObsOperatorFactory {
 public:
  /// \brief Create and return a new linear observation operator.
  ///
  /// The type of the operator is determined by the `name` attribute of \p params. \p params must
  /// be an instance of the subclass of ObsOperatorParametersBase associated with that operator,
  /// otherwise an exception will be thrown.
  static LinearObsOperatorBase * create(const ioda::ObsSpace &, const ObsOperatorParametersBase &);
 
  /// \brief Create and return an instance of the subclass of ObsOperatorParametersBase
  /// storing parameters of linear observation operators of the specified type.
  static std::unique_ptr<ObsOperatorParametersBase> createParameters(
      const std::string &name);
 
  /// \brief Return the names of all operators that can be created by one of the registered makers.
  static std::vector<std::string> getMakerNames() {
    return oops::keys(getMakers());
  }
 
  virtual ~LinearObsOperatorFactory() = default;
 
 protected:
  /// \brief Register a maker able to create linear observation operators of type \p name.
  explicit LinearObsOperatorFactory(const std::string &name);
 
 private:
  virtual LinearObsOperatorBase * make(const ioda::ObsSpace &,
                                       const ObsOperatorParametersBase &) = 0;
 
  virtual std::unique_ptr<ObsOperatorParametersBase> makeParameters() const = 0;
 
  static std::map < std::string, LinearObsOperatorFactory * > & getMakers() {
    static std::map < std::string, LinearObsOperatorFactory * > makers_;
    return makers_;
  }
};
 
// -----------------------------------------------------------------------------
 
template<class T>
class LinearObsOperatorMaker : public LinearObsOperatorFactory {
  /// Defined as T::Parameters_ if T defines a Parameters_ type; otherwise as
  /// GenericObsOperatorParameters.
  typedef oops::TParameters_IfAvailableElseFallbackType_t<T, GenericObsOperatorParameters>
    Parameters_;
 
  LinearObsOperatorBase * make(const ioda::ObsSpace & odb,
                               const ObsOperatorParametersBase & params) override {
    const auto &stronglyTypedParams = dynamic_cast<const Parameters_&>(params);
    return new T(odb,
                 oops::parametersOrConfiguration<oops::HasParameters_<T>::value>(
                   stronglyTypedParams));
  }
 
  std::unique_ptr<ObsOperatorParametersBase> makeParameters() const override {
    return std::make_unique<Parameters_>();
  }
 
 public:
  explicit LinearObsOperatorMaker(const std::string & name) : LinearObsOperatorFactory(name) {}
};
 
// -----------------------------------------------------------------------------
 
}  // namespace ufo
 
#endif  // UFO_LINEAROBSOPERATORBASE_H_