diff --git a/do_imfit_tests b/do_imfit_tests index 1b16ffe8..1a2d8fe5 100755 --- a/do_imfit_tests +++ b/do_imfit_tests @@ -204,11 +204,11 @@ echo -n "*** Diff comparison with archives: first test... " STATUS+=$? echo -n "*** Diff comparison with archives: tiny SDSS cutout image... " -./python/diff_printouts.py tests/imfit_reference/imfit_textout2 temptest/test_dump2 --skip-last=3 +./python/diff_printouts.py tests/imfit_reference/macos_x86-64/imfit_textout2 temptest/test_dump2 --skip-last=3 STATUS+=$? echo -n "*** Diff comparison with archives: tiny SDSS cutout image (take 2)... " -./python/diff_printouts.py tests/imfit_reference/imfit_textout3 temptest/test_dump3 --skip-last=3 +./python/diff_printouts.py tests/imfit_reference/macos_x86-64/imfit_textout3 temptest/test_dump3 --skip-last=3 STATUS+=$? echo -n "*** Diff comparison with archives: tiny SDSS cutout image (N-M simplex fit)... " @@ -226,11 +226,11 @@ then fi echo -n "*** Diff comparison with archives: tiny SDSS cutout image (L-M fit with --loud)... " -./python/diff_printouts.py tests/imfit_reference/imfit_textout3d temptest/test_dump3d --skip-last=3 +./python/diff_printouts.py tests/imfit_reference/macos_x86-64/imfit_textout3d temptest/test_dump3d --skip-last=3 STATUS+=$? echo -n "*** Diff comparison with archives: tiny SDSS cutout image (L-M fit with --loud, fixed parameter)... " -./python/diff_printouts.py tests/imfit_reference/imfit_textout3d2 temptest/test_dump3d2 --skip-last=3 +./python/diff_printouts.py tests/imfit_reference/macos_x86-64/imfit_textout3d2 temptest/test_dump3d2 --skip-last=3 STATUS+=$? echo -n "*** Diff comparison with archives: tiny SDSS cutout image (L-M fit using sigma map)... " @@ -238,16 +238,16 @@ echo -n "*** Diff comparison with archives: tiny SDSS cutout image (L-M fit usin STATUS+=$? echo -n "*** Diff comparison with archives: tiny SDSS cutout image (L-M fit using variance map)... " -./python/diff_printouts.py tests/imfit_reference/imfit_textout3e2 temptest/test_dump3e2 --skip-last=3 +./python/diff_printouts.py tests/imfit_reference/macos_x86-64/imfit_textout3e2 temptest/test_dump3e2 --skip-last=3 STATUS+=$? echo -n "*** Diff comparison with archives: fit to larger (rectangular) SDSS image w/ PSF convolution... " -./python/diff_printouts.py tests/imfit_reference/imfit_textout4 temptest/test_dump4 --skip-last=3 +./python/diff_printouts.py tests/imfit_reference/macos_x86-64/imfit_textout4 temptest/test_dump4 --skip-last=3 STATUS+=$? echo -n "*** Diff comparison with archives: same, but using counts/sec image + --exptime option... " -./python/diff_printouts.py tests/imfit_reference/imfit_textout4b temptest/test_dump4b --skip-last=3 +./python/diff_printouts.py tests/imfit_reference/macos_x86-64/imfit_textout4b temptest/test_dump4b --skip-last=3 STATUS+=$? echo -n "*** Diff comparison with archives: chisquare-only calculation... " @@ -290,22 +290,22 @@ STATUS+=$? echo -n "*** Diff comparison with archives: Poisson-MLR-stat fit to problematic image (+ L-M solver)... " tail -n 18 temptest/test_dump5d > temptest/test_dump5d_tail -./python/diff_printouts.py tests/imfit_reference/imfit_textout5d_tail temptest/test_dump5d_tail --skip-last=3 +./python/diff_printouts.py tests/imfit_reference/macos_x86-64/imfit_textout5d_tail temptest/test_dump5d_tail --skip-last=3 STATUS+=$? echo -n "*** Diff comparison with archives: bootstrap-resampling output... " tail -n 7 temptest/temp_bootstrap_output.dat > temptest/temp_bootstrap_output_tail -./python/diff_printouts.py tests/bootstrap_output_seed10_tail.dat temptest/temp_bootstrap_output_tail +./python/diff_printouts.py tests/imfit_reference/macos_x86-64/bootstrap_output_seed10_tail.dat temptest/temp_bootstrap_output_tail STATUS+=$? echo -n "*** Diff comparison with archives: bootstrap-resampling printout (with image offset)... " tail -n 18 temptest/test_dump5e > temptest/test_dump5e_tail -./python/diff_printouts.py tests/imfit_reference/imfit_textout5e_tail temptest/test_dump5e_tail --skip-last=3 +./python/diff_printouts.py tests/imfit_reference/macos_x86-64/imfit_textout5e_tail temptest/test_dump5e_tail --skip-last=3 STATUS+=$? echo -n "*** Diff comparison with archives: bootstrap-resampling output (with image offset)... " tail -n 7 temptest/temp_bootstrap_output2.dat > temptest/temp_bootstrap_output2_tail -./python/diff_printouts.py tests/bootstrap_output_seed10_2_tail.dat temptest/temp_bootstrap_output2_tail +./python/diff_printouts.py tests/imfit_reference/macos_x86-64/bootstrap_output_seed10_2_tail.dat temptest/temp_bootstrap_output2_tail STATUS+=$? @@ -400,7 +400,7 @@ fi echo -n "*** Diff comparison with archives: fitting with bicubic-interpolated PSF image... " -./python/diff_printouts.py tests/imfit_reference/imfit_textout8a temptest/test_dump8a --skip-last=3 +./python/diff_printouts.py tests/imfit_reference/macos_x86-64/imfit_textout8a temptest/test_dump8a --skip-last=3 STATUS+=$? diff --git a/tests/bootstrap_output_seed10_2_tail.dat b/tests/bootstrap_output_seed10_2_tail.dat deleted file mode 100644 index 17b39115..00000000 --- a/tests/bootstrap_output_seed10_2_tail.dat +++ /dev/null @@ -1,7 +0,0 @@ -# Bootstrap resampling output (5 iterations requested): -# X0_1 Y0_1 PA_1 ell_1 n_1 I_e_1 r_e_1 -32.94726552 34.12621527 18.03158227 0.2333425815 2.321036428 22.04684422 56.50895444 -33.00120115 34.07993250 18.19980483 0.2345082374 2.370800674 20.95444970 58.43736181 -32.96236884 34.12467737 18.19960286 0.2352328103 2.431600282 19.59597362 61.31643093 -33.00436662 34.08148751 17.97633518 0.2363263530 2.355999072 21.08248601 58.39570615 -32.96361379 34.07080346 18.65985530 0.2316460593 2.336338352 21.77805322 56.91927251 diff --git a/tests/bootstrap_output_seed10_tail.dat b/tests/bootstrap_output_seed10_tail.dat deleted file mode 100644 index 5246d5c8..00000000 --- a/tests/bootstrap_output_seed10_tail.dat +++ /dev/null @@ -1,7 +0,0 @@ -# Bootstrap resampling output (5 iterations requested): -# X0_1 Y0_1 PA_1 ell_1 n_1 I_e_1 r_e_1 -32.95308050 34.11944163 18.01701845 0.2328484191 2.415204213 19.68006019 61.34187504 -32.90890710 34.09187402 18.67735761 0.2367431324 2.384060089 20.18405061 60.52670009 -32.93601461 34.02577178 18.31447754 0.2332945768 2.369359259 20.62288895 59.41423925 -32.96227662 34.10840801 18.41373630 0.2352079617 2.404208286 20.09082967 60.45350769 -32.92217909 34.05858506 17.96043527 0.2432690148 2.334375744 21.21359885 58.76094420 diff --git a/tests/imfit_reference/macos_x86-64/bootstrap_output_seed10_2_tail.dat b/tests/imfit_reference/macos_x86-64/bootstrap_output_seed10_2_tail.dat new file mode 100644 index 00000000..183fdfca --- /dev/null +++ b/tests/imfit_reference/macos_x86-64/bootstrap_output_seed10_2_tail.dat @@ -0,0 +1,7 @@ +# Bootstrap resampling output (5 iterations requested): +# X0_1 Y0_1 PA_1 ell_1 n_1 I_e_1 r_e_1 +32.94726551 34.12621527 18.03158181 0.2333425815 2.321036457 22.04684347 56.50895582 +33.00120115 34.07993250 18.19980518 0.2345082359 2.370800655 20.95445017 58.43736083 +32.96236884 34.12467736 18.19960251 0.2352328136 2.431600328 19.59597255 61.31643342 +33.00436662 34.08148751 17.97633490 0.2363263502 2.355999091 21.08248554 58.39570695 +32.96361379 34.07080346 18.65985525 0.2316460605 2.336338391 21.77805224 56.91927431 diff --git a/tests/imfit_reference/macos_x86-64/bootstrap_output_seed10_tail.dat b/tests/imfit_reference/macos_x86-64/bootstrap_output_seed10_tail.dat new file mode 100644 index 00000000..7588483b --- /dev/null +++ b/tests/imfit_reference/macos_x86-64/bootstrap_output_seed10_tail.dat @@ -0,0 +1,7 @@ +# Bootstrap resampling output (5 iterations requested): +# X0_1 Y0_1 PA_1 ell_1 n_1 I_e_1 r_e_1 +32.95308050 34.11944163 18.01701865 0.2328484183 2.415204213 19.68006019 61.34187501 +32.90890710 34.09187402 18.67735766 0.2367431340 2.384060118 20.18404995 60.52670153 +32.93601461 34.02577178 18.31447793 0.2332945762 2.369359219 20.62288987 59.41423744 +32.96227662 34.10840801 18.41373693 0.2352079600 2.404208278 20.09082986 60.45350721 +32.92217910 34.05858507 17.96043499 0.2432690140 2.334375758 21.21359851 58.76094479 diff --git a/tests/imfit_reference/macos_x86-64/imfit_textout2 b/tests/imfit_reference/macos_x86-64/imfit_textout2 new file mode 100644 index 00000000..41e2e6b2 --- /dev/null +++ b/tests/imfit_reference/macos_x86-64/imfit_textout2 @@ -0,0 +1,57 @@ + Image file = tests/ic3478rss_64x64.fits + configuration file = tests/imfit_reference/imfit_config_ic3478_64x64.dat + original sky level = 130.1 ADU + gain = 4.725 e-/ADU + read noise = 4.3 e- +Reading data image ("tests/ic3478rss_64x64.fits") ... +naxis1 [# pixels/row] = 64, naxis2 [# pixels/col] = 64; nPixels_tot = 4096 +* No PSF image supplied -- no image convolution will be done! +* No noise image supplied ... will generate noise image from input data image. +Function: Sersic +7 total parameters +Model Object: 4096 data values (pixels) +ModelObject: mask vector applied to weight vector. (4096 valid pixels remain) +7 free parameters (4089 degrees of freedom) +Estimated memory use: 458752 bytes (448.0 KB) + +Performing fit by minimizing chi^2 (data-based errors): +Calling Levenberg-Marquardt solver ... + mpfit iteration 1: fit statistic = 263212.803059 + mpfit iteration 2: fit statistic = 38837.128651 + mpfit iteration 3: fit statistic = 16396.742094 + mpfit iteration 4: fit statistic = 15914.301369 + mpfit iteration 5: fit statistic = 12632.941400 + mpfit iteration 6: fit statistic = 12632.237160 + mpfit iteration 7: fit statistic = 8036.249910 + mpfit iteration 8: fit statistic = 5002.948200 + mpfit iteration 9: fit statistic = 4569.384236 + mpfit iteration 10: fit statistic = 4555.953001 + mpfit iteration 11: fit statistic = 4555.947657 + mpfit iteration 12: fit statistic = 4555.947654 + +*** mpfit status = 1 -- SUCCESS: Convergence in fit-statistic value. + CHI-SQUARE = 4555.947654 (4089 DOF) + INITIAL CHI^2 = 494330.961892 + NPAR = 7 + NFREE = 7 + NPEGGED = 0 + NITER = 13 + NFEV = 98 + +Reduced Chi^2 = 1.114196 +AIC = 4569.975051, BIC = 4614.172017 + +X0 32.9439 # +/- 0.0128 +Y0 34.0933 # +/- 0.0156 +FUNCTION Sersic +PA 18.2612 # +/- 0.30986 deg (CCW from +y axis) +ell 0.235989 # +/- 0.0023165 +n 2.40028 # +/- 0.018072 +I_e 20.0094 # +/- 0.40396 counts/pixel +r_e 60.7612 # +/- 0.86063 pixels + +Saving best-fit parameters in file "bestfit_parameters_imfit.dat" + +(Elapsed time: 0.035381 sec for fit, 0.037554 sec total) +Done! + diff --git a/tests/imfit_reference/macos_x86-64/imfit_textout3 b/tests/imfit_reference/macos_x86-64/imfit_textout3 new file mode 100644 index 00000000..af33f7a9 --- /dev/null +++ b/tests/imfit_reference/macos_x86-64/imfit_textout3 @@ -0,0 +1,57 @@ + Image file = tests/ic3478rss_64x64.fits + configuration file = tests/imfit_reference/imfit_config_ic3478_64x64b.dat +Value from config file: gain = 4.725000 e-/ADU +Value from config file: read noise = 4.300000 e- +Value from config file: original sky = 130.100000 +Reading data image ("tests/ic3478rss_64x64.fits") ... +naxis1 [# pixels/row] = 64, naxis2 [# pixels/col] = 64; nPixels_tot = 4096 +* No PSF image supplied -- no image convolution will be done! +* No noise image supplied ... will generate noise image from input data image. +Function: Sersic +7 total parameters +Model Object: 4096 data values (pixels) +ModelObject: mask vector applied to weight vector. (4096 valid pixels remain) +7 free parameters (4089 degrees of freedom) +Estimated memory use: 458752 bytes (448.0 KB) + +Performing fit by minimizing chi^2 (data-based errors): +Calling Levenberg-Marquardt solver ... + mpfit iteration 1: fit statistic = 263212.803059 + mpfit iteration 2: fit statistic = 38837.128651 + mpfit iteration 3: fit statistic = 16396.742094 + mpfit iteration 4: fit statistic = 15914.301369 + mpfit iteration 5: fit statistic = 12632.941400 + mpfit iteration 6: fit statistic = 12632.237160 + mpfit iteration 7: fit statistic = 8036.249910 + mpfit iteration 8: fit statistic = 5002.948200 + mpfit iteration 9: fit statistic = 4569.384236 + mpfit iteration 10: fit statistic = 4555.953001 + mpfit iteration 11: fit statistic = 4555.947657 + mpfit iteration 12: fit statistic = 4555.947654 + +*** mpfit status = 1 -- SUCCESS: Convergence in fit-statistic value. + CHI-SQUARE = 4555.947654 (4089 DOF) + INITIAL CHI^2 = 494330.961892 + NPAR = 7 + NFREE = 7 + NPEGGED = 0 + NITER = 13 + NFEV = 98 + +Reduced Chi^2 = 1.114196 +AIC = 4569.975051, BIC = 4614.172017 + +X0 32.9439 # +/- 0.0128 +Y0 34.0933 # +/- 0.0156 +FUNCTION Sersic +PA 18.2612 # +/- 0.30986 deg (CCW from +y axis) +ell 0.235989 # +/- 0.0023165 +n 2.40028 # +/- 0.018072 +I_e 20.0094 # +/- 0.40396 counts/pixel +r_e 60.7612 # +/- 0.86063 pixels + +Saving best-fit parameters in file "bestfit_parameters_imfit.dat" + +(Elapsed time: 0.036661 sec for fit, 0.038524 sec total) +Done! + diff --git a/tests/imfit_reference/macos_x86-64/imfit_textout3d b/tests/imfit_reference/macos_x86-64/imfit_textout3d new file mode 100644 index 00000000..16dabdf1 --- /dev/null +++ b/tests/imfit_reference/macos_x86-64/imfit_textout3d @@ -0,0 +1,153 @@ + Image file = tests/ic3478rss_64x64.fits + configuration file = tests/imfit_reference/imfit_config_ic3478_64x64b.dat +Value from config file: gain = 4.725000 e-/ADU +Value from config file: read noise = 4.300000 e- +Value from config file: original sky = 130.100000 +Reading data image ("tests/ic3478rss_64x64.fits") ... +naxis1 [# pixels/row] = 64, naxis2 [# pixels/col] = 64; nPixels_tot = 4096 +* No PSF image supplied -- no image convolution will be done! +* No noise image supplied ... will generate noise image from input data image. +Function: Sersic +7 total parameters +Model Object: 4096 data values (pixels) +ModelObject: mask vector applied to weight vector. (4096 valid pixels remain) +7 free parameters (4089 degrees of freedom) +Estimated memory use: 458752 bytes (448.0 KB) + +Performing fit by minimizing chi^2 (data-based errors): +Calling Levenberg-Marquardt solver ... + mpfit iteration 1: fit statistic = 263212.803059 +X0 32.3912 +Y0 32.7346 +FUNCTION Sersic +PA 18.0332 # deg (CCW from +y axis) +ell 0 +n 1.81654 +I_e 19.043 # counts/pixel +r_e 30.4484 # pixels + mpfit iteration 2: fit statistic = 38837.128651 +X0 32.8246 +Y0 33.5074 +FUNCTION Sersic +PA 18.0332 # deg (CCW from +y axis) +ell 0 +n 2.20493 +I_e 24.7011 # counts/pixel +r_e 38.3852 # pixels + mpfit iteration 3: fit statistic = 16396.742094 +X0 32.9992 +Y0 34.0253 +FUNCTION Sersic +PA 25.3509 # deg (CCW from +y axis) +ell 0 +n 2.40435 +I_e 19.0899 # counts/pixel +r_e 51.3035 # pixels + mpfit iteration 4: fit statistic = 15914.301369 +X0 32.9958 +Y0 34.0292 +FUNCTION Sersic +PA 90 # deg (CCW from +y axis) +ell 0 +n 2.40632 +I_e 19.1018 # counts/pixel +r_e 51.481 # pixels + mpfit iteration 5: fit statistic = 12632.941400 +X0 32.9480 +Y0 34.0851 +FUNCTION Sersic +PA 90 # deg (CCW from +y axis) +ell 0 +n 2.43403 +I_e 19.2689 # counts/pixel +r_e 53.9773 # pixels + mpfit iteration 6: fit statistic = 12632.237160 +X0 32.9485 +Y0 34.0895 +FUNCTION Sersic +PA 0 # deg (CCW from +y axis) +ell 0 +n 2.43444 +I_e 19.2547 # counts/pixel +r_e 53.9852 # pixels + mpfit iteration 7: fit statistic = 8036.249910 +X0 32.9468 +Y0 34.0973 +FUNCTION Sersic +PA 0.00311914 # deg (CCW from +y axis) +ell 0.18689 +n 2.39424 +I_e 20.3117 # counts/pixel +r_e 56.873 # pixels + mpfit iteration 8: fit statistic = 5002.948200 +X0 32.9462 +Y0 34.0948 +FUNCTION Sersic +PA 22.3975 # deg (CCW from +y axis) +ell 0.194004 +n 2.41769 +I_e 19.5831 # counts/pixel +r_e 59.8389 # pixels + mpfit iteration 9: fit statistic = 4569.384236 +X0 32.9434 +Y0 34.0951 +FUNCTION Sersic +PA 17.2139 # deg (CCW from +y axis) +ell 0.235027 +n 2.40359 +I_e 19.9373 # counts/pixel +r_e 60.7908 # pixels + mpfit iteration 10: fit statistic = 4555.953001 +X0 32.9442 +Y0 34.0933 +FUNCTION Sersic +PA 18.2648 # deg (CCW from +y axis) +ell 0.23583 +n 2.4004 +I_e 20.0071 # counts/pixel +r_e 60.7584 # pixels + mpfit iteration 11: fit statistic = 4555.947657 +X0 32.9440 +Y0 34.0933 +FUNCTION Sersic +PA 18.2612 # deg (CCW from +y axis) +ell 0.235989 +n 2.4003 +I_e 20.009 # counts/pixel +r_e 60.7621 # pixels + mpfit iteration 12: fit statistic = 4555.947654 +X0 32.9439 +Y0 34.0933 +FUNCTION Sersic +PA 18.2612 # deg (CCW from +y axis) +ell 0.235989 +n 2.40028 +I_e 20.0094 # counts/pixel +r_e 60.7612 # pixels + +*** mpfit status = 1 -- SUCCESS: Convergence in fit-statistic value. + CHI-SQUARE = 4555.947654 (4089 DOF) + INITIAL CHI^2 = 494330.961892 + NPAR = 7 + NFREE = 7 + NPEGGED = 0 + NITER = 13 + NFEV = 98 + +Reduced Chi^2 = 1.114196 +AIC = 4569.975051, BIC = 4614.172017 + +X0 32.9439 # +/- 0.0128 +Y0 34.0933 # +/- 0.0156 +FUNCTION Sersic +PA 18.2612 # +/- 0.30986 deg (CCW from +y axis) +ell 0.235989 # +/- 0.0023165 +n 2.40028 # +/- 0.018072 +I_e 20.0094 # +/- 0.40396 counts/pixel +r_e 60.7612 # +/- 0.86063 pixels + +Saving best-fit parameters in file "bestfit_parameters_imfit.dat" + +(Elapsed time: 0.035436 sec for fit, 0.037478 sec total) +Done! + diff --git a/tests/imfit_reference/macos_x86-64/imfit_textout3d2 b/tests/imfit_reference/macos_x86-64/imfit_textout3d2 new file mode 100644 index 00000000..111b1f9c --- /dev/null +++ b/tests/imfit_reference/macos_x86-64/imfit_textout3d2 @@ -0,0 +1,117 @@ + Image file = tests/ic3478rss_64x64.fits + configuration file = tests/imfit_reference/imfit_config_ic3478_64x64c.dat +Value from config file: gain = 4.725000 e-/ADU +Value from config file: read noise = 4.300000 e- +Value from config file: original sky = 130.100000 +Reading data image ("tests/ic3478rss_64x64.fits") ... +naxis1 [# pixels/row] = 64, naxis2 [# pixels/col] = 64; nPixels_tot = 4096 +* No PSF image supplied -- no image convolution will be done! +* No noise image supplied ... will generate noise image from input data image. +Function: Sersic +7 total parameters +Model Object: 4096 data values (pixels) +ModelObject: mask vector applied to weight vector. (4096 valid pixels remain) +6 free parameters (4090 degrees of freedom) +Estimated memory use: 425984 bytes (416.0 KB) + +Performing fit by minimizing chi^2 (data-based errors): +Calling Levenberg-Marquardt solver ... + mpfit iteration 1: fit statistic = 198337.464916 +X0 32.3180 +Y0 32.5876 +FUNCTION Sersic +PA 14.6519 # deg (CCW from +y axis) +ell 0 +n 2.4 +I_e 18.061 # counts/pixel +r_e 32.6624 # pixels + mpfit iteration 2: fit statistic = 15612.629060 +X0 32.9545 +Y0 33.5463 +FUNCTION Sersic +PA 14.6519 # deg (CCW from +y axis) +ell 0 +n 2.4 +I_e 22.9445 # counts/pixel +r_e 46.03 # pixels + mpfit iteration 3: fit statistic = 10937.152322 +X0 32.9686 +Y0 33.8079 +FUNCTION Sersic +PA 14.6519 # deg (CCW from +y axis) +ell 0.0709715 +n 2.4 +I_e 23.0724 # counts/pixel +r_e 47.3492 # pixels + mpfit iteration 4: fit statistic = 7037.449529 +X0 32.9579 +Y0 34.0175 +FUNCTION Sersic +PA 27.4742 # deg (CCW from +y axis) +ell 0.175583 +n 2.4 +I_e 22.6486 # counts/pixel +r_e 51.3717 # pixels + mpfit iteration 5: fit statistic = 5140.174422 +X0 32.9440 +Y0 34.0860 +FUNCTION Sersic +PA 15.1246 # deg (CCW from +y axis) +ell 0.224897 +n 2.4 +I_e 19.8625 # counts/pixel +r_e 59.3537 # pixels + mpfit iteration 6: fit statistic = 4556.680252 +X0 32.9445 +Y0 34.0937 +FUNCTION Sersic +PA 18.4899 # deg (CCW from +y axis) +ell 0.234923 +n 2.4 +I_e 20.0148 # counts/pixel +r_e 60.6989 # pixels + mpfit iteration 7: fit statistic = 4555.947935 +X0 32.9439 +Y0 34.0934 +FUNCTION Sersic +PA 18.2603 # deg (CCW from +y axis) +ell 0.235983 +n 2.4 +I_e 20.0157 # counts/pixel +r_e 60.7479 # pixels + mpfit iteration 8: fit statistic = 4555.947901 +X0 32.9439 +Y0 34.0933 +FUNCTION Sersic +PA 18.2613 # deg (CCW from +y axis) +ell 0.235985 +n 2.4 +I_e 20.0156 # counts/pixel +r_e 60.7483 # pixels + +*** mpfit status = 1 -- SUCCESS: Convergence in fit-statistic value. + CHI-SQUARE = 4555.947901 (4090 DOF) + INITIAL CHI^2 = 439705.788935 + NPAR = 7 + NFREE = 6 + NPEGGED = 0 + NITER = 9 + NFEV = 59 + +Reduced Chi^2 = 1.113924 +AIC = 4567.968444, BIC = 4605.854498 + +X0 32.9439 # +/- 0.0128 +Y0 34.0933 # +/- 0.0156 +FUNCTION Sersic +PA 18.2613 # +/- 0.30982 deg (CCW from +y axis) +ell 0.235985 # +/- 0.0023011 +n 2.4 # +/- 0 +I_e 20.0156 # +/- 0.077391 counts/pixel +r_e 60.7483 # +/- 0.24714 pixels + +Saving best-fit parameters in file "bestfit_parameters_imfit.dat" + +(Elapsed time: 0.020950 sec for fit, 0.022975 sec total) +Done! + diff --git a/tests/imfit_reference/macos_x86-64/imfit_textout3e2 b/tests/imfit_reference/macos_x86-64/imfit_textout3e2 new file mode 100644 index 00000000..50d141af --- /dev/null +++ b/tests/imfit_reference/macos_x86-64/imfit_textout3e2 @@ -0,0 +1,58 @@ + Image file = tests/ic3478rss_64x64.fits + configuration file = tests/imfit_reference/imfit_config_ic3478_64x64c.dat + noise image = tests/ic3478rss_64x64_variance.fits +Value from config file: gain = 4.725000 e-/ADU +Value from config file: read noise = 4.300000 e- +Value from config file: original sky = 130.100000 +Reading data image ("tests/ic3478rss_64x64.fits") ... +naxis1 [# pixels/row] = 64, naxis2 [# pixels/col] = 64; nPixels_tot = 4096 +Reading noise image ("tests/ic3478rss_64x64_variance.fits") ... +* No PSF image supplied -- no image convolution will be done! +Function: Sersic +7 total parameters +Model Object: 4096 data values (pixels) +ModelObject: mask vector applied to weight vector. (4096 valid pixels remain) +6 free parameters (4090 degrees of freedom) +Estimated memory use: 425984 bytes (416.0 KB) + +Performing fit by minimizing chi^2 (user-supplied error image): +Calling Levenberg-Marquardt solver ... + mpfit iteration 1: fit statistic = 198337.457290 + mpfit iteration 2: fit statistic = 15612.628210 + mpfit iteration 3: fit statistic = 15095.395477 + mpfit iteration 4: fit statistic = 12159.620345 + mpfit iteration 5: fit statistic = 10741.693780 + mpfit iteration 6: fit statistic = 8390.025557 + mpfit iteration 7: fit statistic = 6277.909722 + mpfit iteration 8: fit statistic = 5023.975651 + mpfit iteration 9: fit statistic = 4585.082498 + mpfit iteration 10: fit statistic = 4555.947901 + mpfit iteration 11: fit statistic = 4555.947598 + mpfit iteration 12: fit statistic = 4555.947598 + +*** mpfit status = 1 -- SUCCESS: Convergence in fit-statistic value. + CHI-SQUARE = 4555.947598 (4090 DOF) + INITIAL CHI^2 = 439705.768065 + NPAR = 7 + NFREE = 6 + NPEGGED = 0 + NITER = 13 + NFEV = 86 + +Reduced Chi^2 = 1.113924 +AIC = 4567.968141, BIC = 4605.854195 + +X0 32.9439 # +/- 0.0128 +Y0 34.0933 # +/- 0.0156 +FUNCTION Sersic +PA 18.2613 # +/- 0.30981 deg (CCW from +y axis) +ell 0.235985 # +/- 0.0023011 +n 2.4 # +/- 0 +I_e 20.0156 # +/- 0.077391 counts/pixel +r_e 60.7483 # +/- 0.24714 pixels + +Saving best-fit parameters in file "bestfit_parameters_imfit.dat" + +(Elapsed time: 0.030884 sec for fit, 0.032940 sec total) +Done! + diff --git a/tests/imfit_reference/macos_x86-64/imfit_textout4 b/tests/imfit_reference/macos_x86-64/imfit_textout4 new file mode 100644 index 00000000..8e70d615 --- /dev/null +++ b/tests/imfit_reference/macos_x86-64/imfit_textout4 @@ -0,0 +1,72 @@ + Image file = tests/n3073rss_small.fits + configuration file = tests/imfit_reference/imfit_config_n3073.dat + PSF image = tests/psf_moffat_35_n4699z.fits + mask image = tests/n3073rss_small_mask.fits +Value from config file: gain = 4.725000 e-/ADU +Value from config file: read noise = 4.300000 e- +Value from config file: original sky = 154.330000 +Reading data image ("tests/n3073rss_small.fits") ... +naxis1 [# pixels/row] = 150, naxis2 [# pixels/col] = 200; nPixels_tot = 30000 +Reading mask image ("tests/n3073rss_small_mask.fits") ... +Reading PSF image ("tests/psf_moffat_35_n4699z.fits") ... +naxis1 [# pixels/row] = 35, naxis2 [# pixels/col] = 35; nPixels_tot = 1225 +ModelObject::AddMaskVector -- treating zero-valued pixels as good ... +* No noise image supplied ... will generate noise image from input data image. +Function: Sersic +Function: Exponential +11 total parameters +Model Object: 30000 data values (pixels) +ModelObject: mask vector applied to weight vector. (28754 valid pixels remain) +11 free parameters (28743 degrees of freedom) +Estimated memory use: 8747785 bytes (8.3 MB) + +Performing fit by minimizing chi^2 (data-based errors): +Calling Levenberg-Marquardt solver ... + mpfit iteration 1: fit statistic = 49095.070363 + mpfit iteration 2: fit statistic = 37345.298170 + mpfit iteration 3: fit statistic = 37110.187479 + mpfit iteration 4: fit statistic = 37053.537386 + mpfit iteration 5: fit statistic = 36993.750717 + mpfit iteration 6: fit statistic = 36947.464078 + mpfit iteration 7: fit statistic = 36910.889944 + mpfit iteration 8: fit statistic = 36882.958246 + mpfit iteration 9: fit statistic = 36862.501461 + mpfit iteration 10: fit statistic = 36848.786565 + mpfit iteration 11: fit statistic = 36841.049095 + mpfit iteration 12: fit statistic = 36831.936835 + mpfit iteration 13: fit statistic = 36828.584568 + mpfit iteration 14: fit statistic = 36828.578485 + mpfit iteration 15: fit statistic = 36828.577740 + mpfit iteration 16: fit statistic = 36828.577653 + +*** mpfit status = 1 -- SUCCESS: Convergence in fit-statistic value. + CHI-SQUARE = 36828.577653 (28743 DOF) + INITIAL CHI^2 = 235606.146670 + NPAR = 11 + NFREE = 11 + NPEGGED = 0 + NITER = 17 + NFEV = 196 + +Reduced Chi^2 = 1.281306 +AIC = 36850.586838, BIC = 36941.509506 + +X0 100.7073 # +/- 0.0058 +Y0 97.0877 # +/- 0.0052 +FUNCTION Sersic +PA 111.919 # +/- 0.26435 deg (CCW from +y axis) +ell 0.184262 # +/- 0.0015989 +n 3.2805 # +/- 0.034084 +I_e 78.2973 # +/- 2.693 counts/pixel +r_e 26.7975 # +/- 0.6481 pixels +FUNCTION Exponential +PA 51.724 # +/- 0.80649 deg (CCW from +y axis) +ell 0.35226 # +/- 0.014009 +I_0 39.7754 # +/- 1.8008 counts/pixel +h 53.6058 # +/- 1.1641 pixels + +Saving best-fit parameters in file "bestfit_parameters_imfit.dat" + +(Elapsed time: 0.665725 sec for fit, 0.678504 sec total) +Done! + diff --git a/tests/imfit_reference/macos_x86-64/imfit_textout4b b/tests/imfit_reference/macos_x86-64/imfit_textout4b new file mode 100644 index 00000000..29f6fbd8 --- /dev/null +++ b/tests/imfit_reference/macos_x86-64/imfit_textout4b @@ -0,0 +1,70 @@ + Image file = tests/n3073rss_small_cps.fits + configuration file = tests/imfit_reference/imfit_config_n3073.dat + PSF image = tests/psf_moffat_35_n4699z.fits + mask image = tests/n3073rss_small_mask.fits + original sky level = 2.86287 ADU + exposure time = 53.9075 sec +Value from config file: gain = 4.725000 e-/ADU +Value from config file: read noise = 4.300000 e- +Original-sky value in config file ignored (using command-line value) +Reading data image ("tests/n3073rss_small_cps.fits") ... +naxis1 [# pixels/row] = 150, naxis2 [# pixels/col] = 200; nPixels_tot = 30000 +Reading mask image ("tests/n3073rss_small_mask.fits") ... +Reading PSF image ("tests/psf_moffat_35_n4699z.fits") ... +naxis1 [# pixels/row] = 35, naxis2 [# pixels/col] = 35; nPixels_tot = 1225 +ModelObject::AddMaskVector -- treating zero-valued pixels as good ... +* No noise image supplied ... will generate noise image from input data image. +Function: Sersic +Function: Exponential +11 total parameters +Model Object: 30000 data values (pixels) +ModelObject: mask vector applied to weight vector. (28754 valid pixels remain) +11 free parameters (28743 degrees of freedom) +Estimated memory use: 8747785 bytes (8.3 MB) + +Performing fit by minimizing chi^2 (data-based errors): +Calling Levenberg-Marquardt solver ... + mpfit iteration 1: fit statistic = 1377107.029550 + mpfit iteration 2: fit statistic = 1050140.180068 + mpfit iteration 3: fit statistic = 120559.627697 + mpfit iteration 4: fit statistic = 71413.213300 + mpfit iteration 5: fit statistic = 43860.266858 + mpfit iteration 6: fit statistic = 36911.962119 + mpfit iteration 7: fit statistic = 36828.859254 + mpfit iteration 8: fit statistic = 36828.611690 + mpfit iteration 9: fit statistic = 36828.581709 + mpfit iteration 10: fit statistic = 36828.578450 + mpfit iteration 11: fit statistic = 36828.578059 + mpfit iteration 12: fit statistic = 36828.578011 + +*** mpfit status = 1 -- SUCCESS: Convergence in fit-statistic value. + CHI-SQUARE = 36828.578011 (28743 DOF) + INITIAL CHI^2 = 3184994682.091280 + NPAR = 11 + NFREE = 11 + NPEGGED = 0 + NITER = 13 + NFEV = 147 + +Reduced Chi^2 = 1.281306 +AIC = 36850.587197, BIC = 36941.509865 + +X0 100.7073 # +/- 0.0058 +Y0 97.0877 # +/- 0.0052 +FUNCTION Sersic +PA 111.919 # +/- 0.26435 deg (CCW from +y axis) +ell 0.184258 # +/- 0.0015988 +n 3.28046 # +/- 0.034083 +I_e 1.45249 # +/- 0.049947 counts/pixel +r_e 26.797 # +/- 0.64814 pixels +FUNCTION Exponential +PA 51.7247 # +/- 0.80657 deg (CCW from +y axis) +ell 0.352234 # +/- 0.014013 +I_0 0.73783 # +/- 0.0334 counts/pixel +h 53.6064 # +/- 1.1648 pixels + +Saving best-fit parameters in file "bestfit_parameters_imfit.dat" + +(Elapsed time: 0.509414 sec for fit, 0.523494 sec total) +Done! + diff --git a/tests/imfit_reference/macos_x86-64/imfit_textout5d_tail b/tests/imfit_reference/macos_x86-64/imfit_textout5d_tail new file mode 100644 index 00000000..07ff422b --- /dev/null +++ b/tests/imfit_reference/macos_x86-64/imfit_textout5d_tail @@ -0,0 +1,18 @@ + +Reduced Chi^2 equivalent = 0.227929 +AIC = 18123.661998, BIC = 18188.641651 + +X0 492.0205 # +/- 0.0626 +Y0 444.3322 # +/- 0.0641 +FUNCTION Sersic +PA 138.318 # +/- 0.48695 deg (CCW from +y axis) +ell 0.250331 # +/- 0.0037457 +n 1.03804 # +/- 0.0080468 +I_e 7.83608 # +/- 0.06377 counts/pixel +r_e 34.9043 # +/- 0.17459 pixels + +Saving best-fit parameters in file "bestfit_parameters_imfit.dat" + +(Elapsed time: 0.345436 sec for fit, 0.355686 sec total) +Done! + diff --git a/tests/imfit_reference/macos_x86-64/imfit_textout5e_tail b/tests/imfit_reference/macos_x86-64/imfit_textout5e_tail new file mode 100644 index 00000000..d322d072 --- /dev/null +++ b/tests/imfit_reference/macos_x86-64/imfit_textout5e_tail @@ -0,0 +1,18 @@ +Starting bootstrap iterations (L-M solver): +[================> ] 1 (20.0%) [================================> ] 2 (40.0%) [================================================> ] 3 (60.0%) [================================================================> ] 4 (80.0%) [================================================================================] 5 (100.0%) + +Statistics for parameter values from bootstrap resampling (5 successful iterations): +Best-fit Bootstrap [68% conf.int., half-width]; (mean +/- standard deviation) +X0 = 32.9497 +0.0547039, -0.00239716 [32.9473 -- 33.0044, 0.0285506]; (32.9758 +/- 0.0255166) +Y0 = 34.0986 +0.027661, -0.0277509 [34.0708 -- 34.1262, 0.0277059]; (34.0966 +/- 0.026632) +PA = 18.1649 +0.494974, -0.188546 [17.9763 -- 18.6599, 0.34176]; (18.2134 +/- 0.268773) +ell = 0.234412 +0.00191482, -0.00276547 [0.231646 -- 0.236326, 0.00234014]; (0.234211 +/- 0.00179863) +n = 2.35479 +0.0768108, -0.033753 [2.32104 -- 2.4316, 0.0552819]; (2.36315 +/- 0.0426835) +I_e = 21.2168 +0.830046, -1.62082 [19.596 -- 22.0468, 1.22544]; (21.0916 +/- 0.953874) +r_e = 58.0709 +3.24553, -1.56195 [56.509 -- 61.3164, 2.40374]; (58.3155 +/- 1.88679) +Bootstrap-resampling output saved to file "temptest/temp_bootstrap_output2.dat". +Saving best-fit parameters in file "bestfit_parameters_imfit.dat" + +(Elapsed time: 0.035934 sec for fit, 0.054264 for bootstrap, 0.093777 sec total) +Done! + diff --git a/tests/imfit_reference/macos_x86-64/imfit_textout8a b/tests/imfit_reference/macos_x86-64/imfit_textout8a new file mode 100644 index 00000000..88f09661 --- /dev/null +++ b/tests/imfit_reference/macos_x86-64/imfit_textout8a @@ -0,0 +1,50 @@ + Image file = tests/psf_moffat_fwhm2.fits + configuration file = tests/imfit_reference/config_imfit_ptsource.dat + PSF image = tests/psf_moffat_fwhm2_35x35.fits +Reading data image ("tests/psf_moffat_fwhm2.fits") ... +naxis1 [# pixels/row] = 40, naxis2 [# pixels/col] = 40; nPixels_tot = 1600 +Reading PSF image ("tests/psf_moffat_fwhm2_35x35.fits") ... +naxis1 [# pixels/row] = 35, naxis2 [# pixels/col] = 35; nPixels_tot = 1225 +* No noise image supplied ... will generate noise image from input data image. +Function: PointSource +Function: FlatSky +4 total parameters +Model Object: 1600 data values (pixels) +ModelObject: mask vector applied to weight vector. (1600 valid pixels remain) +3 free parameters (1597 degrees of freedom) +Estimated memory use: 1383465 bytes (1.3 MB) + +Performing fit by minimizing chi^2 (data-based errors): +Calling Levenberg-Marquardt solver ... + mpfit iteration 1: fit statistic = 241515.661816 + mpfit iteration 2: fit statistic = 126454.117011 + mpfit iteration 3: fit statistic = 23833.639972 + mpfit iteration 4: fit statistic = 2263.479645 + mpfit iteration 5: fit statistic = 1991.947078 + mpfit iteration 6: fit statistic = 1991.946235 + mpfit iteration 7: fit statistic = 1991.946235 + +*** mpfit status = 1 -- SUCCESS: Convergence in fit-statistic value. + CHI-SQUARE = 1991.946235 (1597 DOF) + INITIAL CHI^2 = 1846834.686862 + NPAR = 4 + NFREE = 3 + NPEGGED = 0 + NITER = 8 + NFEV = 30 + +Reduced Chi^2 = 1.247305 +AIC = 1997.961272, BIC = 2014.079511 + +X0 20.4006 # +/- 0.0021 +Y0 20.3985 # +/- 0.0021 +FUNCTION PointSource +I_tot 263918 # +/- 519.52 counts +FUNCTION FlatSky +I_sky 100 # +/- 0 counts/pixel + +Saving best-fit parameters in file "bestfit_parameters_imfit.dat" + +(Elapsed time: 0.034380 sec for fit, 0.047271 sec total) +Done! +