differentialforms.hh

Go to the documentation of this file.

 
#ifndef differentialforms_hh
#define differentialforms_hh
 
#include "toolbox/array.hh"
#include "integration/quadRule.hh"
#include "geometry/cell3D.hh"
 
namespace concepts
{
 
  class multilevelindex
  {
  public:
    multilevelindex() {};
    ~multilevelindex() {};
    multilevelindex(const Array<int> & Nsize) : L_(Nsize.size()), Nsize_(Nsize)
    {
      N_ = 1;
      for(int l=0; l<L_; ++l) N_ *= Nsize_[l];
    };
    int getL() {return L_; };
    int getN() {return N_; };
    Array<int> getNsize() {return Nsize_; };
    Array<int> operator[] (const int n)
    {
      Array<int> k(L_);
      int idxk=n;
      for(int l=0; l<L_; ++l)
        {
          k[l] = idxk % Nsize_[l];
          idxk /= Nsize_[l];
        }
      return k;
    };
 
  protected:
    int N_;
    int L_;
    Array<int> Nsize_;
  };
 
  Real quadratureWeightTensor(const Array<QuadratureRule1d*> & ArrayQuad1D, const Array<int> & Index);
 
  Array<Real> quadratureWeightTensor(const Array<QuadratureRule1d*> & ArrayQuad1D);
 
  Mapping<Real, 3, 3> jacobianTensor(const Hexahedron3d & Hexa, const Real3d & xi);
 
  Array<Mapping<Real, 3, 3> > jacobianTensor(const Hexahedron3d & Hexa, const Array<QuadratureRule1d*> & ArrayQuad1D);
 
  Real jacobianDeterminant(const Hexahedron3d & Hexa, const Real3d & xi);
 
  Array<Real> jacobianDeterminant(const Hexahedron3d & Hexa, const Array<QuadratureRule1d*> & ArrayQuad1D);
 
}
 
#endif // differentialforms_hh