20 |
|
#include "bsdfrep.h" |
21 |
|
|
22 |
|
#ifndef NEIGH_FACT2 |
23 |
< |
#define NEIGH_FACT2 15. /* empirical neighborhood distance weight */ |
23 |
> |
#define NEIGH_FACT2 0.1 /* empirical neighborhood distance weight */ |
24 |
|
#endif |
25 |
|
/* number of processes to run */ |
26 |
|
int nprocs = 1; |
165 |
|
double d; |
166 |
|
int p[4]; |
167 |
|
int i, j; |
168 |
– |
|
169 |
– |
if ((x0 == x1) & (y0 == y1)) |
170 |
– |
return(0.); |
168 |
|
/* check radius */ |
169 |
|
p[0] = x0; p[1] = y0; p[2] = x1; p[3] = y1; |
170 |
|
for (i = 4; i--; ) { |
559 |
|
double bestd; |
560 |
|
int n; |
561 |
|
|
562 |
+ |
|
563 |
|
if (dsf_list == NULL) |
564 |
|
return; /* XXX should be error? */ |
565 |
|
near_rbf = dsf_list; |
604 |
|
memcpy(mir_rbf, near_rbf, n); |
605 |
|
mir_rbf->ord = near_rbf->ord - 1; /* not used, I think */ |
606 |
|
mir_rbf->next = NULL; |
607 |
+ |
mir_rbf->ejl = NULL; |
608 |
|
rev_rbf_symmetry(mir_rbf, MIRROR_X|MIRROR_Y); |
609 |
|
nprocs = 1; /* compute migration matrix */ |
610 |
< |
if (mig_list != create_migration(mir_rbf, near_rbf)) |
610 |
> |
if (create_migration(mir_rbf, near_rbf) == NULL) |
611 |
|
exit(1); /* XXX should never happen! */ |
612 |
|
/* interpolate normal dist. */ |
613 |
|
rbf = e_advect_rbf(mig_list, norm_vec, 2*near_rbf->nrbf); |