[daip] IMEAN histogram fit fails on my machine, but not on alcor!
Michael Rupen
mrupen at milkyway.gsfc.nasa.goV
Fri Oct 15 10:15:59 EDT 2004
Dear Daip,
in trying out AIPS on my new machine, all seems fine -- except that I can't
get IMEAN to do a histogram fit for the mean/noise estimates. This fails
on several different images. I ran the midnight job but that made no
difference; then I copied the offending file over to alcor, and IMEAN worked
just fine. Apart from the machine architecture the only oddity I know is that
AIPS was compiled here using gnu compiler version 3.4.2 . This specific
problem is basically an annoyance but I worry this may be the sign of more
subtle errors to come...I did run the Y2K (large) test, with acceptable
results in terms of bit errors.
The sample offending file is in /home/alcor/mrupen/v4641sgr.icl ,
as well as on the alcor AIPS disk 1 under user 213.
Thank you much --
Michael
==============================================================================
IMEAN on xeros (at Goddard):
>inp
AIPS 1: IMEAN: Task to print the mean, rms and extrema in an image
AIPS 1: Adverbs Values Comments
AIPS 1: ----------------------------------------------------------------
AIPS 1: DOHIST -1 True (1.0) do histogram plot.
AIPS 1: = 2 => flux on x axis
AIPS 1: USERID 0 User ID. 0=>current user
AIPS 1: 32000=>all users
AIPS 1: INNAME 'V4641SGR 289' Image name (name)
AIPS 1: INCLASS 'ICL001' Image name (class)
AIPS 1: INSEQ 1 Image name (seq. #)
AIPS 1: INDISK 5 Disk drive #
AIPS 1: BLC 35 35 Bottom left corner of image
AIPS 1: *rest 0 0=>entire image
AIPS 1: TRC 480 230 Top right corner of image
AIPS 1: *rest 0 0=>entire image
AIPS 1: NBOXES 0 No. of ranges for histogram.
AIPS 1: PIXRANGE 0 0 Min and max range for hist.
AIPS 1: FUNCTYPE ' ' 'LG' => do log10 plot of #
AIPS 1: samples, else linear
AIPS 1: PIXAVG 0 Estimate of mean noise value
AIPS 1: PIXSTD 0 Estimate of true noise rms
AIPS 1: < 0 => don't do one
AIPS 1: = 0 => 2-passes to get
AIPS 1: DOCAT 1 Put true RMS in header
AIPS 1: LTYPE 3 Type of labeling: 1 border,
AIPS 1: 2 no ticks, 3 - 6 standard,
AIPS 1: 7 - 10 only tick labels
AIPS 1: <0 -> no date/time
AIPS 1: OUTFILE *all ' ' Name of output log file,
AIPS 1: No output to file if blank
AIPS 1: DOTV -1 > 0 Do plot on the TV, else
AIPS 1: make a plot file
AIPS 1: GRCHAN 0 Graphics channel 0 => 1.
>go imean
AIPS 1: Waiting for returned adverbs
IMEAN1: Task IMEAN (release of 31DEC04) begins
IMEAN1: Histogram fit for RMS fails: check PIXAVG, PIXSTD
IMEAN1: Image= V4641SGR 289.ICL001. 1 5 xywind= 35 35 480 230
IMEAN1: Mean and rms found by including all data:
IMEAN1: Mean= 7.5626E-07 Rms= 8.7908E-05 JY/BEAM over 87416 pixels
IMEAN1: Minimum=-3.7414E-04 at 254 76 1 1
IMEAN1: Skypos: RA 18 19 21.65176 DEC -25 24 43.7000
IMEAN1: Maximum= 3.6190E-04 at 124 230 1 1
IMEAN1: Skypos: RA 18 19 22.61123 DEC -25 24 28.2998
IMEAN1: Skypos: IPOL 4860.100 MHZ
IMEAN1: returns adverbs to AIPS
IMEAN1: Appears to have ended successfully
IMEAN1: xeros 31DEC04 TST: Cpu= 0.0 Real= 0
==============================================================================
IMEAN on alcor (after FITTP'ing this data set; using 31DEC04):
>inp
AIPS 2: IMEAN: Task to print the mean, rms and extrema in an image
AIPS 2: Adverbs Values Comments
AIPS 2: ----------------------------------------------------------------
AIPS 2: DOHIST -1 True (1.0) do histogram plot.
AIPS 2: = 2 => flux on x axis
AIPS 2: USERID 0 User ID. 0=>current user
AIPS 2: 32000=>all users
AIPS 2: INNAME 'V4641SGR 289' Image name (name)
AIPS 2: INCLASS 'ICL001' Image name (class)
AIPS 2: INSEQ 1 Image name (seq. #)
AIPS 2: INDISK 1 Disk drive #
AIPS 2: BLC 35 35 Bottom left corner of image
AIPS 2: *rest 0 0=>entire image
AIPS 2: TRC 480 230 Top right corner of image
AIPS 2: *rest 0 0=>entire image
AIPS 2: NBOXES 0 No. of ranges for histogram.
AIPS 2: PIXRANGE 0 0 Min and max range for hist.
AIPS 2: FUNCTYPE ' ' 'LG' => do log10 plot of #
AIPS 2: samples, else linear
AIPS 2: PIXAVG 0 Estimate of mean noise value
AIPS 2: PIXSTD 0 Estimate of true noise rms
AIPS 2: < 0 => don't do one
AIPS 2: = 0 => 2-passes to get
AIPS 2: DOCAT 1 Put true RMS in header
AIPS 2: LTYPE 3 Type of labeling: 1 border,
AIPS 2: 2 no ticks, 3 - 6 standard,
AIPS 2: 7 - 10 only tick labels
AIPS 2: <0 -> no date/time
AIPS 2: OUTFILE *all ' ' Name of output log file,
AIPS 2: No output to file if blank
AIPS 2: DOTV -1 > 0 Do plot on the TV, else
AIPS 2: make a plot file
AIPS 2: GRCHAN 0 Graphics channel 0 => 1.
>go imean
AIPS 2: Waiting for returned adverbs
IMEAN2: Task IMEAN (release of 31DEC04) begins
IMEAN2: Image= V4641SGR 289.ICL001. 1 1 xywind= 35 35 480 230
IMEAN2: Mean and rms found by fitting peak in histogram:
IMEAN2: Mean= 2.6486E-06 Rms= 9.1436E-05 **** from histogram
IMEAN2: Mean and rms found by including all data:
IMEAN2: Mean= 7.5626E-07 Rms= 8.7908E-05 JY/BEAM over 87416 pixels
IMEAN2: Minimum=-3.7414E-04 at 254 76 1 1
IMEAN2: Skypos: RA 18 19 21.65176 DEC -25 24 43.7000
IMEAN2: Maximum= 3.6190E-04 at 124 230 1 1
IMEAN2: Skypos: RA 18 19 22.61123 DEC -25 24 28.2998
IMEAN2: Skypos: IPOL 4860.100 MHZ
IMEAN2: returns adverbs to AIPS
IMEAN2: Appears to have ended successfully
IMEAN2: alcor 31DEC04 TST: Cpu= 0.0 Real= 0
==============================================================================
More information about the Daip
mailing list