View on GitHub

IRAF Community Distribution

IRAF maintained by the community

Home | Installation | Packages | X11IRAF | PyRAF | Forum

segmentation violation in fitparams

Philip Massey wrote on Aug 22, 2016


I'm using fitparams,using 2.16.1 EXPORT on a Mac and I'm getting a segmentation violation. Here are my parameters:

observat= standobs List of observations files

catalogs= nlandolt List of standard catalog files

config = lmi.cfg Configuration file

paramete= lmi.ans Output parameters file

(weighti= uniform) Weighting type (uniform,photometric,equations)

(addscat= yes) Add a scatter term to the weights ?

(toleran= 3.0000000000000E-5) Fit convergence tolerance

(maxiter= 15) Maximum number of fit iterations

(nreject= 0) Number of rejection iterations

(low_rej= 3.) Low sigma rejection factor

(high_re= 3.) High sigma rejection factor

(grow = 0.) Rejection growing radius

(interac= yes) Solve fit interactively ?

(logfile= STDOUT) Output log file

(log_unm= yes) Log any unmatched stars ?

(log_fit= no) Log the fit parameters and statistics ?

(log_res= no) Log the results ?

My "standobs" file is pretty short: here it is:

# FIELD FILTER OTIME AIRMASS XCENTER YCENTER MAG MERR

110_441 U INDEF 1.210 2080.549 1296.255 13.757 0.003

* B INDEF 1.210 2079.600 1294.721 11.168 0.002

* V INDEF 1.210 2080.066 1294.516 10.581 0.001

* R INDEF 1.210 2079.373 1293.331 10.222 0.001

* I INDEF 1.210 2078.205 1291.540 10.228 0.001

110_450 U INDEF 1.210 961.499 2119.321 15.005 0.005

* B INDEF 1.210 960.771 2118.641 12.018 0.002

* V INDEF 1.210 961.161 2117.582 11.105 0.001

* R INDEF 1.210 960.355 2117.109 10.459 0.001

* I INDEF 1.210 958.589 2115.782 10.166 0.001

MARK_A1 U INDEF 1.740 1689.317 1801.150 18.865 0.009

* B INDEF 1.740 1687.136 1799.455 16.288 0.003

* V INDEF 1.730 1685.465 1797.432 15.542 0.002

* R INDEF 1.730 1684.113 1796.474 15.125 0.002

* I INDEF 1.730 1682.638 1796.413 15.058 0.002

MARK_A2 U INDEF 1.740 1900.295 2223.136 17.687 0.005

* B INDEF 1.740 1897.703 2221.268 14.959 0.001

* V INDEF 1.730 1896.460 2219.470 14.189 0.001

* R INDEF 1.730 1894.733 2218.379 13.745 0.001

* I INDEF 1.730 1893.508 2217.972 13.686 0.001

MARK_A3 U INDEF 1.740 1357.699 2192.872 18.781 0.009

* B INDEF 1.740 1355.324 2191.042 15.494 0.002

* V INDEF 1.730 1353.780 2189.044 14.468 0.001

* R INDEF 1.730 1352.493 2188.090 13.823 0.001

* I INDEF 1.730 1351.166 2187.553 13.612 0.001

113_259 U INDEF 1.730 1877.799 972.449 16.455 0.009

* B INDEF 1.730 1875.438 971.193 12.678 0.003

* V INDEF 1.720 1873.533 969.753 11.452 0.002

* R INDEF 1.720 1872.157 967.935 10.765 0.001

* I INDEF 1.720 1870.708 968.842 10.458 0.001

113_260 U INDEF 1.730 1673.974 2527.391 15.422 0.006

* B INDEF 1.730 1671.716 2525.741 12.651 0.003

* V INDEF 1.720 1669.579 2524.443 12.061 0.002

* R INDEF 1.720 1669.083 2523.505 11.734 0.002

* I INDEF 1.720 1666.820 2522.668 11.692 0.002

The fstandobs.dat file is unsurprising too:

# Declare the observations file variables

observations

TU 3 # time of observation in filter U

XU 4 # airmass in filter U

xU 5 # x coordinate in filter U

yU 6 # y coordinate in filter U

mU 7 # instrumental magnitude in filter U

error(mU) 8 # magnitude error in filter U

TB 10 # time of observation in filter B

XB 11 # airmass in filter B

xB 12 # x coordinate in filter B

yB 13 # y coordinate in filter B

mB 14 # instrumental magnitude in filter B

