jacobiantocurl.hh

// -*- tab-width: 8; indent-tabs-mode: nil; c-basic-offset: 2 -*-
// vi: set et ts=8 sw=2 sts=2:
#ifndef DUNE_PDELAB_COMMON_JACOBIANTOCURL_HH
#define DUNE_PDELAB_COMMON_JACOBIANTOCURL_HH
 
#include <cstddef>
 
#include <dune/common/fvector.hh>
 
namespace Dune {
  namespace PDELab {
 
 
    template<typename Jacobian, std::size_t dimR = Jacobian::rows,
             std::size_t dimD = Jacobian::cols>
    class JacobianToCurl;
 
 
    template<typename Jacobian>
    class JacobianToCurl<Jacobian, 1, 2> {
      static_assert
      ( Jacobian::rows == 1 && Jacobian::cols == 2, "This specialization "
        "works only for dimRange == 1 and dimDomain == 2");
 
    public:
      typedef typename Jacobian::block_type CurlField;
      static const std::size_t dimCurl = 2;
      typedef FieldVector<CurlField, dimCurl> Curl;
 
      void operator()(const Jacobian& jacobian, Curl& curl) const {
        curl[0] =  jacobian[0][1];
        curl[1] = -jacobian[0][0];
      }
    };
 
 
    template<typename Jacobian>
    class JacobianToCurl<Jacobian, 2, 2> {
      static_assert
      ( Jacobian::rows == 2 && Jacobian::cols == 2, "This specialization "
        "works only for dimRange == 2 and dimDomain == 2");
 
    public:
      typedef typename Jacobian::block_type CurlField;
      static const std::size_t dimCurl = 1;
      typedef FieldVector<CurlField, dimCurl> Curl;
 
      void operator()(const Jacobian& jacobian, Curl& curl) const {
        curl[0] = jacobian[1][0]-jacobian[0][1];
      }
    };
 
 
    template<typename Jacobian>
    class JacobianToCurl<Jacobian, 3, 3> {
      static_assert
      ( Jacobian::rows == 3 && Jacobian::cols == 3, "This specialization "
        "works only for dimRange == 3 and dimDomain == 3");
 
    public:
      typedef typename Jacobian::block_type CurlField;
      static const std::size_t dimCurl = 3;
      typedef FieldVector<CurlField, dimCurl> Curl;
 
      void operator()(const Jacobian& jacobian, Curl& curl) const {
        for(std::size_t alpha = 0; alpha < 3; ++alpha) {
          std::size_t beta  = (alpha+1)%3;
          std::size_t gamma = (alpha+2)%3;
          curl[alpha] = jacobian[gamma][beta]-jacobian[beta][gamma];
        }
      }
    };
 
  } // namespace PDELab
} //namespace Dune
 
#endif // DUNE_PDELAB_COMMON_JACOBIANTOCURL_HH