--- trunk/test/tracto.diderot 2010/04/26 00:57:58 56 +++ trunk/test/tracto.diderot 2010/05/03 03:46:16 58 @@ -18,26 +18,39 @@ { int stepNum = 0; vec3 pos = seedPoint; + //// principleEvec : tensor[3,3] -> vec3 + //// * : real x vec3 -> vec3 vec3 guide = seedSign * principleEvec(F@seedPoint); real aniso = 0.0; update { if (inside (pos,F)) { + //// @ : field x vec -> tensor[3,3] tensor[3,3] ten = F@pos; + //// CL : tensor[3,3] -> real aniso = CL(ten); // "CL" is a linear anisotropy measure + //// < : real x real -> bool if (aniso < clmin) { // terminate because anisotropy went too low stabilize; } // else vec3 evec = principleEvec(ten); + //// dot : vec3 x vec3 -> real + //// - : vec3 -> vec3 if (dot(evec, guide) < 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 guide = evec; pos = pos + stepSz*evec; + //// + : int x int -> int stepNum = 1 + stepNum; + //// > : real x real -> bool + //// && : bool x bool -> bool if (stepNumMax > 0 && stepNum > stepNumMax) { // terminate because we took too many steps stabilize;
Click to toggle
does not end with </html> tag
does not end with </body> tag
The output has ended thus: tepNumMax) { // terminate because we took too many steps stabilize;