test/interface/GeoVaLs.h

Go to the documentation of this file.

/*
 * (C) Copyright 2017-2018 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 TEST_INTERFACE_GEOVALS_H_
#define TEST_INTERFACE_GEOVALS_H_
 
#include <cmath>
#include <memory>
#include <string>
#include <vector>
 
#define ECKIT_TESTING_SELF_REGISTER_CASES 0
 
#include <boost/noncopyable.hpp>
 
#include "eckit/config/LocalConfiguration.h"
#include "eckit/testing/Test.h"
#include "oops/base/ObsSpaces.h"
#include "oops/base/State.h"
#include "oops/base/Variables.h"
#include "oops/interface/GeoVaLs.h"
#include "oops/mpi/mpi.h"
#include "oops/runs/Test.h"
#include "oops/util/dot_product.h"
#include "oops/util/Logger.h"
#include "test/interface/ObsTestsFixture.h"
#include "test/TestEnvironment.h"
 
namespace test {
 
// -----------------------------------------------------------------------------
/// \brief Tests test-constructor and print method
template <typename OBS> void testConstructor() {
  typedef ObsTestsFixture<OBS>  Test_;
  typedef oops::GeoVaLs<OBS>    GeoVaLs_;
 
  for (std::size_t jj = 0; jj < Test_::obspace().size(); ++jj) {
    eckit::LocalConfiguration gconf(Test_::config(jj), "geovals");
    oops::Variables geovars(gconf, "state variables");
    std::unique_ptr<GeoVaLs_> geovals(new GeoVaLs_(gconf, Test_::obspace()[jj], geovars));
    EXPECT(geovals.get());
    oops::Log::test() << "Testing GeoVaLs: " << *geovals << std::endl;
    geovals.reset();
    EXPECT(!geovals.get());
  }
}
 
// -----------------------------------------------------------------------------
 
template <typename OBS> void testUtils() {
  typedef ObsTestsFixture<OBS>  Test_;
  typedef oops::GeoVaLs<OBS>    GeoVaLs_;
 
  const double tol = 1e-6;
 
  for (std::size_t jj = 0; jj < Test_::obspace().size(); ++jj) {
    eckit::LocalConfiguration gconf(Test_::config(jj), "geovals");
    oops::Variables geovars(gconf, "state variables");
    GeoVaLs_ gval(gconf, Test_::obspace()[jj], geovars);
 
    const double zz = dot_product(gval, gval);
 
    oops::Log::trace() << "Testing copy constructor (=) " << std::endl;
 
    GeoVaLs_ gval2 = gval;
 
    double zz0 = dot_product(gval2, gval2);
 
    EXPECT(zz0 == zz);
 
    oops::Log::trace() << "Testing *= double" << std::endl;
 
    gval2 *= 2.0;
 
    const double zz1 = dot_product(gval2, gval2);
 
    EXPECT(zz1/zz - 4.0 < tol);
 
    oops::Log::trace() << "Testing += GeoVals" << std::endl;
 
    gval2 += gval;
 
    const double zz2 = dot_product(gval2, gval2);
 
    EXPECT(zz2/zz - 9.0  < tol);
 
    oops::Log::trace() << "Testing -= GeoVals" << std::endl;
 
    gval2 -= gval;
 
    const double zz3 = dot_product(gval2, gval2);
 
    EXPECT(zz3/zz - 4.0  < tol);
 
    oops::Log::trace() << "Testing *= GeoVals" << std::endl;
 
    gval2 = gval;
 
    gval2 *= gval;
 
    GeoVaLs_ gval3 = gval2;
 
    gval3 *= gval;
 
    const double zz4 = dot_product(gval2, gval2);
 
    const double zz5 = dot_product(gval3, gval);
 
    EXPECT(zz4/zz5 - 1.0  < tol);
 
    oops::Log::trace() << "Testing random" << std::endl;
 
    gval.random();
    const double zz6 = dot_product(gval, gval);
    EXPECT(zz6 > 0.0);
 
    oops::Log::trace() << "Testing zero" << std::endl;
 
    gval.zero();
    const double zz7 = dot_product(gval, gval);
    EXPECT(zz7 == 0.0);
  }
}
 
// -----------------------------------------------------------------------------
 
template <typename OBS> void testRead() {
  typedef ObsTestsFixture<OBS> Test_;
  typedef oops::GeoVaLs<OBS>   GeoVaLs_;
 
  const double tol = 1.0e-9;
  for (std::size_t jj = 0; jj < Test_::obspace().size(); ++jj) {
    eckit::LocalConfiguration gconf(Test_::config(jj), "geovals");
    oops::Variables geovars(gconf, "state variables");
    GeoVaLs_ gval(gconf, Test_::obspace()[jj], geovars);
 
    const double xx = gconf.getDouble("norm");
    const double zz = sqrt(dot_product(gval, gval));
 
    oops::Log::debug() << "xx: " << std::fixed << std::setprecision(8) << xx << std::endl;
    oops::Log::debug() << "zz: " << std::fixed << std::setprecision(8) << zz << std::endl;
 
    EXPECT(oops::is_close(xx, zz, tol));
  }
}
 
// -----------------------------------------------------------------------------
 
template <typename OBS>
class GeoVaLs : public oops::Test {
  typedef ObsTestsFixture<OBS> Test_;
 public:
  GeoVaLs() {}
  virtual ~GeoVaLs() {}
 private:
  std::string testid() const override {return "test::GeoVaLs<" + OBS::name() + ">";}
 
  void register_tests() const override {
    std::vector<eckit::testing::Test>& ts = eckit::testing::specification();
 
    ts.emplace_back(CASE("interface/GeoVaLs/testConstructor")
      { testConstructor<OBS>(); });
    ts.emplace_back(CASE("interface/GeoVaLs/testUtils")
      { testUtils<OBS>(); });
    ts.emplace_back(CASE("interface/GeoVaLs/testRead")
      { testRead<OBS>(); });
  }
 
  void clear() const override {
    Test_::reset();
  }
};
 
// =============================================================================
 
}  // namespace test
 
#endif  // TEST_INTERFACE_GEOVALS_H_