ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/JSOC/proj/sharp/apps/errorprop.c
(Generate patch)

Comparing proj/sharp/apps/errorprop.c (file contents):
Revision 1.1 by xudong, Thu Aug 23 08:38:33 2012 UTC vs.
Revision 1.2 by xudong, Sat Feb 2 04:53:58 2019 UTC

# Line 32 | Line 32 | int errorprop (float *bTotal, float *bAz
32    // Xudong Oct 18 2011: Fill factor variances covariances are all NaNs. removed
33    // NNB interpolation, just 1 point
34    // Fixed definition of azi for all derivatives
35 <
35 >  // Feb 1 2019, dBt/dXX is actually dBy/dXX, so the sign is wrong; however they are used in pairs (squared) so final result was okay
36 >
37    static double raddeg = M_PI / 180.;
38    double b, inc, azim;
39    double a11, a12, a13, a21, a22, a23, a31, a32, a33;
# Line 90 | Line 91 | int errorprop (float *bTotal, float *bAz
91    dBpdInc = b * (- a11 * cos(inc) * sin(azim) + a12 * cos(inc) * cos(azim) - a13 * sin(inc));
92    dBpdAzim = b * (- a11 * sin(inc) * cos(azim) - a12 * sin(inc) * sin(azim));
93    
94 <  dBtdBtotal = (- a21 * sin(inc) * sin(azim) + a22 * sin(inc) * cos(azim) + a23 * cos(inc));
95 <  dBtdInc = b * (- a21 * cos(inc) * sin(azim) + a22 * cos(inc) * cos(azim) - a23 * sin(inc));
96 <  dBtdAzim = b * (- a21 * sin(inc) * cos(azim) - a22 * sin(inc) * sin(azim));
94 >  dBtdBtotal = (- a21 * sin(inc) * sin(azim) + a22 * sin(inc) * cos(azim) + a23 * cos(inc)) * (-1);
95 >  dBtdInc = b * (- a21 * cos(inc) * sin(azim) + a22 * cos(inc) * cos(azim) - a23 * sin(inc)) * (-1);
96 >  dBtdAzim = b * (- a21 * sin(inc) * cos(azim) - a22 * sin(inc) * sin(azim)) * (-1);
97    
98    dBrdBtotal = (- a31 * sin(inc) * sin(azim) + a32 * sin(inc) * cos(azim) + a33 * cos(inc));
99    dBrdInc = b * (- a31 * cos(inc) * sin(azim) + a32 * cos(inc) * cos(azim) - a33 * sin(inc));

Diff Legend

– Removed lines
+ Added lines
< Changed lines
> Changed lines