Changeset 21363 for trunk/psModules/src/objects/pmSourcePhotometry.c
- Timestamp:
- Feb 5, 2009, 4:31:25 PM (17 years ago)
- File:
-
- 1 edited
-
trunk/psModules/src/objects/pmSourcePhotometry.c (modified) (10 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/psModules/src/objects/pmSourcePhotometry.c
r21183 r21363 3 3 * @author EAM, IfA; GLG, MHPCC 4 4 * 5 * @version $Revision: 1.4 8$ $Name: not supported by cvs2svn $6 * @date $Date: 2009-0 1-27 06:39:38$5 * @version $Revision: 1.49 $ $Name: not supported by cvs2svn $ 6 * @date $Date: 2009-02-06 02:31:25 $ 7 7 * 8 8 * Copyright 2004 Maui High Performance Computing Center, University of Hawaii … … 446 446 assert (Pj != NULL); 447 447 448 const psImage *Wi = Mi-> weight;448 const psImage *Wi = Mi->variance; 449 449 if (!unweighted_sum) { 450 450 assert (Wi != NULL); … … 511 511 assert (Pj != NULL); 512 512 513 const psImage *Wi = Mi-> weight;513 const psImage *Wi = Mi->variance; 514 514 if (!unweighted_sum) { 515 515 assert (Wi != NULL); … … 567 567 const psImage *Pi = Mi->pixels; 568 568 assert (Pi != NULL); 569 const psImage *Wi = Mi-> weight;569 const psImage *Wi = Mi->variance; 570 570 if (!unweighted_sum) { 571 571 assert (Wi != NULL); … … 612 612 # endif 613 613 614 bool pmSourceChisq (pmModel *model, psImage *image, psImage *mask, psImage * weight,614 bool pmSourceChisq (pmModel *model, psImage *image, psImage *mask, psImage *variance, 615 615 psImageMaskType maskVal) 616 616 { … … 618 618 PS_ASSERT_PTR_NON_NULL(image, false); 619 619 PS_ASSERT_PTR_NON_NULL(mask, false); 620 PS_ASSERT_PTR_NON_NULL( weight, false);620 PS_ASSERT_PTR_NON_NULL(variance, false); 621 621 622 622 double dC = 0.0; … … 626 626 if (mask->data.PS_TYPE_IMAGE_MASK_DATA[j][i] & maskVal) 627 627 continue; 628 if ( weight->data.F32[j][i] <= 0)629 continue; 630 dC += PS_SQR (image->data.F32[j][i]) / weight->data.F32[j][i];628 if (variance->data.F32[j][i] <= 0) 629 continue; 630 dC += PS_SQR (image->data.F32[j][i]) / variance->data.F32[j][i]; 631 631 Npix ++; 632 632 } … … 648 648 const psImage *Pi = Mi->modelFlux; 649 649 assert (Pi != NULL); 650 const psImage *Wi = Mi-> weight;650 const psImage *Wi = Mi->variance; 651 651 if (!unweighted_sum) { 652 652 assert (Wi != NULL); … … 706 706 assert (Pj != NULL); 707 707 708 const psImage *Wi = Mi-> weight;708 const psImage *Wi = Mi->variance; 709 709 if (!unweighted_sum) { 710 710 assert (Wi != NULL); … … 770 770 assert (Pj != NULL); 771 771 772 const psImage *Wi = Mi-> weight;772 const psImage *Wi = Mi->variance; 773 773 if (!unweighted_sum) { 774 774 assert (Wi != NULL);
Note:
See TracChangeset
for help on using the changeset viewer.
