funcs.scad.cpp

   1 // -*- C -*-
   2
   3 function dist(a,b) =
   4 #define dsq(i) (a[i]-b[i])*(a[i]-b[i])
   5                sqrt(dsq(0) + dsq(1));
   6 #undef dsq
   7
   8 function circle_point(c, r, alpha) = [ c[0] + r * cos(alpha),
   9                                        c[1] + r * sin(alpha) ];
  10
  11 #define d     (dist(a,c))
  12 #define alpha (atan2(a[1]-c[1],a[0]-c[0]))
  13 #define gamma (asin(r / d))
  14 #define beta  (alpha + 90 - gamma)
  15
  16 function tangent_intersect_beta(c,r,a) =
  17               beta;
  18
  19 function tangent_intersect_b(c,r,a) =
  20               circle_point(c, r, beta);
  21 #undef d
  22 #undef alpha
  23 #undef gamma
  24 #undef beta
  25
  26 function reflect_in_y(p) = [-p[0], p[1]];
  27