AcouSTO  version 2.0

◆ evalRR0_0fs()

static void evalRR0_0fs ( double *  pts,
unsigned  dim,
double *  p,
const double *  c,
const double *  r 
)
static
375 {
376  unsigned i, j;
377 
378  for (i = 0; i < dim - 1; ++i) {
379  p[i] = c[i] - r[i];
380  for (j = i + 1; j < dim; ++j) {
381  p[j] = c[j] - r[j];
382  memcpy(pts, p, sizeof(double) * dim); pts += dim;
383  p[i] = c[i] + r[i];
384  memcpy(pts, p, sizeof(double) * dim); pts += dim;
385  p[j] = c[j] + r[j];
386  memcpy(pts, p, sizeof(double) * dim); pts += dim;
387  p[i] = c[i] - r[i];
388  memcpy(pts, p, sizeof(double) * dim); pts += dim;
389 
390  p[j] = c[j]; /* Done with j -> Restore p[j] */
391  }
392  p[i] = c[i]; /* Done with i -> Restore p[i] */
393  }
394 }