/***** 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" #include "opsymbols.h" using vm::stack; using vm::error; using vm::array; using vm::read; 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 *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: #ifndef NOSYM #include "runpath3d.symbols.h" #endif 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; copyArray2C(A,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 286 "runpath3d.in" // Int size(path3 p); void gen_runpath3d31(stack *Stack) { path3 p=vm::pop(Stack); #line 287 "runpath3d.in" {Stack->push(p.size()); return;} } #line 291 "runpath3d.in" // path3 &(path3 p, path3 q); void gen_runpath3d32(stack *Stack) { path3 q=vm::pop(Stack); path3 p=vm::pop(Stack); #line 292 "runpath3d.in" {Stack->push(camp::concat(p,q)); return;} } #line 296 "runpath3d.in" // triple min(path3 p); void gen_runpath3d33(stack *Stack) { path3 p=vm::pop(Stack); #line 297 "runpath3d.in" {Stack->push(p.min()); return;} } #line 301 "runpath3d.in" // triple max(path3 p); void gen_runpath3d34(stack *Stack) { path3 p=vm::pop(Stack); #line 302 "runpath3d.in" {Stack->push(p.max()); return;} } #line 306 "runpath3d.in" // realarray* mintimes(path3 p); void gen_runpath3d35(stack *Stack) { path3 p=vm::pop(Stack); #line 307 "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 316 "runpath3d.in" // realarray* maxtimes(path3 p); void gen_runpath3d36(stack *Stack) { path3 p=vm::pop(Stack); #line 317 "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 326 "runpath3d.in" // path3 *(realarray2 *t, path3 g); void gen_runpath3d37(stack *Stack) { path3 g=vm::pop(Stack); realarray2 * t=vm::pop(Stack); #line 327 "runpath3d.in" {Stack->push(transformed(*t,g)); return;} } #line 331 "runpath3d.in" // pair minratio(path3 g); void gen_runpath3d38(stack *Stack) { path3 g=vm::pop(Stack); #line 332 "runpath3d.in" {Stack->push(g.ratio(::min)); return;} } #line 336 "runpath3d.in" // pair maxratio(path3 g); void gen_runpath3d39(stack *Stack) { path3 g=vm::pop(Stack); #line 337 "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(), SYM(path3), formal(tripleArray(), SYM(pre), false, false), formal(tripleArray(), SYM(point), false, false), formal(tripleArray(), SYM(post), false, false), formal(booleanArray(), SYM(straight), false, false), formal(primBoolean(), SYM(cyclic), false, false)); #line 56 "runpath3d.in" REGISTER_BLTIN(run::nullPath3,"nullPath3"); #line 61 "runpath3d.in" addFunc(ve, run::gen_runpath3d2, primBoolean(), SYM_EQ, formal(primPath3(), SYM(a), false, false), formal(primPath3(), SYM(b), false, false)); #line 66 "runpath3d.in" addFunc(ve, run::gen_runpath3d3, primBoolean(), SYM_NEQ, formal(primPath3(), SYM(a), false, false), formal(primPath3(), SYM(b), false, false)); #line 71 "runpath3d.in" addFunc(ve, run::gen_runpath3d4, primTriple(), SYM(point), formal(primPath3(), SYM(p), false, false), formal(primInt(), SYM(t), false, false)); #line 76 "runpath3d.in" addFunc(ve, run::gen_runpath3d5, primTriple(), SYM(point), formal(primPath3(), SYM(p), false, false), formal(primReal(), SYM(t), false, false)); #line 81 "runpath3d.in" addFunc(ve, run::gen_runpath3d6, primTriple(), SYM(precontrol), formal(primPath3(), SYM(p), false, false), formal(primInt(), SYM(t), false, false)); #line 86 "runpath3d.in" addFunc(ve, run::gen_runpath3d7, primTriple(), SYM(precontrol), formal(primPath3(), SYM(p), false, false), formal(primReal(), SYM(t), false, false)); #line 91 "runpath3d.in" addFunc(ve, run::gen_runpath3d8, primTriple(), SYM(postcontrol), formal(primPath3(), SYM(p), false, false), formal(primInt(), SYM(t), false, false)); #line 96 "runpath3d.in" addFunc(ve, run::gen_runpath3d9, primTriple(), SYM(postcontrol), formal(primPath3(), SYM(p), false, false), formal(primReal(), SYM(t), false, false)); #line 101 "runpath3d.in" addFunc(ve, run::gen_runpath3d10, primTriple(), SYM(dir), formal(primPath3(), SYM(p), false, false), formal(primInt(), SYM(t), false, false), formal(primInt(), SYM(sign), true, false), formal(primBoolean(), SYM(normalize), true, false)); #line 106 "runpath3d.in" addFunc(ve, run::gen_runpath3d11, primTriple(), SYM(dir), formal(primPath3(), SYM(p), false, false), formal(primReal(), SYM(t), false, false), formal(primBoolean(), SYM(normalize), true, false)); #line 111 "runpath3d.in" addFunc(ve, run::gen_runpath3d12, primTriple(), SYM(accel), formal(primPath3(), SYM(p), false, false), formal(primInt(), SYM(t), false, false), formal(primInt(), SYM(sign), true, false)); #line 116 "runpath3d.in" addFunc(ve, run::gen_runpath3d13, primTriple(), SYM(accel), formal(primPath3(), SYM(p), false, false), formal(primReal(), SYM(t), false, false)); #line 121 "runpath3d.in" addFunc(ve, run::gen_runpath3d14, primReal(), SYM(radius), formal(primPath3(), SYM(p), false, false), formal(primReal(), SYM(t), false, false)); #line 133 "runpath3d.in" addFunc(ve, run::gen_runpath3d15, primReal(), SYM(radius), formal(primTriple(), SYM(z0), false, false), formal(primTriple(), SYM(c0), false, false), formal(primTriple(), SYM(c1), false, false), formal(primTriple(), SYM(z1), false, false), formal(primReal(), SYM(t), false, false)); #line 145 "runpath3d.in" addFunc(ve, run::gen_runpath3d16, primPath3(), SYM(reverse), formal(primPath3(), SYM(p), false, false)); #line 150 "runpath3d.in" addFunc(ve, run::gen_runpath3d17, primPath3(), SYM(subpath), formal(primPath3(), SYM(p), false, false), formal(primInt(), SYM(a), false, false), formal(primInt(), SYM(b), false, false)); #line 155 "runpath3d.in" addFunc(ve, run::gen_runpath3d18, primPath3(), SYM(subpath), formal(primPath3(), SYM(p), false, false), formal(primReal(), SYM(a), false, false), formal(primReal(), SYM(b), false, false)); #line 160 "runpath3d.in" addFunc(ve, run::gen_runpath3d19, primInt(), SYM(length), formal(primPath3(), SYM(p), false, false)); #line 165 "runpath3d.in" addFunc(ve, run::gen_runpath3d20, primBoolean(), SYM(cyclic), formal(primPath3(), SYM(p), false, false)); #line 170 "runpath3d.in" addFunc(ve, run::gen_runpath3d21, primBoolean(), SYM(straight), formal(primPath3(), SYM(p), false, false), formal(primInt(), SYM(t), false, false)); #line 175 "runpath3d.in" addFunc(ve, run::gen_runpath3d22, primPath3(), SYM(unstraighten), formal(primPath3(), SYM(p), false, false)); #line 180 "runpath3d.in" addFunc(ve, run::gen_runpath3d23, primReal(), SYM(straightness), formal(primPath3(), SYM(p), false, false), formal(primInt(), SYM(t), false, false)); #line 191 "runpath3d.in" addFunc(ve, run::gen_runpath3d24, primReal(), SYM(straightness), formal(primTriple(), SYM(z0), false, false), formal(primTriple(), SYM(c0), false, false), formal(primTriple(), SYM(c1), false, false), formal(primTriple(), SYM(z1), false, false)); #line 199 "runpath3d.in" addFunc(ve, run::gen_runpath3d25, primBoolean(), SYM(piecewisestraight), formal(primPath3(), SYM(p), false, false)); #line 204 "runpath3d.in" addFunc(ve, run::gen_runpath3d26, primReal(), SYM(arclength), formal(primPath3(), SYM(p), false, false)); #line 209 "runpath3d.in" addFunc(ve, run::gen_runpath3d27, primReal(), SYM(arctime), formal(primPath3(), SYM(p), false, false), formal(primReal(), SYM(dval), false, false)); #line 214 "runpath3d.in" addFunc(ve, run::gen_runpath3d28, realArray(), SYM(intersect), formal(primPath3(), SYM(p), false, false), formal(primPath3(), SYM(q), false, false), formal(primReal(), SYM(fuzz), true, false)); #line 232 "runpath3d.in" addFunc(ve, run::gen_runpath3d29, realArray2(), SYM(intersections), formal(primPath3(), SYM(p), false, false), formal(primPath3(), SYM(q), false, false), formal(primReal(), SYM(fuzz), true, false)); #line 265 "runpath3d.in" addFunc(ve, run::gen_runpath3d30, realArray2(), SYM(intersections), formal(primPath3(), SYM(p), false, false), formal(tripleArray2(), SYM(p), false, false), formal(primReal(), SYM(fuzz), true, false)); #line 286 "runpath3d.in" addFunc(ve, run::gen_runpath3d31, primInt(), SYM(size), formal(primPath3(), SYM(p), false, false)); #line 291 "runpath3d.in" addFunc(ve, run::gen_runpath3d32, primPath3(), SYM_AMPERSAND, formal(primPath3(), SYM(p), false, false), formal(primPath3(), SYM(q), false, false)); #line 296 "runpath3d.in" addFunc(ve, run::gen_runpath3d33, primTriple(), SYM(min), formal(primPath3(), SYM(p), false, false)); #line 301 "runpath3d.in" addFunc(ve, run::gen_runpath3d34, primTriple(), SYM(max), formal(primPath3(), SYM(p), false, false)); #line 306 "runpath3d.in" addFunc(ve, run::gen_runpath3d35, realArray(), SYM(mintimes), formal(primPath3(), SYM(p), false, false)); #line 316 "runpath3d.in" addFunc(ve, run::gen_runpath3d36, realArray(), SYM(maxtimes), formal(primPath3(), SYM(p), false, false)); #line 326 "runpath3d.in" addFunc(ve, run::gen_runpath3d37, primPath3(), SYM_TIMES, formal(realArray2(), SYM(t), false, false), formal(primPath3(), SYM(g), false, false)); #line 331 "runpath3d.in" addFunc(ve, run::gen_runpath3d38, primPair(), SYM(minratio), formal(primPath3(), SYM(g), false, false)); #line 336 "runpath3d.in" addFunc(ve, run::gen_runpath3d39, primPair(), SYM(maxratio), formal(primPath3(), SYM(g), false, false)); } } // namespace trans