Datasets examples

* fichier :  ssch.dgibi
 SAUT PAGE ;
************************************************************
*
*  test elementaire de l'operateur ssch
*
************************************************************

opti dime 2 ;

pa = 0.0 0.0 ;
pb = 1.0 0.0 ;
pc = 1.0 1.0 ;
pd = 0.0 1.0 ;

su = manu qua4 pa pb pc pd ;
sup = chan su poi1 ;

tabdon = table ;
tabdon . iden = lect 1 101 ;
tabchi = chi1 tabdon COMP '/u/castem/divers/COMPOM'
                     LOGK '/u/castem/divers/COMPOM' ;

tbparm = table ;
tbparm . itmax = 40 ;
tbparm . eps   = 1.d-4 ;
tbparm . sortie = mots 'SOLU' ;

tchi = table ;
tchi . soustype = donnees_chimiques ;
tchi . tot  = manu chpo sup 2 x001 (prog  5.d-5 5.d-4 1.d-3 2.d-5 )
                              x101 (prog  5.d-5 1.d-3 5.d-4 1.d-5 ) ;
tchi . logc = manu chpo sup 2 x001 -4.
                              x101 -4.             ;

tachi = chi2 tabchi tbparm tchi ;

dfi = manu chpo sup 2 x001  (prog 0.d0  1.d-5 2.e-5 0.d0)
                      x101  (prog 0.d0  1.d-5 3.e-5 0.d0) ;
cc = tachi . solu     ;
dcc = manu chpo sup 3 w001 (prog 1.e-6  0.0 0.0  2.e-6)
                      w002 (prog 2.e-6  0.0 0.0  3.e-6)
                      w003 (prog 3.e-6  0.0 0.0  1.e-6) ;
res  = ssch tabchi dfi cc dcc ;


tes = manu chpo sup 3 w001 (prog  1.d-6 6.4432d-6 9.4077d-6   3.d-6)
                      w002 (prog  1.d-6 6.4432d-6 1.94077d-5  3.d-6)
                      w003 (prog -1.d-6 3.5567d-6 1.05923d-5 -3.d-6);

ver  = ( abs(res - tes) )masque superieur 1.d-5 somme ;

si (ver ega 0) ;
  erre 0 ;
sinon ;
  erre 5 ;
finsi ;

fin ;