| | 7 | |
| | 8 | * I worked with Serge to define the cmf translation code for the TRAIL fits. I was getting some rare core dumps when running trail fitting tests. I put it on the shelf for a bit and Bill found the same problem independently. He tracked it down to a thread-unsafe piece of code being used in a threaded context and fixed the bug. Now we can go ahead with TRAIL tests on gpc1 data and finalize the analysis. |
| | 9 | * I implemented Tamas' exact projection cell centers in our skycell generation code (using the C algorithm he kindly sent me). I then explored the differences between RINGS.V3 (our current version) and Tamas' projection cells. I'm writing a more detailed report, but the summary is that the differences are mostly minor. |
| | 10 | * I created a code concept to define the projection cell boundaries (Boundary Tree table) and generated a boundary tree for the RINGS.V3 tessellation. The boundary tree translates between a specified RA & DEC and the primary projection cell based on the lines of constant RA & DEC at the mid points of the projection cell centers. Because RINGS.V3 has uniform spacing of the zones in DEC and uniform spacing of the projection cell centers in a zone, the translation is a two step evaluation of a linear equation, and thus very fast. Evaluation of 10^8^ points only takes about 0.6 seconds on my desktop machine. (Tamas' projections cells also have linear spacing of RA along a zone, but the zone centers are not uniformly spaced. I have not (yet?) implemented a boundary tree sufficiently generic for that tessellation as it requires an initial search by bisection and a bit more information to define each zone.) |
| | 11 | * I added code to DVO to set the secfilt mean flux (psf or kron) based on the primary projection cell using the boundary tree code above. This is a two-step selection: |
| | 12 | * determine the primary projection cell for the object's RA & DEC |
| | 13 | * if multiple stack detections come from the primary projection, use the one closest to its skycell center. |
| | 14 | * if no stack detections come from the primary projection, use the non-primary detection which is closest to its skycell center. |
| | 15 | * I also worked on the text and figures for Ken's NSF proposal response. |
