51 |
|
if (i) spinvector(tvec, tvec, outvec, phinc); |
52 |
|
if (tvec[2] > 0 ^ output_orient > 0) |
53 |
|
continue; |
54 |
< |
sum += eval_rbfrep(rbf, tvec) * output_orient*tvec[2]; |
54 |
> |
sum += eval_rbfrep(rbf, tvec) * COSF(tvec[2]); |
55 |
|
++n; |
56 |
|
} |
57 |
|
if (n < 2) /* should never happen! */ |
64 |
|
est_DSFrad(const RBFNODE *rbf, const FVECT outvec) |
65 |
|
{ |
66 |
|
const double rad_epsilon = 0.03; |
67 |
< |
const double DSFtarget = 0.60653066 * eval_rbfrep(rbf,outvec) |
68 |
< |
* output_orient*outvec[2]; |
67 |
> |
const double DSFtarget = 0.60653066 * eval_rbfrep(rbf,outvec) * |
68 |
> |
COSF(outvec[2]); |
69 |
|
double inside_rad = rad_epsilon; |
70 |
|
double outside_rad = 0.5; |
71 |
|
double DSFinside = eval_DSFsurround(rbf, outvec, inside_rad); |