/***** Autogenerated from runpath3d.in; changes will be overwritten *****/ #line 1 "runtimebase.in" /***** * runtimebase.in * Andy Hammerlindl 2009/07/28 * * Common declarations needed for all code-generating .in files. * *****/ #line 1 "runpath3d.in" /***** * runpath3.in * * Runtime functions for path3 operations. * *****/ #line 1 "runtimebase.in" #include "stack.h" #include "types.h" #include "builtin.h" #include "entry.h" #include "errormsg.h" #include "array.h" #include "triple.h" #include "callable.h" using vm::stack; using vm::error; using vm::array; using vm::callable; using types::formal; using types::function; using camp::triple; #define PRIMITIVE(name,Name,asyName) using types::prim##Name; #include #undef PRIMITIVE typedef double real; void unused(void *); namespace run { array *copyArray(array *a); array *copyArray2(array *a); array *copyArray3(array *a); double *copyArrayC(const array *a, size_t dim=0, GCPlacement placement=NoGC); double *copyArray2C(const array *a, bool square=true, size_t dim2=0, GCPlacement placement=NoGC); triple *copyTripleArrayC(const array *a, size_t dim=0); triple *copyTripleArray2C(const array *a, bool square=true, size_t dim2=0); double *copyTripleArray2Components(array *a, bool square=true, size_t dim2=0, GCPlacement placement=NoGC); } function *realRealFunction(); #define CURRENTPEN processData().currentpen #line 17 "runpath3d.in" #include "path3.h" #include "array.h" #include "drawsurface.h" using namespace camp; using namespace vm; typedef array boolarray; typedef array realarray; typedef array realarray2; typedef array triplearray; typedef array triplearray2; using types::booleanArray; using types::realArray; using types::realArray2; using types::tripleArray; using types::tripleArray2; // Autogenerated routines: namespace run { #line 39 "runpath3d.in" // path3 path3(triplearray *pre, triplearray *point, triplearray *post, boolarray *straight, bool cyclic); void gen_runpath3d0(stack *Stack) { bool cyclic=vm::pop(Stack); boolarray * straight=vm::pop(Stack); triplearray * post=vm::pop(Stack); triplearray * point=vm::pop(Stack); triplearray * pre=vm::pop(Stack); #line 41 "runpath3d.in" size_t n=checkArrays(pre,point); checkEqual(n,checkArray(post)); checkEqual(n,checkArray(straight)); mem::vector nodes(n); for(size_t i=0; i < n; ++i) { nodes[i].pre=read(pre,i); nodes[i].point=read(point,i); nodes[i].post=read(post,i); nodes[i].straight=read(straight,i); } {Stack->push(path3(nodes,(Int) n,cyclic)); return;} } #line 56 "runpath3d.in" void nullPath3(stack *Stack) { #line 57 "runpath3d.in" {Stack->push(nullpath3); return;} } #line 61 "runpath3d.in" // bool ==(path3 a, path3 b); void gen_runpath3d2(stack *Stack) { path3 b=vm::pop(Stack); path3 a=vm::pop(Stack); #line 62 "runpath3d.in" {Stack->push(a == b); return;} } #line 66 "runpath3d.in" // bool !=(path3 a, path3 b); void gen_runpath3d3(stack *Stack) { path3 b=vm::pop(Stack); path3 a=vm::pop(Stack); #line 67 "runpath3d.in" {Stack->push(!(a == b)); return;} } #line 71 "runpath3d.in" // triple point(path3 p, Int t); void gen_runpath3d4(stack *Stack) { Int t=vm::pop(Stack); path3 p=vm::pop(Stack); #line 72 "runpath3d.in" {Stack->push(p.point((Int) t)); return;} } #line 76 "runpath3d.in" // triple point(path3 p, real t); void gen_runpath3d5(stack *Stack) { real t=vm::pop(Stack); path3 p=vm::pop(Stack); #line 77 "runpath3d.in" {Stack->push(p.point(t)); return;} } #line 81 "runpath3d.in" // triple precontrol(path3 p, Int t); void gen_runpath3d6(stack *Stack) { Int t=vm::pop(Stack); path3 p=vm::pop(Stack); #line 82 "runpath3d.in" {Stack->push(p.precontrol((Int) t)); return;} } #line 86 "runpath3d.in" // triple precontrol(path3 p, real t); void gen_runpath3d7(stack *Stack) { real t=vm::pop(Stack); path3 p=vm::pop(Stack); #line 87 "runpath3d.in" {Stack->push(p.precontrol(t)); return;} } #line 91 "runpath3d.in" // triple postcontrol(path3 p, Int t); void gen_runpath3d8(stack *Stack) { Int t=vm::pop(Stack); path3 p=vm::pop(Stack); #line 92 "runpath3d.in" {Stack->push(p.postcontrol((Int) t)); return;} } #line 96 "runpath3d.in" // triple postcontrol(path3 p, real t); void gen_runpath3d9(stack *Stack) { real t=vm::pop(Stack); path3 p=vm::pop(Stack); #line 97 "runpath3d.in" {Stack->push(p.postcontrol(t)); return;} } #line 101 "runpath3d.in" // triple dir(path3 p, Int t, Int sign=0, bool normalize=true); void gen_runpath3d10(stack *Stack) { bool normalize=vm::pop(Stack,true); Int sign=vm::pop(Stack,0); Int t=vm::pop(Stack); path3 p=vm::pop(Stack); #line 102 "runpath3d.in" {Stack->push(p.dir(t,sign,normalize)); return;} } #line 106 "runpath3d.in" // triple dir(path3 p, real t, bool normalize=true); void gen_runpath3d11(stack *Stack) { bool normalize=vm::pop(Stack,true); real t=vm::pop(Stack); path3 p=vm::pop(Stack); #line 107 "runpath3d.in" {Stack->push(p.dir(t,normalize)); return;} } #line 111 "runpath3d.in" // triple accel(path3 p, Int t, Int sign=0); void gen_runpath3d12(stack *Stack) { Int sign=vm::pop(Stack,0); Int t=vm::pop(Stack); path3 p=vm::pop(Stack); #line 112 "runpath3d.in" {Stack->push(p.accel(t,sign)); return;} } #line 116 "runpath3d.in" // triple accel(path3 p, real t); void gen_runpath3d13(stack *Stack) { real t=vm::pop(Stack); path3 p=vm::pop(Stack); #line 117 "runpath3d.in" {Stack->push(p.accel(t)); return;} } #line 121 "runpath3d.in" // real radius(path3 p, real t); void gen_runpath3d14(stack *Stack) { real t=vm::pop(Stack); path3 p=vm::pop(Stack); #line 122 "runpath3d.in" triple v=p.dir(t,false); triple a=p.accel(t); real d=dot(a,v); real v2=v.abs2(); real a2=a.abs2(); real denom=v2*a2-d*d; real r=v2*sqrt(v2); {Stack->push(denom > 0 ? r/sqrt(denom) : 0.0); return;} } #line 133 "runpath3d.in" // real radius(triple z0, triple c0, triple c1, triple z1, real t); void gen_runpath3d15(stack *Stack) { real t=vm::pop(Stack); triple z1=vm::pop(Stack); triple c1=vm::pop(Stack); triple c0=vm::pop(Stack); triple z0=vm::pop(Stack); #line 134 "runpath3d.in" triple v=(3.0*(z1-z0)+9.0*(c0-c1))*t*t+(6.0*(z0+c1)-12.0*c0)*t+3.0*(c0-z0); triple a=6.0*(z1-z0+3.0*(c0-c1))*t+6.0*(z0+c1)-12.0*c0; real d=dot(a,v); real v2=v.abs2(); real a2=a.abs2(); real denom=v2*a2-d*d; real r=v2*sqrt(v2); {Stack->push(denom > 0 ? r/sqrt(denom) : 0.0); return;} } #line 145 "runpath3d.in" // path3 reverse(path3 p); void gen_runpath3d16(stack *Stack) { path3 p=vm::pop(Stack); #line 146 "runpath3d.in" {Stack->push(p.reverse()); return;} } #line 150 "runpath3d.in" // path3 subpath(path3 p, Int a, Int b); void gen_runpath3d17(stack *Stack) { Int b=vm::pop(Stack); Int a=vm::pop(Stack); path3 p=vm::pop(Stack); #line 151 "runpath3d.in" {Stack->push(p.subpath((Int) a, (Int) b)); return;} } #line 155 "runpath3d.in" // path3 subpath(path3 p, real a, real b); void gen_runpath3d18(stack *Stack) { real b=vm::pop(Stack); real a=vm::pop(Stack); path3 p=vm::pop(Stack); #line 156 "runpath3d.in" {Stack->push(p.subpath(a,b)); return;} } #line 160 "runpath3d.in" // Int length(path3 p); void gen_runpath3d19(stack *Stack) { path3 p=vm::pop(Stack); #line 161 "runpath3d.in" {Stack->push(p.length()); return;} } #line 165 "runpath3d.in" // bool cyclic(path3 p); void gen_runpath3d20(stack *Stack) { path3 p=vm::pop(Stack); #line 166 "runpath3d.in" {Stack->push(p.cyclic()); return;} } #line 170 "runpath3d.in" // bool straight(path3 p, Int t); void gen_runpath3d21(stack *Stack) { Int t=vm::pop(Stack); path3 p=vm::pop(Stack); #line 171 "runpath3d.in" {Stack->push(p.straight(t)); return;} } #line 175 "runpath3d.in" // path3 unstraighten(path3 p); void gen_runpath3d22(stack *Stack) { path3 p=vm::pop(Stack); #line 176 "runpath3d.in" {Stack->push(p.unstraighten()); return;} } // Return the maximum perpendicular deviation of segment i of path3 g // from a straight line. #line 182 "runpath3d.in" // real straightness(path3 p, Int t); void gen_runpath3d23(stack *Stack) { Int t=vm::pop(Stack); path3 p=vm::pop(Stack); #line 183 "runpath3d.in" if(p.straight(t)) {Stack->push(0); return;} triple z0=p.point(t); triple u=unit(p.point(t+1)-z0); {Stack->push(::max(length(perp(p.postcontrol(t)-z0,u)), length(perp(p.precontrol(t+1)-z0,u)))); return;} } // Return the maximum perpendicular deviation of z0..controls c0 and c1..z1 // from a straight line. #line 193 "runpath3d.in" // real straightness(triple z0, triple c0, triple c1, triple z1); void gen_runpath3d24(stack *Stack) { triple z1=vm::pop(Stack); triple c1=vm::pop(Stack); triple c0=vm::pop(Stack); triple z0=vm::pop(Stack); #line 194 "runpath3d.in" triple u=unit(z1-z0); {Stack->push(::max(length(perp(c0-z0,u)),length(perp(c1-z0,u)))); return;} } #line 199 "runpath3d.in" // bool piecewisestraight(path3 p); void gen_runpath3d25(stack *Stack) { path3 p=vm::pop(Stack); #line 200 "runpath3d.in" {Stack->push(p.piecewisestraight()); return;} } #line 204 "runpath3d.in" // real arclength(path3 p); void gen_runpath3d26(stack *Stack) { path3 p=vm::pop(Stack); #line 205 "runpath3d.in" {Stack->push(p.arclength()); return;} } #line 209 "runpath3d.in" // real arctime(path3 p, real dval); void gen_runpath3d27(stack *Stack) { real dval=vm::pop(Stack); path3 p=vm::pop(Stack); #line 210 "runpath3d.in" {Stack->push(p.arctime(dval)); return;} } #line 214 "runpath3d.in" // realarray* intersect(path3 p, path3 q, real fuzz=-1); void gen_runpath3d28(stack *Stack) { real fuzz=vm::pop(Stack,-1); path3 q=vm::pop(Stack); path3 p=vm::pop(Stack); #line 215 "runpath3d.in" bool exact=fuzz <= 0.0; if(fuzz < 0) fuzz=BigFuzz*::max(::max(length(p.max()),length(p.min())), ::max(length(q.max()),length(q.min()))); std::vector S,T; real s,t; if(intersections(s,t,S,T,p,q,fuzz,true,exact)) { array *V=new array(2); (*V)[0]=s; (*V)[1]=t; {Stack->push(V); return;} } else {Stack->push(new array(0)); return;} } #line 232 "runpath3d.in" // realarray2* intersections(path3 p, path3 q, real fuzz=-1); void gen_runpath3d29(stack *Stack) { real fuzz=vm::pop(Stack,-1); path3 q=vm::pop(Stack); path3 p=vm::pop(Stack); #line 233 "runpath3d.in" bool exact=fuzz <= 0.0; if(fuzz < 0) fuzz=BigFuzz*::max(::max(length(p.max()),length(p.min())), ::max(length(q.max()),length(q.min()))); bool single=!exact; real s,t; std::vector S,T; bool found=intersections(s,t,S,T,p,q,fuzz,single,exact); if(!found) {Stack->push(new array(0)); return;} array *V; if(single) { V=new array(1); array *Vi=new array(2); (*V)[0]=Vi; (*Vi)[0]=s; (*Vi)[1]=t; } else { size_t n=S.size(); V=new array(n); for(size_t i=0; i < n; ++i) { array *Vi=new array(2); (*V)[i]=Vi; (*Vi)[0]=S[i]; (*Vi)[1]=T[i]; } } stable_sort(V->begin(),V->end(),run::compare2()); {Stack->push(V); return;} } #line 265 "runpath3d.in" // realarray2* intersections(path3 p, triplearray2 *P, real fuzz=-1); void gen_runpath3d30(stack *Stack) { real fuzz=vm::pop(Stack,-1); triplearray2 * P=vm::pop(Stack); path3 p=vm::pop(Stack); #line 266 "runpath3d.in" triple *A=copyTripleArray2C(P,true,4); if(fuzz <= 0) fuzz=BigFuzz*::max(::max(length(p.max()),length(p.min())), norm(A,16)); std::vector T,U,V; intersections(T,U,V,p,A,fuzz); delete[] A; size_t n=T.size(); array *W=new array(n); for(size_t i=0; i < n; ++i) { array *Wi=new array(3); (*W)[i]=Wi; (*Wi)[0]=T[i]; (*Wi)[1]=U[i]; (*Wi)[2]=V[i]; } {Stack->push(W); return;} // Sorting will done in asy. } #line 285 "runpath3d.in" // Int size(path3 p); void gen_runpath3d31(stack *Stack) { path3 p=vm::pop(Stack); #line 286 "runpath3d.in" {Stack->push(p.size()); return;} } #line 290 "runpath3d.in" // path3 &(path3 p, path3 q); void gen_runpath3d32(stack *Stack) { path3 q=vm::pop(Stack); path3 p=vm::pop(Stack); #line 291 "runpath3d.in" {Stack->push(camp::concat(p,q)); return;} } #line 295 "runpath3d.in" // triple min(path3 p); void gen_runpath3d33(stack *Stack) { path3 p=vm::pop(Stack); #line 296 "runpath3d.in" {Stack->push(p.min()); return;} } #line 300 "runpath3d.in" // triple max(path3 p); void gen_runpath3d34(stack *Stack) { path3 p=vm::pop(Stack); #line 301 "runpath3d.in" {Stack->push(p.max()); return;} } #line 305 "runpath3d.in" // realarray* mintimes(path3 p); void gen_runpath3d35(stack *Stack) { path3 p=vm::pop(Stack); #line 306 "runpath3d.in" array *V=new array(3); triple v=p.mintimes(); (*V)[0]=v.getx(); (*V)[1]=v.gety(); (*V)[2]=v.getz(); {Stack->push(V); return;} } #line 315 "runpath3d.in" // realarray* maxtimes(path3 p); void gen_runpath3d36(stack *Stack) { path3 p=vm::pop(Stack); #line 316 "runpath3d.in" array *V=new array(3); triple v=p.maxtimes(); (*V)[0]=v.getx(); (*V)[1]=v.gety(); (*V)[2]=v.getz(); {Stack->push(V); return;} } #line 325 "runpath3d.in" // path3 *(realarray2 *t, path3 g); void gen_runpath3d37(stack *Stack) { path3 g=vm::pop(Stack); realarray2 * t=vm::pop(Stack); #line 326 "runpath3d.in" {Stack->push(transformed(*t,g)); return;} } #line 330 "runpath3d.in" // pair minratio(path3 g); void gen_runpath3d38(stack *Stack) { path3 g=vm::pop(Stack); #line 331 "runpath3d.in" {Stack->push(g.ratio(::min)); return;} } #line 335 "runpath3d.in" // pair maxratio(path3 g); void gen_runpath3d39(stack *Stack) { path3 g=vm::pop(Stack); #line 336 "runpath3d.in" {Stack->push(g.ratio(::max)); return;} } } // namespace run namespace trans { void gen_runpath3d_venv(venv &ve) { #line 39 "runpath3d.in" addFunc(ve, run::gen_runpath3d0, primPath3(), "path3", formal(tripleArray(), "pre", false, false), formal(tripleArray(), "point", false, false), formal(tripleArray(), "post", false, false), formal(booleanArray(), "straight", false, false), formal(primBoolean(), "cyclic", false, false)); #line 56 "runpath3d.in" REGISTER_BLTIN(run::nullPath3,"nullPath3"); #line 61 "runpath3d.in" addFunc(ve, run::gen_runpath3d2, primBoolean(), "==", formal(primPath3(), "a", false, false), formal(primPath3(), "b", false, false)); #line 66 "runpath3d.in" addFunc(ve, run::gen_runpath3d3, primBoolean(), "!=", formal(primPath3(), "a", false, false), formal(primPath3(), "b", false, false)); #line 71 "runpath3d.in" addFunc(ve, run::gen_runpath3d4, primTriple(), "point", formal(primPath3(), "p", false, false), formal(primInt(), "t", false, false)); #line 76 "runpath3d.in" addFunc(ve, run::gen_runpath3d5, primTriple(), "point", formal(primPath3(), "p", false, false), formal(primReal(), "t", false, false)); #line 81 "runpath3d.in" addFunc(ve, run::gen_runpath3d6, primTriple(), "precontrol", formal(primPath3(), "p", false, false), formal(primInt(), "t", false, false)); #line 86 "runpath3d.in" addFunc(ve, run::gen_runpath3d7, primTriple(), "precontrol", formal(primPath3(), "p", false, false), formal(primReal(), "t", false, false)); #line 91 "runpath3d.in" addFunc(ve, run::gen_runpath3d8, primTriple(), "postcontrol", formal(primPath3(), "p", false, false), formal(primInt(), "t", false, false)); #line 96 "runpath3d.in" addFunc(ve, run::gen_runpath3d9, primTriple(), "postcontrol", formal(primPath3(), "p", false, false), formal(primReal(), "t", false, false)); #line 101 "runpath3d.in" addFunc(ve, run::gen_runpath3d10, primTriple(), "dir", formal(primPath3(), "p", false, false), formal(primInt(), "t", false, false), formal(primInt(), "sign", true, false), formal(primBoolean(), "normalize", true, false)); #line 106 "runpath3d.in" addFunc(ve, run::gen_runpath3d11, primTriple(), "dir", formal(primPath3(), "p", false, false), formal(primReal(), "t", false, false), formal(primBoolean(), "normalize", true, false)); #line 111 "runpath3d.in" addFunc(ve, run::gen_runpath3d12, primTriple(), "accel", formal(primPath3(), "p", false, false), formal(primInt(), "t", false, false), formal(primInt(), "sign", true, false)); #line 116 "runpath3d.in" addFunc(ve, run::gen_runpath3d13, primTriple(), "accel", formal(primPath3(), "p", false, false), formal(primReal(), "t", false, false)); #line 121 "runpath3d.in" addFunc(ve, run::gen_runpath3d14, primReal(), "radius", formal(primPath3(), "p", false, false), formal(primReal(), "t", false, false)); #line 133 "runpath3d.in" addFunc(ve, run::gen_runpath3d15, primReal(), "radius", formal(primTriple(), "z0", false, false), formal(primTriple(), "c0", false, false), formal(primTriple(), "c1", false, false), formal(primTriple(), "z1", false, false), formal(primReal(), "t", false, false)); #line 145 "runpath3d.in" addFunc(ve, run::gen_runpath3d16, primPath3(), "reverse", formal(primPath3(), "p", false, false)); #line 150 "runpath3d.in" addFunc(ve, run::gen_runpath3d17, primPath3(), "subpath", formal(primPath3(), "p", false, false), formal(primInt(), "a", false, false), formal(primInt(), "b", false, false)); #line 155 "runpath3d.in" addFunc(ve, run::gen_runpath3d18, primPath3(), "subpath", formal(primPath3(), "p", false, false), formal(primReal(), "a", false, false), formal(primReal(), "b", false, false)); #line 160 "runpath3d.in" addFunc(ve, run::gen_runpath3d19, primInt(), "length", formal(primPath3(), "p", false, false)); #line 165 "runpath3d.in" addFunc(ve, run::gen_runpath3d20, primBoolean(), "cyclic", formal(primPath3(), "p", false, false)); #line 170 "runpath3d.in" addFunc(ve, run::gen_runpath3d21, primBoolean(), "straight", formal(primPath3(), "p", false, false), formal(primInt(), "t", false, false)); #line 175 "runpath3d.in" addFunc(ve, run::gen_runpath3d22, primPath3(), "unstraighten", formal(primPath3(), "p", false, false)); #line 180 "runpath3d.in" addFunc(ve, run::gen_runpath3d23, primReal(), "straightness", formal(primPath3(), "p", false, false), formal(primInt(), "t", false, false)); #line 191 "runpath3d.in" addFunc(ve, run::gen_runpath3d24, primReal(), "straightness", formal(primTriple(), "z0", false, false), formal(primTriple(), "c0", false, false), formal(primTriple(), "c1", false, false), formal(primTriple(), "z1", false, false)); #line 199 "runpath3d.in" addFunc(ve, run::gen_runpath3d25, primBoolean(), "piecewisestraight", formal(primPath3(), "p", false, false)); #line 204 "runpath3d.in" addFunc(ve, run::gen_runpath3d26, primReal(), "arclength", formal(primPath3(), "p", false, false)); #line 209 "runpath3d.in" addFunc(ve, run::gen_runpath3d27, primReal(), "arctime", formal(primPath3(), "p", false, false), formal(primReal(), "dval", false, false)); #line 214 "runpath3d.in" addFunc(ve, run::gen_runpath3d28, realArray(), "intersect", formal(primPath3(), "p", false, false), formal(primPath3(), "q", false, false), formal(primReal(), "fuzz", true, false)); #line 232 "runpath3d.in" addFunc(ve, run::gen_runpath3d29, realArray2(), "intersections", formal(primPath3(), "p", false, false), formal(primPath3(), "q", false, false), formal(primReal(), "fuzz", true, false)); #line 265 "runpath3d.in" addFunc(ve, run::gen_runpath3d30, realArray2(), "intersections", formal(primPath3(), "p", false, false), formal(tripleArray2(), "p", false, false), formal(primReal(), "fuzz", true, false)); #line 285 "runpath3d.in" addFunc(ve, run::gen_runpath3d31, primInt(), "size", formal(primPath3(), "p", false, false)); #line 290 "runpath3d.in" addFunc(ve, run::gen_runpath3d32, primPath3(), "&", formal(primPath3(), "p", false, false), formal(primPath3(), "q", false, false)); #line 295 "runpath3d.in" addFunc(ve, run::gen_runpath3d33, primTriple(), "min", formal(primPath3(), "p", false, false)); #line 300 "runpath3d.in" addFunc(ve, run::gen_runpath3d34, primTriple(), "max", formal(primPath3(), "p", false, false)); #line 305 "runpath3d.in" addFunc(ve, run::gen_runpath3d35, realArray(), "mintimes", formal(primPath3(), "p", false, false)); #line 315 "runpath3d.in" addFunc(ve, run::gen_runpath3d36, realArray(), "maxtimes", formal(primPath3(), "p", false, false)); #line 325 "runpath3d.in" addFunc(ve, run::gen_runpath3d37, primPath3(), "*", formal(realArray2(), "t", false, false), formal(primPath3(), "g", false, false)); #line 330 "runpath3d.in" addFunc(ve, run::gen_runpath3d38, primPair(), "minratio", formal(primPath3(), "g", false, false)); #line 335 "runpath3d.in" addFunc(ve, run::gen_runpath3d39, primPair(), "maxratio", formal(primPath3(), "g", false, false)); } } // namespace trans