This loads a font easier to read for people with dyslexia.
This renders the document in high contrast mode.
This renders the document as white on black
This can help those with trouble processing rapid screen movements.

[atnf-data-reduction] diverging CLEAN solutions

From: <Ilana.Feain_at_email.protected>
Date: Wed, 05 Dec 2007 16:15:52 +1100

Hi

I wonder whether anyone has noticed the following diverging behaviour of
CLEAN or can offer an educated guess on what's going on:

Using a single dirty map and single dirty beam (the same one every time)
I run a straightforward CLEAN as follows:

   Task: clean
   map = p_148.ch0.imap
   beam = p_148.ch0.beam/
   model =
   out = p_148.ch0.cmp
   gain =
   options = positive
   cutoff = 2e-3
   niters = 1000
   region =
   phat =
   minpatch =
   speed =
   mode =
   clip =

What comes out usually makes sense and is identically reproducible in
most cases. But once every 8 times (or so) I run the task above, it
diverges and fails producing nonsense results. An example of this is
attached to this email where the first 3 runs of clean produce the
identical result but the fourth time produces the divergence.

This problem is non reproducible in that clean tends to fail every so
often, but not in a periodic way.

The problem occurs when I run clean over and over both with the same map
and beam, and also when I run it over and over in a loop to clean
different pointings.

thanks,
Ilana
clean% inp
  Task: clean
  map = p_148.ch0.imap
  beam = p_148.ch0.beam/
  model =
  out = p_148.ch0.cmp
  gain =
  options = positive
  cutoff = 2e-3
  niters = 1000
  region =
  phat =
  minpatch =
  speed =
  mode =
  clip =

clean% go
Clean: version 1.0 07-Apr-04
Begin iterating
 Clark Iterations: 18
 Residual min,max,rms: -2.623E-02 3.476E-02 4.954E-03
 Total CLEANed flux: 9.487E-02
 Clark Iterations: 183
 Residual min,max,rms: -2.068E-02 1.592E-02 4.352E-03
 Total CLEANed flux: 4.040E-01
 Clark Iterations: 688
 Residual min,max,rms: -1.459E-02 1.144E-02 3.455E-03
 Total CLEANed flux: 1.028E+00
 Clark Iterations: 1000
 Residual min,max,rms: -1.366E-02 9.226E-03 3.123E-03
 Total CLEANed flux: 1.316E+00
 Stopping -- Maximum iterations performed

clean% rm -rf p_148.ch0.cmp/

clean% go
Clean: version 1.0 07-Apr-04
Begin iterating
 Clark Iterations: 18
 Residual min,max,rms: -2.623E-02 3.476E-02 4.954E-03
 Total CLEANed flux: 9.487E-02
 Clark Iterations: 183
 Residual min,max,rms: -2.068E-02 1.592E-02 4.352E-03
 Total CLEANed flux: 4.040E-01
 Clark Iterations: 688
 Residual min,max,rms: -1.459E-02 1.144E-02 3.455E-03
 Total CLEANed flux: 1.028E+00
 Clark Iterations: 1000
 Residual min,max,rms: -1.366E-02 9.226E-03 3.123E-03
 Total CLEANed flux: 1.316E+00
 Stopping -- Maximum iterations performed

clean% rm -rf p_148.ch0.cmp/

clean% go
Clean: version 1.0 07-Apr-04
Begin iterating
 Clark Iterations: 18
 Residual min,max,rms: -2.623E-02 3.476E-02 4.954E-03
 Total CLEANed flux: 9.487E-02
 Clark Iterations: 183
 Residual min,max,rms: -2.068E-02 1.592E-02 4.352E-03
 Total CLEANed flux: 4.040E-01
 Clark Iterations: 688
 Residual min,max,rms: -1.459E-02 1.144E-02 3.455E-03
 Total CLEANed flux: 1.028E+00
 Clark Iterations: 1000
 Residual min,max,rms: -1.366E-02 9.226E-03 3.123E-03
 Total CLEANed flux: 1.316E+00
 Stopping -- Maximum iterations performed

clean% rm -rf p_148.ch0.cmp/

clean% go
Clean: version 1.0 07-Apr-04
Begin iterating
 Clark Iterations: 18
 Residual min,max,rms: -2.623E-02 3.476E-02 4.954E-03
 Total CLEANed flux: 9.487E-02
 Clark Iterations: 183
 Residual min,max,rms: -2.068E-02 1.592E-02 4.352E-03
 Total CLEANed flux: 4.040E-01
 Clark Iterations: 688
 Residual min,max,rms: -1.459E-02 1.144E-02 3.455E-03
 Total CLEANed flux: 1.028E+00
 Clark Iterations: 1000
 Residual min,max,rms: -1.366E-02 9.226E-03 3.123E-03
 Total CLEANed flux: 1.316E+00
 Stopping -- Maximum iterations performed

clean% rm -rf p_148.ch0.cmp/

