JmmcModelFittingFeedBack < Jmmc/Software

Les donnees suivantes sont partagees et [[http://jmmc.fr/oidata][mises en ligne actuellement sur le serveur d'application du JMMC]].
Chaque fichier fits partage dispose sur cette page d'une section qui peut recueillir les commentaires et retours pour completer les [[JmmcModelFitting][informations de R&D sur le logiciel]].

Si cela aide a la lisibilite, il ne faut pas hesiter a apposer la signature ( deja formatee sur le bas du formulaire de saisie ) comme par ex: 

-- Main.GuillaumeMella - 18 Mar 2008

%TOC%

---+ VLTI/MIDI observations of RY Sgr

---++ [[http://jmmc.fr/oidata/#data_VLTI/MIDI observations of RY Sgr_tabvis0525A_e.fits.gz][tabvis0525A_e.fits.gz]]
-- Main.MichelTallon - 19 Mar 2008

Erreur obtenue par Olivier Chesneau:

<verbatim>
ERROR (_lpf_optimal_step) Floating point interrupt (SIGFPE)
  LINE: 1426  FILE: /home/users/mella/ysDev/yoga/LITpro/yorick/LITpro_fitter.i
</verbatim>

L'examen du fichier données montre plusieurs problèmes:

   * Une seule base à (ucoord,vcoord)= (0.,0.). Donc quelque soit la longueur d'onde, la base est toujours nulle. Il est donc impossible d'ajuster quoi que ce soit. Néanmoins, le fitter devrait renvoyer une erreur autre que "Floating point interrupt".

   * La première table (WAVELENGTH) contient un spectre de 231 longueurs d'onde de 7.54422 à 12.745375. Puisqu'il s'agit d'observations en bande N, l'unité semble être des microns, mais le standard oidata est de fournir des mètres !

   * Une seule table OI_VIS. Alors que les données VISAMP et VISAMPERR semblent correctes, VISPHI et VISPHIERR contiennent chacun la même suite de nombres de 0.0 à 230.0. Cela ne gène pas le fitter si l'on sélectionne uniquement les données VISAMP pour l'ajustement.

---++ [[http://jmmc.fr/oidata/#data_VLTI/MIDI observations of RY Sgr_tabvis0525B_g.fits.gz][tabvis0525B_g.fits.gz]]
---++ [[http://jmmc.fr/oidata/#data_VLTI/MIDI observations of RY Sgr_tabvis0624A_a.fits.gz][tabvis0624A_a.fits.gz]]
---++ [[http://jmmc.fr/oidata/#data_VLTI/MIDI observations of RY Sgr_tabvis0624B_b.fits.gz][tabvis0624B_b.fits.gz]]
---++ [[http://jmmc.fr/oidata/#data_VLTI/MIDI observations of RY Sgr_tabvis0624C_c.fits.gz][tabvis0624C_c.fits.gz]]
---++ [[http://jmmc.fr/oidata/#data_VLTI/MIDI observations of RY Sgr_tabvis0625_d.fits.gz][tabvis0625_d.fits.gz]]
---++ [[http://jmmc.fr/oidata/#data_VLTI/MIDI observations of RY Sgr_tabvis0627_f.fits.gz][tabvis0627_f.fits.gz]]
---+ High spatial resolution observations of HM Sge (Sacuto+, 2007)

---++ [[http://jmmc.fr/oidata/#data_High spatial resolution observations of HM Sge (Sacuto+, 2007)_fig1.fits][fig1.fits]]
---++ [[http://jmmc.fr/oidata/#data_High spatial resolution observations of HM Sge (Sacuto+, 2007)_fig2.fits][fig2.fits]]
---++ [[http://jmmc.fr/oidata/#data_High spatial resolution observations of HM Sge (Sacuto+, 2007)_fig3.fits][fig3.fits]]
---++ [[http://jmmc.fr/oidata/#data_High spatial resolution observations of HM Sge (Sacuto+, 2007)_fig4.fits][fig4.fits]]
---++ [[http://jmmc.fr/oidata/#data_High spatial resolution observations of HM Sge (Sacuto+, 2007)_fig5.fits][fig5.fits]]
---++ [[http://jmmc.fr/oidata/#data_High spatial resolution observations of HM Sge (Sacuto+, 2007)_fig6.fits][fig6.fits]]
---++ [[http://jmmc.fr/oidata/#data_High spatial resolution observations of HM Sge (Sacuto+, 2007)_fig7.fits][fig7.fits]]
---++ [[http://jmmc.fr/oidata/#data_High spatial resolution observations of HM Sge (Sacuto+, 2007)_fig8.fits][fig8.fits]]
---+ Other

---++ [[http://jmmc.fr/oidata/#data_Other_2007-06-29.fits][2007-06-29.fits]]
-- Main.IsaTallonBosc - 19 Mar 2008<br />
<verbatim>
sur Target MIR0003, fit effectué avec 2 disques :
Valeurs initiales:
     name     value  prev_val    vmin    vmax   scale  fixed  units
    "i1":       0.5       0.5       0       1    AUTO      0  
    "i2":       0.5       0.5       0       1    AUTO      0  
    "r1":         0         0       0      20    AUTO      0  mas
    "r2":         0         0       0      20    AUTO      0  mas
    "x1":         0         0     -40      40    AUTO      1  mas
    "x2":         0         0     -40      40    AUTO      0  mas
    "y1":         0         0     -40      40    AUTO      1  mas
    "y2":         0         0     -40      40    AUTO      0  mas
Number of iterations  44 + 200 =244  (Max Number of iterations 200) 
Final values for fitted parameters and standard deviation:

    i1 =   0.443251 +/-  0.0311645
    i2 =   0.556749 +/-  0.0370634
    x2 =   -18.0838 +/-  0.0804189  mas
    r1 =         20 +/-    0.55042  mas
    y2 =   -26.3237 +/-   0.133774  mas
    r2 =         20 +/-   0.425638  mas

Initial Chi2 = 1.28836e+06    -  Final Chi2 = 1153.77 
Initial reduced Chi2 = 4667.95   -  Final reduced Chi2 = 4.18034 
Number of degrees of freedom = 276
Confidence Level = 1.40617e-107

--- Covariance matrix ---
         i1        i2        x2        r1        y2        r2   
  i1   9.7e-04   6.0e-04  -5.4e-04   9.1e-03   2.1e-03  -7.0e-03
  i2   6.0e-04   1.4e-03   5.4e-04  -9.1e-03  -2.1e-03   7.0e-03
  x2  -5.4e-04   5.4e-04   6.5e-03  -1.8e-02  -6.6e-03   1.3e-02
  r1   9.1e-03  -9.1e-03  -1.8e-02   3.0e-01   6.9e-02  -2.3e-01
  y2   2.1e-03  -2.1e-03  -6.6e-03   6.9e-02   1.8e-02  -5.4e-02
  r2  -7.0e-03   7.0e-03   1.3e-02  -2.3e-01  -5.4e-02   1.8e-01

--- Correlation matrix ---
      i1     i2     x2     r1     y2     r2   
  i1      1   0.52  -0.21   0.53    0.5  -0.53
  i2   0.52      1   0.18  -0.45  -0.42   0.45
  x2  -0.21   0.18      1  -0.41  -0.61   0.39
  r1   0.53  -0.45  -0.41      1   0.94     -1
  y2    0.5  -0.42  -0.61   0.94      1  -0.94
  r2  -0.53   0.45   0.39     -1  -0.94      1
  
Les rayons étant aux valeurs des bornes :
> set_parameter,wh1,"r1","r2",vmax=40
> go_fit,wh1

Number of iterations  10   (Max Number of iterations 200) 
Final values for fitted parameters and standard deviation:

    i1 =   0.584265 +/-  0.0354177
    i2 =   0.415737 +/-   0.025628
    x2 =   -20.0322 +/-   0.175057  mas
    r1 =    24.9394 +/-    0.20153  mas
    y2 =   -24.9518 +/-   0.110816  mas
    r2 =    17.4562 +/-   0.118216  mas

Initial Chi2 = 1153.77    -  Final Chi2 = 542.2 
Initial reduced Chi2 = 4.18034   -  Final reduced Chi2 = 1.96449 
Number of degrees of freedom = 276
Confidence Level = 1.59928e-19

--- Covariance matrix ---
         i1        i2        x2        r1        y2        r2   
  i1   1.3e-03   8.2e-04  -9.4e-04   1.3e-03   6.2e-04  -7.6e-04
  i2   8.2e-04   6.6e-04   9.4e-04  -1.3e-03  -6.2e-04   7.6e-04
  x2  -9.4e-04   9.4e-04   3.1e-02  -2.8e-02  -1.8e-02   1.6e-02
  r1   1.3e-03  -1.3e-03  -2.8e-02   4.1e-02   1.8e-02  -2.2e-02
  y2   6.2e-04  -6.2e-04  -1.8e-02   1.8e-02   1.2e-02  -1.1e-02
  r2  -7.6e-04   7.6e-04   1.6e-02  -2.2e-02  -1.1e-02   1.4e-02

--- Correlation matrix ---
      i1     i2     x2     r1     y2     r2   
  i1      1    0.9  -0.15   0.18   0.16  -0.18
  i2    0.9      1   0.21  -0.25  -0.22   0.25
  x2  -0.15   0.21      1  -0.78  -0.95   0.76
  r1   0.18  -0.25  -0.78      1   0.82  -0.93
  y2   0.16  -0.22  -0.95   0.82      1  -0.81
  r2  -0.18   0.25   0.76  -0.93  -0.81      1

Visualisation du fit:
</verbatim>

---+ Simulated

---++ [[http://jmmc.fr/oidata/#data_Simulated_ASPRO-STAR_1-AMBER-04-DEC-2006.oifits][ASPRO-STAR_1-AMBER-04-DEC-2006.oifits]]
---++ [[http://jmmc.fr/oidata/#data_Simulated_binary_punct.oifits][binary_punct.oifits]]
---++ [[http://jmmc.fr/oidata/#data_Simulated_gaussian_noise.oifits][gaussian_noise.oifits]]
---++ [[http://jmmc.fr/oidata/#data_Simulated_unif_disk.oifits][unif_disk.oifits]]
---++ [[http://jmmc.fr/oidata/#data_Simulated_unif_ellipse.oifits][unif_ellipse.oifits]]
---+ %NOP%BeautyContest

---++ [[http://jmmc.fr/oidata/#data_BeautyContest_2004-BSC1948I.fits][2004-BSC1948I.fits]]
---++ [[http://jmmc.fr/oidata/#data_BeautyContest_2004-FKV1137.fits][2004-FKV1137.fits]]
---++ [[http://jmmc.fr/oidata/#data_BeautyContest_2004-data1.fits][2004-data1.fits]]
---++ [[http://jmmc.fr/oidata/#data_BeautyContest_2004-data2.fits][2004-data2.fits]]
---++ [[http://jmmc.fr/oidata/#data_BeautyContest_2006-03-03.fits][2006-03-03.fits]]
---++ [[http://jmmc.fr/oidata/#data_BeautyContest_2006-03-04.fits][2006-03-04.fits]]
---++ [[http://jmmc.fr/oidata/#data_BeautyContest_2006-03-05.fits][2006-03-05.fits]]
---++ [[http://jmmc.fr/oidata/#data_BeautyContest_2006-03-06.fits][2006-03-06.fits]]
This topic: Jmmc/Software > WebHome > JmmcModelFitting > JmmcModelFittingFeedBack
Topic revision: r3 - 2008-03-19 - IsaTallonBosc