EqualOrderOperators.h

Go to the documentation of this file.

//==============================================================================
//==============================================================================
 
#ifndef EQUAL_ORDER_OPERATORS_H
#define EQUAL_ORDER_OPERATORS_H
 
class Vec3;
 
#include "FiniteElement.h"
#include "MatVec.h"
 
 
namespace WeakOperators
{
  enum ConvectionForm
  {
    CONVECTIVE = 0,   
    CONSERVATIVE = 1, 
    SKEWSYMMETRIC = 2 
  };
}
 
 
class EqualOrderOperators
{
public:
  class Weak {
  public:
    static void Advection(Matrix& EM, const FiniteElement& fe,
                          const Vec3& AC, double scale=1.0,
                          WeakOperators::ConvectionForm form = WeakOperators::CONVECTIVE,
                          int basis=1);
 
    static void Convection(Matrix& EM, const FiniteElement& fe,
                           const Vec3& U, const Tensor& dUdX, double scale,
                           WeakOperators::ConvectionForm form=WeakOperators::CONVECTIVE,
                           int basis=1);
 
    static void Divergence(Matrix& EM, const FiniteElement& fe,
                           double scale=1.0, int basis=1, int tbasis=1);
 
    static void Divergence(Vector& EV, const FiniteElement& fe,
                           const Vec3& D, double scale=1.0, int basis=1);
 
    static void Gradient(Matrix& EM, const FiniteElement& fe,
                         double scale=1.0, int basis=1, int tbasis=1);
 
    static void Laplacian(Matrix& EM, const FiniteElement& fe,
                          double scale=1.0, bool stress=false, int basis=1);
 
    static void Stress(Matrix& EM, const FiniteElement& fe,
                       double scale=1.0, int basis=1);
 
    static void LaplacianCoeff(Matrix& EM, const Matrix& K, const FiniteElement& fe,
                               double scale=1.0, int basis=1);
 
    static void ItgConstraint(Matrix& EM, const FiniteElement& fe,
                              double scale = 1.0, int basis = 1);
 
    static void Mass(Matrix& EM, const FiniteElement& fe,
                     double scale=1.0, int basis=1);
 
    static void Source(Vector& EV, const FiniteElement& fe,
                       double scale=1.0, int cmp=1, int basis=1);
 
    static void Source(Vector& EV, const FiniteElement& fe,
                       const Vec3& f, double scale=1.0, int basis=1);
  };
 
  class Residual {
  public:
    static void Advection(Vector& EV, const FiniteElement& fe,
                          const Vec3& AC, const Tensor& g,
                          double scale = 1.0, int basis=1);
 
    static void Convection(Vector& EV, const FiniteElement& fe,
                           const Vec3& U, const Tensor& dUdX, const Vec3& UC,
                           double scale,
                           WeakOperators::ConvectionForm form=WeakOperators::CONVECTIVE,
                           int basis=1);
 
    static void Divergence(Vector& EV, const FiniteElement& fe,
                           const Tensor& dUdX, double scale=1.0,
                           size_t basis=1);
 
    static void Gradient(Vector& EV, const FiniteElement& fe,
                         double scale=1.0, int basis=1);
 
    static void Laplacian(Vector& EV, const FiniteElement& fe,
                          const Vec3& dUdX, double scale = 1.0,
                          int basis = 1);
 
    static void Laplacian(Vector& EV, const FiniteElement& fe,
                          const Tensor& dUdX, double scale = 1.0,
                          bool stress = false, int basis = 1);
  };
};
 
#endif