Mercurial > hg > fxanalyse

comparison FXAnalyse.c @ 221:c14f5a95850d

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Fix N3 measurement
author Daniele Nicolodi <daniele.nicolodi@obspm.fr>
date Thu, 18 Sep 2014 18:01:56 +0200
parents 28a4ff70fe48
children db051dbf8d6e
comparison
equal deleted inserted replaced
220:28a4ff70fe48 221:c14f5a95850d
1212 f_rep_minus = f_rep_minus / nobs; 1212 f_rep_minus = f_rep_minus / nobs;
1213 f_beat_minus = f_beat_minus / nobs; 1213 f_beat_minus = f_beat_minus / nobs;
1214 nobs = 0; 1214 nobs = 0;
1215 1215
1216 // check delta frep 1216 // check delta frep
1217 double delta_f_rep_m = f_rep_plus - f_rep_minus; 1217 double delta_f_rep_m = f_rep_minus - f_rep_plus;
1218 double delta_f_rep = Sign1 * Ndiv * 2.0 * delta_f_lock_3 / N1; 1218 double delta_f_rep = Sign1 * Ndiv * 2.0 * delta_f_lock_3 / N1;
1219 double delta = delta_f_rep_m - delta_f_rep; 1219 double delta = delta_f_rep_m - delta_f_rep;
1220 1220
1221 logmsg("delta frep: measured=%e expected=%e difference=%e rel=%e", 1221 logmsg("delta frep: measured=%e expected=%e difference=%e rel=%e",
1222 delta_f_rep_m, delta_f_rep, delta, delta / delta_f_rep); 1222 delta_f_rep_m, delta_f_rep, delta, delta / delta_f_rep);
1223 1223
1224 // compute N3 1224 // compute N3
1225 double delta_f_beat = f_beat_plus - f_beat_minus + 2.0 * Sign1 * Sign3 * N3/N1 * Ndiv * delta_f_lock_3; 1225 double delta_f_beat = f_beat_minus - f_beat_plus + 2.0 * Sign1 * Sign3 * N3/N1 * Ndiv * delta_f_lock_3;
1226 double delta_f_beat_expected = delta_f_rep * N3; 1226 double delta_f_beat_expected = delta_f_rep * N3;
1227 1227
1228 logmsg("delta fbeat: measured=%e expected=%e difference=%e", 1228 logmsg("delta fbeat: measured=%e expected=%e difference=%e",
1229 delta_f_beat, delta_f_beat_expected, delta_f_beat - delta_f_beat_expected); 1229 delta_f_beat, delta_f_beat_expected, delta_f_beat - delta_f_beat_expected);
1230 1230