meshpatch.h

Go to the documentation of this file.

#ifndef MESHPATCH
#define MESHPATCH

#include <point.h>
#include <print.h>

#include <bernstein.h>

typedef point3<double> pt3;

class meshpatch
{
public:

  uintc M;
  uintc N;

protected:
  bernsteinPoly * const vbu;
  bernsteinPoly * const vbv;

  double * const vbuval;
  double * const vbvval;

public:

  meshpatch(uintc M_, uintc N_);

  ~meshpatch();

  template< typename T, typename W >
  void eval
  ( 
    T & p, 
    double const u, 
    double const v,
    W const & pij 
  ) const
  {
    for (uint i=0; i<=M; ++i)
      vbuval[i] = vbu[i].eval(u);

    for (uint i=0; i<=N; ++i)
      vbvval[i] = vbv[i].eval(v);

    p = T();
    T p2;
    uint k;
    for (uint i=0; i<=M; ++i)
    {
      T p2;
      for (k=0; k<=N; ++k)
        p2 += pij(i,k)*vbvval[k];

      p += p2*vbuval[i];
    }
  }

};



#endif