| 324 |
|
if (normalize(np->v) == 0.0) { |
| 325 |
|
if (fabs(np->u_alpha - np->v_alpha) > 0.001) |
| 326 |
|
objerror(np->mp, WARNING, "illegal orientation vector"); |
| 327 |
< |
getperpendicular(np->u, np->pnorm, 1); /* punting */ |
| 327 |
> |
getperpendicular(np->u, np->pnorm, 0); /* punting */ |
| 328 |
|
fcross(np->v, np->pnorm, np->u); |
| 329 |
|
np->u_alpha = np->v_alpha = sqrt( 0.5 * |
| 330 |
|
(np->u_alpha*np->u_alpha + np->v_alpha*np->v_alpha) ); |
| 361 |
|
nstarget = 1; |
| 362 |
|
} |
| 363 |
|
scolorblack(scol); |
| 364 |
< |
dimlist[ndims++] = (int)(size_t)np->mp; |
| 364 |
> |
dimlist[ndims_inc] = (int)(size_t)np->mp; |
| 365 |
|
maxiter = MAXITER*nstarget; |
| 366 |
|
for (nstaken = ntrials = 0; (nstaken < nstarget) & |
| 367 |
|
(ntrials < maxiter); ntrials++) { |
| 428 |
|
} else |
| 429 |
|
nstarget = 1; |
| 430 |
|
} |
| 431 |
< |
dimlist[ndims++] = (int)(size_t)np->mp; |
| 431 |
> |
dimlist[ndims_inc] = (int)(size_t)np->mp; |
| 432 |
|
maxiter = MAXITER*nstarget; |
| 433 |
|
for (nstaken = ntrials = 0; (nstaken < nstarget) & |
| 434 |
|
(ntrials < maxiter); ntrials++) { |
