This had been giving me jokes for two months. Had been testing mpf1t3xy2 function all that time. Finally noticed that I was calling the function wrongly at 2am last night.
bembo:tully$ hg diff -r 608 tfall.pro
diff -r 2bfbd2cb20a2 tfall.pro
--- a/tfall.pro Wed Jun 10 00:41:50 2009 +0100
+++ b/tfall.pro Thu Jun 11 15:12:22 2009 +0100
@@ -7,15 +7,6 @@
sp = ['NGC_3390', 'PGC_44931', 'ESO_443G42', 'NGC_5746', 'IC_4767', $
'NGC_6722', 'ESO_185G53', 'IC_4937', 'IC_5096', 'ESO_240G11', $
'NGC_1886', 'NGC_4703', 'NGC_7123', 'IC_5176']
-
- ; Old galaxy order (FIX: this gives different answers for log mass fits,
- ; and possibly others. This may be a bug in mpf1t3xy)
- ;s0 = ['NGC_0128', 'NGC_1381', 'NGC_1596', 'NGC_2310', 'ESO_311G12', $
- ; 'NGC_3203', 'NGC_4469', 'NGC_4710', 'NGC_6771', 'ESO_597G36', $
- ; 'NGC_1032', 'NGC_3957', 'NGC_5084', 'ESO_151G04']
- ;sp = ['ESO_185G53', 'ESO_240G11', 'ESO_443G42', 'IC_4767', 'IC_4937', $
- ; 'IC_5096', 'IC_5176', 'NGC_1886', 'NGC_3390', 'NGC_4703', $
- ; 'NGC_5746', 'NGC_6722', 'NGC_7123', 'PGC_44931']
; Fix distance and apparent magnitude errors
fixderror = 0.0
@@ -113,32 +104,32 @@
print, "log v_obs -- log mass (stellar)"
logv_logmass = mpf1t3xy2(logvsp, logmasssp, logvs0, logmasss0, $
- e_logvsp, e_logmasssp, e_logvs0, e_logmasssp, /quiet, x0 = x0, $
+ e_logvsp, e_logmasssp, e_logvs0, e_logmasss0, /quiet, x0 = x0, $
guess = vmassguess, reduce = reduce, /latex)
print, "log v_drift -- log mass (stellar)"
logvdrift_logmass = mpf1t3xy2(logvdriftsp, logmasssp, logvdrifts0, $
- logmasss0, e_logvdriftsp, e_logmasssp, e_logvdrifts0, e_logmasssp, $
+ logmasss0, e_logvdriftsp, e_logmasssp, e_logvdrifts0, e_logmasss0, $
/quiet, x0 = x0, guess = vmassguess, reduce = reduce, /latex)
print, "log v_c -- log mass (stellar)"
logvc_logmass = mpf1t3xy2(logvcsp, logmasssp, logvcs0, logmasss0, $
- e_logvcsp, e_logmasssp, e_logvcs0, e_logmasssp, /quiet, x0 = x0, $
+ e_logvcsp, e_logmasssp, e_logvcs0, e_logmasss0, /quiet, x0 = x0, $
guess = vmassguess, reduce = reduce, /latex)
print, "log v_obs -- log mass (dynamical)"
logv_logmassdyn = mpf1t3xy2(logvsp, logmassspdyn, logvs0, logmasss0dyn, $
- e_logvsp, e_logmassspdyn, e_logvs0, e_logmassspdyn, /quiet, x0 = x0, $
+ e_logvsp, e_logmassspdyn, e_logvs0, e_logmasss0dyn, /quiet, x0 = x0, $
guess = vmassguess, reduce = reduce, /latex)
print, "log v_drift -- log mass (dynamical)"
logvdrift_logmassdyn = mpf1t3xy2(logvdriftsp, logmassspdyn, logvdrifts0, $
logmasss0dyn, e_logvdriftsp, e_logmassspdyn, e_logvdrifts0, $
- e_logmassspdyn, /quiet, x0 = x0, guess = vmassguess, reduce = reduce, /latex)
+ e_logmasss0dyn, /quiet, x0 = x0, guess = vmassguess, reduce = reduce, /latex)
print, "log v_c -- log mass (dynamical)"
logvc_logmassdyn = mpf1t3xy2(logvcsp, logmassspdyn, logvcs0, logmasss0dyn, $
- e_logvcsp, e_logmassspdyn, e_logvcs0, e_logmassspdyn, /quiet, x0 = x0, $
+ e_logvcsp, e_logmassspdyn, e_logvcs0, e_logmasss0dyn, /quiet, x0 = x0, $
guess = vmassguess, reduce = reduce, /latex)
― caek, Thursday, 11 June 2009 14:14 (fourteen years ago) link
one month passes...
two months pass...
two months pass...
four weeks pass...
one month passes...
three weeks pass...
three weeks pass...