corrupting star at 567.25 1094.32 corrupting star at 1532.74 1743.83 start with_hsm for imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=100, xmax=147, ymin=322, ymax=369), array= array([[ 1.57293516e-06, 4.83428551e-07, 1.01336786e-06, ..., -1.52099460e-06, 1.62390904e-06, -3.06897618e-07], [ -4.01491747e-07, -6.06370918e-07, -1.84310323e-07, ..., -3.39801630e-07, 4.16551785e-07, 1.06846937e-06], [ -4.68112916e-07, 2.85567381e-07, -4.94964070e-07, ..., -1.98679504e-06, -9.59963018e-07, 1.69272738e-07], ..., [ 4.87750185e-07, 5.51780090e-07, 4.79130222e-07, ..., -2.40008035e-06, 1.37110010e-06, 1.22177823e-06], [ -2.55823608e-07, -8.19296361e-07, 2.35463608e-06, ..., -1.41771136e-06, -1.60905051e-06, 4.51424427e-07], [ -2.71851150e-07, 1.25770748e-06, 1.75458837e-07, ..., -1.10357121e-06, 1.51166330e-06, -2.85283708e-09]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=100, xmax=147, ymin=322, ymax=369), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(123.12,345.43) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=100, xmax=147, ymin=322, ymax=369), moments_status=0, observed_shape=galsim.Shear((0.23000160481118787-0.16999831959428038j)), moments_sigma=4.9998674392700195, moments_amp=0.9999805688858032, moments_centroid=galsim.PositionD(x=123.11974024161614, y=345.42992834522084), moments_rho4=1.9999460194128353, moments_n_iter=8) got hsm parameters start with_hsm for imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=100, xmax=147, ymin=322, ymax=369), array= array([[ 1.57293516e-06, 4.83428551e-07, 1.01336786e-06, ..., -1.52099460e-06, 1.62390904e-06, -3.06897618e-07], [ -4.01491747e-07, -6.06370918e-07, -1.84310323e-07, ..., -3.39801630e-07, 4.16551785e-07, 1.06846937e-06], [ -4.68112916e-07, 2.85567381e-07, -4.94964070e-07, ..., -1.98679504e-06, -9.59963018e-07, 1.69272738e-07], ..., [ 4.87750185e-07, 5.51780090e-07, 4.79130222e-07, ..., -2.40008035e-06, 1.37110010e-06, 1.22177823e-06], [ -2.55823608e-07, -8.19296361e-07, 2.35463608e-06, ..., -1.41771136e-06, -1.60905051e-06, 4.51424427e-07], [ -2.71851150e-07, 1.25770748e-06, 1.75458837e-07, ..., -1.10357121e-06, 1.51166330e-06, -2.85283708e-09]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=100, xmax=147, ymin=322, ymax=369), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(123.12,345.43) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=100, xmax=147, ymin=322, ymax=369), moments_status=0, observed_shape=galsim.Shear((0.23000160481118787-0.16999831959428038j)), moments_sigma=4.9998674392700195, moments_amp=0.9999805688858032, moments_centroid=galsim.PositionD(x=123.11974024161614, y=345.42992834522084), moments_rho4=1.9999460194128353, moments_n_iter=8) got hsm parameters imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=100, xmax=147, ymin=322, ymax=369), array= array([[ 9.56835177e-19, 4.41491191e-18, 1.90391963e-17, ..., 6.30861288e-18, 1.38961659e-18, 2.86086748e-19], [ 4.50840689e-18, 2.08021399e-17, 8.97087102e-17, ..., 2.97248644e-17, 6.54758205e-18, 1.34798087e-18], [ 1.98541264e-17, 9.16084794e-17, 3.95059290e-16, ..., 1.30902377e-16, 2.88342489e-17, 5.93623933e-18], ..., [ 1.61987789e-17, 7.47424255e-17, 3.22324836e-16, ..., 1.06801914e-16, 2.35255675e-17, 4.84331702e-18], [ 3.64371530e-18, 1.68123847e-17, 7.25029961e-17, ..., 2.40237718e-17, 5.29178623e-18, 1.08944443e-18], [ 7.66034439e-19, 3.53454234e-18, 1.52426259e-17, ..., 5.05062420e-18, 1.11251574e-18, 2.29038738e-19]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=100, xmax=147, ymin=322, ymax=369), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(123.12,345.43) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=100, xmax=147, ymin=322, ymax=369), moments_status=0, observed_shape=galsim.Shear((-5.913742412433365e-11+3.0522600003024536e-09j)), moments_sigma=3.846153736114502, moments_amp=0.7143693566322327, moments_centroid=galsim.PositionD(x=123.12000001547273, y=345.42999999212833), moments_rho4=1.9999999871822245, moments_n_iter=6) start with_hsm for imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=100, xmax=147, ymin=322, ymax=369), array= array([[ 1.57293516e-06, 4.83428551e-07, 1.01336786e-06, ..., -1.52099460e-06, 1.62390904e-06, -3.06897618e-07], [ -4.01491747e-07, -6.06370918e-07, -1.84310323e-07, ..., -3.39801630e-07, 4.16551785e-07, 1.06846937e-06], [ -4.68112916e-07, 2.85567381e-07, -4.94964070e-07, ..., -1.98679504e-06, -9.59963018e-07, 1.69272738e-07], ..., [ 4.87750185e-07, 5.51780090e-07, 4.79130222e-07, ..., -2.40008035e-06, 1.37110010e-06, 1.22177823e-06], [ -2.55823608e-07, -8.19296361e-07, 2.35463608e-06, ..., -1.41771136e-06, -1.60905051e-06, 4.51424427e-07], [ -2.71851150e-07, 1.25770748e-06, 1.75458837e-07, ..., -1.10357121e-06, 1.51166330e-06, -2.85283708e-09]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=100, xmax=147, ymin=322, ymax=369), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(123.12,345.43) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=100, xmax=147, ymin=322, ymax=369), moments_status=0, observed_shape=galsim.Shear((0.23000160481118787-0.16999831959428038j)), moments_sigma=4.9998674392700195, moments_amp=0.9999805688858032, moments_centroid=galsim.PositionD(x=123.11974024161614, y=345.42992834522084), moments_rho4=1.9999460194128353, moments_n_iter=8) got hsm parameters imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=100, xmax=147, ymin=322, ymax=369), array= array([[ 1.47229563e-17, 3.84441123e-17, 9.77012204e-17, ..., 1.89917956e-10, 1.46552076e-10, 1.10066019e-10], [ 9.67435968e-17, 2.48899573e-16, 6.23248748e-16, ..., 6.40787590e-10, 4.87199836e-10, 3.60525165e-10], [ 5.94398753e-16, 1.50676986e-15, 3.71750309e-15, ..., 2.02157846e-09, 1.51443558e-09, 1.10419651e-09], ..., [ 1.20617749e-09, 1.61721181e-09, 2.11036677e-09, ..., 1.46359808e-15, 5.79920790e-16, 2.23640706e-16], [ 3.85762394e-10, 5.09615572e-10, 6.55240695e-10, ..., 2.40354043e-16, 9.38351106e-17, 3.56545251e-17], [ 1.15360485e-10, 1.50157498e-10, 1.90227056e-10, ..., 3.69070619e-17, 1.41968016e-17, 5.31504530e-18]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=100, xmax=147, ymin=322, ymax=369), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(123.12,345.43) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=100, xmax=147, ymin=322, ymax=369), moments_status=0, observed_shape=galsim.Shear((0.22999950307438866-0.1699977333434136j)), moments_sigma=4.9998555183410645, moments_amp=0.999978244304657, moments_centroid=galsim.PositionD(x=123.11974364918903, y=345.4299275928385), moments_rho4=1.9999711387393857, moments_n_iter=8) start with_hsm for imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=323, xmax=370, ymin=544, ymax=591), array= array([[ -4.64985504e-07, -5.38259087e-07, 1.27860790e-06, ..., 9.59986892e-07, 5.81386814e-07, -1.69146315e-08], [ 5.36656898e-07, 6.98122562e-07, -2.27826763e-06, ..., -4.54603338e-07, -4.67363577e-07, 7.88441866e-07], [ -4.12713433e-07, 8.72180976e-07, -1.06762320e-06, ..., 9.32305170e-07, 9.44102894e-08, -2.08292727e-06], ..., [ 5.33941829e-07, 1.06433686e-06, 1.69580323e-06, ..., 8.69430949e-07, -1.51942322e-06, -1.98380121e-07], [ 8.97153313e-07, -1.10840438e-06, 6.98605405e-08, ..., 1.42423278e-06, 3.54525412e-07, -1.99148940e-06], [ -6.01974818e-07, 1.06053415e-06, 1.12591874e-06, ..., -1.30754131e-07, 3.08745797e-07, -9.19329295e-07]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=323, xmax=370, ymin=544, ymax=591), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(345.98,567.45) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=323, xmax=370, ymin=544, ymax=591), moments_status=0, observed_shape=galsim.Shear((0.22998069831116813-0.16999026883163862j)), moments_sigma=4.999975204467773, moments_amp=0.999998152256012, moments_centroid=galsim.PositionD(x=345.980138171224, y=567.4498072462465), moments_rho4=1.999932797336089, moments_n_iter=8) got hsm parameters start with_hsm for imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=323, xmax=370, ymin=544, ymax=591), array= array([[ -4.64985504e-07, -5.38259087e-07, 1.27860790e-06, ..., 9.59986892e-07, 5.81386814e-07, -1.69146315e-08], [ 5.36656898e-07, 6.98122562e-07, -2.27826763e-06, ..., -4.54603338e-07, -4.67363577e-07, 7.88441866e-07], [ -4.12713433e-07, 8.72180976e-07, -1.06762320e-06, ..., 9.32305170e-07, 9.44102894e-08, -2.08292727e-06], ..., [ 5.33941829e-07, 1.06433686e-06, 1.69580323e-06, ..., 8.69430949e-07, -1.51942322e-06, -1.98380121e-07], [ 8.97153313e-07, -1.10840438e-06, 6.98605405e-08, ..., 1.42423278e-06, 3.54525412e-07, -1.99148940e-06], [ -6.01974818e-07, 1.06053415e-06, 1.12591874e-06, ..., -1.30754131e-07, 3.08745797e-07, -9.19329295e-07]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=323, xmax=370, ymin=544, ymax=591), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(345.98,567.45) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=323, xmax=370, ymin=544, ymax=591), moments_status=0, observed_shape=galsim.Shear((0.22998069831116813-0.16999026883163862j)), moments_sigma=4.999975204467773, moments_amp=0.999998152256012, moments_centroid=galsim.PositionD(x=345.980138171224, y=567.4498072462465), moments_rho4=1.999932797336089, moments_n_iter=8) got hsm parameters imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=323, xmax=370, ymin=544, ymax=591), array= array([[ 1.15294895e-18, 5.26968727e-18, 2.25113366e-17, ..., 4.96534601e-18, 1.08342897e-18, 2.20949532e-19], [ 5.43980399e-18, 2.48632542e-17, 1.06212201e-16, ..., 2.34273233e-17, 5.11179710e-18, 1.04247639e-18], [ 2.39882229e-17, 1.09640957e-16, 4.68370203e-16, ..., 1.03308843e-16, 2.25417917e-17, 4.59706949e-18], ..., [ 2.07432942e-17, 9.48096286e-17, 4.05012933e-16, ..., 8.93340768e-17, 1.94925244e-17, 3.97521578e-18], [ 4.67226062e-18, 2.13551085e-17, 9.12259116e-17, ..., 2.01217838e-17, 4.39053448e-18, 8.95385480e-19], [ 9.83599019e-19, 4.49565342e-18, 1.92047764e-17, ..., 4.23601508e-18, 9.24290376e-19, 1.88495527e-19]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=323, xmax=370, ymin=544, ymax=591), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(345.98,567.45) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=323, xmax=370, ymin=544, ymax=591), moments_status=0, observed_shape=galsim.Shear((4.510695636450457e-11+1.0815878193071171e-09j)), moments_sigma=3.846153736114502, moments_amp=0.714364230632782, moments_centroid=galsim.PositionD(x=345.9800000142484, y=567.4499999910925), moments_rho4=1.9999999975422822, moments_n_iter=6) start with_hsm for imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=323, xmax=370, ymin=544, ymax=591), array= array([[ -4.64985504e-07, -5.38259087e-07, 1.27860790e-06, ..., 9.59986892e-07, 5.81386814e-07, -1.69146315e-08], [ 5.36656898e-07, 6.98122562e-07, -2.27826763e-06, ..., -4.54603338e-07, -4.67363577e-07, 7.88441866e-07], [ -4.12713433e-07, 8.72180976e-07, -1.06762320e-06, ..., 9.32305170e-07, 9.44102894e-08, -2.08292727e-06], ..., [ 5.33941829e-07, 1.06433686e-06, 1.69580323e-06, ..., 8.69430949e-07, -1.51942322e-06, -1.98380121e-07], [ 8.97153313e-07, -1.10840438e-06, 6.98605405e-08, ..., 1.42423278e-06, 3.54525412e-07, -1.99148940e-06], [ -6.01974818e-07, 1.06053415e-06, 1.12591874e-06, ..., -1.30754131e-07, 3.08745797e-07, -9.19329295e-07]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=323, xmax=370, ymin=544, ymax=591), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(345.98,567.45) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=323, xmax=370, ymin=544, ymax=591), moments_status=0, observed_shape=galsim.Shear((0.22998069831116813-0.16999026883163862j)), moments_sigma=4.999975204467773, moments_amp=0.999998152256012, moments_centroid=galsim.PositionD(x=345.980138171224, y=567.4498072462465), moments_rho4=1.999932797336089, moments_n_iter=8) got hsm parameters imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=323, xmax=370, ymin=544, ymax=591), array= array([[ 1.62732516e-17, 4.23425021e-17, 1.07229481e-16, ..., 1.79137066e-10, 1.37747480e-10, 1.03090189e-10], [ 1.06833294e-16, 2.73890274e-16, 6.83411169e-16, ..., 6.03905093e-10, 4.57545946e-10, 3.37393086e-10], [ 6.55797948e-16, 1.65656256e-15, 4.07268772e-15, ..., 1.90363258e-09, 1.42107470e-09, 1.03249076e-09], ..., [ 1.28756261e-09, 1.72036863e-09, 2.23722751e-09, ..., 1.33743580e-15, 5.28107202e-16, 2.02958266e-16], [ 4.11530948e-10, 5.41780620e-10, 6.94192481e-10, ..., 2.19511897e-16, 8.54033741e-17, 3.23390661e-17], [ 1.22989438e-10, 1.59535218e-10, 2.01409903e-10, ..., 3.36879297e-17, 1.29139495e-17, 4.81813888e-18]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=323, xmax=370, ymin=544, ymax=591), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(345.98,567.45) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=323, xmax=370, ymin=544, ymax=591), moments_status=0, observed_shape=galsim.Shear((0.22997845008461243-0.1699896373640655j)), moments_sigma=4.99996280670166, moments_amp=0.9999957084655762, moments_centroid=galsim.PositionD(x=345.9801429698122, y=567.4498061911613), moments_rho4=1.999969451271802, moments_n_iter=8) start with_hsm for imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=1072, xmax=1119, ymin=901, ymax=948), array= array([[ -1.50366873e-06, -1.14600471e-06, 1.25665451e-06, ..., -1.32513458e-06, -7.00539260e-07, 8.74302089e-07], [ -8.53509903e-07, 7.94479604e-07, -1.04325943e-06, ..., -1.24426458e-06, 8.49906655e-07, 4.03015946e-07], [ -1.13993985e-06, 1.78198775e-06, -3.41910862e-07, ..., 1.77908319e-07, -1.18304445e-06, 3.23419158e-07], ..., [ -1.69720124e-07, -7.90254887e-07, -3.06249355e-08, ..., 5.63925653e-07, -1.12936516e-06, -4.91445121e-07], [ -1.56005638e-06, 4.69508507e-07, -1.35530729e-06, ..., 1.76694653e-06, -1.54620722e-07, -8.53312031e-07], [ -8.67083600e-07, 4.62170050e-07, -1.95648198e-07, ..., -3.77680294e-07, -2.61958689e-07, -3.48743470e-07]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=1072, xmax=1119, ymin=901, ymax=948), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(1094.94,924.29) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=1072, xmax=1119, ymin=901, ymax=948), moments_status=0, observed_shape=galsim.Shear((0.229992433331852-0.17002166569114985j)), moments_sigma=4.999765396118164, moments_amp=0.9999370574951172, moments_centroid=galsim.PositionD(x=1094.9398594882296, y=924.2902140730297), moments_rho4=2.0000078873293408, moments_n_iter=8) got hsm parameters start with_hsm for imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=1072, xmax=1119, ymin=901, ymax=948), array= array([[ -1.50366873e-06, -1.14600471e-06, 1.25665451e-06, ..., -1.32513458e-06, -7.00539260e-07, 8.74302089e-07], [ -8.53509903e-07, 7.94479604e-07, -1.04325943e-06, ..., -1.24426458e-06, 8.49906655e-07, 4.03015946e-07], [ -1.13993985e-06, 1.78198775e-06, -3.41910862e-07, ..., 1.77908319e-07, -1.18304445e-06, 3.23419158e-07], ..., [ -1.69720124e-07, -7.90254887e-07, -3.06249355e-08, ..., 5.63925653e-07, -1.12936516e-06, -4.91445121e-07], [ -1.56005638e-06, 4.69508507e-07, -1.35530729e-06, ..., 1.76694653e-06, -1.54620722e-07, -8.53312031e-07], [ -8.67083600e-07, 4.62170050e-07, -1.95648198e-07, ..., -3.77680294e-07, -2.61958689e-07, -3.48743470e-07]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=1072, xmax=1119, ymin=901, ymax=948), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(1094.94,924.29) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=1072, xmax=1119, ymin=901, ymax=948), moments_status=0, observed_shape=galsim.Shear((0.229992433331852-0.17002166569114985j)), moments_sigma=4.999765396118164, moments_amp=0.9999370574951172, moments_centroid=galsim.PositionD(x=1094.9398594882296, y=924.2902140730297), moments_rho4=2.0000078873293408, moments_n_iter=8) got hsm parameters imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=1072, xmax=1119, ymin=901, ymax=948), array= array([[ 1.57958544e-18, 7.20018289e-18, 3.06750732e-17, ..., 6.02334269e-18, 1.31073282e-18, 2.66583018e-19], [ 7.37257064e-18, 3.36061943e-17, 1.43173106e-16, ..., 2.81134007e-17, 6.11772565e-18, 1.24425180e-18], [ 3.21615101e-17, 1.46600968e-16, 6.24566890e-16, ..., 1.22639651e-16, 2.66874770e-17, 5.42782435e-18], ..., [ 1.74675068e-17, 7.96216794e-17, 3.39213746e-16, ..., 6.66078453e-17, 1.44944592e-17, 2.94795114e-18], [ 3.89209034e-18, 1.77412143e-17, 7.55832248e-17, ..., 1.48414863e-17, 3.22963902e-18, 6.56859213e-19], [ 8.10544046e-19, 3.69468162e-18, 1.57405224e-17, ..., 3.09080130e-18, 6.72585778e-19, 1.36793663e-19]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=1072, xmax=1119, ymin=901, ymax=948), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(1094.94,924.29) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=1072, xmax=1119, ymin=901, ymax=948), moments_status=0, observed_shape=galsim.Shear((-1.4565343375849693e-09-2.610427873861454e-09j)), moments_sigma=3.846153736114502, moments_amp=0.7143539786338806, moments_centroid=galsim.PositionD(x=1094.9400000126518, y=924.2900000043493), moments_rho4=1.9999999881742239, moments_n_iter=6) start with_hsm for imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=1072, xmax=1119, ymin=901, ymax=948), array= array([[ -1.50366873e-06, -1.14600471e-06, 1.25665451e-06, ..., -1.32513458e-06, -7.00539260e-07, 8.74302089e-07], [ -8.53509903e-07, 7.94479604e-07, -1.04325943e-06, ..., -1.24426458e-06, 8.49906655e-07, 4.03015946e-07], [ -1.13993985e-06, 1.78198775e-06, -3.41910862e-07, ..., 1.77908319e-07, -1.18304445e-06, 3.23419158e-07], ..., [ -1.69720124e-07, -7.90254887e-07, -3.06249355e-08, ..., 5.63925653e-07, -1.12936516e-06, -4.91445121e-07], [ -1.56005638e-06, 4.69508507e-07, -1.35530729e-06, ..., 1.76694653e-06, -1.54620722e-07, -8.53312031e-07], [ -8.67083600e-07, 4.62170050e-07, -1.95648198e-07, ..., -3.77680294e-07, -2.61958689e-07, -3.48743470e-07]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=1072, xmax=1119, ymin=901, ymax=948), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(1094.94,924.29) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=1072, xmax=1119, ymin=901, ymax=948), moments_status=0, observed_shape=galsim.Shear((0.229992433331852-0.17002166569114985j)), moments_sigma=4.999765396118164, moments_amp=0.9999370574951172, moments_centroid=galsim.PositionD(x=1094.9398594882296, y=924.2902140730297), moments_rho4=2.0000078873293408, moments_n_iter=8) got hsm parameters imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=1072, xmax=1119, ymin=901, ymax=948), array= array([[ 2.28303735e-17, 5.92081258e-17, 1.49445914e-16, ..., 2.16100679e-10, 1.65604849e-10, 1.23516322e-10], [ 1.48237396e-16, 3.78784373e-16, 9.42021340e-16, ..., 7.20392246e-10, 5.43941114e-10, 3.99732247e-10], [ 8.99974095e-16, 2.26584511e-15, 5.55219678e-15, ..., 2.24548091e-09, 1.67054515e-09, 1.20959909e-09], ..., [ 1.09068277e-09, 1.45222823e-09, 1.88193994e-09, ..., 9.65846365e-16, 3.80008309e-16, 1.45516608e-16], [ 3.44649975e-10, 4.52148041e-10, 5.77321135e-10, ..., 1.56695312e-16, 6.07445112e-17, 2.29188305e-17], [ 1.01832286e-10, 1.31629693e-10, 1.65598604e-10, ..., 2.37701218e-17, 9.07921712e-18, 3.37520192e-18]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=1072, xmax=1119, ymin=901, ymax=948), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(1094.94,924.29) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=1072, xmax=1119, ymin=901, ymax=948), moments_status=0, observed_shape=galsim.Shear((0.22999012972995844-0.17002101290668628j)), moments_sigma=4.999752521514893, moments_amp=0.9999344348907471, moments_centroid=galsim.PositionD(x=1094.9398646818527, y=924.2902129245838), moments_rho4=1.9999688817256012, moments_n_iter=8) start with_hsm for imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=901, xmax=948, ymin=1510, ymax=1557), array= array([[ 1.22655206e-06, -1.42661361e-06, 8.50180811e-07, ..., 1.01519606e-06, -6.49568221e-07, -1.39539134e-06], [ -8.76452361e-07, 8.22429968e-07, -1.82626911e-06, ..., 1.45980687e-06, -1.06748121e-06, 5.19500674e-08], [ 6.82359428e-07, 6.66908988e-07, 8.72911301e-07, ..., 7.77042999e-07, 1.26813586e-06, -2.53787363e-07], ..., [ -5.39905542e-08, 7.49546530e-07, 4.23453258e-07, ..., 2.20813604e-07, -4.12538640e-07, 9.22425102e-07], [ 1.45666706e-06, 7.74921460e-09, -1.05382276e-06, ..., 8.90651997e-10, 3.04599666e-07, 9.49413618e-07], [ -1.26131816e-09, 1.45115553e-06, 1.08020208e-06, ..., -1.93169990e-06, 9.34230684e-07, -6.01052079e-07]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=901, xmax=948, ymin=1510, ymax=1557), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(924.15,1532.92) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=901, xmax=948, ymin=1510, ymax=1557), moments_status=0, observed_shape=galsim.Shear((0.22997834151042335-0.17000010236816496j)), moments_sigma=4.99981164932251, moments_amp=0.9999456405639648, moments_centroid=galsim.PositionD(x=924.1498952345294, y=1532.9199801983234), moments_rho4=1.9999540838151564, moments_n_iter=8) got hsm parameters start with_hsm for imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=901, xmax=948, ymin=1510, ymax=1557), array= array([[ 1.22655206e-06, -1.42661361e-06, 8.50180811e-07, ..., 1.01519606e-06, -6.49568221e-07, -1.39539134e-06], [ -8.76452361e-07, 8.22429968e-07, -1.82626911e-06, ..., 1.45980687e-06, -1.06748121e-06, 5.19500674e-08], [ 6.82359428e-07, 6.66908988e-07, 8.72911301e-07, ..., 7.77042999e-07, 1.26813586e-06, -2.53787363e-07], ..., [ -5.39905542e-08, 7.49546530e-07, 4.23453258e-07, ..., 2.20813604e-07, -4.12538640e-07, 9.22425102e-07], [ 1.45666706e-06, 7.74921460e-09, -1.05382276e-06, ..., 8.90651997e-10, 3.04599666e-07, 9.49413618e-07], [ -1.26131816e-09, 1.45115553e-06, 1.08020208e-06, ..., -1.93169990e-06, 9.34230684e-07, -6.01052079e-07]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=901, xmax=948, ymin=1510, ymax=1557), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(924.15,1532.92) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=901, xmax=948, ymin=1510, ymax=1557), moments_status=0, observed_shape=galsim.Shear((0.22997834151042335-0.17000010236816496j)), moments_sigma=4.99981164932251, moments_amp=0.9999456405639648, moments_centroid=galsim.PositionD(x=924.1498952345294, y=1532.9199801983234), moments_rho4=1.9999540838151564, moments_n_iter=8) got hsm parameters imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=901, xmax=948, ymin=1510, ymax=1557), array= array([[ 2.02977715e-18, 9.38456175e-18, 4.05528472e-17, ..., 1.46615150e-17, 3.23609062e-18, 6.67581335e-19], [ 9.23977950e-18, 4.27196071e-17, 1.84601236e-16, ..., 6.67409023e-17, 1.47310577e-17, 3.03890730e-18], [ 3.93112193e-17, 1.81753246e-16, 7.85397536e-16, ..., 2.83953370e-16, 6.26742045e-17, 1.29292217e-17], ..., [ 7.28316813e-18, 3.36733250e-17, 1.45510175e-16, ..., 5.26078829e-17, 1.16116165e-17, 2.39538994e-18], [ 1.58274072e-18, 7.31771367e-18, 3.16215233e-17, ..., 1.14324748e-17, 2.52337670e-18, 5.20553807e-19], [ 3.21470329e-19, 1.48630025e-18, 6.42264469e-18, ..., 2.32204901e-18, 5.12522814e-19, 1.05729638e-19]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=901, xmax=948, ymin=1510, ymax=1557), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(924.15,1532.92) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=901, xmax=948, ymin=1510, ymax=1557), moments_status=0, observed_shape=galsim.Shear((-3.1837090741504426e-09+1.221435397269488e-09j)), moments_sigma=3.846153736114502, moments_amp=0.7143598198890686, moments_centroid=galsim.PositionD(x=924.1500000061219, y=1532.919999994482), moments_rho4=1.999999991270585, moments_n_iter=6) start with_hsm for imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=901, xmax=948, ymin=1510, ymax=1557), array= array([[ 1.22655206e-06, -1.42661361e-06, 8.50180811e-07, ..., 1.01519606e-06, -6.49568221e-07, -1.39539134e-06], [ -8.76452361e-07, 8.22429968e-07, -1.82626911e-06, ..., 1.45980687e-06, -1.06748121e-06, 5.19500674e-08], [ 6.82359428e-07, 6.66908988e-07, 8.72911301e-07, ..., 7.77042999e-07, 1.26813586e-06, -2.53787363e-07], ..., [ -5.39905542e-08, 7.49546530e-07, 4.23453258e-07, ..., 2.20813604e-07, -4.12538640e-07, 9.22425102e-07], [ 1.45666706e-06, 7.74921460e-09, -1.05382276e-06, ..., 8.90651997e-10, 3.04599666e-07, 9.49413618e-07], [ -1.26131816e-09, 1.45115553e-06, 1.08020208e-06, ..., -1.93169990e-06, 9.34230684e-07, -6.01052079e-07]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=901, xmax=948, ymin=1510, ymax=1557), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(924.15,1532.92) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=901, xmax=948, ymin=1510, ymax=1557), moments_status=0, observed_shape=galsim.Shear((0.22997834151042335-0.17000010236816496j)), moments_sigma=4.99981164932251, moments_amp=0.9999456405639648, moments_centroid=galsim.PositionD(x=924.1498952345294, y=1532.9199801983234), moments_rho4=1.9999540838151564, moments_n_iter=8) got hsm parameters imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=901, xmax=948, ymin=1510, ymax=1557), array= array([[ 3.76554916e-17, 9.76695758e-17, 2.46561232e-16, ..., 3.58850033e-10, 2.75041406e-10, 2.05171838e-10], [ 2.39199431e-16, 6.11305266e-16, 1.52051467e-15, ..., 1.17046761e-09, 8.83916995e-10, 6.49677978e-10], [ 1.42076161e-15, 3.57755377e-15, 8.76769987e-15, ..., 3.56972718e-09, 2.65615929e-09, 1.92356886e-09], ..., [ 6.74042822e-10, 8.97703745e-10, 1.16362520e-09, ..., 6.03848253e-16, 2.37644331e-16, 9.10252112e-17], [ 2.08417492e-10, 2.73493311e-10, 3.49296120e-10, ..., 9.58713977e-17, 3.71754059e-17, 1.40299758e-17], [ 6.02573547e-11, 7.79093526e-11, 9.80400761e-11, ..., 1.42324448e-17, 5.43767400e-18, 2.02200103e-18]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=901, xmax=948, ymin=1510, ymax=1557), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(924.15,1532.92) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=901, xmax=948, ymin=1510, ymax=1557), moments_status=0, observed_shape=galsim.Shear((0.22997627579613808-0.16999951888399664j)), moments_sigma=4.999800205230713, moments_amp=0.9999433159828186, moments_centroid=galsim.PositionD(x=924.1498983627727, y=1532.91997950834), moments_rho4=1.9999714557718207, moments_n_iter=8) start with_hsm for imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=1720, xmax=1767, ymin=866, ymax=913), array= array([[ -5.37911546e-07, 1.03728951e-06, 3.91244612e-07, ..., -7.42700934e-07, -1.16977367e-06, -2.59595254e-07], [ -8.73514125e-07, 8.34733783e-07, -1.94831369e-08, ..., 1.04750166e-06, -2.09089421e-06, 9.26570635e-07], [ -1.04949578e-07, 1.87344984e-08, -5.22215998e-07, ..., -1.56790168e-06, 4.86066597e-07, -5.73833233e-07], ..., [ -1.04590720e-06, 1.03001787e-06, -1.26062423e-06, ..., -7.01691647e-07, 1.46049240e-06, 2.37872541e-07], [ -1.98352586e-06, 2.16189301e-06, 1.11704185e-06, ..., 9.22360300e-07, -1.08025311e-07, -3.53845451e-08], [ 1.40458994e-07, 4.57331964e-07, 8.61098215e-07, ..., -2.45549427e-06, 7.74880448e-07, -5.09472670e-07]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=1720, xmax=1767, ymin=866, ymax=913), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(1743.11,888.83) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=1720, xmax=1767, ymin=866, ymax=913), moments_status=0, observed_shape=galsim.Shear((0.22997952813520628-0.17001402130371981j)), moments_sigma=5.00012731552124, moments_amp=1.0000312328338623, moments_centroid=galsim.PositionD(x=1743.1098299311714, y=888.830096802052), moments_rho4=1.9999653539793665, moments_n_iter=8) got hsm parameters start with_hsm for imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=1720, xmax=1767, ymin=866, ymax=913), array= array([[ -5.37911546e-07, 1.03728951e-06, 3.91244612e-07, ..., -7.42700934e-07, -1.16977367e-06, -2.59595254e-07], [ -8.73514125e-07, 8.34733783e-07, -1.94831369e-08, ..., 1.04750166e-06, -2.09089421e-06, 9.26570635e-07], [ -1.04949578e-07, 1.87344984e-08, -5.22215998e-07, ..., -1.56790168e-06, 4.86066597e-07, -5.73833233e-07], ..., [ -1.04590720e-06, 1.03001787e-06, -1.26062423e-06, ..., -7.01691647e-07, 1.46049240e-06, 2.37872541e-07], [ -1.98352586e-06, 2.16189301e-06, 1.11704185e-06, ..., 9.22360300e-07, -1.08025311e-07, -3.53845451e-08], [ 1.40458994e-07, 4.57331964e-07, 8.61098215e-07, ..., -2.45549427e-06, 7.74880448e-07, -5.09472670e-07]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=1720, xmax=1767, ymin=866, ymax=913), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(1743.11,888.83) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=1720, xmax=1767, ymin=866, ymax=913), moments_status=0, observed_shape=galsim.Shear((0.22997952813520628-0.17001402130371981j)), moments_sigma=5.00012731552124, moments_amp=1.0000312328338623, moments_centroid=galsim.PositionD(x=1743.1098299311714, y=888.830096802052), moments_rho4=1.9999653539793665, moments_n_iter=8) got hsm parameters imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=1720, xmax=1767, ymin=866, ymax=913), array= array([[ 2.48344084e-18, 1.14510450e-17, 4.93489766e-17, ..., 1.58832492e-17, 3.49628524e-18, 7.19309902e-19], [ 1.12363387e-17, 5.18103022e-17, 2.23279658e-16, ..., 7.18638237e-17, 1.58189574e-17, 3.25452082e-18], [ 4.75157392e-17, 2.19093158e-16, 9.44195250e-16, ..., 3.03894616e-16, 6.68945197e-17, 1.37625756e-17], ..., [ 6.77680000e-18, 3.12475517e-17, 1.34663217e-16, ..., 4.33421234e-17, 9.54064411e-18, 1.96284914e-18], [ 1.46376671e-18, 6.74936904e-18, 2.90868169e-17, ..., 9.36175720e-18, 2.06074803e-18, 4.23968993e-19], [ 2.95502209e-19, 1.36254876e-18, 5.87198684e-18, ..., 1.88993234e-18, 4.16019581e-19, 8.55899871e-20]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=1720, xmax=1767, ymin=866, ymax=913), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(1743.11,888.83) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=1720, xmax=1767, ymin=866, ymax=913), moments_status=0, observed_shape=galsim.Shear((1.9743705281793922e-10-1.1816683187504395e-09j)), moments_sigma=3.846153736114502, moments_amp=0.7143626809120178, moments_centroid=galsim.PositionD(x=1743.1099999958167, y=888.8300000162733), moments_rho4=1.9999999911962614, moments_n_iter=6) start with_hsm for imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=1720, xmax=1767, ymin=866, ymax=913), array= array([[ -5.37911546e-07, 1.03728951e-06, 3.91244612e-07, ..., -7.42700934e-07, -1.16977367e-06, -2.59595254e-07], [ -8.73514125e-07, 8.34733783e-07, -1.94831369e-08, ..., 1.04750166e-06, -2.09089421e-06, 9.26570635e-07], [ -1.04949578e-07, 1.87344984e-08, -5.22215998e-07, ..., -1.56790168e-06, 4.86066597e-07, -5.73833233e-07], ..., [ -1.04590720e-06, 1.03001787e-06, -1.26062423e-06, ..., -7.01691647e-07, 1.46049240e-06, 2.37872541e-07], [ -1.98352586e-06, 2.16189301e-06, 1.11704185e-06, ..., 9.22360300e-07, -1.08025311e-07, -3.53845451e-08], [ 1.40458994e-07, 4.57331964e-07, 8.61098215e-07, ..., -2.45549427e-06, 7.74880448e-07, -5.09472670e-07]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=1720, xmax=1767, ymin=866, ymax=913), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(1743.11,888.83) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=1720, xmax=1767, ymin=866, ymax=913), moments_status=0, observed_shape=galsim.Shear((0.22997952813520628-0.17001402130371981j)), moments_sigma=5.00012731552124, moments_amp=1.0000312328338623, moments_centroid=galsim.PositionD(x=1743.1098299311714, y=888.830096802052), moments_rho4=1.9999653539793665, moments_n_iter=8) got hsm parameters imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=1720, xmax=1767, ymin=866, ymax=913), array= array([[ 4.64873784e-17, 1.20276270e-16, 3.02873094e-16, ..., 3.97095051e-10, 3.03637088e-10, 2.25970131e-10], [ 2.93296011e-16, 7.47684294e-16, 1.85509554e-15, ..., 1.28644539e-09, 9.69212266e-10, 7.10693560e-10], [ 1.73025400e-15, 4.34599476e-15, 1.06244012e-14, ..., 3.89691390e-09, 2.89278135e-09, 2.09000017e-09], ..., [ 6.17019047e-10, 8.19726897e-10, 1.05992648e-09, ..., 4.96040599e-16, 1.94761608e-16, 7.44261116e-17], [ 1.89555413e-10, 2.48127019e-10, 3.16116994e-10, ..., 7.82493338e-17, 3.02715017e-17, 1.13978573e-17], [ 5.44510687e-11, 7.02281983e-11, 8.81561352e-11, ..., 1.15418819e-17, 4.39943826e-18, 1.63212599e-18]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=1720, xmax=1767, ymin=866, ymax=913), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(1743.11,888.83) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=1720, xmax=1767, ymin=866, ymax=913), moments_status=0, observed_shape=galsim.Shear((0.2299774250773424-0.17001341791742625j)), moments_sigma=5.000115394592285, moments_amp=1.0000288486480713, moments_centroid=galsim.PositionD(x=1743.1098334505755, y=888.8300960084571), moments_rho4=1.9999710124831105, moments_n_iter=8) start with_hsm for imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=865, xmax=912, ymin=1010, ymax=1057), array= array([[ 1.86102497e-07, -2.02780100e-07, -2.36136088e-06, ..., 1.07062846e-07, 8.24924939e-07, 7.30341299e-07], [ -1.83401448e-06, -4.21695347e-07, 1.00835814e-06, ..., 1.79026301e-07, 1.33499654e-07, -6.96232689e-07], [ 6.34924959e-07, 1.12090866e-07, 1.31596710e-06, ..., 4.33459036e-07, 9.78215695e-08, 1.06373204e-08], ..., [ -4.43534617e-07, 2.16792714e-07, -1.25495069e-06, ..., -1.43113311e-06, 8.94063987e-07, 1.57721004e-07], [ -2.08316820e-07, -2.97888789e-07, -9.34537582e-07, ..., 2.77731885e-07, -7.33344962e-07, 4.84322072e-07], [ -2.09632844e-06, 1.93295890e-07, 4.53361828e-08, ..., 4.46277426e-07, 8.89419255e-07, -6.07596178e-07]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=865, xmax=912, ymin=1010, ymax=1057), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(888.39,1033.19) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=865, xmax=912, ymin=1010, ymax=1057), moments_status=0, observed_shape=galsim.Shear((0.22998844396853807-0.17003052331984353j)), moments_sigma=4.999911785125732, moments_amp=0.9999915957450867, moments_centroid=galsim.PositionD(x=888.3901571587685, y=1033.1900065377993), moments_rho4=1.9999689146688417, moments_n_iter=8) got hsm parameters start with_hsm for imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=865, xmax=912, ymin=1010, ymax=1057), array= array([[ 1.86102497e-07, -2.02780100e-07, -2.36136088e-06, ..., 1.07062846e-07, 8.24924939e-07, 7.30341299e-07], [ -1.83401448e-06, -4.21695347e-07, 1.00835814e-06, ..., 1.79026301e-07, 1.33499654e-07, -6.96232689e-07], [ 6.34924959e-07, 1.12090866e-07, 1.31596710e-06, ..., 4.33459036e-07, 9.78215695e-08, 1.06373204e-08], ..., [ -4.43534617e-07, 2.16792714e-07, -1.25495069e-06, ..., -1.43113311e-06, 8.94063987e-07, 1.57721004e-07], [ -2.08316820e-07, -2.97888789e-07, -9.34537582e-07, ..., 2.77731885e-07, -7.33344962e-07, 4.84322072e-07], [ -2.09632844e-06, 1.93295890e-07, 4.53361828e-08, ..., 4.46277426e-07, 8.89419255e-07, -6.07596178e-07]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=865, xmax=912, ymin=1010, ymax=1057), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(888.39,1033.19) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=865, xmax=912, ymin=1010, ymax=1057), moments_status=0, observed_shape=galsim.Shear((0.22998844396853807-0.17003052331984353j)), moments_sigma=4.999911785125732, moments_amp=0.9999915957450867, moments_centroid=galsim.PositionD(x=888.3901571587685, y=1033.1900065377993), moments_rho4=1.9999689146688417, moments_n_iter=8) got hsm parameters imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=865, xmax=912, ymin=1010, ymax=1057), array= array([[ 9.13568461e-19, 4.29291973e-18, 1.88541130e-17, ..., 1.36943464e-17, 3.07205615e-18, 6.44107965e-19], [ 4.23527020e-18, 1.99018183e-17, 8.74069776e-17, ..., 6.34864892e-17, 1.42419392e-17, 2.98606080e-18], [ 1.83511293e-17, 8.62331884e-17, 3.78728310e-16, ..., 2.75082499e-16, 6.17093281e-17, 1.29383927e-17], ..., [ 7.45275008e-18, 3.50209709e-17, 1.53808928e-16, ..., 1.11716353e-16, 2.50613574e-17, 5.25453291e-18], [ 1.64942554e-18, 7.75076175e-18, 3.40406395e-17, ..., 2.47248074e-17, 5.54652229e-18, 1.16292113e-18], [ 3.41185507e-19, 1.60325366e-18, 7.04134436e-18, ..., 5.11435392e-18, 1.14730427e-18, 2.40551541e-19]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=865, xmax=912, ymin=1010, ymax=1057), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(888.39,1033.19) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=865, xmax=912, ymin=1010, ymax=1057), moments_status=0, observed_shape=galsim.Shear((6.654236051062412e-10+2.702282619893026e-09j)), moments_sigma=3.846153736114502, moments_amp=0.7143665552139282, moments_centroid=galsim.PositionD(x=888.3900000097059, y=1033.189999999904), moments_rho4=1.999999992508359, moments_n_iter=6) start with_hsm for imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=865, xmax=912, ymin=1010, ymax=1057), array= array([[ 1.86102497e-07, -2.02780100e-07, -2.36136088e-06, ..., 1.07062846e-07, 8.24924939e-07, 7.30341299e-07], [ -1.83401448e-06, -4.21695347e-07, 1.00835814e-06, ..., 1.79026301e-07, 1.33499654e-07, -6.96232689e-07], [ 6.34924959e-07, 1.12090866e-07, 1.31596710e-06, ..., 4.33459036e-07, 9.78215695e-08, 1.06373204e-08], ..., [ -4.43534617e-07, 2.16792714e-07, -1.25495069e-06, ..., -1.43113311e-06, 8.94063987e-07, 1.57721004e-07], [ -2.08316820e-07, -2.97888789e-07, -9.34537582e-07, ..., 2.77731885e-07, -7.33344962e-07, 4.84322072e-07], [ -2.09632844e-06, 1.93295890e-07, 4.53361828e-08, ..., 4.46277426e-07, 8.89419255e-07, -6.07596178e-07]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=865, xmax=912, ymin=1010, ymax=1057), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(888.39,1033.19) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=865, xmax=912, ymin=1010, ymax=1057), moments_status=0, observed_shape=galsim.Shear((0.22998844396853807-0.17003052331984353j)), moments_sigma=4.999911785125732, moments_amp=0.9999915957450867, moments_centroid=galsim.PositionD(x=888.3901571587685, y=1033.1900065377993), moments_rho4=1.9999689146688417, moments_n_iter=8) got hsm parameters imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=865, xmax=912, ymin=1010, ymax=1057), array= array([[ 1.78633491e-17, 4.68244015e-17, 1.19458354e-16, ..., 2.73809697e-10, 2.12094453e-10, 1.59898622e-10], [ 1.15968485e-16, 2.99512756e-16, 7.52879677e-16, ..., 9.12631637e-10, 6.96533498e-10, 5.17396403e-10], [ 7.03957277e-16, 1.79137960e-15, 4.43673910e-15, ..., 2.84427437e-09, 2.13886797e-09, 1.56542168e-09], ..., [ 8.50910953e-10, 1.14515197e-09, 1.49995083e-09, ..., 1.22040412e-15, 4.85348608e-16, 1.87861955e-16], [ 2.68891603e-10, 3.56551455e-10, 4.60152860e-10, ..., 1.98000556e-16, 7.75858489e-17, 2.95892233e-17], [ 7.94510430e-11, 1.03803161e-10, 1.31994693e-10, ..., 3.00371283e-17, 1.15968629e-17, 4.35770203e-18]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=865, xmax=912, ymin=1010, ymax=1057), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(888.39,1033.19) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=865, xmax=912, ymin=1010, ymax=1057), moments_status=0, observed_shape=galsim.Shear((0.22998648880592482-0.17002998228987015j)), moments_sigma=4.999900817871094, moments_amp=0.9999894499778748, moments_centroid=galsim.PositionD(x=888.3901581391518, y=1033.1900063154787), moments_rho4=1.9999728878179026, moments_n_iter=8) start with_hsm for imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=1386, xmax=1433, ymin=100, ymax=147), array= array([[ 6.79682955e-07, 1.15071759e-06, -4.70900972e-07, ..., 7.57913881e-07, -1.55136479e-06, -3.53039326e-07], [ 2.38386217e-08, -1.17431443e-06, 6.36037896e-07, ..., -3.79154898e-07, -2.48615635e-07, 1.47407411e-06], [ 6.88235559e-07, -4.41434395e-07, -3.82128604e-07, ..., -3.33914443e-07, 4.39471748e-08, 2.89131037e-07], ..., [ -2.67208947e-07, 1.41121552e-06, 3.05258254e-07, ..., 7.14137400e-07, -8.23862479e-07, -5.45802152e-07], [ 1.87825393e-07, -1.74170270e-06, -1.01417974e-07, ..., -7.47040900e-08, -6.94858898e-07, -1.44421804e-07], [ -1.40851830e-06, 6.48677826e-07, 1.15912587e-06, ..., 2.23111806e-06, 2.72399376e-07, 2.64195052e-07]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=1386, xmax=1433, ymin=100, ymax=147), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(1409.31,123.11) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=1386, xmax=1433, ymin=100, ymax=147), moments_status=0, observed_shape=galsim.Shear((0.22998119832422834-0.1700170269405822j)), moments_sigma=4.999868392944336, moments_amp=0.9999514818191528, moments_centroid=galsim.PositionD(x=1409.3099308144467, y=123.10999759606494), moments_rho4=1.9999471805643343, moments_n_iter=8) got hsm parameters start with_hsm for imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=1386, xmax=1433, ymin=100, ymax=147), array= array([[ 6.79682955e-07, 1.15071759e-06, -4.70900972e-07, ..., 7.57913881e-07, -1.55136479e-06, -3.53039326e-07], [ 2.38386217e-08, -1.17431443e-06, 6.36037896e-07, ..., -3.79154898e-07, -2.48615635e-07, 1.47407411e-06], [ 6.88235559e-07, -4.41434395e-07, -3.82128604e-07, ..., -3.33914443e-07, 4.39471748e-08, 2.89131037e-07], ..., [ -2.67208947e-07, 1.41121552e-06, 3.05258254e-07, ..., 7.14137400e-07, -8.23862479e-07, -5.45802152e-07], [ 1.87825393e-07, -1.74170270e-06, -1.01417974e-07, ..., -7.47040900e-08, -6.94858898e-07, -1.44421804e-07], [ -1.40851830e-06, 6.48677826e-07, 1.15912587e-06, ..., 2.23111806e-06, 2.72399376e-07, 2.64195052e-07]], dtype=float32), wcs=galsim.PixelScale(0.26))Reading image file data/simple_image.fits Reading star catalog data/simple_cat.fits. Removing objects with flag (col flag) != 0 Removing objects with use (col use) == 0 Processing catalog data/simple_cat.fits with 7 stars Read a total of 7 stars from 1 image weight = galsim.Image(bounds=galsim.BoundsI(xmin=1386, xmax=1433, ymin=100, ymax=147), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(1409.31,123.11) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=1386, xmax=1433, ymin=100, ymax=147), moments_status=0, observed_shape=galsim.Shear((0.22998119832422834-0.1700170269405822j)), moments_sigma=4.999868392944336, moments_amp=0.9999514818191528, moments_centroid=galsim.PositionD(x=1409.3099308144467, y=123.10999759606494), moments_rho4=1.9999471805643343, moments_n_iter=8) got hsm parameters imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=1386, xmax=1433, ymin=100, ymax=147), array= array([[ 1.17474660e-18, 5.49044027e-18, 2.39834616e-17, ..., 1.38055749e-17, 3.08030459e-18, 6.42354136e-19], [ 5.41670911e-18, 2.53162001e-17, 1.10586775e-16, ..., 6.36569480e-17, 1.42031593e-17, 2.96186908e-18], [ 2.33436374e-17, 1.09101702e-16, 4.76580441e-16, ..., 2.74333484e-16, 6.12093801e-17, 1.27643547e-17], ..., [ 7.51328233e-18, 3.51150048e-17, 1.53390123e-16, ..., 8.82957931e-17, 1.97005884e-17, 4.10828019e-18], [ 1.65385427e-18, 7.72965706e-18, 3.37648542e-17, ..., 1.94360279e-17, 4.33657336e-18, 9.04331335e-19], [ 3.40256506e-19, 1.59026465e-18, 6.94662887e-18, ..., 3.99868083e-18, 8.92186979e-19, 1.86053031e-19]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=1386, xmax=1433, ymin=100, ymax=147), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(1409.31,123.11) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=1386, xmax=1433, ymin=100, ymax=147), moments_status=0, observed_shape=galsim.Shear((-1.1230435470466205e-09-1.7155336040630687e-09j)), moments_sigma=3.846153736114502, moments_amp=0.7143504023551941, moments_centroid=galsim.PositionD(x=1409.3099999985277, y=123.10999999779726), moments_rho4=1.9999999930382648, moments_n_iter=6) start with_hsm for imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=1386, xmax=1433, ymin=100, ymax=147), array= array([[ 6.79682955e-07, 1.15071759e-06, -4.70900972e-07, ..., 7.57913881e-07, -1.55136479e-06, -3.53039326e-07], [ 2.38386217e-08, -1.17431443e-06, 6.36037896e-07, ..., -3.79154898e-07, -2.48615635e-07, 1.47407411e-06], [ 6.88235559e-07, -4.41434395e-07, -3.82128604e-07, ..., -3.33914443e-07, 4.39471748e-08, 2.89131037e-07], ..., [ -2.67208947e-07, 1.41121552e-06, 3.05258254e-07, ..., 7.14137400e-07, -8.23862479e-07, -5.45802152e-07], [ 1.87825393e-07, -1.74170270e-06, -1.01417974e-07, ..., -7.47040900e-08, -6.94858898e-07, -1.44421804e-07], [ -1.40851830e-06, 6.48677826e-07, 1.15912587e-06, ..., 2.23111806e-06, 2.72399376e-07, 2.64195052e-07]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=1386, xmax=1433, ymin=100, ymax=147), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(1409.31,123.11) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=1386, xmax=1433, ymin=100, ymax=147), moments_status=0, observed_shape=galsim.Shear((0.22998119832422834-0.1700170269405822j)), moments_sigma=4.999868392944336, moments_amp=0.9999514818191528, moments_centroid=galsim.PositionD(x=1409.3099308144467, y=123.10999759606494), moments_rho4=1.9999471805643343, moments_n_iter=8) got hsm parameters imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=1386, xmax=1433, ymin=100, ymax=147), array= array([[ 2.25019828e-17, 5.87849534e-17, 1.49467315e-16, ..., 2.96212221e-10, 2.28670041e-10, 1.71810718e-10], [ 1.45122036e-16, 3.73546930e-16, 9.35818623e-16, ..., 9.80856840e-10, 7.46067930e-10, 5.52313695e-10], [ 8.75135542e-16, 2.21949208e-15, 5.47856167e-15, ..., 3.03695025e-09, 2.27602426e-09, 1.66016301e-09], ..., [ 7.96903987e-10, 1.06891085e-09, 1.39544154e-09, ..., 9.83654491e-16, 3.89887201e-16, 1.50407693e-16], [ 2.50174020e-10, 3.30631439e-10, 4.25285668e-10, ..., 1.58550962e-16, 6.19201135e-17, 2.35358163e-17], [ 7.34357991e-11, 9.56260487e-11, 1.21193486e-10, ..., 2.38959674e-17, 9.19504557e-18, 3.44364098e-18]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=1386, xmax=1433, ymin=100, ymax=147), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(1409.31,123.11) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=1386, xmax=1433, ymin=100, ymax=147), moments_status=0, observed_shape=galsim.Shear((0.22997922384394187-0.17001646742819265j)), moments_sigma=4.999857425689697, moments_amp=0.9999492764472961, moments_centroid=galsim.PositionD(x=1409.309932470826, y=123.10999723787525), moments_rho4=1.9999725819747456, moments_n_iter=8) mean = [ 1.29997811 0.22998814 -0.17001087] start with_hsm for imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=100, xmax=147, ymin=322, ymax=369), array= array([[ 1.57293516e-06, 4.83428551e-07, 1.01336786e-06, ..., -1.52099460e-06, 1.62390904e-06, -3.06897618e-07], [ -4.01491747e-07, -6.06370918e-07, -1.84310323e-07, ..., -3.39801630e-07, 4.16551785e-07, 1.06846937e-06], [ -4.68112916e-07, 2.85567381e-07, -4.94964070e-07, ..., -1.98679504e-06, -9.59963018e-07, 1.69272738e-07], ..., [ 4.87750185e-07, 5.51780090e-07, 4.79130222e-07, ..., -2.40008035e-06, 1.37110010e-06, 1.22177823e-06], [ -2.55823608e-07, -8.19296361e-07, 2.35463608e-06, ..., -1.41771136e-06, -1.60905051e-06, 4.51424427e-07], [ -2.71851150e-07, 1.25770748e-06, 1.75458837e-07, ..., -1.10357121e-06, 1.51166330e-06, -2.85283708e-09]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=100, xmax=147, ymin=322, ymax=369), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(123.12,345.43) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=100, xmax=147, ymin=322, ymax=369), moments_status=0, observed_shape=galsim.Shear((0.23000160481118787-0.16999831959428038j)), moments_sigma=4.9998674392700195, moments_amp=0.9999805688858032, moments_centroid=galsim.PositionD(x=123.11974024161614, y=345.42992834522084), moments_rho4=1.9999460194128353, moments_n_iter=8) got hsm parameters start with_hsm for imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=100, xmax=147, ymin=322, ymax=369), array= array([[ 1.57293516e-06, 4.83428551e-07, 1.01336786e-06, ..., -1.52099460e-06, 1.62390904e-06, -3.06897618e-07], [ -4.01491747e-07, -6.06370918e-07, -1.84310323e-07, ..., -3.39801630e-07, 4.16551785e-07, 1.06846937e-06], [ -4.68112916e-07, 2.85567381e-07, -4.94964070e-07, ..., -1.98679504e-06, -9.59963018e-07, 1.69272738e-07], ..., [ 4.87750185e-07, 5.51780090e-07, 4.79130222e-07, ..., -2.40008035e-06, 1.37110010e-06, 1.22177823e-06], [ -2.55823608e-07, -8.19296361e-07, 2.35463608e-06, ..., -1.41771136e-06, -1.60905051e-06, 4.51424427e-07], [ -2.71851150e-07, 1.25770748e-06, 1.75458837e-07, ..., -1.10357121e-06, 1.51166330e-06, -2.85283708e-09]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=100, xmax=147, ymin=322, ymax=369), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(123.12,345.43) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=100, xmax=147, ymin=322, ymax=369), moments_status=0, observed_shape=galsim.Shear((0.23000160481118787-0.16999831959428038j)), moments_sigma=4.9998674392700195, moments_amp=0.9999805688858032, moments_centroid=galsim.PositionD(x=123.11974024161614, y=345.42992834522084), moments_rho4=1.9999460194128353, moments_n_iter=8) got hsm parameters imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=100, xmax=147, ymin=322, ymax=369), array= array([[ 9.56835177e-19, 4.41491191e-18, 1.90391963e-17, ..., 6.30861288e-18, 1.38961659e-18, 2.86086748e-19], [ 4.50840689e-18, 2.08021399e-17, 8.97087102e-17, ..., 2.97248644e-17, 6.54758205e-18, 1.34798087e-18], [ 1.98541264e-17, 9.16084794e-17, 3.95059290e-16, ..., 1.30902377e-16, 2.88342489e-17, 5.93623933e-18], ..., [ 1.61987789e-17, 7.47424255e-17, 3.22324836e-16, ..., 1.06801914e-16, 2.35255675e-17, 4.84331702e-18], [ 3.64371530e-18, 1.68123847e-17, 7.25029961e-17, ..., 2.40237718e-17, 5.29178623e-18, 1.08944443e-18], [ 7.66034439e-19, 3.53454234e-18, 1.52426259e-17, ..., 5.05062420e-18, 1.11251574e-18, 2.29038738e-19]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=100, xmax=147, ymin=322, ymax=369), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(123.12,345.43) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=100, xmax=147, ymin=322, ymax=369), moments_status=0, observed_shape=galsim.Shear((-5.913742412433365e-11+3.0522600003024536e-09j)), moments_sigma=3.846153736114502, moments_amp=0.7143693566322327, moments_centroid=galsim.PositionD(x=123.12000001547273, y=345.42999999212833), moments_rho4=1.9999999871822245, moments_n_iter=6) start with_hsm for imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=323, xmax=370, ymin=544, ymax=591), array= array([[ -4.64985504e-07, -5.38259087e-07, 1.27860790e-06, ..., 9.59986892e-07, 5.81386814e-07, -1.69146315e-08], [ 5.36656898e-07, 6.98122562e-07, -2.27826763e-06, ..., -4.54603338e-07, -4.67363577e-07, 7.88441866e-07], [ -4.12713433e-07, 8.72180976e-07, -1.06762320e-06, ..., 9.32305170e-07, 9.44102894e-08, -2.08292727e-06], ..., [ 5.33941829e-07, 1.06433686e-06, 1.69580323e-06, ..., 8.69430949e-07, -1.51942322e-06, -1.98380121e-07], [ 8.97153313e-07, -1.10840438e-06, 6.98605405e-08, ..., 1.42423278e-06, 3.54525412e-07, -1.99148940e-06], [ -6.01974818e-07, 1.06053415e-06, 1.12591874e-06, ..., -1.30754131e-07, 3.08745797e-07, -9.19329295e-07]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=323, xmax=370, ymin=544, ymax=591), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(345.98,567.45) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=323, xmax=370, ymin=544, ymax=591), moments_status=0, observed_shape=galsim.Shear((0.22998069831116813-0.16999026883163862j)), moments_sigma=4.999975204467773, moments_amp=0.999998152256012, moments_centroid=galsim.PositionD(x=345.980138171224, y=567.4498072462465), moments_rho4=1.999932797336089, moments_n_iter=8) got hsm parameters start with_hsm for imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=323, xmax=370, ymin=544, ymax=591), array= array([[ -4.64985504e-07, -5.38259087e-07, 1.27860790e-06, ..., 9.59986892e-07, 5.81386814e-07, -1.69146315e-08], [ 5.36656898e-07, 6.98122562e-07, -2.27826763e-06, ..., -4.54603338e-07, -4.67363577e-07, 7.88441866e-07], [ -4.12713433e-07, 8.72180976e-07, -1.06762320e-06, ..., 9.32305170e-07, 9.44102894e-08, -2.08292727e-06], ..., [ 5.33941829e-07, 1.06433686e-06, 1.69580323e-06, ..., 8.69430949e-07, -1.51942322e-06, -1.98380121e-07], [ 8.97153313e-07, -1.10840438e-06, 6.98605405e-08, ..., 1.42423278e-06, 3.54525412e-07, -1.99148940e-06], [ -6.01974818e-07, 1.06053415e-06, 1.12591874e-06, ..., -1.30754131e-07, 3.08745797e-07, -9.19329295e-07]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=323, xmax=370, ymin=544, ymax=591), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(345.98,567.45) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=323, xmax=370, ymin=544, ymax=591), moments_status=0, observed_shape=galsim.Shear((0.22998069831116813-0.16999026883163862j)), moments_sigma=4.999975204467773, moments_amp=0.999998152256012, moments_centroid=galsim.PositionD(x=345.980138171224, y=567.4498072462465), moments_rho4=1.999932797336089, moments_n_iter=8) got hsm parameters imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=323, xmax=370, ymin=544, ymax=591), array= array([[ 1.15294895e-18, 5.26968727e-18, 2.25113366e-17, ..., 4.96534601e-18, 1.08342897e-18, 2.20949532e-19], [ 5.43980399e-18, 2.48632542e-17, 1.06212201e-16, ..., 2.34273233e-17, 5.11179710e-18, 1.04247639e-18], [ 2.39882229e-17, 1.09640957e-16, 4.68370203e-16, ..., 1.03308843e-16, 2.25417917e-17, 4.59706949e-18], ..., [ 2.07432942e-17, 9.48096286e-17, 4.05012933e-16, ..., 8.93340768e-17, 1.94925244e-17, 3.97521578e-18], [ 4.67226062e-18, 2.13551085e-17, 9.12259116e-17, ..., 2.01217838e-17, 4.39053448e-18, 8.95385480e-19], [ 9.83599019e-19, 4.49565342e-18, 1.92047764e-17, ..., 4.23601508e-18, 9.24290376e-19, 1.88495527e-19]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=323, xmax=370, ymin=544, ymax=591), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(345.98,567.45) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=323, xmax=370, ymin=544, ymax=591), moments_status=0, observed_shape=galsim.Shear((4.510695636450457e-11+1.0815878193071171e-09j)), moments_sigma=3.846153736114502, moments_amp=0.714364230632782, moments_centroid=galsim.PositionD(x=345.9800000142484, y=567.4499999910925), moments_rho4=1.9999999975422822, moments_n_iter=6) start with_hsm for imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=1072, xmax=1119, ymin=901, ymax=948), array= array([[ -1.50366873e-06, -1.14600471e-06, 1.25665451e-06, ..., -1.32513458e-06, -7.00539260e-07, 8.74302089e-07], [ -8.53509903e-07, 7.94479604e-07, -1.04325943e-06, ..., -1.24426458e-06, 8.49906655e-07, 4.03015946e-07], [ -1.13993985e-06, 1.78198775e-06, -3.41910862e-07, ..., 1.77908319e-07, -1.18304445e-06, 3.23419158e-07], ..., [ -1.69720124e-07, -7.90254887e-07, -3.06249355e-08, ..., 5.63925653e-07, -1.12936516e-06, -4.91445121e-07], [ -1.56005638e-06, 4.69508507e-07, -1.35530729e-06, ..., 1.76694653e-06, -1.54620722e-07, -8.53312031e-07], [ -8.67083600e-07, 4.62170050e-07, -1.95648198e-07, ..., -3.77680294e-07, -2.61958689e-07, -3.48743470e-07]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=1072, xmax=1119, ymin=901, ymax=948), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(1094.94,924.29) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=1072, xmax=1119, ymin=901, ymax=948), moments_status=0, observed_shape=galsim.Shear((0.229992433331852-0.17002166569114985j)), moments_sigma=4.999765396118164, moments_amp=0.9999370574951172, moments_centroid=galsim.PositionD(x=1094.9398594882296, y=924.2902140730297), moments_rho4=2.0000078873293408, moments_n_iter=8) got hsm parameters start with_hsm for imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=1072, xmax=1119, ymin=901, ymax=948), array= array([[ -1.50366873e-06, -1.14600471e-06, 1.25665451e-06, ..., -1.32513458e-06, -7.00539260e-07, 8.74302089e-07], [ -8.53509903e-07, 7.94479604e-07, -1.04325943e-06, ..., -1.24426458e-06, 8.49906655e-07, 4.03015946e-07], [ -1.13993985e-06, 1.78198775e-06, -3.41910862e-07, ..., 1.77908319e-07, -1.18304445e-06, 3.23419158e-07], ..., [ -1.69720124e-07, -7.90254887e-07, -3.06249355e-08, ..., 5.63925653e-07, -1.12936516e-06, -4.91445121e-07], [ -1.56005638e-06, 4.69508507e-07, -1.35530729e-06, ..., 1.76694653e-06, -1.54620722e-07, -8.53312031e-07], [ -8.67083600e-07, 4.62170050e-07, -1.95648198e-07, ..., -3.77680294e-07, -2.61958689e-07, -3.48743470e-07]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=1072, xmax=1119, ymin=901, ymax=948), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(1094.94,924.29) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=1072, xmax=1119, ymin=901, ymax=948), moments_status=0, observed_shape=galsim.Shear((0.229992433331852-0.17002166569114985j)), moments_sigma=4.999765396118164, moments_amp=0.9999370574951172, moments_centroid=galsim.PositionD(x=1094.9398594882296, y=924.2902140730297), moments_rho4=2.0000078873293408, moments_n_iter=8) got hsm parameters imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=1072, xmax=1119, ymin=901, ymax=948), array= array([[ 1.57958544e-18, 7.20018289e-18, 3.06750732e-17, ..., 6.02334269e-18, 1.31073282e-18, 2.66583018e-19], [ 7.37257064e-18, 3.36061943e-17, 1.43173106e-16, ..., 2.81134007e-17, 6.11772565e-18, 1.24425180e-18], [ 3.21615101e-17, 1.46600968e-16, 6.24566890e-16, ..., 1.22639651e-16, 2.66874770e-17, 5.42782435e-18], ..., [ 1.74675068e-17, 7.96216794e-17, 3.39213746e-16, ..., 6.66078453e-17, 1.44944592e-17, 2.94795114e-18], [ 3.89209034e-18, 1.77412143e-17, 7.55832248e-17, ..., 1.48414863e-17, 3.22963902e-18, 6.56859213e-19], [ 8.10544046e-19, 3.69468162e-18, 1.57405224e-17, ..., 3.09080130e-18, 6.72585778e-19, 1.36793663e-19]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=1072, xmax=1119, ymin=901, ymax=948), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(1094.94,924.29) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=1072, xmax=1119, ymin=901, ymax=948), moments_status=0, observed_shape=galsim.Shear((-1.4565343375849693e-09-2.610427873861454e-09j)), moments_sigma=3.846153736114502, moments_amp=0.7143539786338806, moments_centroid=galsim.PositionD(x=1094.9400000126518, y=924.2900000043493), moments_rho4=1.9999999881742239, moments_n_iter=6) start with_hsm for imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=901, xmax=948, ymin=1510, ymax=1557), array= array([[ 1.22655206e-06, -1.42661361e-06, 8.50180811e-07, ..., 1.01519606e-06, -6.49568221e-07, -1.39539134e-06], [ -8.76452361e-07, 8.22429968e-07, -1.82626911e-06, ..., 1.45980687e-06, -1.06748121e-06, 5.19500674e-08], [ 6.82359428e-07, 6.66908988e-07, 8.72911301e-07, ..., 7.77042999e-07, 1.26813586e-06, -2.53787363e-07], ..., [ -5.39905542e-08, 7.49546530e-07, 4.23453258e-07, ..., 2.20813604e-07, -4.12538640e-07, 9.22425102e-07], [ 1.45666706e-06, 7.74921460e-09, -1.05382276e-06, ..., 8.90651997e-10, 3.04599666e-07, 9.49413618e-07], [ -1.26131816e-09, 1.45115553e-06, 1.08020208e-06, ..., -1.93169990e-06, 9.34230684e-07, -6.01052079e-07]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=901, xmax=948, ymin=1510, ymax=1557), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(924.15,1532.92) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=901, xmax=948, ymin=1510, ymax=1557), moments_status=0, observed_shape=galsim.Shear((0.22997834151042335-0.17000010236816496j)), moments_sigma=4.99981164932251, moments_amp=0.9999456405639648, moments_centroid=galsim.PositionD(x=924.1498952345294, y=1532.9199801983234), moments_rho4=1.9999540838151564, moments_n_iter=8) got hsm parameters start with_hsm for imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=901, xmax=948, ymin=1510, ymax=1557), array= array([[ 1.22655206e-06, -1.42661361e-06, 8.50180811e-07, ..., 1.01519606e-06, -6.49568221e-07, -1.39539134e-06], [ -8.76452361e-07, 8.22429968e-07, -1.82626911e-06, ..., 1.45980687e-06, -1.06748121e-06, 5.19500674e-08], [ 6.82359428e-07, 6.66908988e-07, 8.72911301e-07, ..., 7.77042999e-07, 1.26813586e-06, -2.53787363e-07], ..., [ -5.39905542e-08, 7.49546530e-07, 4.23453258e-07, ..., 2.20813604e-07, -4.12538640e-07, 9.22425102e-07], [ 1.45666706e-06, 7.74921460e-09, -1.05382276e-06, ..., 8.90651997e-10, 3.04599666e-07, 9.49413618e-07], [ -1.26131816e-09, 1.45115553e-06, 1.08020208e-06, ..., -1.93169990e-06, 9.34230684e-07, -6.01052079e-07]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=901, xmax=948, ymin=1510, ymax=1557), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(924.15,1532.92) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=901, xmax=948, ymin=1510, ymax=1557), moments_status=0, observed_shape=galsim.Shear((0.22997834151042335-0.17000010236816496j)), moments_sigma=4.99981164932251, moments_amp=0.9999456405639648, moments_centroid=galsim.PositionD(x=924.1498952345294, y=1532.9199801983234), moments_rho4=1.9999540838151564, moments_n_iter=8) got hsm parameters imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=901, xmax=948, ymin=1510, ymax=1557), array= array([[ 2.02977715e-18, 9.38456175e-18, 4.05528472e-17, ..., 1.46615150e-17, 3.23609062e-18, 6.67581335e-19], [ 9.23977950e-18, 4.27196071e-17, 1.84601236e-16, ..., 6.67409023e-17, 1.47310577e-17, 3.03890730e-18], [ 3.93112193e-17, 1.81753246e-16, 7.85397536e-16, ..., 2.83953370e-16, 6.26742045e-17, 1.29292217e-17], ..., [ 7.28316813e-18, 3.36733250e-17, 1.45510175e-16, ..., 5.26078829e-17, 1.16116165e-17, 2.39538994e-18], [ 1.58274072e-18, 7.31771367e-18, 3.16215233e-17, ..., 1.14324748e-17, 2.52337670e-18, 5.20553807e-19], [ 3.21470329e-19, 1.48630025e-18, 6.42264469e-18, ..., 2.32204901e-18, 5.12522814e-19, 1.05729638e-19]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=901, xmax=948, ymin=1510, ymax=1557), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(924.15,1532.92) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=901, xmax=948, ymin=1510, ymax=1557), moments_status=0, observed_shape=galsim.Shear((-3.1837090741504426e-09+1.221435397269488e-09j)), moments_sigma=3.846153736114502, moments_amp=0.7143598198890686, moments_centroid=galsim.PositionD(x=924.1500000061219, y=1532.919999994482), moments_rho4=1.999999991270585, moments_n_iter=6) start with_hsm for imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=1720, xmax=1767, ymin=866, ymax=913), array= array([[ -5.37911546e-07, 1.03728951e-06, 3.91244612e-07, ..., -7.42700934e-07, -1.16977367e-06, -2.59595254e-07], [ -8.73514125e-07, 8.34733783e-07, -1.94831369e-08, ..., 1.04750166e-06, -2.09089421e-06, 9.26570635e-07], [ -1.04949578e-07, 1.87344984e-08, -5.22215998e-07, ..., -1.56790168e-06, 4.86066597e-07, -5.73833233e-07], ..., [ -1.04590720e-06, 1.03001787e-06, -1.26062423e-06, ..., -7.01691647e-07, 1.46049240e-06, 2.37872541e-07], [ -1.98352586e-06, 2.16189301e-06, 1.11704185e-06, ..., 9.22360300e-07, -1.08025311e-07, -3.53845451e-08], [ 1.40458994e-07, 4.57331964e-07, 8.61098215e-07, ..., -2.45549427e-06, 7.74880448e-07, -5.09472670e-07]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=1720, xmax=1767, ymin=866, ymax=913), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(1743.11,888.83) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=1720, xmax=1767, ymin=866, ymax=913), moments_status=0, observed_shape=galsim.Shear((0.22997952813520628-0.17001402130371981j)), moments_sigma=5.00012731552124, moments_amp=1.0000312328338623, moments_centroid=galsim.PositionD(x=1743.1098299311714, y=888.830096802052), moments_rho4=1.9999653539793665, moments_n_iter=8) got hsm parameters start with_hsm for imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=1720, xmax=1767, ymin=866, ymax=913), array= array([[ -5.37911546e-07, 1.03728951e-06, 3.91244612e-07, ..., -7.42700934e-07, -1.16977367e-06, -2.59595254e-07], [ -8.73514125e-07, 8.34733783e-07, -1.94831369e-08, ..., 1.04750166e-06, -2.09089421e-06, 9.26570635e-07], [ -1.04949578e-07, 1.87344984e-08, -5.22215998e-07, ..., -1.56790168e-06, 4.86066597e-07, -5.73833233e-07], ..., [ -1.04590720e-06, 1.03001787e-06, -1.26062423e-06, ..., -7.01691647e-07, 1.46049240e-06, 2.37872541e-07], [ -1.98352586e-06, 2.16189301e-06, 1.11704185e-06, ..., 9.22360300e-07, -1.08025311e-07, -3.53845451e-08], [ 1.40458994e-07, 4.57331964e-07, 8.61098215e-07, ..., -2.45549427e-06, 7.74880448e-07, -5.09472670e-07]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=1720, xmax=1767, ymin=866, ymax=913), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(1743.11,888.83) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=1720, xmax=1767, ymin=866, ymax=913), moments_status=0, observed_shape=galsim.Shear((0.22997952813520628-0.17001402130371981j)), moments_sigma=5.00012731552124, moments_amp=1.0000312328338623, moments_centroid=galsim.PositionD(x=1743.1098299311714, y=888.830096802052), moments_rho4=1.9999653539793665, moments_n_iter=8) got hsm parameters imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=1720, xmax=1767, ymin=866, ymax=913), array= array([[ 2.48344084e-18, 1.14510450e-17, 4.93489766e-17, ..., 1.58832492e-17, 3.49628524e-18, 7.19309902e-19], [ 1.12363387e-17, 5.18103022e-17, 2.23279658e-16, ..., 7.18638237e-17, 1.58189574e-17, 3.25452082e-18], [ 4.75157392e-17, 2.19093158e-16, 9.44195250e-16, ..., 3.03894616e-16, 6.68945197e-17, 1.37625756e-17], ..., [ 6.77680000e-18, 3.12475517e-17, 1.34663217e-16, ..., 4.33421234e-17, 9.54064411e-18, 1.96284914e-18], [ 1.46376671e-18, 6.74936904e-18, 2.90868169e-17, ..., 9.36175720e-18, 2.06074803e-18, 4.23968993e-19], [ 2.95502209e-19, 1.36254876e-18, 5.87198684e-18, ..., 1.88993234e-18, 4.16019581e-19, 8.55899871e-20]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=1720, xmax=1767, ymin=866, ymax=913), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(1743.11,888.83) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=1720, xmax=1767, ymin=866, ymax=913), moments_status=0, observed_shape=galsim.Shear((1.9743705281793922e-10-1.1816683187504395e-09j)), moments_sigma=3.846153736114502, moments_amp=0.7143626809120178, moments_centroid=galsim.PositionD(x=1743.1099999958167, y=888.8300000162733), moments_rho4=1.9999999911962614, moments_n_iter=6) start with_hsm for imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=865, xmax=912, ymin=1010, ymax=1057), array= array([[ 1.86102497e-07, -2.02780100e-07, -2.36136088e-06, ..., 1.07062846e-07, 8.24924939e-07, 7.30341299e-07], [ -1.83401448e-06, -4.21695347e-07, 1.00835814e-06, ..., 1.79026301e-07, 1.33499654e-07, -6.96232689e-07], [ 6.34924959e-07, 1.12090866e-07, 1.31596710e-06, ..., 4.33459036e-07, 9.78215695e-08, 1.06373204e-08], ..., [ -4.43534617e-07, 2.16792714e-07, -1.25495069e-06, ..., -1.43113311e-06, 8.94063987e-07, 1.57721004e-07], [ -2.08316820e-07, -2.97888789e-07, -9.34537582e-07, ..., 2.77731885e-07, -7.33344962e-07, 4.84322072e-07], [ -2.09632844e-06, 1.93295890e-07, 4.53361828e-08, ..., 4.46277426e-07, 8.89419255e-07, -6.07596178e-07]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=865, xmax=912, ymin=1010, ymax=1057), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(888.39,1033.19) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=865, xmax=912, ymin=1010, ymax=1057), moments_status=0, observed_shape=galsim.Shear((0.22998844396853807-0.17003052331984353j)), moments_sigma=4.999911785125732, moments_amp=0.9999915957450867, moments_centroid=galsim.PositionD(x=888.3901571587685, y=1033.1900065377993), moments_rho4=1.9999689146688417, moments_n_iter=8) got hsm parameters start with_hsm for imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=865, xmax=912, ymin=1010, ymax=1057), array= array([[ 1.86102497e-07, -2.02780100e-07, -2.36136088e-06, ..., 1.07062846e-07, 8.24924939e-07, 7.30341299e-07], [ -1.83401448e-06, -4.21695347e-07, 1.00835814e-06, ..., 1.79026301e-07, 1.33499654e-07, -6.96232689e-07], [ 6.34924959e-07, 1.12090866e-07, 1.31596710e-06, ..., 4.33459036e-07, 9.78215695e-08, 1.06373204e-08], ..., [ -4.43534617e-07, 2.16792714e-07, -1.25495069e-06, ..., -1.43113311e-06, 8.94063987e-07, 1.57721004e-07], [ -2.08316820e-07, -2.97888789e-07, -9.34537582e-07, ..., 2.77731885e-07, -7.33344962e-07, 4.84322072e-07], [ -2.09632844e-06, 1.93295890e-07, 4.53361828e-08, ..., 4.46277426e-07, 8.89419255e-07, -6.07596178e-07]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=865, xmax=912, ymin=1010, ymax=1057), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(888.39,1033.19) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=865, xmax=912, ymin=1010, ymax=1057), moments_status=0, observed_shape=galsim.Shear((0.22998844396853807-0.17003052331984353j)), moments_sigma=4.999911785125732, moments_amp=0.9999915957450867, moments_centroid=galsim.PositionD(x=888.3901571587685, y=1033.1900065377993), moments_rho4=1.9999689146688417, moments_n_iter=8) got hsm parameters imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=865, xmax=912, ymin=1010, ymax=1057), array= array([[ 9.13568461e-19, 4.29291973e-18, 1.88541130e-17, ..., 1.36943464e-17, 3.07205615e-18, 6.44107965e-19], [ 4.23527020e-18, 1.99018183e-17, 8.74069776e-17, ..., 6.34864892e-17, 1.42419392e-17, 2.98606080e-18], [ 1.83511293e-17, 8.62331884e-17, 3.78728310e-16, ..., 2.75082499e-16, 6.17093281e-17, 1.29383927e-17], ..., [ 7.45275008e-18, 3.50209709e-17, 1.53808928e-16, ..., 1.11716353e-16, 2.50613574e-17, 5.25453291e-18], [ 1.64942554e-18, 7.75076175e-18, 3.40406395e-17, ..., 2.47248074e-17, 5.54652229e-18, 1.16292113e-18], [ 3.41185507e-19, 1.60325366e-18, 7.04134436e-18, ..., 5.11435392e-18, 1.14730427e-18, 2.40551541e-19]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=865, xmax=912, ymin=1010, ymax=1057), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(888.39,1033.19) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=865, xmax=912, ymin=1010, ymax=1057), moments_status=0, observed_shape=galsim.Shear((6.654236051062412e-10+2.702282619893026e-09j)), moments_sigma=3.846153736114502, moments_amp=0.7143665552139282, moments_centroid=galsim.PositionD(x=888.3900000097059, y=1033.189999999904), moments_rho4=1.999999992508359, moments_n_iter=6) start with_hsm for imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=1386, xmax=1433, ymin=100, ymax=147), array= array([[ 6.79682955e-07, 1.15071759e-06, -4.70900972e-07, ..., 7.57913881e-07, -1.55136479e-06, -3.53039326e-07], [ 2.38386217e-08, -1.17431443e-06, 6.36037896e-07, ..., -3.79154898e-07, -2.48615635e-07, 1.47407411e-06], [ 6.88235559e-07, -4.41434395e-07, -3.82128604e-07, ..., -3.33914443e-07, 4.39471748e-08, 2.89131037e-07], ..., [ -2.67208947e-07, 1.41121552e-06, 3.05258254e-07, ..., 7.14137400e-07, -8.23862479e-07, -5.45802152e-07], [ 1.87825393e-07, -1.74170270e-06, -1.01417974e-07, ..., -7.47040900e-08, -6.94858898e-07, -1.44421804e-07], [ -1.40851830e-06, 6.48677826e-07, 1.15912587e-06, ..., 2.23111806e-06, 2.72399376e-07, 2.64195052e-07]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=1386, xmax=1433, ymin=100, ymax=147), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(1409.31,123.11) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=1386, xmax=1433, ymin=100, ymax=147), moments_status=0, observed_shape=galsim.Shear((0.22998119832422834-0.1700170269405822j)), moments_sigma=4.999868392944336, moments_amp=0.9999514818191528, moments_centroid=galsim.PositionD(x=1409.3099308144467, y=123.10999759606494), moments_rho4=1.9999471805643343, moments_n_iter=8)Iteration 1: Fitting 7 stars got hsm parameters start with_hsm for imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=1386, xmax=1433, ymin=100, ymax=147), array= array([[ 6.79682955e-07, 1.15071759e-06, -4.70900972e-07, ..., 7.57913881e-07, -1.55136479e-06, -3.53039326e-07], [ 2.38386217e-08, -1.17431443e-06, 6.36037896e-07, ..., -3.79154898e-07, -2.48615635e-07, 1.47407411e-06], [ 6.88235559e-07, -4.41434395e-07, -3.82128604e-07, ..., -3.33914443e-07, 4.39471748e-08, 2.89131037e-07], ..., [ -2.67208947e-07, 1.41121552e-06, 3.05258254e-07, ..., 7.14137400e-07, -8.23862479e-07, -5.45802152e-07], [ 1.87825393e-07, -1.74170270e-06, -1.01417974e-07, ..., -7.47040900e-08, -6.94858898e-07, -1.44421804e-07], [ -1.40851830e-06, 6.48677826e-07, 1.15912587e-06, ..., 2.23111806e-06, 2.72399376e-07, 2.64195052e-07]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=1386, xmax=1433, ymin=100, ymax=147), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(1409.31,123.11) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=1386, xmax=1433, ymin=100, ymax=147), moments_status=0, observed_shape=galsim.Shear((0.22998119832422834-0.1700170269405822j)), moments_sigma=4.999868392944336, moments_amp=0.9999514818191528, moments_centroid=galsim.PositionD(x=1409.3099308144467, y=123.10999759606494), moments_rho4=1.9999471805643343, moments_n_iter=8) got hsm parameters imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=1386, xmax=1433, ymin=100, ymax=147), array= array([[ 1.17474660e-18, 5.49044027e-18, 2.39834616e-17, ..., 1.38055749e-17, 3.08030459e-18, 6.42354136e-19], [ 5.41670911e-18, 2.53162001e-17, 1.10586775e-16, ..., 6.36569480e-17, 1.42031593e-17, 2.96186908e-18], [ 2.33436374e-17, 1.09101702e-16, 4.76580441e-16, ..., 2.74333484e-16, 6.12093801e-17, 1.27643547e-17], ..., [ 7.51328233e-18, 3.51150048e-17, 1.53390123e-16, ..., 8.82957931e-17, 1.97005884e-17, 4.10828019e-18], [ 1.65385427e-18, 7.72965706e-18, 3.37648542e-17, ..., 1.94360279e-17, 4.33657336e-18, 9.04331335e-19], [ 3.40256506e-19, 1.59026465e-18, 6.94662887e-18, ..., 3.99868083e-18, 8.92186979e-19, 1.86053031e-19]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=1386, xmax=1433, ymin=100, ymax=147), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(1409.31,123.11) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=1386, xmax=1433, ymin=100, ymax=147), moments_status=0, observed_shape=galsim.Shear((-1.1230435470466205e-09-1.7155336040630687e-09j)), moments_sigma=3.846153736114502, moments_amp=0.7143504023551941, moments_centroid=galsim.PositionD(x=1409.3099999985277, y=123.10999999779726), moments_rho4=1.9999999930382648, moments_n_iter=6) start with_hsm for imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=100, xmax=147, ymin=322, ymax=369), array= array([[ 1.57293516e-06, 4.83428551e-07, 1.01336786e-06, ..., -1.52099460e-06, 1.62390904e-06, -3.06897618e-07], [ -4.01491747e-07, -6.06370918e-07, -1.84310323e-07, ..., -3.39801630e-07, 4.16551785e-07, 1.06846937e-06], [ -4.68112916e-07, 2.85567381e-07, -4.94964070e-07, ..., -1.98679504e-06, -9.59963018e-07, 1.69272738e-07], ..., [ 4.87750185e-07, 5.51780090e-07, 4.79130222e-07, ..., -2.40008035e-06, 1.37110010e-06, 1.22177823e-06], [ -2.55823608e-07, -8.19296361e-07, 2.35463608e-06, ..., -1.41771136e-06, -1.60905051e-06, 4.51424427e-07], [ -2.71851150e-07, 1.25770748e-06, 1.75458837e-07, ..., -1.10357121e-06, 1.51166330e-06, -2.85283708e-09]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=100, xmax=147, ymin=322, ymax=369), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(123.12,345.43) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=100, xmax=147, ymin=322, ymax=369), moments_status=0, observed_shape=galsim.Shear((0.23000160481118787-0.16999831959428038j)), moments_sigma=4.9998674392700195, moments_amp=0.9999805688858032, moments_centroid=galsim.PositionD(x=123.11974024161614, y=345.42992834522084), moments_rho4=1.9999460194128353, moments_n_iter=8) got hsm parameters imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=100, xmax=147, ymin=322, ymax=369), array= array([[ 1.47255024e-17, 3.84519011e-17, 9.77238720e-17, ..., 1.90109053e-10, 1.46700527e-10, 1.10178207e-10], [ 9.67560839e-17, 2.48938907e-16, 6.23365004e-16, ..., 6.41381837e-10, 4.87654639e-10, 3.60863672e-10], [ 5.94451005e-16, 1.50694488e-15, 3.71803757e-15, ..., 2.02329953e-09, 1.51573276e-09, 1.10514753e-09], ..., [ 1.20722488e-09, 1.61860758e-09, 2.11217621e-09, ..., 1.46378622e-15, 5.79978970e-16, 2.23656747e-16], [ 3.86127269e-10, 5.10094522e-10, 6.55852261e-10, ..., 2.40395124e-16, 9.38484381e-17, 3.56585419e-17], [ 1.15478842e-10, 1.50310528e-10, 1.90419555e-10, ..., 3.69150392e-17, 1.41994486e-17, 5.31587579e-18]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=100, xmax=147, ymin=322, ymax=369), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(123.12,345.43) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=100, xmax=147, ymin=322, ymax=369), moments_status=0, observed_shape=galsim.Shear((0.22998603559062844-0.17001027777899747j)), moments_sigma=4.999904155731201, moments_amp=0.9999879002571106, moments_centroid=galsim.PositionD(x=123.11974366658731, y=345.4299275851176), moments_rho4=1.9999711398422115, moments_n_iter=8) start with_hsm for imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=323, xmax=370, ymin=544, ymax=591), array= array([[ -4.64985504e-07, -5.38259087e-07, 1.27860790e-06, ..., 9.59986892e-07, 5.81386814e-07, -1.69146315e-08], [ 5.36656898e-07, 6.98122562e-07, -2.27826763e-06, ..., -4.54603338e-07, -4.67363577e-07, 7.88441866e-07], [ -4.12713433e-07, 8.72180976e-07, -1.06762320e-06, ..., 9.32305170e-07, 9.44102894e-08, -2.08292727e-06], ..., [ 5.33941829e-07, 1.06433686e-06, 1.69580323e-06, ..., 8.69430949e-07, -1.51942322e-06, -1.98380121e-07], [ 8.97153313e-07, -1.10840438e-06, 6.98605405e-08, ..., 1.42423278e-06, 3.54525412e-07, -1.99148940e-06], [ -6.01974818e-07, 1.06053415e-06, 1.12591874e-06, ..., -1.30754131e-07, 3.08745797e-07, -9.19329295e-07]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=323, xmax=370, ymin=544, ymax=591), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(345.98,567.45) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=323, xmax=370, ymin=544, ymax=591), moments_status=0, observed_shape=galsim.Shear((0.22998069831116813-0.16999026883163862j)), moments_sigma=4.999975204467773, moments_amp=0.999998152256012, moments_centroid=galsim.PositionD(x=345.980138171224, y=567.4498072462465), moments_rho4=1.999932797336089, moments_n_iter=8) got hsm parameters imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=323, xmax=370, ymin=544, ymax=591), array= array([[ 1.62344188e-17, 4.22444778e-17, 1.06988791e-16, ..., 1.79149792e-10, 1.37762399e-10, 1.03105108e-10], [ 1.06592512e-16, 2.73291845e-16, 6.81964543e-16, ..., 6.03965267e-10, 4.57607480e-10, 3.37449985e-10], [ 6.54404526e-16, 1.65315304e-15, 4.06457357e-15, ..., 1.90386973e-09, 1.42129819e-09, 1.03268605e-09], ..., [ 1.28776068e-09, 1.72057679e-09, 2.23742291e-09, ..., 1.33466272e-15, 5.26976942e-16, 2.02510159e-16], [ 4.11584711e-10, 5.41832357e-10, 6.94233837e-10, ..., 2.19028850e-16, 8.52095359e-17, 3.22634056e-17], [ 1.23002233e-10, 1.59545835e-10, 2.01415593e-10, ..., 3.36094005e-17, 1.28829236e-17, 4.80621549e-18]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=323, xmax=370, ymin=544, ymax=591), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(345.98,567.45) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=323, xmax=370, ymin=544, ymax=591), moments_status=0, observed_shape=galsim.Shear((0.22998588902696587-0.1700102325261405j)), moments_sigma=4.999903202056885, moments_amp=0.9999837875366211, moments_centroid=galsim.PositionD(x=345.9801429637442, y=567.449806211252), moments_rho4=1.9999694516540425, moments_n_iter=8) start with_hsm for imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=1072, xmax=1119, ymin=901, ymax=948), array= array([[ -1.50366873e-06, -1.14600471e-06, 1.25665451e-06, ..., -1.32513458e-06, -7.00539260e-07, 8.74302089e-07], [ -8.53509903e-07, 7.94479604e-07, -1.04325943e-06, ..., -1.24426458e-06, 8.49906655e-07, 4.03015946e-07], [ -1.13993985e-06, 1.78198775e-06, -3.41910862e-07, ..., 1.77908319e-07, -1.18304445e-06, 3.23419158e-07], ..., [ -1.69720124e-07, -7.90254887e-07, -3.06249355e-08, ..., 5.63925653e-07, -1.12936516e-06, -4.91445121e-07], [ -1.56005638e-06, 4.69508507e-07, -1.35530729e-06, ..., 1.76694653e-06, -1.54620722e-07, -8.53312031e-07], [ -8.67083600e-07, 4.62170050e-07, -1.95648198e-07, ..., -3.77680294e-07, -2.61958689e-07, -3.48743470e-07]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=1072, xmax=1119, ymin=901, ymax=948), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(1094.94,924.29) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=1072, xmax=1119, ymin=901, ymax=948), moments_status=0, observed_shape=galsim.Shear((0.229992433331852-0.17002166569114985j)), moments_sigma=4.999765396118164, moments_amp=0.9999370574951172, moments_centroid=galsim.PositionD(x=1094.9398594882296, y=924.2902140730297), moments_rho4=2.0000078873293408, moments_n_iter=8) got hsm parameters imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=1072, xmax=1119, ymin=901, ymax=948), array= array([[ 2.28946305e-17, 5.93698958e-17, 1.49842153e-16, ..., 2.16265672e-10, 1.65730096e-10, 1.23609067e-10], [ 1.48630883e-16, 3.79759360e-16, 9.44371856e-16, ..., 7.20889515e-10, 5.44313761e-10, 4.00004668e-10], [ 9.02223497e-16, 2.27133071e-15, 5.56521187e-15, ..., 2.24687802e-09, 1.67157888e-09, 1.21034560e-09], ..., [ 1.09138776e-09, 1.45317325e-09, 1.88317584e-09, ..., 9.68255115e-16, 3.80986288e-16, 1.45902921e-16], [ 3.44896389e-10, 4.52474141e-10, 5.77742132e-10, ..., 1.57110464e-16, 6.09103907e-17, 2.29833242e-17], [ 1.01912596e-10, 1.31734582e-10, 1.65732178e-10, ..., 2.38369133e-17, 9.10548590e-18, 3.38525485e-18]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=1072, xmax=1119, ymin=901, ymax=948), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(1094.94,924.29) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=1072, xmax=1119, ymin=901, ymax=948), moments_status=0, observed_shape=galsim.Shear((0.22998585307249675-0.17001022978037547j)), moments_sigma=4.999902725219727, moments_amp=0.9999644756317139, moments_centroid=galsim.PositionD(x=1094.9398646691948, y=924.2902129330915), moments_rho4=1.999968859635792, moments_n_iter=8) start with_hsm for imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=901, xmax=948, ymin=1510, ymax=1557), array= array([[ 1.22655206e-06, -1.42661361e-06, 8.50180811e-07, ..., 1.01519606e-06, -6.49568221e-07, -1.39539134e-06], [ -8.76452361e-07, 8.22429968e-07, -1.82626911e-06, ..., 1.45980687e-06, -1.06748121e-06, 5.19500674e-08], [ 6.82359428e-07, 6.66908988e-07, 8.72911301e-07, ..., 7.77042999e-07, 1.26813586e-06, -2.53787363e-07], ..., [ -5.39905542e-08, 7.49546530e-07, 4.23453258e-07, ..., 2.20813604e-07, -4.12538640e-07, 9.22425102e-07], [ 1.45666706e-06, 7.74921460e-09, -1.05382276e-06, ..., 8.90651997e-10, 3.04599666e-07, 9.49413618e-07], [ -1.26131816e-09, 1.45115553e-06, 1.08020208e-06, ..., -1.93169990e-06, 9.34230684e-07, -6.01052079e-07]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=901, xmax=948, ymin=1510, ymax=1557), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(924.15,1532.92) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=901, xmax=948, ymin=1510, ymax=1557), moments_status=0, observed_shape=galsim.Shear((0.22997834151042335-0.17000010236816496j)), moments_sigma=4.99981164932251, moments_amp=0.9999456405639648, moments_centroid=galsim.PositionD(x=924.1498952345294, y=1532.9199801983234), moments_rho4=1.9999540838151564, moments_n_iter=8) got hsm parameters imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=901, xmax=948, ymin=1510, ymax=1557), array= array([[ 3.76693849e-17, 9.77035630e-17, 2.46642203e-16, ..., 3.59135860e-10, 2.75271861e-10, 2.05352499e-10], [ 2.39284902e-16, 6.11510566e-16, 1.52099493e-15, ..., 1.17136234e-09, 8.84628815e-10, 6.50228593e-10], [ 1.42125416e-15, 3.57871547e-15, 8.77036718e-15, ..., 3.57234486e-09, 2.65821454e-09, 1.92513760e-09], ..., [ 6.74597544e-10, 8.98405350e-10, 1.16448806e-09, ..., 6.04048470e-16, 2.37728135e-16, 9.10593307e-17], [ 2.08595863e-10, 2.73715911e-10, 3.49566293e-10, ..., 9.59043790e-17, 3.71889551e-17, 1.40353955e-17], [ 6.03109646e-11, 7.79753623e-11, 9.81191378e-11, ..., 1.42375320e-17, 5.43972624e-18, 2.02280753e-18]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=901, xmax=948, ymin=1510, ymax=1557), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(924.15,1532.92) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=901, xmax=948, ymin=1510, ymax=1557), moments_status=0, observed_shape=galsim.Shear((0.22998607154510722-0.170010280524766j)), moments_sigma=4.999904155731201, moments_amp=0.9999641180038452, moments_centroid=galsim.PositionD(x=924.149898370813, y=1532.9199794973422), moments_rho4=1.9999714361993695, moments_n_iter=8) start with_hsm for imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=1720, xmax=1767, ymin=866, ymax=913), array= array([[ -5.37911546e-07, 1.03728951e-06, 3.91244612e-07, ..., -7.42700934e-07, -1.16977367e-06, -2.59595254e-07], [ -8.73514125e-07, 8.34733783e-07, -1.94831369e-08, ..., 1.04750166e-06, -2.09089421e-06, 9.26570635e-07], [ -1.04949578e-07, 1.87344984e-08, -5.22215998e-07, ..., -1.56790168e-06, 4.86066597e-07, -5.73833233e-07], ..., [ -1.04590720e-06, 1.03001787e-06, -1.26062423e-06, ..., -7.01691647e-07, 1.46049240e-06, 2.37872541e-07], [ -1.98352586e-06, 2.16189301e-06, 1.11704185e-06, ..., 9.22360300e-07, -1.08025311e-07, -3.53845451e-08], [ 1.40458994e-07, 4.57331964e-07, 8.61098215e-07, ..., -2.45549427e-06, 7.74880448e-07, -5.09472670e-07]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=1720, xmax=1767, ymin=866, ymax=913), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(1743.11,888.83) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=1720, xmax=1767, ymin=866, ymax=913), moments_status=0, observed_shape=galsim.Shear((0.22997952813520628-0.17001402130371981j)), moments_sigma=5.00012731552124, moments_amp=1.0000312328338623, moments_centroid=galsim.PositionD(x=1743.1098299311714, y=888.830096802052), moments_rho4=1.9999653539793665, moments_n_iter=8) got hsm parameters imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=1720, xmax=1767, ymin=866, ymax=913), array= array([[ 4.63605087e-17, 1.19955588e-16, 3.02084110e-16, ..., 3.96448124e-10, 3.03137626e-10, 2.25594460e-10], [ 2.92546757e-16, 7.45820610e-16, 1.85058297e-15, ..., 1.28451505e-09, 9.67741665e-10, 7.09601988e-10], [ 1.72612429e-15, 4.33588669e-15, 1.06003183e-14, ..., 3.89154131e-09, 2.88874147e-09, 2.08704032e-09], ..., [ 6.16054097e-10, 8.18458856e-10, 1.05830278e-09, ..., 4.94781961e-16, 1.94255360e-16, 7.42279258e-17], [ 1.89234406e-10, 2.47710769e-10, 3.15591137e-10, ..., 7.80370396e-17, 3.01874734e-17, 1.13654815e-17], [ 5.43514331e-11, 7.01007308e-11, 8.79972831e-11, ..., 1.15084654e-17, 4.38641926e-18, 1.62718917e-18]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=1720, xmax=1767, ymin=866, ymax=913), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(1743.11,888.83) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=1720, xmax=1767, ymin=866, ymax=913), moments_status=0, observed_shape=galsim.Shear((0.22998603559062844-0.17001027777899747j)), moments_sigma=4.999903678894043, moments_amp=0.9999865293502808, moments_centroid=galsim.PositionD(x=1743.1098334115925, y=888.8300960272085), moments_rho4=1.9999710392770154, moments_n_iter=8) start with_hsm for imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=865, xmax=912, ymin=1010, ymax=1057), array= array([[ 1.86102497e-07, -2.02780100e-07, -2.36136088e-06, ..., 1.07062846e-07, 8.24924939e-07, 7.30341299e-07], [ -1.83401448e-06, -4.21695347e-07, 1.00835814e-06, ..., 1.79026301e-07, 1.33499654e-07, -6.96232689e-07], [ 6.34924959e-07, 1.12090866e-07, 1.31596710e-06, ..., 4.33459036e-07, 9.78215695e-08, 1.06373204e-08], ..., [ -4.43534617e-07, 2.16792714e-07, -1.25495069e-06, ..., -1.43113311e-06, 8.94063987e-07, 1.57721004e-07], [ -2.08316820e-07, -2.97888789e-07, -9.34537582e-07, ..., 2.77731885e-07, -7.33344962e-07, 4.84322072e-07], [ -2.09632844e-06, 1.93295890e-07, 4.53361828e-08, ..., 4.46277426e-07, 8.89419255e-07, -6.07596178e-07]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=865, xmax=912, ymin=1010, ymax=1057), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(888.39,1033.19) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=865, xmax=912, ymin=1010, ymax=1057), moments_status=0, observed_shape=galsim.Shear((0.22998844396853807-0.17003052331984353j)), moments_sigma=4.999911785125732, moments_amp=0.9999915957450867, moments_centroid=galsim.PositionD(x=888.3901571587685, y=1033.1900065377993), moments_rho4=1.9999689146688417, moments_n_iter=8) got hsm parameters imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=865, xmax=912, ymin=1010, ymax=1057), array= array([[ 1.78882010e-17, 4.68869860e-17, 1.19611547e-16, ..., 2.73657375e-10, 2.11969803e-10, 1.59799701e-10], [ 1.16122367e-16, 2.99894371e-16, 7.53799607e-16, ..., 9.12141196e-10, 6.96138536e-10, 5.17087984e-10], [ 7.04846767e-16, 1.79355182e-15, 4.44189541e-15, ..., 2.84280199e-09, 2.13770135e-09, 1.56452529e-09], ..., [ 8.50416515e-10, 1.14451915e-09, 1.49916546e-09, ..., 1.22188124e-15, 4.85961172e-16, 1.88108852e-16], [ 2.68729317e-10, 3.56347063e-10, 4.59903254e-10, ..., 1.98252217e-16, 7.76886178e-17, 2.96300132e-17], [ 7.94013744e-11, 1.03741606e-10, 1.31920738e-10, ..., 3.00771540e-17, 1.16129590e-17, 4.36399398e-18]], dtype=float32), wcs=galsim.PixelScale(0.26)) Total chisq = 0.00 / 16107 dof Iteration 2: Fitting 7 stars weight = galsim.Image(bounds=galsim.BoundsI(xmin=865, xmax=912, ymin=1010, ymax=1057), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(888.39,1033.19) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=865, xmax=912, ymin=1010, ymax=1057), moments_status=0, observed_shape=galsim.Shear((0.2299862001315683-0.17001032440475575j)), moments_sigma=4.999904632568359, moments_amp=0.9999902248382568, moments_centroid=galsim.PositionD(x=888.3901581229189, y=1033.1900063268417), moments_rho4=1.9999728864836148, moments_n_iter=8) start with_hsm for imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=1386, xmax=1433, ymin=100, ymax=147), array= array([[ 6.79682955e-07, 1.15071759e-06, -4.70900972e-07, ..., 7.57913881e-07, -1.55136479e-06, -3.53039326e-07], [ 2.38386217e-08, -1.17431443e-06, 6.36037896e-07, ..., -3.79154898e-07, -2.48615635e-07, 1.47407411e-06], [ 6.88235559e-07, -4.41434395e-07, -3.82128604e-07, ..., -3.33914443e-07, 4.39471748e-08, 2.89131037e-07], ..., [ -2.67208947e-07, 1.41121552e-06, 3.05258254e-07, ..., 7.14137400e-07, -8.23862479e-07, -5.45802152e-07], [ 1.87825393e-07, -1.74170270e-06, -1.01417974e-07, ..., -7.47040900e-08, -6.94858898e-07, -1.44421804e-07], [ -1.40851830e-06, 6.48677826e-07, 1.15912587e-06, ..., 2.23111806e-06, 2.72399376e-07, 2.64195052e-07]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=1386, xmax=1433, ymin=100, ymax=147), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(1409.31,123.11) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=1386, xmax=1433, ymin=100, ymax=147), moments_status=0, observed_shape=galsim.Shear((0.22998119832422834-0.1700170269405822j)), moments_sigma=4.999868392944336, moments_amp=0.9999514818191528, moments_centroid=galsim.PositionD(x=1409.3099308144467, y=123.10999759606494), moments_rho4=1.9999471805643343, moments_n_iter=8) got hsm parameters imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=1386, xmax=1433, ymin=100, ymax=147), array= array([[ 2.25210228e-17, 5.88322152e-17, 1.49581347e-16, ..., 2.96205532e-10, 2.28665906e-10, 1.71808581e-10], [ 1.45240356e-16, 3.73836059e-16, 9.36505249e-16, ..., 9.80839743e-10, 7.46058937e-10, 5.52310531e-10], [ 8.75822644e-16, 2.22114464e-15, 5.48242456e-15, ..., 3.03691472e-09, 2.27601160e-09, 1.66016578e-09], ..., [ 7.96898658e-10, 1.06889708e-09, 1.39541634e-09, ..., 9.84384845e-16, 3.90192503e-16, 1.50531718e-16], [ 2.50170690e-10, 3.30625249e-10, 4.25275870e-10, ..., 1.58673345e-16, 6.19704657e-17, 2.35559515e-17], [ 7.34343766e-11, 9.56237520e-11, 1.21190141e-10, ..., 2.39151266e-17, 9.20280535e-18, 3.44669576e-18]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=1386, xmax=1433, ymin=100, ymax=147), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(1409.31,123.11) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=1386, xmax=1433, ymin=100, ymax=147), moments_status=0, observed_shape=galsim.Shear((0.22998616417708262-0.17001032165898455j)), moments_sigma=4.999904632568359, moments_amp=0.9999587535858154, moments_centroid=galsim.PositionD(x=1409.3099324909167, y=123.10999723657686), moments_rho4=1.999972576339034, moments_n_iter=8) start with_hsm for imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=100, xmax=147, ymin=322, ymax=369), array= array([[ 1.57293516e-06, 4.83428551e-07, 1.01336786e-06, ..., -1.52099460e-06, 1.62390904e-06, -3.06897618e-07], [ -4.01491747e-07, -6.06370918e-07, -1.84310323e-07, ..., -3.39801630e-07, 4.16551785e-07, 1.06846937e-06], [ -4.68112916e-07, 2.85567381e-07, -4.94964070e-07, ..., -1.98679504e-06, -9.59963018e-07, 1.69272738e-07], ..., [ 4.87750185e-07, 5.51780090e-07, 4.79130222e-07, ..., -2.40008035e-06, 1.37110010e-06, 1.22177823e-06], [ -2.55823608e-07, -8.19296361e-07, 2.35463608e-06, ..., -1.41771136e-06, -1.60905051e-06, 4.51424427e-07], [ -2.71851150e-07, 1.25770748e-06, 1.75458837e-07, ..., -1.10357121e-06, 1.51166330e-06, -2.85283708e-09]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=100, xmax=147, ymin=322, ymax=369), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(123.12,345.43) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=100, xmax=147, ymin=322, ymax=369), moments_status=0, observed_shape=galsim.Shear((0.23000160481118787-0.16999831959428038j)), moments_sigma=4.9998674392700195, moments_amp=0.9999805688858032, moments_centroid=galsim.PositionD(x=123.11974024161614, y=345.42992834522084), moments_rho4=1.9999460194128353, moments_n_iter=8) got hsm parameters imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=100, xmax=147, ymin=322, ymax=369), array= array([[ 1.47255321e-17, 3.84519772e-17, 9.77240506e-17, ..., 1.90108693e-10, 1.46700249e-10, 1.10177985e-10], [ 9.67562758e-17, 2.48939383e-16, 6.23366115e-16, ..., 6.41380671e-10, 4.87653640e-10, 3.60862923e-10], [ 5.94452222e-16, 1.50694774e-15, 3.71804477e-15, ..., 2.02329575e-09, 1.51572987e-09, 1.10514531e-09], ..., [ 1.20722710e-09, 1.61861047e-09, 2.11217954e-09, ..., 1.46378389e-15, 5.79977964e-16, 2.23656350e-16], [ 3.86127963e-10, 5.10095410e-10, 6.55853316e-10, ..., 2.40394727e-16, 9.38482793e-17, 3.56584790e-17], [ 1.15479050e-10, 1.50310792e-10, 1.90419874e-10, ..., 3.69149796e-17, 1.41994254e-17, 5.31586628e-18]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=100, xmax=147, ymin=322, ymax=369), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(123.12,345.43) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=100, xmax=147, ymin=322, ymax=369), moments_status=0, observed_shape=galsim.Shear((0.2299860535678676-0.17001027915188166j)), moments_sigma=4.999904155731201, moments_amp=0.9999879002571106, moments_centroid=galsim.PositionD(x=123.11974023234252, y=345.42992835577724), moments_rho4=1.999971140946323, moments_n_iter=8) start with_hsm for imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=323, xmax=370, ymin=544, ymax=591), array= array([[ -4.64985504e-07, -5.38259087e-07, 1.27860790e-06, ..., 9.59986892e-07, 5.81386814e-07, -1.69146315e-08], [ 5.36656898e-07, 6.98122562e-07, -2.27826763e-06, ..., -4.54603338e-07, -4.67363577e-07, 7.88441866e-07], [ -4.12713433e-07, 8.72180976e-07, -1.06762320e-06, ..., 9.32305170e-07, 9.44102894e-08, -2.08292727e-06], ..., [ 5.33941829e-07, 1.06433686e-06, 1.69580323e-06, ..., 8.69430949e-07, -1.51942322e-06, -1.98380121e-07], [ 8.97153313e-07, -1.10840438e-06, 6.98605405e-08, ..., 1.42423278e-06, 3.54525412e-07, -1.99148940e-06], [ -6.01974818e-07, 1.06053415e-06, 1.12591874e-06, ..., -1.30754131e-07, 3.08745797e-07, -9.19329295e-07]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=323, xmax=370, ymin=544, ymax=591), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(345.98,567.45) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=323, xmax=370, ymin=544, ymax=591), moments_status=0, observed_shape=galsim.Shear((0.22998069831116813-0.16999026883163862j)), moments_sigma=4.999975204467773, moments_amp=0.999998152256012, moments_centroid=galsim.PositionD(x=345.980138171224, y=567.4498072462465), moments_rho4=1.999932797336089, moments_n_iter=8) got hsm parameters imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=323, xmax=370, ymin=544, ymax=591), array= array([[ 1.62344651e-17, 4.22445936e-17, 1.06989069e-16, ..., 1.79149334e-10, 1.37762024e-10, 1.03104823e-10], [ 1.06592810e-16, 2.73292586e-16, 6.81966290e-16, ..., 6.03963712e-10, 4.57606258e-10, 3.37449041e-10], [ 6.54406326e-16, 1.65315749e-15, 4.06458374e-15, ..., 1.90386484e-09, 1.42129442e-09, 1.03268316e-09], ..., [ 1.28776390e-09, 1.72058090e-09, 2.23742824e-09, ..., 1.33465944e-15, 5.26975619e-16, 2.02509630e-16], [ 4.11585765e-10, 5.41833689e-10, 6.94235391e-10, ..., 2.19028321e-16, 8.52093175e-17, 3.22633195e-17], [ 1.23002539e-10, 1.59546210e-10, 2.01416051e-10, ..., 3.36093177e-17, 1.28828905e-17, 4.80620267e-18]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=323, xmax=370, ymin=544, ymax=591), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(345.98,567.45) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=323, xmax=370, ymin=544, ymax=591), moments_status=0, observed_shape=galsim.Shear((0.22998588902696587-0.1700102325261405j)), moments_sigma=4.999903202056885, moments_amp=0.9999837875366211, moments_centroid=galsim.PositionD(x=345.98013818032746, y=567.4498072281253), moments_rho4=1.9999694506086794, moments_n_iter=8) start with_hsm for imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=1072, xmax=1119, ymin=901, ymax=948), array= array([[ -1.50366873e-06, -1.14600471e-06, 1.25665451e-06, ..., -1.32513458e-06, -7.00539260e-07, 8.74302089e-07], [ -8.53509903e-07, 7.94479604e-07, -1.04325943e-06, ..., -1.24426458e-06, 8.49906655e-07, 4.03015946e-07], [ -1.13993985e-06, 1.78198775e-06, -3.41910862e-07, ..., 1.77908319e-07, -1.18304445e-06, 3.23419158e-07], ..., [ -1.69720124e-07, -7.90254887e-07, -3.06249355e-08, ..., 5.63925653e-07, -1.12936516e-06, -4.91445121e-07], [ -1.56005638e-06, 4.69508507e-07, -1.35530729e-06, ..., 1.76694653e-06, -1.54620722e-07, -8.53312031e-07], [ -8.67083600e-07, 4.62170050e-07, -1.95648198e-07, ..., -3.77680294e-07, -2.61958689e-07, -3.48743470e-07]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=1072, xmax=1119, ymin=901, ymax=948), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(1094.94,924.29) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=1072, xmax=1119, ymin=901, ymax=948), moments_status=0, observed_shape=galsim.Shear((0.229992433331852-0.17002166569114985j)), moments_sigma=4.999765396118164, moments_amp=0.9999370574951172, moments_centroid=galsim.PositionD(x=1094.9398594882296, y=924.2902140730297), moments_rho4=2.0000078873293408, moments_n_iter=8) got hsm parameters imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=1072, xmax=1119, ymin=901, ymax=948), array= array([[ 2.28946984e-17, 5.93700613e-17, 1.49842564e-16, ..., 2.16265061e-10, 1.65729611e-10, 1.23608693e-10], [ 1.48631333e-16, 3.79760419e-16, 9.44374397e-16, ..., 7.20887516e-10, 5.44312206e-10, 4.00003447e-10], [ 9.02226143e-16, 2.27133706e-15, 5.56522669e-15, ..., 2.24687180e-09, 1.67157410e-09, 1.21034194e-09], ..., [ 1.09139076e-09, 1.45317702e-09, 1.88318050e-09, ..., 9.68252679e-16, 3.80985255e-16, 1.45902510e-16], [ 3.44897333e-10, 4.52475335e-10, 5.77743575e-10, ..., 1.57110067e-16, 6.09102319e-17, 2.29832613e-17], [ 1.01912874e-10, 1.31734915e-10, 1.65732594e-10, ..., 2.38368538e-17, 9.10546191e-18, 3.38524534e-18]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=1072, xmax=1119, ymin=901, ymax=948), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(1094.94,924.29) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=1072, xmax=1119, ymin=901, ymax=948), moments_status=0, observed_shape=galsim.Shear((0.22998585169961452-0.17001021264548202j)), moments_sigma=4.999902725219727, moments_amp=0.9999644756317139, moments_centroid=galsim.PositionD(x=1094.939859555632, y=924.2902140497815), moments_rho4=1.9999688613199207, moments_n_iter=8) start with_hsm for imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=901, xmax=948, ymin=1510, ymax=1557), array= array([[ 1.22655206e-06, -1.42661361e-06, 8.50180811e-07, ..., 1.01519606e-06, -6.49568221e-07, -1.39539134e-06], [ -8.76452361e-07, 8.22429968e-07, -1.82626911e-06, ..., 1.45980687e-06, -1.06748121e-06, 5.19500674e-08], [ 6.82359428e-07, 6.66908988e-07, 8.72911301e-07, ..., 7.77042999e-07, 1.26813586e-06, -2.53787363e-07], ..., [ -5.39905542e-08, 7.49546530e-07, 4.23453258e-07, ..., 2.20813604e-07, -4.12538640e-07, 9.22425102e-07], [ 1.45666706e-06, 7.74921460e-09, -1.05382276e-06, ..., 8.90651997e-10, 3.04599666e-07, 9.49413618e-07], [ -1.26131816e-09, 1.45115553e-06, 1.08020208e-06, ..., -1.93169990e-06, 9.34230684e-07, -6.01052079e-07]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=901, xmax=948, ymin=1510, ymax=1557), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(924.15,1532.92) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=901, xmax=948, ymin=1510, ymax=1557), moments_status=0, observed_shape=galsim.Shear((0.22997834151042335-0.17000010236816496j)), moments_sigma=4.99981164932251, moments_amp=0.9999456405639648, moments_centroid=galsim.PositionD(x=924.1498952345294, y=1532.9199801983234), moments_rho4=1.9999540838151564, moments_n_iter=8) got hsm parameters imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=901, xmax=948, ymin=1510, ymax=1557), array= array([[ 3.76694544e-17, 9.77037351e-17, 2.46642627e-16, ..., 3.59135222e-10, 2.75271361e-10, 2.05352124e-10], [ 2.39285352e-16, 6.11511678e-16, 1.52099758e-15, ..., 1.17136034e-09, 8.84627260e-10, 6.50227372e-10], [ 1.42125670e-15, 3.57872183e-15, 8.77038242e-15, ..., 3.57233865e-09, 2.65820965e-09, 1.92513405e-09], ..., [ 6.74598655e-10, 8.98406793e-10, 1.16448984e-09, ..., 6.04047570e-16, 2.37727764e-16, 9.10591851e-17], [ 2.08596210e-10, 2.73716355e-10, 3.49566820e-10, ..., 9.59042334e-17, 3.71888988e-17, 1.40353740e-17], [ 6.03110617e-11, 7.79754872e-11, 9.81192835e-11, ..., 1.42375105e-17, 5.43971797e-18, 2.02280422e-18]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=901, xmax=948, ymin=1510, ymax=1557), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(924.15,1532.92) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=901, xmax=948, ymin=1510, ymax=1557), moments_status=0, observed_shape=galsim.Shear((0.22998607154510722-0.170010280524766j)), moments_sigma=4.999904155731201, moments_amp=0.9999641180038452, moments_centroid=galsim.PositionD(x=924.1498952189185, y=1532.9199802228143), moments_rho4=1.9999714362554022, moments_n_iter=8) start with_hsm for imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=1720, xmax=1767, ymin=866, ymax=913), array= array([[ -5.37911546e-07, 1.03728951e-06, 3.91244612e-07, ..., -7.42700934e-07, -1.16977367e-06, -2.59595254e-07], [ -8.73514125e-07, 8.34733783e-07, -1.94831369e-08, ..., 1.04750166e-06, -2.09089421e-06, 9.26570635e-07], [ -1.04949578e-07, 1.87344984e-08, -5.22215998e-07, ..., -1.56790168e-06, 4.86066597e-07, -5.73833233e-07], ..., [ -1.04590720e-06, 1.03001787e-06, -1.26062423e-06, ..., -7.01691647e-07, 1.46049240e-06, 2.37872541e-07], [ -1.98352586e-06, 2.16189301e-06, 1.11704185e-06, ..., 9.22360300e-07, -1.08025311e-07, -3.53845451e-08], [ 1.40458994e-07, 4.57331964e-07, 8.61098215e-07, ..., -2.45549427e-06, 7.74880448e-07, -5.09472670e-07]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=1720, xmax=1767, ymin=866, ymax=913), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(1743.11,888.83) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=1720, xmax=1767, ymin=866, ymax=913), moments_status=0, observed_shape=galsim.Shear((0.22997952813520628-0.17001402130371981j)), moments_sigma=5.00012731552124, moments_amp=1.0000312328338623, moments_centroid=galsim.PositionD(x=1743.1098299311714, y=888.830096802052), moments_rho4=1.9999653539793665, moments_n_iter=8) got hsm parameters imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=1720, xmax=1767, ymin=866, ymax=913), array= array([[ 4.63606047e-17, 1.19955827e-16, 3.02084692e-16, ..., 3.96447375e-10, 3.03137015e-10, 2.25593988e-10], [ 2.92547366e-16, 7.45822039e-16, 1.85058657e-15, ..., 1.28451261e-09, 9.67739777e-10, 7.09600489e-10], [ 1.72612779e-15, 4.33589516e-15, 1.06003378e-14, ..., 3.89153376e-09, 2.88873570e-09, 2.08703610e-09], ..., [ 6.16055262e-10, 8.18460300e-10, 1.05830467e-09, ..., 4.94781114e-16, 1.94255029e-16, 7.42277868e-17], [ 1.89234767e-10, 2.47711213e-10, 3.15591664e-10, ..., 7.80369138e-17, 3.01874204e-17, 1.13654616e-17], [ 5.43515337e-11, 7.01008557e-11, 8.79974357e-11, ..., 1.15084464e-17, 4.38641182e-18, 1.62718617e-18]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=1720, xmax=1767, ymin=866, ymax=913), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(1743.11,888.83) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=1720, xmax=1767, ymin=866, ymax=913), moments_status=0, observed_shape=galsim.Shear((0.22998603559062844-0.17001027777899747j)), moments_sigma=4.999903678894043, moments_amp=0.9999865293502808, moments_centroid=galsim.PositionD(x=1743.109829944854, y=888.8300967922208), moments_rho4=1.999971036622578, moments_n_iter=8) start with_hsm for imported galsim: 1.5 /home/esheldon/miniconda3/lib/python3.5/site-packages/galsim/__init__.py image = galsim.Image(bounds=galsim.BoundsI(xmin=865, xmax=912, ymin=1010, ymax=1057), array= array([[ 1.86102497e-07, -2.02780100e-07, -2.36136088e-06, ..., 1.07062846e-07, 8.24924939e-07, 7.30341299e-07], [ -1.83401448e-06, -4.21695347e-07, 1.00835814e-06, ..., 1.79026301e-07, 1.33499654e-07, -6.96232689e-07], [ 6.34924959e-07, 1.12090866e-07, 1.31596710e-06, ..., 4.33459036e-07, 9.78215695e-08, 1.06373204e-08], ..., [ -4.43534617e-07, 2.16792714e-07, -1.25495069e-06, ..., -1.43113311e-06, 8.94063987e-07, 1.57721004e-07], [ -2.08316820e-07, -2.97888789e-07, -9.34537582e-07, ..., 2.77731885e-07, -7.33344962e-07, 4.84322072e-07], [ -2.09632844e-06, 1.93295890e-07, 4.53361828e-08, ..., 4.46277426e-07, 8.89419255e-07, -6.07596178e-07]], dtype=float32), wcs=galsim.PixelScale(0.26)) weight = galsim.Image(bounds=galsim.BoundsI(xmin=865, xmax=912, ymin=1010, ymax=1057), array= array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.]]), wcs=None) image_pos = galsim.PositionD(888.39,1033.19) mom = galsim.hsm.ShapeData(image_bounds=galsim.BoundsI(xmin=865, xmax=912, ymin=1010, ymax=1057), moments_status=0, observed_shape=galsim.Shear((0.22998844396853807-0.17003052331984353j)), moments_sigma=4.999911785125732, moments_amp=0.9999915957450867, moments_centroid=galsim.PositionD(x=888.3901571587685, y=1033.1900065377993), moments_rho4=1.9999689146688417, moments_n_iter=8) Total chisq = 0.00 / 16107 dof Writing PSF to file output/simple_psf.fits Wrote the basic PSF information to extname psf Wrote the PSF stars to extname psf_stars Wrote the PSF WCS to extname psf_wcs Reading PSF from file output/simple_psf.fits Reading image file data/simple_image.fits Reading star catalog data/simple_cat.fits. Removing objects with flag (col flag) != 0 Removing objects with use (col use) == 0 Processing catalog data/simple_cat.fits with 7 stars Read a total of 7 stars from 1 image Building SimplePSF