QPMS
/
qpms
1
0
Fork 0
Code Issues Pull Requests Projects Releases Wiki Activity

Fix "zero" ray in plane wave test 

Former-commit-id: cd30bf7acb9c68bc325fd60524a3a21df26629bc
This commit is contained in:
Marek Nečada 2018-05-01 20:13:39 +03:00
parent 7eefd445e0
commit 0a7da94b4f
1 changed files with 3 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
qpms/tests/ray_planewaves.c
View File

@ -1,4 +1,4 @@
// c99 -o raypwtest1 -I .. ray_planewaves.c -lgsl -lm ../legendre.c ../vswf.c ../bessel.c -lblas // c99 -o ../../tests/pwraytests/raypwtest1 -I .. ray_planewaves.c -lgsl -lm ../legendre.c ../vswf.c ../bessel.c -lblas
#include <gsl/gsl_rng.h> #include <gsl/gsl_rng.h>
#include <gsl/gsl_math.h> #include <gsl/gsl_math.h>
#include <gsl/gsl_randist.h> #include <gsl/gsl_randist.h>
@ -57,9 +57,9 @@ int main(int argc, char **argv) {
cart3_t rays[nrays]; cart3_t rays[nrays];
double relerrs[nrays]; double relerrs[nrays];
memset(rays, 0, nrays * sizeof(cart3_t)); memset(rays, 0, nrays * sizeof(cart3_t));
rays[1].x = rays[2].y = rays[3].z = sigma; rays[0].x = rays[1].y = rays[2].z = sigma;
double relerrthreshold = 1e-11; double relerrthreshold = 1e-11;
for (unsigned i = 4; i < nrays; ++i) { for (unsigned i = 3; i < nrays; ++i) {
cart3_t *w = rays+i; cart3_t *w = rays+i;
w->x = gsl_ran_gaussian(rng, sigma); w->x = gsl_ran_gaussian(rng, sigma);
w->y = gsl_ran_gaussian(rng, sigma); w->y = gsl_ran_gaussian(rng, sigma);