Iterative Deconvolution Algorithm in C - idac

The Diagnostic File


The diagnostic file is very important for checking your input into the idac program. It includes not only the start date and time of the program but also the stop date and time and the stopping criteria computed from the `N0' parameters. The mean value for each frame read in is included as well as the volume normalization values.

A sample is included below. It is used to verify the correct reads of the input files as well as giving the normalization parameters. Highlighted in red are the inputs from the FITS header of the convolution file. Highlighted in blue is the computed value of the Global Stopping Criterion at the end of the file. The algorithm will terminate when this value has been reached (see the Error Metric page).

*************************************************
******  IDAC Version 2.7 12 August 2000  **********
*************************************************


Program Start Date and Time: Mon Dec 11 13:37:41 2000

COMMAND LINE:
        /u/christou/PROGS/idac/bin/idac_v27 conv obj psf wt 1000 10 1 1 0 9 1 

# -------------- CONTROL BITS ----------------
# 0 : PSF/BL constraint    on
# 1 : PSF[saa] contraint   off
# 2 : Fourier constraint  off
# 3 : Object/PSF fixed    off
# 4 : Sky fitting      off
# --------------------------------------------
# conv = 64 x 64 x 10, Type 7
# ALPHA = 1.000000e+00
# PA = 0.000000e+00
# Stopping Criteria Frame 0 = 6.253000e-01
# Stopping Criteria Frame 1 = 0.000000e+00
# Stopping Criteria Frame 2 = 0.000000e+00
# Stopping Criteria Frame 3 = 0.000000e+00
# Stopping Criteria Frame 4 = 0.000000e+00
# Stopping Criteria Frame 5 = 0.000000e+00
# Stopping Criteria Frame 6 = 0.000000e+00
# Stopping Criteria Frame 7 = 0.000000e+00
# Stopping Criteria Frame 8 = 0.000000e+00
# Stopping Criteria Frame 9 = 0.000000e+00
# Reading and writing 10 frames = 40960 pixels
# R0X = 2.400000e+01
# R0Y = 2.400000e+01
# conv_sup = 64 x 64 x 1
# input conv_sup = 64 x 64 x 1
# output will be 64 x 64 x 10
# making 10 copies of frame 1 data
# mean value of conv_sup[,,0] = 1.000000e+00
# mean value of conv_sup[,,1] = 1.000000e+00
# mean value of conv_sup[,,2] = 1.000000e+00
# mean value of conv_sup[,,3] = 1.000000e+00
# mean value of conv_sup[,,4] = 1.000000e+00
# mean value of conv_sup[,,5] = 1.000000e+00
# mean value of conv_sup[,,6] = 1.000000e+00
# mean value of conv_sup[,,7] = 1.000000e+00
# mean value of conv_sup[,,8] = 1.000000e+00
# mean value of conv_sup[,,9] = 1.000000e+00
# conv_sup read in correctly...
# Convolution Data: conv = [64,64,10]
# total frames to process = 10
# begin frame = 0
# end frame = 9
# --------------------------------------------
# sky = 64 x 64 x 1
# input sky = 64 x 64 x 1
# output will be 64 x 64 x 10
# making 10 copies of frame 1 data
# mean value of sky[,,0] = 0.000000e+00
# mean value of sky[,,1] = 0.000000e+00
# mean value of sky[,,2] = 0.000000e+00
# mean value of sky[,,3] = 0.000000e+00
# mean value of sky[,,4] = 0.000000e+00
# mean value of sky[,,5] = 0.000000e+00
# mean value of sky[,,6] = 0.000000e+00
# mean value of sky[,,7] = 0.000000e+00
# mean value of sky[,,8] = 0.000000e+00
# mean value of sky[,,9] = 0.000000e+00
# sky read in correctly...
# CONV: 
        0 : Normalization Factor = 4.463885e+03
        1 : Normalization Factor = 4.427188e+03
        2 : Normalization Factor = 4.481470e+03
        3 : Normalization Factor = 4.442515e+03
        4 : Normalization Factor = 4.474363e+03
        5 : Normalization Factor = 4.369575e+03
        6 : Normalization Factor = 4.447770e+03
        7 : Normalization Factor = 4.445880e+03
        8 : Normalization Factor = 4.497589e+03
        9 : Normalization Factor = 4.510292e+03
# SKY:
        0 : Normalization Factor = 4.463885e+03
        1 : Normalization Factor = 4.427188e+03
        2 : Normalization Factor = 4.481470e+03
        3 : Normalization Factor = 4.442515e+03
        4 : Normalization Factor = 4.474363e+03
        5 : Normalization Factor = 4.369575e+03
        6 : Normalization Factor = 4.447770e+03
        7 : Normalization Factor = 4.445880e+03
        8 : Normalization Factor = 4.497589e+03
        9 : Normalization Factor = 4.510292e+03
# We ALWAYS read and subtract the Sky...
# --------------------------------------------
# obj = 64 x 64 x 1
# input obj = 64 x 64 x 1
# output will be 64 x 64 x 1
# mean value of obj[,,0] = 1.087904e+00
# obj read in correctly...
# obj_sup = 64 x 64 x 1
# input obj_sup = 64 x 64 x 1
# output will be 64 x 64 x 1
# mean value of obj_sup[,,0] = 1.000000e+00
# obj_sup read in correctly...
Iteration 0
        0 : Normalization Factor = 4.792589e+03
# --------------------------------------------
# psf = 64 x 64 x 1
# input psf = 64 x 64 x 1
# output will be 64 x 64 x 10
# making 10 copies of frame 1 data
# mean value of psf[,,0] = 2.421114e-04
# mean value of psf[,,1] = 2.421114e-04
# mean value of psf[,,2] = 2.421114e-04
# mean value of psf[,,3] = 2.421114e-04
# mean value of psf[,,4] = 2.421114e-04
# mean value of psf[,,5] = 2.421114e-04
# mean value of psf[,,6] = 2.421114e-04
# mean value of psf[,,7] = 2.421114e-04
# mean value of psf[,,8] = 2.421114e-04
# mean value of psf[,,9] = 2.421114e-04
# psf read in correctly...
# psf_sup = 64 x 64 x 1
# input psf_sup = 64 x 64 x 1
# output will be 64 x 64 x 10
# making 10 copies of frame 1 data
# mean value of psf_sup[,,0] = 1.000000e+00
# mean value of psf_sup[,,1] = 1.000000e+00
# mean value of psf_sup[,,2] = 1.000000e+00
# mean value of psf_sup[,,3] = 1.000000e+00
# mean value of psf_sup[,,4] = 1.000000e+00
# mean value of psf_sup[,,5] = 1.000000e+00
# mean value of psf_sup[,,6] = 1.000000e+00
# mean value of psf_sup[,,7] = 1.000000e+00
# mean value of psf_sup[,,8] = 1.000000e+00
# mean value of psf_sup[,,9] = 1.000000e+00
# psf_sup read in correctly...
# --------------------------------------------
# Not Reading SAA image...
# **** PSF normalization ****
Iteration 0
        0 : Normalization Factor = 9.916882e-01
        1 : Normalization Factor = 9.916882e-01
        2 : Normalization Factor = 9.916882e-01
        3 : Normalization Factor = 9.916882e-01
        4 : Normalization Factor = 9.916882e-01
        5 : Normalization Factor = 9.916882e-01
        6 : Normalization Factor = 9.916882e-01
        7 : Normalization Factor = 9.916882e-01
        8 : Normalization Factor = 9.916882e-01
        9 : Normalization Factor = 9.916882e-01
# --------------------------------------------
# wt = 64 x 64 x 1
# input wt = 64 x 64 x 1
# output will be 64 x 64 x 10
# making 10 copies of frame 1 data
# mean value of wt[,,0] = 1.000000e+00
# mean value of wt[,,1] = 1.000000e+00
# mean value of wt[,,2] = 1.000000e+00
# mean value of wt[,,3] = 1.000000e+00
# mean value of wt[,,4] = 1.000000e+00
# mean value of wt[,,5] = 1.000000e+00
# mean value of wt[,,6] = 1.000000e+00
# mean value of wt[,,7] = 1.000000e+00
# mean value of wt[,,8] = 1.000000e+00
# mean value of wt[,,9] = 1.000000e+00
# wt read in correctly...
# --------------------------------------------
# conv_wf = 64 x 64 x 1
# input conv_wf = 64 x 64 x 1
# output will be 64 x 64 x 10
# making 10 copies of frame 1 data
# mean value of conv_wf[,,0] = 1.000000e+00
# mean value of conv_wf[,,1] = 1.000000e+00
# mean value of conv_wf[,,2] = 1.000000e+00
# mean value of conv_wf[,,3] = 1.000000e+00
# mean value of conv_wf[,,4] = 1.000000e+00
# mean value of conv_wf[,,5] = 1.000000e+00
# mean value of conv_wf[,,6] = 1.000000e+00
# mean value of conv_wf[,,7] = 1.000000e+00
# mean value of conv_wf[,,8] = 1.000000e+00
# mean value of conv_wf[,,9] = 1.000000e+00
# conv_wf read in correctly...
# obj: 4096, psf: 40960 sky: 0
# --------------------------------------------
# ---------------INPUT COMPLETE---------------
# --------------------------------------------
Global Stopping Criterion = 8.065590e-04
# total allocated global bytes = 2605576
# total allocated function bytes = 3047468
# Stopping Criterion Met.  Quitting...
Program End Date and Time: Mon Dec 11 13:39:32 2000



Written by J. C. Christou
Updated:  11 December 2000 (JCC)