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) + |