Simulated data from Gilles Duvert sent for testing yoga and the GUI
developped at LAOG by Guillaume Mella.

Uniform disk:
    i = 1
    x = 0
    y = 0
    r = 1.5

--------------------------------------------------------------------------
Original mail (2006.12.04) from Gilles:
<<

Bonjour,

Je viens de tester la version de developpement GUI+yoga, et c'est assez
sympa je trouve, qu'en pensez-vous?
Maintenant, j'ai un problme: Je lui ai fait manger un fichier OI-FITS
qui est une simulation VLTI+AMBER sortie d'ASPRO, avec un disque
stellaire de 3mas.... et... quand je fais tourner yoga, il ne converge
pas? cad, si je rentre comme diamtre de dŽpart zŽro, il sort zŽro, si
je rentre 2 mas il sort 2 mas, etc.. Seule le premier paramtre (flux)
est fittŽ ?
Est-ce que vous pouvez regarder ce qui se passe? Je vous joins le xml
qui encapsule ˆ la fois le fichier fits et le setup que j'ai employŽ, il
n'y a plus qu'ˆ le loader... La version de demo d'aujourd'hui permet la
sauvegarde du fichier oi-fits contenu dans le xml.

>>

--------------------------------------------------------------------------

Actually, the data seems quite good. There is no problem when using 
directly yoga (2006.12.05)

<<
Isa et moi venons de regarder ce problme. Cela semble marcher
parfaitement ˆ la main. Pour chercher l'erreur, une possibilitŽ pourrait
tre d'interrompre le batch de yorick (il semble que ce soit possible!) et 
jeter un oeil aux donnŽes lues (plot_radial, world, ...).
Isa et Michel

1/ Extraction du fichier oifits du xml de Gilles.

(fichier ASPRO-STAR_1-AMBER-04-DEC-2006.oifits)

2/ Copie du fichier settings.i listŽ lors de l'exŽcution via l'interface. 
Juste le nom du fichier est changŽ par le nom du fichier extrait du xml.

_LPP_MODEL_SETTINGS = h_new(
  target = h_new(
    TG1 = h_new(
      ident = "STAR_1",
      files = YOGA_SITE+"data/Simul_data/ASPRO-STAR_1-AMBER-04-DEC-2006.oifits",
    /* Indicates files corresponding to given ids 
     * - id0 file is /home/gildas/ASPRO-STAR_1-AMBER-04-DEC-2006.fits
     */
      OBJ1 = h_new(
        name = "disk",
        params = h_new(weight = "i1", x = "x1", y = "y1", radius = "r1")
      ),
      residuals = h_new(
        NORMALIZE = TRUE
      )
    )
  ),
  param = h_new(
    i1 = h_new(value=0.5, vmin=0.0, vmax=1.0, units=""),
    x1 = h_new(value=0.0, vmin=-40.0, vmax=40.0, units="mas"),
    y1 = h_new(value=0.0, vmin=-40.0, vmax=40.0, units="mas"),
    r1 = h_new(value=0.0, vmin=0.0, vmax=20.0, units="mas")
  ),
  fitter = "standard"
);



3/ Session yoga. Le rŽsultat est tout ˆ fait satisfaisant.

> world = load_model(LIT_MODELS+"settings.i");
 --- Loading model description in "settings.i" and data...
skipping empty IMAGE extension in unit 1
loading OI_ARRAY extension data from unit 2
loading OI_TARGET extension data from unit 3
loading OI_WAVELENGTH extension data from unit 4
loading OI_VIS extension data from unit 5
loading OI_VIS2 extension data from unit 6
loading OI_T3 extension data from unit 7
FIXME: OIFITS_ARRAY not yet implemented
WARNING: all spectral bandwidth(s) are zero!
 --- Finalizing data fields in world...
 --- Finalizing configuration for fitting...
 --- Ready to fit.
> show_parameters, world
     name     value  prev_val    vmin    vmax   scale  fixed  units
    "i1":       0.5       0.5       0       1    AUTO      0  
    "r1":         0         0       0      20    AUTO      0  mas
    "x1":         0         0     -40      40    AUTO      0  mas
    "y1":         0         0     -40      40    AUTO      0  mas
> go_fit, world;
 Chi2:30192.5  
 no progress

Number of iterations  33   (Max Number of iterations 200)
 
Final values for fitted parameters and standard deviation:
    i1 =          1 +/-  0.0741249
    x1 = -2.17996e-21 +/- 0.000247825
    y1 = 5.15966e-22 +/- 0.000284445
    r1 =        1.5 +/- 0.00620989

First Chi2 = 30192.5    -  Last Chi2 = 1.59607e-10 
Confidence Level = 0

--- Covariance matrix ---
         i1        x1        y1        r1   
  i1    0.0055   6.5e-27  -9.9e-26   1.2e-12
  x1   6.5e-27   6.1e-08  -4.9e-08   2.5e-18
  y1  -9.9e-26  -4.9e-08   8.1e-08  -1.6e-18
  r1   1.2e-12   2.5e-18  -1.6e-18   3.9e-05

--- Correlation matrix ---
         i1       x1        y1        r1   
  i1         1  3.5e-22  -4.7e-21   2.6e-09
  x1   3.5e-22        1     -0.69   1.6e-12
  y1  -4.7e-21    -0.69         1  -9.2e-13
  r1   2.6e-09  1.6e-12  -9.2e-13         1
> 
--------------------------------------------------------------------------