eic-opticks/ana/ellipse_intersect_circle.cc

0001 //  gcc ellipse_intersect_circle.cc -lm -lstdc++ -o /tmp/ellipse_intersect_circle  && /tmp/ellipse_intersect_circle
0002 //  head -1 ellipse_intersect_circle.cc | perl -pe 's,//,,' - | sh
0003
0004 #include "ellipse_intersect_circle.hh"
0005
0006 int main(int argc, char** argv)
0007 {
0008      double e_cx = 0. ;
0009      double e_cy = -5. ;
0010      double e_ax = 254. ;
0011      double e_ay = 190. ;
0012
0013      double c_cx = 207. ;   // torus_x
0014      double c_cy = -210. ;   // torus_z
0015      double c_r = 80. ;  // torus_r
0016
0017      int n = 1000000 ;
0018      bool verbose = false ;
0019
0020      Ellipse_Intersect_Circle ec = Ellipse_Intersect_Circle::make( e_cx, e_cy, e_ax, e_ay, c_cx, c_cy, c_r, n, verbose );
0021
0022      if(ec.intersect.count > 0)
0023      printf(" (%10.4f, %10.4f) \n", ec.intersect.p[0].x, ec.intersect.p[0].y );
0024
0025      return 0 ;
0026 }