Changeset 34078 for branches/eam_branches/ipp-20120601/psphot/src
- Timestamp:
- Jun 26, 2012, 10:54:57 AM (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/eam_branches/ipp-20120601/psphot/src/psphotSourceFits.c
r32744 r34078 224 224 if (source->type == PM_SOURCE_TYPE_DEFECT) return false; 225 225 if (source->type == PM_SOURCE_TYPE_SATURATED) return false; 226 227 # define TEST_X -420.0228 # define TEST_Y 300.0229 230 if ((fabs(source->peak->xf - TEST_X) < 5) && (fabs(source->peak->yf - TEST_Y) < 5)) {231 fprintf (stderr, "test galaxy\n");232 }233 234 # undef TEST_X235 # undef TEST_Y236 226 237 227 // set the radius based on the footprint (also sets the mask pixels) … … 513 503 } 514 504 515 # define TEST_X -540.0516 # define TEST_Y 540.0517 518 if ((fabs(source->peak->xf - TEST_X) < 5) && (fabs(source->peak->yf - TEST_Y) < 5)) {519 psTraceSetLevel("psModules.objects.pmPCM_MinimizeChisq", 5);520 }521 522 505 float t1, t2, t4, t5; 523 506 if (TIMING) { psTimerStart ("psphotFitPCM"); } … … 569 552 } 570 553 571 if ((fabs(source->peak->xf - TEST_X) < 5) && (fabs(source->peak->yf - TEST_Y) < 5)) {572 psTraceSetLevel("psModules.objects.pmPCM_MinimizeChisq", 0);573 }574 575 554 // psTraceSetLevel("psLib.math.psMinimizeLMChi2", 0); 576 555 psFree (pcm); … … 578 557 return model; 579 558 } 580 581 # undef TEST_X582 # undef TEST_Y583 559 584 560 // note that these should be 1/2n of the standard sersic index … … 603 579 float xMin = NAN; 604 580 float chiSquare[N_INDEX_GUESS]; 605 606 # define TEST_X -540.0607 # define TEST_Y 540.0608 609 if ((fabs(source->peak->xf - TEST_X) < 5) && (fabs(source->peak->yf - TEST_Y) < 5)) {610 psTraceSetLevel("psModules.objects.pmPCM_MinimizeChisq", 5);611 }612 581 613 582 for (int i = 0; i < N_INDEX_GUESS; i++) { … … 635 604 assert (iMin >= 0); 636 605 637 if ((fabs(source->peak->xf - TEST_X) < 5) && (fabs(source->peak->yf - TEST_Y) < 5)) {638 psTraceSetLevel("psModules.objects.pmPCM_MinimizeChisq", 0);639 }640 641 606 model->flags = PM_MODEL_STATUS_NONE; // do not attempt to handle failures here, let the next iteration deal with it 642 607 model->params->data.F32[PM_PAR_7] = 0.5/indexGuess[iMin]; … … 666 631 float xMin = NAN; 667 632 float chiSquare[N_INDEX_GUESS]; 668 669 if ((fabs(source->peak->xf - TEST_X) < 5) && (fabs(source->peak->yf - TEST_Y) < 5)) {670 psTraceSetLevel("psModules.objects.pmPCM_MinimizeChisq", 5);671 }672 633 673 634 for (int i = 0; i < N_INDEX_GUESS; i++) { … … 701 662 } 702 663 assert (iMin >= 0); 703 704 if ((fabs(source->peak->xf - TEST_X) < 5) && (fabs(source->peak->yf - TEST_Y) < 5)) {705 psTraceSetLevel("psModules.objects.pmPCM_MinimizeChisq", 0);706 }707 664 708 665 model->flags = PM_MODEL_STATUS_NONE; // do not attempt to handle failures here, let the next iteration deal with it
Note:
See TracChangeset
for help on using the changeset viewer.