error(m 15 # magnitude error in filter B

TV 17 # time of observation in filter V

XV 18 # airmass in filter V

xV 19 # x coordinate in filter V

yV 20 # y coordinate in filter V

mV 21 # instrumental magnitude in filter V

error(mV) 22 # magnitude error in filter V

TR 24 # time of observation in filter R

XR 25 # airmass in filter R

xR 26 # x coordinate in filter R

yR 27 # y coordinate in filter R

mR 28 # instrumental magnitude in filter R

error(mR) 29 # magnitude error in filter R

TI 31 # time of observation in filter I

XI 32 # airmass in filter I

xI 33 # x coordinate in filter I

yI 34 # y coordinate in filter I

mI 35 # instrumental magnitude in filter I

error(mI) 36 # magnitude error in filter I

Finally, my lmi.cfg file is:

# Declare the new Landolt UBVRI standards catalog variables

catalog

V 4 # the V magnitude

BV 5 # the (B-V) color

UB 6 # the (U- color

VR 7 # the (V-R) color

RI 8 # the (R-I) color

VI 9 # the (V-I) color

error(V) 12 # the V magnitude error

error(BV) 13 # the (B-V) color error

error(U 14 # the (U- color error

error(VR) 15 # the (V-R) color error

error(RI) 16 # the (R-I) color error

error(VI) 17 # the (V-I) color error

# Declare the observations file variables

observations

TU 3 # time of observation in filter U

XU 4 # airmass in filter U

xU 5 # x coordinate in filter U

yU 6 # y coordinate in filter U

mU 7 # instrumental magnitude in filter U

error(mU) 8 # magnitude error in filter U

TB 10 # time of observation in filter B

XB 11 # airmass in filter B

xB 12 # x coordinate in filter B

yB 13 # y coordinate in filter B

mB 14 # instrumental magnitude in filter B

error(m 15 # magnitude error in filter B

TV 17 # time of observation in filter V

XV 18 # airmass in filter V

xV 19 # x coordinate in filter V

yV 20 # y coordinate in filter V

mV 21 # instrumental magnitude in filter V

error(mV) 22 # magnitude error in filter V

TR 24 # time of observation in filter R

XR 25 # airmass in filter R

xR 26 # x coordinate in filter R

yR 27 # y coordinate in filter R

mR 28 # instrumental magnitude in filter R

error(mR) 29 # magnitude error in filter R

TI 31 # time of observation in filter I

XI 32 # airmass in filter I

xI 33 # x coordinate in filter I

yI 34 # y coordinate in filter I

mI 35 # instrumental magnitude in filter I

error(mI) 36 # magnitude error in filter I

# Sample transformation section for the new Landolt UBVRI system

transformation

fit u1=0.0, u2=0.65, u3=0.0

const u4=0.0

UFIT : mU = (UB + BV + V) + u1 + u2 * XU + u3 * UB + u4 * UB * XU

fit b1=0.0, b2=0.35, b3=0.0

const b4=0.0

BFIT : mB = (BV + V) + b1 + b2 * XB + b3 * BV + b4 * BV * XB

fit v1=0.0, v2=0.17, v3=0.0

const v4=0.0

VFIT : mV = V + v1 + v2 * XV + v3 * BV + v4 * BV * XV

fit r1=0.0, r2=0.08, r3=0.0

const r4=0.0

RFIT : mR = (V - VR) + r1 + r2 * XR + r3 * VR + r4 * VR * XR

fit i1=0.0, i2=0.03, i3=0.0

const i4=0.0

IFIT : mI = (V - VI) + i1 + i2 * XI + i3 * VI + i4 * VI * XI

So, I've been staring at this all morning and I can't see that there's any problems. Thoughts?

Thanks!

---phil

(If it's easlier, I've also made copies of these files available by download:

http://www2.lowell.edu/users/massey/lmi.cfg

http://www2.lowell.edu/users/massey/standobs

http://www2.lowell.edu/users/massey/fstandobs

Mike Fitzpatrick wrote on Aug 22, 2016


Hi Phil,

There is a know, but undiagnosed, 64-bit issue with FITPARAMS that we should rule out first. To do this, try installing the 32-bit binary from

ftp://iraf.noao.edu/iraf/v216/support/macosx/x_photcal.e

Download in binary mode and replace (after backing up) your existing $iraf/noao/bin.macintel/x_photcal.e binary, be sure the execute mode is also set.

If that works you can continue to use the 32-bit binary and I'll grab your test files to see if I can reproduce it here (not always true in the past reports). Hope this helps.

Philip Massey wrote on Aug 22, 2016


Mike---that solved it! MANY THANKS!!!

---phil

Last post on Aug 22, 2016