diff options
author | Norbert Preining <norbert@preining.info> | 2020-08-07 03:02:55 +0000 |
---|---|---|
committer | Norbert Preining <norbert@preining.info> | 2020-08-07 03:02:55 +0000 |
commit | 3710c23789342d9d88783cbe9ad23eb5263a7c33 (patch) | |
tree | 1c40b2bfb9f6b652f7877ad0991f5d6b44391154 /graphics/asymptote/knot.cc | |
parent | d99fea72e367e8d4d2809b43b3b0206b1bb526b8 (diff) |
CTAN sync 202008070302
Diffstat (limited to 'graphics/asymptote/knot.cc')
-rw-r--r-- | graphics/asymptote/knot.cc | 12 |
1 files changed, 6 insertions, 6 deletions
diff --git a/graphics/asymptote/knot.cc b/graphics/asymptote/knot.cc index 199be8249b..09ea2fe3b2 100644 --- a/graphics/asymptote/knot.cc +++ b/graphics/asymptote/knot.cc @@ -71,7 +71,7 @@ double velocity(double theta, double phi, tension t) double denom = t.val * (3.0 + c*ct + d*cf); - double r = denom != 0.0 ? (2.0 + a*(st - b*sf)*(sf - b*st)*(ct-cf)) / denom + double r = denom != 0.0 ? (2.0 + a*(st - b*sf)*(sf - b*st)*(ct-cf)) / denom : VELOCITY_BOUND; //cerr << " velocity(" << theta << "," << phi <<")= " << r << endl; @@ -282,7 +282,7 @@ struct eqnprop : public knotprop<eqn> { // Ax=0), there is no need to solve for theta; we can just use zeros for the // thetas. In fact, our general solving method may not work in this case. // A common example of this is -// +// // a{curl 1}..{curl 1}b // // which arises when solving a one-length path a..b or in a larger path a @@ -320,7 +320,7 @@ weqn scale(weqn q) { /* Recalculate the equations in the form: * theta[j] + post * theta[j+1] = aug + w * theta[0] - * + * * Used as the first step in solve cyclic equations. */ cvector<weqn> recalc(cvector<eqn>& e) @@ -375,7 +375,7 @@ double solveForTheta0(cvector<weqn>& we) } // After the iteration we have - // + // // theta[n] = a + b*theta[0] + c*theta[n] // // where theta[n]=theta[0], so @@ -601,7 +601,7 @@ void encodeStraight(protopath& p, Int k, knotlist& l) pair b=l.back().z; double bt=l.back().tin.val; pair step=(b-a)/3.0; - + if (at==1.0 && bt==1.0) { p.straight(k)=true; p.post(k)=a+step; @@ -804,7 +804,7 @@ path solveSimple(cvector<pair>& z) { // The two specifiers used: an open spec and a curl spec for the ends. spec open; - + // curlSpec curl; // curlSpec curly(2.0); // dirSpec E(0); |