diff options
author | ultramn <dchapm2@umbc.edu> | 2017-03-09 17:40:19 -0800 |
---|---|---|
committer | ultramn <dchapm2@umbc.edu> | 2017-03-09 17:40:19 -0800 |
commit | a956d29b35aa9c8b7cfd869683afb382da65b3ea (patch) | |
tree | a39afb19d7278007fa8bfff9458eb1b9697c7df5 | |
parent | 3e39224bdd65b8a49de44f007add854dc626a209 (diff) | |
download | libsurvive-a956d29b35aa9c8b7cfd869683afb382da65b3ea.tar.gz libsurvive-a956d29b35aa9c8b7cfd869683afb382da65b3ea.tar.bz2 |
Added a loop to check the z values
-rw-r--r-- | dave/AffineSolve.c | 4 |
1 files changed, 3 insertions, 1 deletions
diff --git a/dave/AffineSolve.c b/dave/AffineSolve.c index 1c1d0f1..d036858 100644 --- a/dave/AffineSolve.c +++ b/dave/AffineSolve.c @@ -261,6 +261,7 @@ printf("rhat %f %f (len %f)\n", rhat[0][0], rhat[1][0], rhat_len); //------------------------- FLOAT ydist2; FLOAT bestBestErr = 9999.0; + FLOAT bestYdist = 0; for (ydist2=ydist-0.1; ydist2<ydist+0.1; ydist2+=0.0001) { FLOAT x2[3][1] = { {M[0][0]*ydist2}, {0.0}, {M[1][0]*ydist2} }; @@ -309,9 +310,10 @@ printf("rhat %f %f (len %f)\n", rhat[0][0], rhat[1][0], rhat_len); memcpy(y,y2,3*sizeof(FLOAT)); memcpy(z,z2,3*sizeof(FLOAT)); bestBestErr = bestErr; + bestYdist = ydist2; } } - + ydist = bestYdist; /* for (i=0; i<nPoints; i++) { |