Télécharger la plaquette Cast3M


fr en

Télécharger ssch.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  ssch.dgibi
  2.  SAUT PAGE ;
  3. ************************************************************
  4. *
  5. *  test elementaire de l'operateur ssch
  6. *
  7. ************************************************************
  8.  
  9. * repertoire des fichiers "divers"
  10. DIVERS = VENV 'CASTEM_DIVERS';
  11. *
  13.  
  14. pa = 0.0 0.0 ;
  15. pb = 1.0 0.0 ;
  16. pc = 1.0 1.0 ;
  17. pd = 0.0 1.0 ;
  18.  
  19. su = manu qua4 pa pb pc pd ;
  20. sup = chan su poi1 ;
  21.  
  22. tabdon = table ;
  23. tabdon . iden = lect 1 101 ;
  24. tabchi = chi1 tabdon COMP ('CHAINE' DIVERS '/COMPOM')   
  25.                      LOGK ('CHAINE' DIVERS '/COMPOM') ;
  26.  
  27. tbparm = table ;
  28. tbparm . itmax = 40 ;
  29. tbparm . eps   = 1.d-4 ;
  30. tbparm . sortie = mots 'SOLU' ;
  31.  
  32. tchi = table ;
  33. tchi . soustype = donnees_chimiques ;
  34. tchi . tot  = manu chpo sup 2 x001 (prog  5.d-5 5.d-4 1.d-3 2.d-5 )
  35.                               x101 (prog  5.d-5 1.d-3 5.d-4 1.d-5 ) ;
  36. tchi . logc = manu chpo sup 2 x001 -4.
  37.                               x101 -4.             ;
  38.  
  39. tachi = chi2 tabchi tbparm tchi ;
  40.  
  41. dfi = manu chpo sup 2 x001  (prog 0.d0  1.d-5 2.e-5 0.d0)
  42.                       x101  (prog 0.d0  1.d-5 3.e-5 0.d0) ;
  43. cc = tachi . solu     ;
  44. dcc = manu chpo sup 3 w001 (prog 1.e-6  0.0 0.0  2.e-6)
  45.                       w002 (prog 2.e-6  0.0 0.0  3.e-6)
  46.                       w003 (prog 3.e-6  0.0 0.0  1.e-6) ;
  47. res  = ssch tabchi dfi cc dcc ;
  48.  
  49.  
  50. tes = manu chpo sup 3 w001 (prog  1.d-6 6.4432d-6 9.4077d-6   3.d-6)
  51.                       w002 (prog  1.d-6 6.4432d-6 1.94077d-5  3.d-6)
  52.                       w003 (prog -1.d-6 3.5567d-6 1.05923d-5 -3.d-6);
  53.  
  54. ver  = ( abs(res - tes) )masque superieur 1.d-5 somme ;
  55.  
  56. si (ver ega 0) ;
  57.   erre 0 ;
  58. sinon ;
  59.   erre 5 ;
  60. finsi ;
  61.  
  62. fin ;
  63.  
  64.  
  65.  
  66.  
  67.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales