| 553 |
|
static void |
| 554 |
|
check_normal_incidence(void) |
| 555 |
|
{ |
| 556 |
< |
static const FVECT norm_vec = {.0, .0, 1.}; |
| 556 |
> |
static FVECT norm_vec = {.0, .0, 1.}; |
| 557 |
|
const int saved_nprocs = nprocs; |
| 558 |
|
RBFNODE *near_rbf, *mir_rbf, *rbf; |
| 559 |
|
double bestd; |
| 560 |
|
int n; |
| 561 |
|
|
| 562 |
– |
|
| 562 |
|
if (dsf_list == NULL) |
| 563 |
|
return; /* XXX should be error? */ |
| 564 |
|
near_rbf = dsf_list; |
| 608 |
|
nprocs = 1; /* compute migration matrix */ |
| 609 |
|
if (create_migration(mir_rbf, near_rbf) == NULL) |
| 610 |
|
exit(1); /* XXX should never happen! */ |
| 611 |
< |
/* interpolate normal dist. */ |
| 611 |
> |
norm_vec[2] = input_orient; /* interpolate normal dist. */ |
| 612 |
|
rbf = e_advect_rbf(mig_list, norm_vec, 2*near_rbf->nrbf); |
| 613 |
|
nprocs = saved_nprocs; /* final clean-up */ |
| 614 |
|
free(mir_rbf); |