ngsolve 6.2.2501.post21.dev1__cp312-cp312-macosx_10_15_universal2.whl → 6.2.2501.post37.dev1__cp312-cp312-macosx_10_15_universal2.whl

netgen/include/hdivfe_utils.hpp

@@ -1,716 +0,0 @@
-#ifndef HDIVFE_UTILS
-#define HDIVFE_UTILS
-
-
-#include "hcurlfe_utils.hpp"
-
-
-namespace ngfem
-{
-  
-
-  // rotated gradient
-  template <int DIM> class DuRot;
-
-  template <> class DuRot<2>
-  {
-
-  public:
-    const AutoDiff<2> u;
-
-    DuRot (const AutoDiff<2> au)
-      : u(au) { ; }
-
-    Vec<2> Value () const
-    {
-      Vec<2> val;
-      val(0) = u.DValue(1);
-      val(1) = -u.DValue(0);
-      return val;
-    }
-
-    /*
-    Vec<DIM_CURL> CurlValue () const
-    {
-      return Vec<DIM> (0.0);
-    }
-    */
-  };
-
-
-
-
-  /*
-  template <int DIM, typename SCAL>
-  class Class_uDvDw_Cyclic
-  {
-  public:
-    const AutoDiff<DIM,SCAL> u, v, w;
-    Class_uDvDw_Cyclic (const AutoDiff<DIM,SCAL> au, 
-                        const AutoDiff<DIM,SCAL> av,
-                        const AutoDiff<DIM,SCAL> aw)
-      : u(au), v(av), w(aw) { ; }
-  };
-
-  template <int DIM, typename SCAL>
-  INLINE Class_uDvDw_Cyclic<DIM,SCAL> 
-  uDvDw_Cyclic (AutoDiff<DIM,SCAL> u, AutoDiff<DIM,SCAL> v, AutoDiff<DIM,SCAL> w)
-  { return Class_uDvDw_Cyclic<DIM,SCAL> (u,v,w); }
-  */
-
-  template <int DIM, typename SCAL>
-  class uDvDw_Cyclic
-  {
-  public:
-    const AutoDiff<DIM,SCAL> u, v, w;
-    uDvDw_Cyclic (const AutoDiff<DIM,SCAL> au, 
-                  const AutoDiff<DIM,SCAL> av,
-                  const AutoDiff<DIM,SCAL> aw)
-      : u(au), v(av), w(aw) { ; }
-  };
-
-  
-  /*
-  template <int DIM, typename SCAL>
-  INLINE Class_uDvDw_Cyclic<DIM,SCAL> 
-  uDvDw_Cyclic (AutoDiff<DIM,SCAL> u, AutoDiff<DIM,SCAL> v, AutoDiff<DIM,SCAL> w)
-  { return Class_uDvDw_Cyclic<DIM,SCAL> (u,v,w); }
-  */
-
-
-
-  template <int DIM, typename SCAL>
-  class z_times_uDvDw_Cyclic
-  {
-  public:
-    const AutoDiff<DIM,SCAL> u, v, w, z;
-    z_times_uDvDw_Cyclic (const AutoDiff<DIM,SCAL> au, 
-                          const AutoDiff<DIM,SCAL> av,
-                          const AutoDiff<DIM,SCAL> aw,
-                          const AutoDiff<DIM,SCAL> az)
-      : u(au), v(av), w(aw), z(az) { ; }
-  };
-
-  /*
-  template <int DIM, typename SCAL>
-  INLINE auto
-  z_times_uDvDw_Cyclic (AutoDiff<DIM,SCAL> u, AutoDiff<DIM,SCAL> v,
-                        AutoDiff<DIM,SCAL> w, AutoDiff<DIM,SCAL> z)
-  { return Class_z_times_uDvDw_Cyclic<DIM,SCAL> (u,v,w,z); }
-  */
-
-  
-  /*
-  template <int DIM, typename SCAL>
-  INLINE auto
-  z_times_uDvDw_Cyclic (AutoDiff<DIM,SCAL> u, AutoDiff<DIM,SCAL> v,
-                AutoDiff<DIM,SCAL> w, AutoDiff<DIM,SCAL> z)
-  { return Class_z_times_uDvDw_Cyclic<DIM,SCAL> (u,v,w,z); }
-  */
-
-
-  
-  template <int DIM, typename SCAL>
-  class Du_Cross_Dv
-  {
-  public:
-    const AutoDiff<DIM,SCAL> u, v;
-    Du_Cross_Dv (const AutoDiff<DIM,SCAL> au, 
-                 const AutoDiff<DIM,SCAL> av)
-      : u(au), v(av) { ; }
-  };
-
-  /*
-  template <int DIM, typename SCAL>
-  INLINE Class_Du_Cross_Dv<DIM,SCAL> 
-  Du_Cross_Dv (AutoDiff<DIM,SCAL> u, AutoDiff<DIM,SCAL> v)
-  { return Class_Du_Cross_Dv<DIM,SCAL> (u,v); }
-  */
-
-  
-  /*
-  template <int DIM, typename SCAL>
-  INLINE Class_Du_Cross_Dv<DIM,SCAL>
-  Du_Cross_Dv (AutoDiff<DIM,SCAL> u, AutoDiff<DIM,SCAL> v)
-  { return Class_Du_Cross_Dv<DIM,SCAL> (u,v); }
-  */
-
-  template <int DIM, typename SCAL>
-  class wDu_Cross_Dv
-  {
-  public:
-    const AutoDiff<DIM,SCAL> u, v, w;
-    wDu_Cross_Dv (const AutoDiff<DIM,SCAL> au, 
-                  const AutoDiff<DIM,SCAL> av,
-                  const AutoDiff<DIM,SCAL> aw)
-      : u(au), v(av), w(aw) { ; }
-  };
-
-  /*
-  template <int DIM, typename SCAL>
-  INLINE Class_wDu_Cross_Dv<DIM,SCAL> 
-  wDu_Cross_Dv(AutoDiff<DIM,SCAL> u, AutoDiff<DIM,SCAL> v, AutoDiff<DIM,SCAL> w)
-  { return Class_wDu_Cross_Dv<DIM,SCAL> (u,v,w); }
-  */
-  
-  /*
-  template <int DIM, typename SCAL>
-  INLINE Class_wDu_Cross_Dv<DIM,SCAL>
-  wDu_Cross_Dv(AutoDiff<DIM,SCAL> u, AutoDiff<DIM,SCAL> v, AutoDiff<DIM,SCAL> w)
-  { return Class_wDu_Cross_Dv<DIM,SCAL> (u,v,w); }
-  */
-
-  
-  template <int DIM, typename SCAL>
-  class uDvDw_minus_DuvDw
-  {
-  public:
-    const AutoDiff<DIM,SCAL> u, v, w;
-    uDvDw_minus_DuvDw (const AutoDiff<DIM,SCAL> au, 
-                       const AutoDiff<DIM,SCAL> av,
-                       const AutoDiff<DIM,SCAL> aw)
-      : u(au), v(av), w(aw) { ; }
-  };
-
-  /*
-  template <int DIM, typename SCAL>
-  INLINE Class_uDvDw_minus_DuvDw<DIM,SCAL> 
-  uDvDw_minus_DuvDw (AutoDiff<DIM,SCAL> u, AutoDiff<DIM,SCAL> v, AutoDiff<DIM,SCAL> w)
-  { return Class_uDvDw_minus_DuvDw<DIM,SCAL> (u,v,w); }
-  */
-  
-  /*
-  template <int DIM, typename SCAL>
-  INLINE Class_uDvDw_minus_DuvDw<DIM,SCAL>
-  uDvDw_minus_DuvDw (AutoDiff<DIM,SCAL> u, AutoDiff<DIM,SCAL> v, AutoDiff<DIM,SCAL> w)
-  { return Class_uDvDw_minus_DuvDw<DIM,SCAL> (u,v,w); }
-  */
-
-  template <int DIM, typename SCAL>
-  class curl_uDvw_minus_Duvw
-  {
-  public:
-    const AutoDiff<DIM,SCAL> u, v, w;
-    curl_uDvw_minus_Duvw (const AutoDiff<DIM,SCAL> au, 
-                          const AutoDiff<DIM,SCAL> av,
-                          const AutoDiff<DIM,SCAL> aw)
-      : u(au), v(av), w(aw) { ; }
-  };
-
-  /*
-  template <int DIM, typename SCAL>
-  INLINE Class_curl_uDvw_minus_Duvw<DIM,SCAL> 
-  curl_uDvw_minus_Duvw (AutoDiff<DIM,SCAL> u, AutoDiff<DIM,SCAL> v, AutoDiff<DIM,SCAL> w)
-  { return Class_curl_uDvw_minus_Duvw<DIM,SCAL> (u,v,w); }
-  */
-  
-  /*
-  template <int DIM, typename SCAL>
-  INLINE Class_curl_uDvw_minus_Duvw<DIM,SCAL> 
-  curl_uDvw_minus_Duvw (AutoDiff<DIM,SCAL> u, AutoDiff<DIM,SCAL> v, AutoDiff<DIM,SCAL> w)
-  { return Class_curl_uDvw_minus_Duvw<DIM,SCAL> (u,v,w); }
-  */
-
-
-#ifdef OLDTHDIV2SHAPE
-  template <int DIM, typename SCAL = double> class THDiv2Shape
-  {
-  public:
-    INLINE operator Vec<DIM,SCAL> () { return 0.0; }
-  };
-
-
-  template <typename SCAL> class THDiv2Shape<2,SCAL>
-  {
-    Vec<2,SCAL> data;
-  public:
-    INLINE THDiv2Shape (Du<2,SCAL> uv)
-    {
-      data = Vec<2,SCAL> (uv.u.DValue(1), -uv.u.DValue(0));
-    }
-    
-    INLINE THDiv2Shape (uDv<2,SCAL> uv)
-    {
-      data = Vec<2,SCAL> (-uv.u.Value()*uv.v.DValue(1), 
-                          uv.u.Value()*uv.v.DValue(0));
-    }
-
-    INLINE THDiv2Shape (const uDv_minus_vDu<2,SCAL> & uv) 
-    { 
-      data(0) = -uv.u.Value() * uv.v.DValue(1) + uv.u.DValue(1) * uv.v.Value();
-      data(1) =  uv.u.Value() * uv.v.DValue(0) - uv.u.DValue(0) * uv.v.Value();
-    }
-
-    INLINE THDiv2Shape (const wuDv_minus_wvDu<2,SCAL> & uv) 
-    { 
-      data[0] = -uv.u.Value() * uv.v.DValue(1) + uv.u.DValue(1) * uv.v.Value();
-      data[1] =  uv.u.Value() * uv.v.DValue(0) - uv.u.DValue(0) * uv.v.Value();
-      data[0] *= uv.w.Value();
-      data[1] *= uv.w.Value();
-    }
-    
-    INLINE operator Vec<2,SCAL> () const { return data; }
-  };
-
-
-  template <typename SCAL> class THDiv2Shape<3,SCAL>
-  {
-    Vec<3,SCAL> data;
-  public:
-
-    INLINE THDiv2Shape (const uDvDw_Cyclic<3,SCAL> & uvw) 
-    { 
-      /*
-      AutoDiff<3,SCAL> hv =
-        uvw.u.Value() * Cross (uvw.v, uvw.w) +
-        uvw.v.Value() * Cross (uvw.w, uvw.u) +
-        uvw.w.Value() * Cross (uvw.u, uvw.v);
-
-      for (int i = 0; i < 3; i++)
-        data[i] = hv.DValue(i);
-      */
-      AutoDiff<3,SCAL> p1 = Cross (uvw.v, uvw.w);
-      AutoDiff<3,SCAL> p2 = Cross (uvw.w, uvw.u);
-      AutoDiff<3,SCAL> p3 = Cross (uvw.u, uvw.v);
-
-      for (int i = 0; i < 3; i++)
-        data[i] =
-          uvw.u.Value() * p1.DValue(i) + 
-          uvw.v.Value() * p2.DValue(i) + 
-          uvw.w.Value() * p3.DValue(i);
-    }
-
-
-    INLINE THDiv2Shape (const z_times_uDvDw_Cyclic<3,SCAL> & uvwz) 
-    { 
-      /*
-      AutoDiff<3,SCAL> hv =
-        uvw.u.Value() * Cross (uvw.v, uvw.w) +
-        uvw.v.Value() * Cross (uvw.w, uvw.u) +
-        uvw.w.Value() * Cross (uvw.u, uvw.v);
-
-      for (int i = 0; i < 3; i++)
-        data[i] = hv.DValue(i);
-      */
-      AutoDiff<3,SCAL> p1 = Cross (uvwz.v, uvwz.w);
-      AutoDiff<3,SCAL> p2 = Cross (uvwz.w, uvwz.u);
-      AutoDiff<3,SCAL> p3 = Cross (uvwz.u, uvwz.v);
-
-      for (int i = 0; i < 3; i++)
-        data[i] =
-          uvwz.z.Value() * uvwz.u.Value() * p1.DValue(i) + 
-          uvwz.z.Value() * uvwz.v.Value() * p2.DValue(i) + 
-          uvwz.z.Value() * uvwz.w.Value() * p3.DValue(i);
-    }
-    
-    INLINE THDiv2Shape (const Du_Cross_Dv<3,SCAL> & uv) 
-    { 
-      AutoDiff<3,SCAL> hv = Cross (uv.u, uv.v);
-      for (int i = 0; i < 3; i++)
-        data[i] = hv.DValue(i);
-    }
-
-    INLINE THDiv2Shape (const wDu_Cross_Dv<3,SCAL> & uvw) 
-    { 
-      AutoDiff<3,SCAL> hv = Cross (uvw.u, uvw.v);
-      for (int i = 0; i < 3; i++)
-        data[i] = uvw.w.Value() * hv.DValue(i);
-    }
-
-
-    INLINE THDiv2Shape (const uDvDw_minus_DuvDw<3,SCAL> & uvw) 
-    { 
-      /*
-      AutoDiff<3,SCAL> hv =
-        uvw.u.Value() * Cross (uvw.v, uvw.w) +
-        uvw.v.Value() * Cross (uvw.w, uvw.u);
-
-      for (int i = 0; i < 3; i++)
-        data[i] = hv.DValue(i);
-      */
-      AutoDiff<3,SCAL> p1 = Cross (uvw.v, uvw.w);
-      AutoDiff<3,SCAL> p2 = Cross (uvw.w, uvw.u);
-
-      for (int i = 0; i < 3; i++)
-        data[i] =
-          uvw.u.Value() * p1.DValue(i) + 
-          uvw.v.Value() * p2.DValue(i);
-    }
-
-    INLINE THDiv2Shape (const curl_uDvw_minus_Duvw<3,SCAL> & uvw) 
-    { 
-      AutoDiff<3,SCAL> hv = Cross (uvw.u*uvw.w, uvw.v) - Cross (uvw.v*uvw.w, uvw.u);
-      for (int i = 0; i < 3; i++)
-        data[i] = hv.DValue(i);
-    }
-
-    INLINE operator Vec<3,SCAL> () const { return data; }
-  };
-#endif
-
-
-
-  /////////////////////
-
-
-  template <typename SCAL>
-  INLINE auto HDiv2ShapeNew (Du<2,SCAL> uv)
-  {
-    return Vec<2,SCAL> (-uv.u.DValue(0), -uv.u.DValue(1));  // signs to fit historic inconsistency
-  }
-
-  
-  template <typename SCAL>
-  INLINE auto HDiv2ShapeNew (uDv<2,SCAL> uv)
-  {
-    return Vec<2,SCAL> (uv.u.Value()*uv.v.DValue(0), 
-                        uv.u.Value()*uv.v.DValue(1));
-  }
-  
-  template <typename SCAL>
-  INLINE auto HDiv2ShapeNew (const uDv_minus_vDu<2,SCAL> & uv) 
-  {
-    Vec<2,SCAL> data;
-    data(0) = uv.u.Value() * uv.v.DValue(0) - uv.u.DValue(0) * uv.v.Value();
-    data(1) = uv.u.Value() * uv.v.DValue(1) - uv.u.DValue(1) * uv.v.Value();
-    return data;
-  }
-
-  template <typename SCAL>
-  INLINE auto HDiv2ShapeNew (const wuDv_minus_wvDu<2,SCAL> & uv) 
-  {
-    Vec<2,SCAL> data;    
-    data[0] =  uv.u.Value() * uv.v.DValue(0) - uv.u.DValue(0) * uv.v.Value();
-    data[1] =  uv.u.Value() * uv.v.DValue(1) - uv.u.DValue(1) * uv.v.Value();
-    data[0] *= uv.w.Value();
-    data[1] *= uv.w.Value();
-    return data;
-  }
-  
-  template <typename SCAL>  
-  INLINE auto HDiv2ShapeNew (Du<3,SCAL> uv)
-  {
-    return Vec<3,SCAL> (-uv.u.DValue(0), -uv.u.DValue(1), -uv.u.DValue(2));
-  }
-  
-  template <typename SCAL>
-  INLINE auto HDiv2ShapeNew (uDv<3,SCAL> uv)
-  {
-    return Vec<3,SCAL> (uv.u.Value()*uv.v.DValue(0), 
-                        uv.u.Value()*uv.v.DValue(1), 
-                        uv.u.Value()*uv.v.DValue(2));
-  }
-  
-  template <typename SCAL>
-  INLINE auto HDiv2ShapeNew (const uDv_minus_vDu<3,SCAL> & uv) 
-  {
-    Vec<3,SCAL> data;
-    data(0) = uv.u.Value() * uv.v.DValue(0) - uv.u.DValue(0) * uv.v.Value();
-    data(1) = uv.u.Value() * uv.v.DValue(1) - uv.u.DValue(1) * uv.v.Value();
-    data(2) = uv.u.Value() * uv.v.DValue(2) - uv.u.DValue(2) * uv.v.Value();
-    return data;
-  }
-  
-  template <typename SCAL>
-  INLINE auto HDiv2ShapeNew (const wuDv_minus_wvDu<3,SCAL> & uv) 
-  {
-    Vec<3,SCAL> data;
-    data[0] =  uv.u.Value() * uv.v.DValue(0) - uv.u.DValue(0) * uv.v.Value();
-    data[1] =  uv.u.Value() * uv.v.DValue(1) - uv.u.DValue(1) * uv.v.Value();
-    data[2] =  uv.u.Value() * uv.v.DValue(2) - uv.u.DValue(2) * uv.v.Value();
-    data[0] *= uv.w.Value();
-    data[1] *= uv.w.Value();
-    data[2] *= uv.w.Value();
-    return data;
-  }
-    
-  template <typename SCAL>
-  INLINE auto HDiv2ShapeNew (const uDvDw_Cyclic<3,SCAL> & uvw) 
-    { 
-      AutoDiff<3,SCAL> p1 = Cross (uvw.v, uvw.w);
-      AutoDiff<3,SCAL> p2 = Cross (uvw.w, uvw.u);
-      AutoDiff<3,SCAL> p3 = Cross (uvw.u, uvw.v);
-
-      Vec<3,SCAL> data;
-      for (int i = 0; i < 3; i++)
-        data[i] =
-          uvw.u.Value() * p1.DValue(i) + 
-          uvw.v.Value() * p2.DValue(i) + 
-          uvw.w.Value() * p3.DValue(i);
-      return data;
-    }
-
-  template <typename SCAL>
-  INLINE auto HDiv2ShapeNew (const z_times_uDvDw_Cyclic<3,SCAL> & uvwz) 
-    { 
-      AutoDiff<3,SCAL> p1 = Cross (uvwz.v, uvwz.w);
-      AutoDiff<3,SCAL> p2 = Cross (uvwz.w, uvwz.u);
-      AutoDiff<3,SCAL> p3 = Cross (uvwz.u, uvwz.v);
-
-      Vec<3,SCAL> data;
-      for (int i = 0; i < 3; i++)
-        data[i] =
-          uvwz.z.Value() * uvwz.u.Value() * p1.DValue(i) + 
-          uvwz.z.Value() * uvwz.v.Value() * p2.DValue(i) + 
-          uvwz.z.Value() *  uvwz.w.Value() * p3.DValue(i);
-      return data;
-    }
-
-  
-  template <typename SCAL>
-  INLINE auto HDiv2ShapeNew (const Du_Cross_Dv<3,SCAL> & uv) 
-  { 
-    AutoDiff<3,SCAL> hv = Cross (uv.u, uv.v);
-    Vec<3,SCAL> data;
-    for (int i = 0; i < 3; i++)
-      data[i] = hv.DValue(i);
-    return data;
-  }
-  
-  template <typename SCAL>
-  INLINE auto HDiv2ShapeNew (const wDu_Cross_Dv<3,SCAL> & uvw) 
-  { 
-    AutoDiff<3,SCAL> hv = Cross (uvw.u, uvw.v);
-    Vec<3,SCAL> data;
-    for (int i = 0; i < 3; i++)
-      data[i] = uvw.w.Value() * hv.DValue(i);
-    return data;
-  }
-  
-
-  template <typename SCAL>
-  INLINE auto HDiv2ShapeNew (const uDvDw_minus_DuvDw<3,SCAL> & uvw) 
-  { 
-    /*
-      AutoDiff<3,SCAL> hv =
-      uvw.u.Value() * Cross (uvw.v, uvw.w) +
-        uvw.v.Value() * Cross (uvw.w, uvw.u);
-
-      for (int i = 0; i < 3; i++)
-        data[i] = hv.DValue(i);
-    */
-    AutoDiff<3,SCAL> p1 = Cross (uvw.v, uvw.w);
-    AutoDiff<3,SCAL> p2 = Cross (uvw.w, uvw.u);
-
-    Vec<3,SCAL> data;
-    for (int i = 0; i < 3; i++)
-      data[i] =
-        uvw.u.Value() * p1.DValue(i) + 
-        uvw.v.Value() * p2.DValue(i);
-    return data;
-  }
-  
-  template <typename SCAL>
-  INLINE auto HDiv2ShapeNew (const curl_uDvw_minus_Duvw<3,SCAL> & uvw) 
-  { 
-    AutoDiff<3,SCAL> hv = Cross (uvw.u*uvw.w, uvw.v) - Cross (uvw.v*uvw.w, uvw.u);
-    Vec<3,SCAL> data;
-    for (int i = 0; i < 3; i++)
-      data[i] = hv.DValue(i);
-    return data;
-  }
-
-
-  //////////////////////////
-
-
-
-  
-
-  
-
-  /*
-  // 2D 
-  template <int DIM>
-  class HDivShapeElement
-  {
-    double * data;
-  public:
-    HDivShapeElement (double * adata) : data(adata) { ; }
-
-    void operator= (THDiv2Shape<DIM> hd2vec)
-    {
-      Vec<DIM> v = hd2vec;
-      for (int j = 0; j < DIM; j++)
-        data[j] = v(j);
-    }
-  };
-
-
-  template <int DIM>
-  class HDivEvaluateShapeElement
-  {
-    const double * coefs;
-    Vec<DIM> & sum;
-  public:
-    HDivEvaluateShapeElement (const double * acoefs, Vec<DIM> & asum)
-      : coefs(acoefs), sum(asum) { ; }
-
-
-    void operator= (THDiv2Shape<DIM> hd2vec)
-    {
-      sum += *coefs * Vec<DIM> (hd2vec);
-    }
-  };
-  */
-
-
-
-
-
-
-  template <int DIM, typename SCAL = double> class THDiv2DivShape
-  {
-  public:
-    INLINE operator SCAL () const { return SCAL(0.0); }
-  };
-
-  template <typename SCAL> class THDiv2DivShape<2,SCAL>
-  {
-    SCAL data;
-  public:
-    INLINE THDiv2DivShape (Du<2,SCAL> uv)
-    {
-      data = SCAL(0.0);
-    }
-    
-    INLINE THDiv2DivShape (uDv<2,SCAL> uv)
-    {
-      AutoDiff<1,SCAL> hd = Cross (uv.u, uv.v);
-      data = -hd.DValue(0);
-    }
-
-    INLINE THDiv2DivShape (const uDv_minus_vDu<2,SCAL> & uv) 
-    { 
-      data = -2*uv.u.DValue(0) * uv.v.DValue(1) 
-        + 2*uv.u.DValue(1) * uv.v.DValue(0);
-    }
-
-    INLINE THDiv2DivShape (const wuDv_minus_wvDu<2,SCAL> & uv) 
-    { 
-      AutoDiff<1,SCAL> hd = Cross (uv.u*uv.w, uv.v) + Cross(uv.u, uv.v*uv.w);
-      data = -hd.DValue(0);
-    }
-    
-    INLINE operator SCAL () const { return data; }
-    INLINE SCAL Get() const { return data; }
-  };
-
-
-  template <typename SCAL> class THDiv2DivShape<3,SCAL>
-  {
-    SCAL data;
-  public:
-
-    INLINE THDiv2DivShape (const uDvDw_Cyclic<3,SCAL> & uvw) 
-    { 
-      data = 
-        Dot (uvw.u, Cross (uvw.v, uvw.w)) +
-        Dot (uvw.v, Cross (uvw.w, uvw.u)) +
-        Dot (uvw.w, Cross (uvw.u, uvw.v));
-    }
-
-    INLINE THDiv2DivShape (const z_times_uDvDw_Cyclic<3,SCAL> & uvwz) 
-    { 
-      data = 
-        Dot (uvwz.z * uvwz.u, Cross (uvwz.v, uvwz.w)) +
-        Dot (uvwz.z * uvwz.v, Cross (uvwz.w, uvwz.u)) +
-        Dot (uvwz.z * uvwz.w, Cross (uvwz.u, uvwz.v));
-    }
-
-    INLINE THDiv2DivShape (const Du_Cross_Dv<3,SCAL> & uv) 
-    { 
-      data = 0.0;
-    }
-
-    INLINE THDiv2DivShape (const wDu_Cross_Dv<3,SCAL> & uvw) 
-    { 
-      data = Dot (uvw.w, Cross (uvw.u, uvw.v));
-    }
-
-    INLINE THDiv2DivShape (const uDvDw_minus_DuvDw<3,SCAL> & uvw) 
-    { 
-      data = 
-        Dot (uvw.u, Cross (uvw.v, uvw.w)) +
-        Dot (uvw.v, Cross (uvw.w, uvw.u));
-    }
-
-    INLINE THDiv2DivShape (const curl_uDvw_minus_Duvw<3,SCAL> & uvw) 
-    { 
-      data = SCAL(0.0);
-    }
-
-    INLINE operator SCAL () const { return data; }
-    INLINE SCAL Get() const { return data; }    
-  };
-
-
-
-
-  /*
-  template <int DIM, typename SCAL = double>
-  class HDivDivShapeElement
-  {
-    SCAL * data;
-  public:
-    HDivDivShapeElement (SCAL * adata) : data(adata) { ; }
-
-
-    void operator= (const THDiv2DivShape<DIM,SCAL> & hd2dshape)
-    {
-      *data = hd2dshape;
-    }
-  };
-
-  template <int DIM>
-  class HDivShapeAssign
-  {
-    double * dshape;
-  public:
-    HDivShapeAssign (FlatMatrixFixWidth<DIM> mat)
-    { dshape = &mat(0,0); }
-    
-    HDivShapeElement<DIM> operator[] (int i) const
-    { return HDivShapeElement<DIM> (dshape + i*DIM); }
-  };
-
-
-  template <int DIM>
-  class HDivDivShapeAssign
-  {
-    SliceVector<> dshape;
-  public:
-    HDivDivShapeAssign (SliceVector<>  mat)
-      : dshape(mat) { ; }
-
-    HDivDivShapeElement<DIM> operator[] (int i) const
-    { return HDivDivShapeElement<DIM> (&dshape(i)); }
-  };
-
-  template <int DIM>
-  class HDivEvaluateShape
-  {
-    const double * coefs;
-    Vec<DIM> sum;
-  public:
-    HDivEvaluateShape (FlatVector<> acoefs)
-    { coefs = &acoefs(0); sum = 0.0; }
-
-    HDivEvaluateShapeElement<DIM> operator[] (int i) 
-    { return HDivEvaluateShapeElement<DIM> (coefs+i, sum); }
-
-    Vec<DIM> Sum() { return sum; }
-  };
-  */
-
-
-
-
-
-  
-
-
-}
-
-
-#endif