addstar doesn't load dM, dRA, dDec from GPC1 pt 2 .cmf files

[Michael Wood-Vasey]

I ran addstar on the 2008/10/06 GPC1 pt 2 release .cmf files and successfully put all of the 8 FPAs and all their chips into a DVO database. But the database did not have dM, dRA, or dDec loaded in the measures table so relphot and relastro both fail. I'm trying to trace through how dM gets loaded. Things still seem to be fine as of the LoadST

---
[wwoodvas@serenity gpc1] ls
gpc1
o4732g0168o.34857/
o4732g0169o.34858/
o4732g0170o.34859/
o4732g0171o.34860/
o4732g0172o.34861/
o4732g0179o.34868/
o4732g0190o.34879/
o4732g0212o.34901/
[wwoodvas@serenity gpc1] # Clear the directory and start again
[wwoodvas@serenity gpc1] rm -rf gpc1
[wwoodvas@serenity gpc1] foreach cmf ( */*.cmf )
foreach? addstar -accept-astrom -D CATDIR gpc1 -D CAMERA gpc1 $cmf
foreach? end
no data in db gpc1/Photcodes.dat
directory gpc1 does not exist, creating...
SUCCESS: elapsed time 0.7603 sec for 320 stars ( 0 matches), 320 average, 320 measure
SUCCESS: elapsed time 0.3961 sec for 448 stars ( 0 matches), 768 average, 768 measure
SUCCESS: elapsed time 0.3962 sec for 379 stars ( 0 matches), 1147 average, 1147 measure
SUCCESS: elapsed time 0.3992 sec for 359 stars ( 0 matches), 1506 average, 1506 measure
[...]
SUCCESS: elapsed time 0.6230 sec for 160 stars ( 22 matches), 16016 average, 16452 measure
SUCCESS: elapsed time 0.6211 sec for 69 stars ( 2 matches), 16083 average, 16521 measure

---

Then we can look at things in DVO:

---

[wwoodvas@serenity gpc1] dvo -D CATDIR gpc1
no file /Users/wwoodvas/.dvorc

Welcome to DVO - desktop virtual observatory

dvo: region 360.5 0.25 1.5
dvo: images
current style (0): -x 2 -c black -pt 2 -lt 0 -lw 0.500000 -sz 1.000000
reading images (mode DVO_MODE_MEF)
plotting 477 images
dvo: cgrid
current style (0): -x 2 -c black -pt 2 -lt 0 -lw 0.500000 -sz 1.000000
dvo: mextract -region ra, dec, i, i:err
dvo: czplot ra dec i 17 12 -c red -pt 7

---

And here's where we see the first problem. 'i:err' is NaN for every value:

---

dvo: limits i i:err
dvo: limits
limits: 6.623180 29.526000 nan nan [-a] [-n device]

---

Here we've asked the "limits" function to automatically figure out and set the limits for "i" and "i:err" (for example because we were going to plot them). It finds that the min and max values for "i:err" are NaN, i.e., all of the values of "i:err" are NaN. This holds true for all of the filters (griy). A bit of further investigation reveals that all of the ra:err and dec:err values are 0, but I don't think the inputs to these values from the .cmf value are reasonable in the first place since they tend to be in the 1000s or 0 (are the units on X_PSF_SIG thousandths of a pixel or something?).

I don't know what's going on here? "PSF_INST_MAG" and "PSF_INST_MAG_SIG" are both defined and reasonable in the .cmf files. I can only assume that they're not getting set right by 'addstar' in the load into the DVO database.

Even more details follow as the gdb debug of the relastro command which ideally shouldn't seg fault even if things aren't defined in the DVO database:

[wwoodvas@serenity GPC1] gdb which relastro
GNU gdb 6.3.50-20050815 (Apple version gdb-768) (Tue Oct 2 04:11:19 UTC 2007)
Copyright 2004 Free Software Foundation, Inc.
GDB is free software, covered by the GNU General Public License, and you are
welcome to change it and/or distribute copies of it under certain conditions.
Type "show copying" to see the conditions.
There is absolutely no warranty for GDB. Type "show warranty" for details.
This GDB was configured as "powerpc-apple-darwin"...Warning: /Volumes/data/PS1/data/
GPC1/src: No such file or directory.

Reading symbols for shared libraries .............. done

(gdb) run -D CATDIR gpc1 -D CAMERA gpc1 -region 359.5 361.5 -1.25 1.75 -update-simp
le -v
Starting program: /Volumes/data/PS1/code/ipp/default.darwin/bin/relastro -D CATDIR g
pc1 -D CAMERA gpc1 -region 359.5 361.5 -1.25 1.75 -update-simple -v
Reading symbols for shared libraries +++++++++++++.....warning: Could not find objec
t file "/sw/src/fink.build/gcc43-4.3.1-1000/darwin_objdir/powerpc-apple-darwin9.4.0/
libgcc/_muldi3_s.o" - no debug information available for "../../../gcc-4.3.1/libgcc/
../gcc/libgcc2.c".
[...]
image 472 has 22 measures (2008/09/23,00:00:00, GPC1.y.XY72)
image 473 has 31 measures (2008/09/23,00:00:00, GPC1.y.XY73)image 474 has 6 measures (2008/09/23,00:00:00, GPC1.y.XY74)
image 475 has 22 measures (2008/09/23,00:00:00, GPC1.y.XY75)image 476 has 2 measures (2008/09/23,00:00:00, GPC1.y.XY76)

Program received signal EXC_BAD_ACCESS, Could not access memory.
Reason: KERN_INVALID_ADDRESS at address: 0xc075bdb7
0x0001141c in CoordsGetCenter (fit=0x1d05ca0, tol=0.001, xo=0xbfffa330, yo=0xbfffa338) at /Volumes/data/PS1/src/ipp/Ohana/src/relastro/src/mkpolyterm.c:31
31 R = (xfit[1][0]*yfit[0][1] - xfit[0][1]*yfit[1][0]);
(gdb) p xfit
$1 = (double ) 0x1d7b250
(gdb) p yfit
$2 = (double ) 0x1d6f6d0
(gdb) p yfia[0]t
No symbol "yfia" in current context.
(gdb) p yfit[0]
$3 = (double *) 0x1d6f6e0
(gdb) p xfit[0]
$4 = (double *) 0x1d7b260
(gdb) p xfit[1][0]
Cannot access memory at address 0xc075bdb7
(gdb) p yfit[0][1]
$5 = 6.2747415538852984e-311
(gdb) p fit[0].xfit
$6 = (double ) 0x1d7b250
(gdb) p fit[0].xfia[0]t
A syntax error in expression, near `t'.
(gdb) p fit[0].xfit[0]
$7 = (double *) 0x1d7b260
(gdb) p fit[0].xfit[0][1]
$8 = 1.674254679072588e-311
(gdb) p fit[0].xfit[1][0]
Cannot access memory at address 0xc075bdb7
(gdb) p fit[0].xfit[1]
$9 = (double *) 0xc075bdb7
(gdb) p *fit[0].xfit[1]
Cannot access memory at address 0xc075bdb7
(gdb) p *fit[0].xfit[0]
$10 = 1.2704902607295261
(gdb) bt
#0 0x0001141c in CoordsGetCenter (fit=0x1d05ca0, tol=0.001, xo=0xbfffa330, yo=0xbfffa338) at /Volumes/data/PS1/src/ipp/Ohana/src/relastro/src/mkpolyterm.c:31
#1 0x0000f574 in fit_apply_coords (fit=0x1d05ca0, coords=0x7ab000) at /Volumes/data/PS1/src/ipp/Ohana/src/relastro/src/fitpoly.c:283
#2 0x00004808 in FitSimple (raw=0x825a00, ref=0x82e800, Nmatch=15, coords=0x7ab000) at /Volumes/data/PS1/src/ipp/Ohana/src/relastro/src/FitSimple.c:28
#3 0x0000b840 in UpdateSimple (catalog=0x800000, Ncatalog=4) at /Volumes/data/PS1/src/ipp/Ohana/src/relastro/src/UpdateSimple.c:24
#4 0x00014148 in main (argc=14, argv=0xbfffa658) at /Volumes/data/PS1/src/ipp/Ohana/src/relastro/src/relastro.c:47