hp2D/shortestDist.hh

Go to the documentation of this file.

/* weight functions for weighted regularization
 */

#ifndef shortestDist2D_hh
#define shortestDist2D_hh

#include <cmath>
#include "basics/vectorsMatrices.hh"
#include "hp2D/singularSet.hh"

namespace hp2D {

  // forward declarations
  template<class F>
  class Quad;
  using concepts::Real;

  // ********************************************************* TrivialWeight **

  class TrivialWeight {
  public:
    Real operator()(const SingularSet&, concepts::Real2d,
        const Quad<Real>*) const { return 1.0; }
  };

  std::ostream& operator<<(std::ostream& os, const TrivialWeight& p);

  // ********************************************************** ShortestDist **

  class ShortestDist {
  public:
    Real operator()(const SingularSet& singularities,
        concepts::Real2d point, const Quad<Real>* elm) const;
  };

  std::ostream& operator<<(std::ostream& os, const ShortestDist& p);

  // ********************************************************** DistancePost **

  template<typename DistClass, typename Function>
  class DistancePost {
  public:
    Real operator()(const SingularSet& singularities,
        concepts::Real2d point, const Quad<Real>* elm) const {
      return Function::post(dist_(singularities, point, elm));
    }
  private:
    DistClass dist_;
  };

  template<typename DistClass, typename Function>
  inline std::ostream& operator<<(std::ostream& os,
          const DistancePost<DistClass, Function>& p) {
    return os << Function() << '(' << DistClass() << ')';
  }


  // ********************************************************** Postprocess4 **

  class Postprocess4 {
  public:
    static Real post(Real x) { return std::pow(x, 0.4); }
  };

  std::ostream& operator<<(std::ostream& os, const Postprocess4& p);

  // ********************************************************** Postprocess7 **

  class Postprocess7 {
  public:
    static Real post(Real x) { return std::pow(x, 0.7); }
  };

  std::ostream& operator<<(std::ostream& os, const Postprocess7& p);

  // ********************************************************** Postprocess8 **

  class Postprocess8 {
  public:
    static Real post(Real x) { return std::pow(x, 0.8); }
  };

  std::ostream& operator<<(std::ostream& os, const Postprocess8& p);

  // ********************************************************** Postprocess9 **

  class Postprocess9 {
  public:
    static Real post(Real x) { return std::pow(x, 0.9); }
  };

  std::ostream& operator<<(std::ostream& os, const Postprocess9& p);

  // ******************************************************* PostprocessSqrt **

  class PostprocessSqrt {
  public:
    static Real post(Real x) { return std::sqrt(x); }
  };

  std::ostream& operator<<(std::ostream& os, const PostprocessSqrt& p);

} // namespace hp2D

#endif // shortestDist2D_hh