u_smspl3_curve1d.hpp

/* u_smspl3_curve1d.hpp */
/* $Id$ */
#ifndef __u_smspl3_curve1d_hpp
#define __u_smspl3_curve1d_hpp

#include <vector>
#include <mth/u2d_wpoint.hpp>
#include <mth/u2d_box.hpp>
#include <mth/u_curve1d.hpp>

class USmspl3Curve1d
{
public:
  
  struct Node
  {
    Node(double ix = 0., double iw = 0.):
      x(ix), w(iw) {/* Nothing to do. */};
    
    double x; 
    double w; 
  };
  typedef std::vector<Node> Nodes;
    
  USmspl3Curve1d(double alpha = 0);
  
  USmspl3Curve1d(const Nodes &nodes, double alpha = 0.);
  
  virtual ~USmspl3Curve1d();
  
  
  virtual void approximate(const Nodes &nodes, double alpha = 0.);
  
  virtual double operator()(double p) const;
    
  virtual bool ready() const
  {return m_ready;}
  
  
  virtual void addNode(const Node &node, int id)
  {printf("USmspl3Curve1d::addNode: not implemented.\n");}
  
  virtual bool removeNode(size_t id)
  {printf("USmspl3Curve1d::removeNode: not implemented.\n"); return false;}
  
protected:
  
  struct Coeffs
  {double a,b,c,d,p;};
  typedef std::vector<Coeffs> CoeffsAr;
  
  double m_alpha;

  CoeffsAr m_coeffs;
  
  U2dBox m_bound;
  
  bool m_ready;
  
  // spline linear system matricies and vectors elements
  inline double Am(size_t i, size_t j, const Nodes &nodes) const;
  inline double Hm(size_t i, size_t j, const Nodes &nodes) const;
  inline double Km(size_t i, size_t j, const Nodes &nodes) const;
  inline double bv(size_t i, const Nodes &nodes) const;
  
  bool calc_coeffs(const Nodes &nodes, 
           double alpha,
           CoeffsAr &coeffs);
  
  inline double get_s(double p) const;

private:
  
  double m_h_value;
};

#endif /* u_smspl3_curve1d.hpp */