CompatibleOperators.h

Go to the documentation of this file.

//==============================================================================
//==============================================================================
 
#ifndef COMPATIBLE_OPERATORS_H_
#define COMPATIBLE_OPERATORS_H_
 
class FiniteElement;
class Tensor;
class Vec3;
 
#include "MatVec.h"
#include "EqualOrderOperators.h"
 
#include <array>
 
class CompatibleOperators
{
public:
  class Weak {
  public:
    static void Advection(std::vector<Matrix>& EM, const FiniteElement& fe,
                          const Vec3& AC,
                          const std::array<std::array<int,3>,3>& idx,
                          double scale = 1.0,
                          WeakOperators::ConvectionForm cnvForm = WeakOperators::CONVECTIVE);
 
    static void Convection(std::vector<Matrix>& EM, const FiniteElement& fe,
                           const Vec3& U, const Tensor& dUdX,
                           const std::array<std::array<int,3>,3>& idx,
                           double scale,
                           WeakOperators::ConvectionForm form = WeakOperators::CONVECTIVE);
 
    static void Gradient(std::vector<Matrix>& EM,
                         const FiniteElement& fe,
                         const std::array<int,3>& idx,
                         double scale = 1.0);
 
    static void ItgConstraint(std::vector<Matrix>& EM,
                              const FiniteElement& fe,
                              const std::array<int,3>& idx);
 
    static void Laplacian(std::vector<Matrix>& EM, const FiniteElement& fe,
                          const std::array<std::array<int,3>,3>& idx,
                          double scale = 1.0, bool stress = false);
 
    static void Mass(std::vector<Matrix>& EM,
                     const FiniteElement& fe,
                     const std::array<std::array<int,3>,3>& idx,
                     double scale = 1.0);
 
    static void Source(Vectors& EV, const FiniteElement& fe,
                       const std::array<int,3>& idx, double scale = 1.0);
 
    static void Source(Vectors& EV, const FiniteElement& fe,
                       const Vec3& f, const std::array<int,3>& idx,
                       double scale = 1.0);
  };
 
  class Residual
  {
  public:
    static void Convection(Vectors& EV, const FiniteElement& fe,
                           const Vec3& U, const Tensor& dUdX,
                           const Vec3& UC,
                           const std::array<int,3>& idx,
                           double scale,
                           WeakOperators::ConvectionForm form = WeakOperators::CONVECTIVE);
 
    static void Gradient(Vectors& EV, const FiniteElement& fe,
                         const std::array<int,3>& idx, double scale = 1.0);
 
    static void Laplacian(Vectors& EV, const FiniteElement& fe,
                          const Tensor& dUdX,
                          const std::array<int,3>& idx,
                          double scale = 1.0,
                          bool stress = false);
  };
};
 
#endif