Welcome to iraf.net Saturday, January 20 2018 @ 10:56 PM GMT


 Forum Index > Help Desk > Applications New Topic Post Reply
 segmentation violation in fitparams
   
massey
 08/22/2016 07:14PM (Read 772 times)  
+++++
Active Member
sad

Status: offline


Registered: 02/10/2006
Posts: 158
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(mCool 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-Cool 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(UCool 14 # the (U-Cool 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(mCool 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

 
Profile Email Website
 Quote
fitz
 08/22/2016 07:21PM  
AAAAA
Admin

Status: offline


Registered: 09/30/2005
Posts: 3991
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.

 
Profile Email
 Quote
massey
 08/22/2016 08:18PM  
+++++
Active Member
happy

Status: offline


Registered: 02/10/2006
Posts: 158
Mike---that solved it! MANY THANKS!!!

---phil

 
Profile Email Website
 Quote
   
Content generated in: 0.09 seconds
New Topic Post Reply

Normal Topic Normal Topic
Sticky Topic Sticky Topic
Locked Topic Locked Topic
New Post New Post
Sticky Topic W/ New Post Sticky Topic W/ New Post
Locked Topic W/ New Post Locked Topic W/ New Post
View Anonymous Posts 
Anonymous users can post 
Filtered HTML Allowed 
Censored Content 
dog allergies remedies cialis 20 mg chilblain remedies


Privacy Policy
Terms of Use

User Functions

Login