478 |
|
double t, full_dist; |
479 |
|
/* get relative position */ |
480 |
|
t = Acos(DOT(invec, mig->rbfv[0]->invec)); |
481 |
< |
if (t < M_PI/grid_res) { /* near first DSF */ |
481 |
> |
if (t <= .001) { /* near first DSF */ |
482 |
|
n = sizeof(RBFNODE) + sizeof(RBFVAL)*(mig->rbfv[0]->nrbf-1); |
483 |
|
rbf = (RBFNODE *)malloc(n); |
484 |
|
if (rbf == NULL) |
488 |
|
return(rbf); |
489 |
|
} |
490 |
|
full_dist = acos(DOT(mig->rbfv[0]->invec, mig->rbfv[1]->invec)); |
491 |
< |
if (t > full_dist-M_PI/grid_res) { /* near second DSF */ |
491 |
> |
if (t >= full_dist-.001) { /* near second DSF */ |
492 |
|
n = sizeof(RBFNODE) + sizeof(RBFVAL)*(mig->rbfv[1]->nrbf-1); |
493 |
|
rbf = (RBFNODE *)malloc(n); |
494 |
|
if (rbf == NULL) |