source: trunk/examples/berrada/MFLBC_RD_perfomanceMeasurement/mflbc_rd.i @ 408

#"""""""""""""""""""""""""""""""""""""""""""""""""""""""""""
#"""""""""""""""""  M.Berrada (02-2005)  """""""""""""""""""
#""""""""""" fichier de D INSTRUCTIONS mufrlbc2.i """"""""""
#""""""""""""""""""   PE LBC MultiFreq   """""""""""""""""""
#"""""""""""""""""""""""""""""""""""""""""""""""""""""""""""
# Application acoustique marine equation parabolique standard de Tappert avec condition locale d'impedence.
# Cas multifrequences où les frequences sont sur l'axe de temps.
# Ce fichier permet de faire des experiences jumelles differentes de celles dans la thèse.
# MFLBC_RD où RD est Range Dependent.
# La difference avec MFLBC au niveau du graphe de dependences est que les modules du sous-espace 1D sont passés
# à l'espace 2D. Cela permet de faire un test interessant sur l'algorithme parallèle : article HPCC.

echo on
#-----------------------------------------------------------
#-------------->  INITIALISATION DU MODELE  <---------------
xset_deltazr 0.26367188 0.48804295 |1.46
xset_c0    1520
xset_tab_freq  400 4025 450 475 500
xset_zs    93.00
xset_fac   10.00
xcoef_OBS  1  |1:NZ Obs  2:NZ/2 Obs ...
xset_init
#-----------------------------------------------------------
# etat vrai
savestate gama 1 i 0 A 0 ./data_out/gamtr.dat
savestate gama 2 i 0 A 0 ./data_out/gamti.dat
#-----------------------------------------------------------
print_time ON
set_modeltime 0
FORWARD
set_modeltime 0
print_time OFF
#-----------------------------------------------------------
#-----> sauvegarde de la trajectoire true
savestate psi 1 ij 1 A 0 ./data_out/psitr.dat
savestate psi 2 ij 1 A 0 ./data_out/psiti.dat

xnorme

#-----------------------------------------------------------
#-----> recuperation des sorties obtenues comme observations
#-----> (experience jumelle)

OUTOOBS psi 0 1 2 3 4 5

#-----------------------------------------------------------
#-----> regularisation
#set_bcoef gama 0.03
#xset_gama  0.0 0.00
#OUTOEBX gama 0 0

print_cost ON
#-----------------------------------------------------------
#---------->  PERTURBATIONS DES PARAMETRES   <--------------
xset_gama  -0.3 0.00
#xgama0 300
savestate gama 1 i 0 A 0 ./data_out/gamir.dat
savestate gama 2 i 0 A 0 ./data_out/gamii.dat
#-------------->  FIN DES PERTURBATIONS   <-----------------
print_time ON
set_modeltime 0
FORWARD
set_modeltime 0
print_time OFF
#-----------------------------------------------------------
#-----> sauvegarde de la trajectoire  init
savestate psi 1 ij 1 A 0 ./data_out/psinr.dat
savestate psi 2 ij 1 A 0 ./data_out/psini.dat

#BACKWARD # pour tester le parallélisme on ajoute un backward ici, après avoir 
          # chargé les observations et on peut effecer tout ce qu'il y a après.

GOTO FINTESTOF
print_cost off
#testof 0.1 10 10.0 10 0
testof 0.1 10 10.0 10 1
FINTESTOF

#""""""""""""""""""""""""""""""""""""""""""""""""""""""""""
#-------------> choix du run
#-------> 1. RUNYAO
#-------> 2. RUNM2QN1
#-------> 3. RUNM1QN3
#-----------------------------------------------------------
GOTO 3
#------

# le RUN standard de YAO ..................................
1
setepsi  gama  0.001
set_nbiter 50

RUN
goto FINRUN

# le RUN sous M1QN3 .......................................
3
setm_impres  5
setm_io      6
setm_mode    0
set_nbiter   50
setm_nsim    50
setm_dxmin   1.0e-20
setm_epsg    1.0e-20
setm_ddf1    1.0

runm
goto FINRUN

FINRUN
xnorme
#-----------------------------------------------------------
#-----> sauvegarde de la trajectoire calculee
savestate psi 1 ij 1 A 0 ./data_out/psyr.dat
savestate psi 2 ij 1 A 0 ./data_out/psyi.dat

savestate gama 1 i 0 A 0 ./data_out/gamyr.dat
savestate gama 2 i 0 A 0 ./data_out/gamyi.dat
#============================= FIN =========================