contour_interpolation.h

/* Copyright (C) 2004-2009 MBSim Development Team
 *
 * This library is free software; you can redistribute it and/or 
 * modify it under the terms of the GNU Lesser General Public 
 * License as published by the Free Software Foundation; either 
 * version 2.1 of the License, or (at your option) any later version. 
 *  
 * This library is distributed in the hope that it will be useful, 
 * but WITHOUT ANY WARRANTY; without even the implied warranty of 
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU 
 * Lesser General Public License for more details. 
 *  
 * You should have received a copy of the GNU Lesser General Public 
 * License along with this library; if not, write to the Free Software 
 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301 USA
 *
 * Contact: martin.o.foerg@googlemail.com
 */
 
#ifndef _CONTOUR_INTERPOLATION_H_
#define _CONTOUR_INTERPOLATION_H_
 
#include "mbsim/contours/contour.h"
 
namespace MBSim {
 
  class Point;
 
  class ContourInterpolation : public Contour {
    public:
      ContourInterpolation(const std::string &name, int parameters_, int nPoints_);
      void setPoint(Point *point_, int n);
      std::vector<Point*> getPoints()   const {return iPoints;}
      Point* getPoint(const int n) const {return iPoints[n];}
      int getNPoints() const {return numberOfPoints;}
      int getNZeta() const {return contourParameters;}
       virtual bool testInsideBounds(const fmatvec::Vec2 &zeta) = 0;
       virtual double computePointWeight(const fmatvec::Vec2 &zeta, int i) = 0;
       virtual double computePointWeight(const fmatvec::Vec2 &zeta, int i, int diff) = 0;
//       fmatvec::VecV computePointWeights(const fmatvec::Vec2 &zeta);
//
//       fmatvec::Vec3 getWrOC(const fmatvec::Vec2& s);// {ContourPointData cp; cp.type=EXTINTERPOL;cp.alpha=s; return getWrOC(cp);};
//       fmatvec::Vec3 computeWvC (const fmatvec::Vec2& s);// {ContourPointData cp; cp.type=EXTINTERPOL;cp.alpha=s; return computeWvC (cp);};
//       fmatvec::Mat3xV computeWt  (const fmatvec::Vec2& s);// {ContourPointData cp; cp.type=EXTINTERPOL;cp.alpha=s; return computeWt  (cp);};
//       fmatvec::Vec3 computeWn  (const fmatvec::Vec2& s);// {ContourPointData cp; cp.type=EXTINTERPOL;cp.alpha=s; return computeWn  (cp);};
//
//       fmatvec::Vec3 getWrOC(const ContourPointData &cp);
//       fmatvec::Vec3 computeWvC (const ContourPointData &cp);
//
//       fmatvec::Mat3xV computeWt  (const ContourPointData &cp);
//       virtual fmatvec::Vec3 computeWn  (const ContourPointData &cp) = 0;
 
    protected:
       std::vector<Point*> iPoints;
       int contourParameters;
       int numberOfPoints;
  };
}
 
#endif /* _CONTOUR_INTERPOLATION_H_ */