--- trunk/test/tracto.diderot 2010/05/27 15:16:36 91 +++ trunk/test/tracto.diderot 2010/05/27 15:16:48 92 @@ -20,7 +20,7 @@ vec3 pos = seedPoint; //// principleEvec : tensor[3,3] -> vec3 //// * : real x vec3 -> vec3 - vec3 guide = seedSign * principleEvec(F@seedPoint); + vec3 guide = real(seedSign) * principleEvec(F@seedPoint); real aniso = 0.0; update @@ -31,7 +31,7 @@ //// CL : tensor[3,3] -> real aniso = CL(ten); // "CL" is a linear anisotropy measure //// < : real x real -> bool - if (aniso < clmin) { + if (aniso < CLmin) { // terminate because anisotropy went too low stabilize; } @@ -39,12 +39,12 @@ vec3 evec = principleEvec(ten); //// dot : vec3 x vec3 -> real //// - : vec3 -> vec3 - if (dot(evec, guide) < 0) { evec = -evec; } // fix eigenvector sign + if (dot(evec, guide) < 0.0) { evec = -evec; } // fix eigenvector sign //// * : real x real -> real //// * : real x vec3 -> vec3 //// + : vec3 x vec3 -> vec3 evec = principleEvec(F @ (pos + 0.5*stepSz*evec)); - if (dot(evec, guide) < 0) { evec = -evec; } // fix eigenvector sign + if (dot(evec, guide) < 0.0) { evec = -evec; } // fix eigenvector sign guide = evec; pos = pos + stepSz*evec; //// + : int x int -> int