| 82 |
|
int |
| 83 |
|
use_symmetry(FVECT vec) |
| 84 |
|
{ |
| 85 |
< |
double phi = get_phi360(vec); |
| 85 |
> |
const double phi = get_phi360(vec); |
| 86 |
|
|
| 87 |
|
switch (inp_coverage) { |
| 88 |
|
case INP_QUAD1|INP_QUAD2|INP_QUAD3|INP_QUAD4: |
| 189 |
|
int pos[2]; |
| 190 |
|
int n; |
| 191 |
|
|
| 192 |
< |
for (n = rbf->nrbf; n-- > 0; ) { |
| 192 |
> |
for (n = ((-.01 > phi) | (phi > .01))*rbf->nrbf; n-- > 0; ) { |
| 193 |
|
ovec_from_pos(outvec, rbf->rbfa[n].gx, rbf->rbfa[n].gy); |
| 194 |
|
spinvector(outvec, outvec, vnorm, phi); |
| 195 |
|
pos_from_vec(pos, outvec); |
| 485 |
|
rbfh.invec[0] = getflt(ifp); |
| 486 |
|
rbfh.invec[1] = getflt(ifp); |
| 487 |
|
rbfh.invec[2] = getflt(ifp); |
| 488 |
+ |
if (normalize(rbfh.invec) == 0) { |
| 489 |
+ |
fprintf(stderr, "%s: zero incident vector\n", progname); |
| 490 |
+ |
return(0); |
| 491 |
+ |
} |
| 492 |
|
rbfh.vtotal = getflt(ifp); |
| 493 |
|
rbfh.nrbf = getint(4, ifp); |
| 494 |
|
newrbf = (RBFNODE *)malloc(sizeof(RBFNODE) + |