void oldCode(const GlobalPoint & P1, const GlobalPoint & P2) {
typedef TkRotation<double> Rotation;
typedef Basic2DVector<double> Point2D;
GlobalVector aX = GlobalVector( P1.x(), P1.y(), 0.).unit();
GlobalVector aY( -aX.y(), aX.x(), 0.);
GlobalVector aZ( 0., 0., 1.);
TkRotation<double > theRotation = Rotation(aX,aY,aZ);
PointUV p1(Point2D(P1.x(),P1.y()), &theRotation);
PointUV p2(Point2D(P2.x(),P2.y()), &theRotation);
std::cout << "\nold for " << P1 <<", " << P2 << std::endl;
std::cout << theRotation << std::endl;
std::cout << p1.u() << " " << p1.v() << std::endl;
std::cout << p2.u() << " " << p2.v() << std::endl;
std::cout << p1.unmap() << std::endl;
std::cout << p2.unmap() << std::endl;
}
virtual GlobalVector inTesla ( const GlobalPoint& ) const {
return GlobalVector(0,0,4);
}
