numerics/corner_singularity_biharmonic.h

// -*- c++ -*-

// +--------------------------------------------------------------------+
// | This file is part of MathTL - the Mathematical Template Library    |
// |                                                                    |
// | Copyright (c) 2002-2009                                            |
// | Thorsten Raasch, Manuel Werner                                     |
// +--------------------------------------------------------------------+

#ifndef _MATHTL_CORNER_SINGULARITY_BIHARMONIC_H
#define _MATHTL_CORNER_SINGULARITY_BIHARMONIC_H

#include <geometry/point.h>
#include <utils/function.h>

namespace MathTL
{
  class CornerSingularityBiharmonic
    : public Function<2>
  {
  public:
    CornerSingularityBiharmonic(const Point<2>& x0,
                      const double theta0,
                      const double omega,
                      const double r0 = 0.01,
                      const double r1 = 0.99);

    virtual ~CornerSingularityBiharmonic() {}

    double value(const Point<2>& x, const unsigned int component = 0) const;

    void vector_value(const Point<2>& p, Vector<double>& values) const;
    
  protected:
    Point<2> x0;

    double theta0;

    double omega;

    double r0, r1;

    double zeta(const double r) const;
  };

  class CornerSingularityBiharmonicRHS
    : public Function<2>
  {
  public:
    CornerSingularityBiharmonicRHS(const Point<2>& x0,
                         const double theta0,
                         const double omega,
                         const double r0 = 0.01,
                         const double r1 = 0.99);
    
    virtual ~CornerSingularityBiharmonicRHS() {}

    double value(const Point<2>& x, const unsigned int component = 0) const;

    void vector_value(const Point<2>& p, Vector<double>& values) const;
    
  protected:
    Point<2> x0;

    double theta0;

    double omega;

    double r0, r1;

    double zeta(const double r) const;

    double zeta_prime(const double r) const;

    double zeta_primeprime(const double r) const;
 
    double zeta_third_der(const double r) const;

    double zeta_fourth_der(const double r) const;
 };

}

// implementations of inline functions
#include "numerics/corner_singularity_biharmonic.cpp"

#endif