PiolaOperators.h

Go to the documentation of this file.

//==============================================================================
//==============================================================================
 
#ifndef PIOLA_OPERATORS_H_
#define PIOLA_OPERATORS_H_
 
#include "MatVec.h"
#include "EqualOrderOperators.h"
 
#include <array>
 
class FiniteElement;
class Tensor;
class Vec3;
 
 
class PiolaOperators
{
public:
  class Weak
  {
  public:
    static void Advection(Matrices& 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(Matrices& 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(Matrices& EM,
                         const FiniteElement& fe,
                         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(Matrices& EM,
                          const FiniteElement& fe,
                          const std::array<std::array<int,3>,3>& idx,
                          double scale, bool stress);
 
    static void Mass(Matrices& EM,
                     const FiniteElement& fe,
                     const std::array<std::array<int,3>,3>& idx,
                     double scale);
 
    static void Source(Vectors& EV, const FiniteElement& fe,
                       const Vec3& f, const std::array<int,3>& idx,
                       double scale);
  };
 
  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);
  };
 
  static void Copy(Matrices& EM,
                   const FiniteElement& fe,
                   const std::array<std::array<int,3>,3>& idx,
                   const Matrix& A);
 
  static void Copy(Vectors& EV,
                   const FiniteElement& fe,
                   const std::array<int,3>& idx,
                   const RealArray& V);
};
 
#endif