clean% go
Clean: version 1.0 07-Apr-04
Begin iterating
 Clark Iterations: 18
 Residual min,max,rms: -1.091E+00 1.189E+00 2.955E-01
 Total CLEANed flux: 9.487E-02
 Clark Iterations: 24
 Residual min,max,rms: -8.468E+00 1.270E+01 2.343E+00
 Total CLEANed flux: 7.628E-01
 Clark Iterations: 30
 Residual min,max,rms: -1.066E+02 1.141E+02 2.617E+01
 Total CLEANed flux: 6.871E+00
 Clark Iterations: 36
 Residual min,max,rms: -8.379E+02 1.235E+03 2.284E+02
 Total CLEANed flux: 6.982E+01
 Clark Iterations: 42
 Residual min,max,rms: -9.795E+03 1.175E+04 2.555E+03
 Total CLEANed flux: 6.594E+02
 Clark Iterations: 48
 Residual min,max,rms: -8.926E+04 1.232E+05 2.237E+04
 Total CLEANed flux: 6.871E+03
 Clark Iterations: 54
 Residual min,max,rms: -1.012E+06 1.247E+06 2.615E+05
 Total CLEANed flux: 6.688E+04
 Clark Iterations: 60
 Residual min,max,rms: -9.817E+06 1.324E+07 2.551E+06
 Total CLEANed flux: 7.172E+05
 Clark Iterations: 67
 Residual min,max,rms: -1.339E+08 1.577E+08 3.240E+07
 Total CLEANed flux: 8.543E+06
 Clark Iterations: 74
 Residual min,max,rms: -1.377E+09 1.925E+09 3.690E+08
 Total CLEANed flux: 1.057E+08
 Clark Iterations: 81
 Residual min,max,rms: -2.127E+10 2.397E+10 5.045E+09
 Total CLEANed flux: 1.279E+09
 Clark Iterations: 88
 Residual min,max,rms: -2.130E+11 2.976E+11 5.731E+10
 Total CLEANed flux: 1.639E+10
 Clark Iterations: 94
 Residual min,max,rms: -2.852E+12 3.258E+12 6.835E+11
 Total CLEANed flux: 1.765E+11
 Clark Iterations: 101
 Residual min,max,rms: -2.857E+13 4.056E+13 7.786E+12
 Total CLEANed flux: 2.233E+12
 Clark Iterations: 108
 Residual min,max,rms: -4.475E+14 5.086E+14 1.065E+14
 Total CLEANed flux: 2.699E+13
 Clark Iterations: 115
 Residual min,max,rms: -4.452E+15 6.301E+15 1.209E+15
 Total CLEANed flux: 3.464E+14
 Clark Iterations: 122
 Residual min,max,rms: -6.949E+16 7.897E+16 1.654E+16
 Total CLEANed flux: 4.191E+15
 Clark Iterations: 129
 Residual min,max,rms: -6.913E+17 9.783E+17 INF
 Total CLEANed flux: 5.378E+16
 Clark Iterations: 136
 Residual min,max,rms: -1.079E+19 1.226E+19 INF
 Total CLEANed flux: 6.506E+17
 Clark Iterations: 143
 Residual min,max,rms: -1.073E+20 1.519E+20 INF
 Total CLEANed flux: 8.350E+18
 Clark Iterations: 144
 Residual min,max,rms: -3.057E+20 3.361E+20 INF
 Total CLEANed flux: 2.354E+19
 Clark Iterations: 145
 Residual min,max,rms: -7.865E+20 7.958E+20 INF
 Total CLEANed flux: 5.715E+19
 Clark Iterations: 146
 Residual min,max,rms: -2.219E+21 2.090E+21 INF
 Total CLEANed flux: 1.367E+20
 Clark Iterations: 147
 Residual min,max,rms: -4.591E+21 4.831E+21 INF
 Total CLEANed flux: 3.458E+20
 Clark Iterations: 148
 Residual min,max,rms: -1.151E+22 1.290E+22 INF
 Total CLEANed flux: 8.289E+20
 Clark Iterations: 149
 Residual min,max,rms: -2.984E+22 3.023E+22 INF
 Total CLEANed flux: 2.119E+21
 Clark Iterations: 150
 Residual min,max,rms: -8.375E+22 7.920E+22 INF
 Total CLEANed flux: 5.142E+21
 Clark Iterations: 151
 Residual min,max,rms: -1.755E+23 1.831E+23 INF
 Total CLEANed flux: 1.306E+22
 Clark Iterations: 152
 Residual min,max,rms: -4.364E+23 4.891E+23 INF
 Total CLEANed flux: 3.137E+22
 Clark Iterations: 153
 Residual min,max,rms: -1.131E+24 1.146E+24 INF
 Total CLEANed flux: 8.028E+22
 Clark Iterations: 154
 Residual min,max,rms: -3.174E+24 3.002E+24 INF
 Total CLEANed flux: 1.949E+23
 Clark Iterations: 155
 Residual min,max,rms: -6.653E+24 6.939E+24 INF
 Total CLEANed flux: 4.951E+23
 Clark Iterations: 156
 Residual min,max,rms: -1.654E+25 1.854E+25 INF
 Total CLEANed flux: 1.189E+24
 Clark Iterations: 157
 Residual min,max,rms: -4.288E+25 4.344E+25 INF
 Total CLEANed flux: 3.043E+24
 Clark Iterations: 158
 Residual min,max,rms: -1.203E+26 1.138E+26 INF
 Total CLEANed flux: 7.387E+24
 Clark Iterations: 159
 Residual min,max,rms: -2.522E+26 2.630E+26 INF
 Total CLEANed flux: 1.877E+25
 Clark Iterations: 160
 Residual min,max,rms: -6.270E+26 7.028E+26 INF
 Total CLEANed flux: 4.507E+25
 Clark Iterations: 161
 Residual min,max,rms: -1.626E+27 1.647E+27 INF
 Total CLEANed flux: 1.154E+26
 Clark Iterations: 162
 Residual min,max,rms: -4.561E+27 4.314E+27 INF
 Total CLEANed flux: 2.800E+26
 Clark Iterations: 163
 Residual min,max,rms: -9.560E+27 9.971E+27 INF
 Total CLEANed flux: 7.114E+26
 Clark Iterations: 164
 Residual min,max,rms: -2.377E+28 2.664E+28 INF
 Total CLEANed flux: 1.708E+27
 Clark Iterations: 165
 Residual min,max,rms: -6.162E+28 6.241E+28 INF
 Total CLEANed flux: 4.373E+27
 Clark Iterations: 166
 Residual min,max,rms: -1.729E+29 1.635E+29 INF
 Total CLEANed flux: 1.061E+28
 Clark Iterations: 167
 Residual min,max,rms: -3.624E+29 3.780E+29 INF
 Total CLEANed flux: 2.696E+28
 Clark Iterations: 168
 Residual min,max,rms: -9.010E+29 1.010E+30 INF
 Total CLEANed flux: 6.476E+28
 Clark Iterations: 169
 Residual min,max,rms: -2.336E+30 2.366E+30 INF
 Total CLEANed flux: 1.657E+29
 Clark Iterations: 170
 Residual min,max,rms: -6.553E+30 6.198E+30 INF
 Total CLEANed flux: 4.023E+29
 Clark Iterations: 171
 Residual min,max,rms: -1.374E+31 1.433E+31 INF
 Total CLEANed flux: 1.022E+30
 Clark Iterations: 172
 Residual min,max,rms: -3.415E+31 3.828E+31 INF
 Total CLEANed flux: 2.455E+30
 Clark Iterations: 173
 Residual min,max,rms: -8.854E+31 8.968E+31 INF
 Total CLEANed flux: 6.283E+30
 Clark Iterations: 174
 Residual min,max,rms: -2.484E+32 2.349E+32 INF
 Total CLEANed flux: 1.525E+31
 Clark Iterations: 175
 Residual min,max,rms: -5.207E+32 5.431E+32 INF
 Total CLEANed flux: 3.875E+31
 Clark Iterations: 176
 Residual min,max,rms: -1.295E+33 1.451E+33 INF
 Total CLEANed flux: 9.305E+31
 Clark Iterations: 177
 Residual min,max,rms: -3.356E+33 3.399E+33 INF
 Total CLEANed flux: 2.382E+32
 Clark Iterations: 178
 Residual min,max,rms: -9.416E+33 8.906E+33 INF
 Total CLEANed flux: 5.781E+32
 Clark Iterations: 179
 Residual min,max,rms: -1.974E+34 2.059E+34 INF
 Total CLEANed flux: 1.469E+33
 Clark Iterations: 180
 Residual min,max,rms: -4.907E+34 5.500E+34 INF
 Total CLEANed flux: 3.527E+33
 Clark Iterations: 181
 Residual min,max,rms: -1.272E+35 1.289E+35 INF
 Total CLEANed flux: 9.028E+33
 Clark Iterations: 182
 Residual min,max,rms: -3.569E+35 3.376E+35 INF
 Total CLEANed flux: 2.191E+34
 Clark Iterations: 183
 Residual min,max,rms: -7.482E+35 7.803E+35 INF
 Total CLEANed flux: 5.567E+34
 Clark Iterations: 184
 Residual min,max,rms: -1.860E+36 2.085E+36 INF
 Total CLEANed flux: 1.337E+35
 Clark Iterations: 185
 Residual min,max,rms: -4.822E+36 4.885E+36 INF
 Total CLEANed flux: 3.422E+35
 Clark Iterations: 186
 Residual min,max,rms: -1.353E+37 1.280E+37 INF
 Total CLEANed flux: 8.307E+35
 Clark Iterations: 187
 Residual min,max,rms: -2.836E+37 2.958E+37 INF
 Total CLEANed flux: 2.110E+36
 Clark Iterations: 188
 Residual min,max,rms: -7.051E+37 7.903E+37 INF
 Total CLEANed flux: 5.068E+36
 Clark Iterations: 189
 Residual min,max,rms: -1.828E+38 1.852E+38 INF
 Total CLEANed flux: 1.297E+37
 Clark Iterations: 190
 Residual min,max,rms: -INF INF INF
 Total CLEANed flux: 3.149E+37
### Fatal Error: All pixel values are identical
Received on 2007-12-05 16:16:19