#ifndef jbt_bend_h
#define jbt_bend_h

#include "warp.h"


class jbt_bend : public jbt_warp
{
public:
    

    jbt_bend( jFlt theta = 90, jFlt ymin = -1, jFlt ymax = 1, jFlt y0 = 0,
              jbt_blobtree* obj = 0 );

    jbt_blobtree::Type type() const { return jbt_blobtree::jbtBend; }


    void precompute( const jVec3& );
    jVec3 warp( const jVec3& ) const;
    jVec3 unwarp( const jVec3& ) const;

private:

    jFlt ymin,ymax, y0, theta, thetamin, thetamax, k, ik;
    jFlt Ctheta, Stheta;
    jFlt a,b;
};

#endif