Accueil | Cast3M

   1 : $$$$ GANE     NOTICE  CHAT      11/09/12    21:16:22     7124           
   2 :                                              DATE     11/09/12
   3 : 
   4 :   Operateur GA(us-)NE(wton)               Voir aussi : MOCA, LEVM
   5 :     -------------------------                            AJUSTE  EXCE  
   6 :     TAB2=GANE TAB1 ('AMOR' FLOT1);
   7 :     
   8 :     CHPO1 RIGI1=GANE TAB1 'MATR' ('AMOR' FLOT1);
   9 : 
  10 :     Objet :
  11 :     _______
  12 : 
  13 :     L'operateur GANE construit la matrice et le second membre de la
  14 :     methode de Gauss-Newton ou de Levenberg-Marquardt.
  15 : 
  16 :     - En donnant le mot clef 'MATR', ces objets sont retournes dans
  17 :       CHPO1 (type 'CHPOINT') and RIGI1 (type 'RIGIDITE').
  18 : 
  19 :     - Sans ce mot clef, le syteme lineaire est resolu et le resultat
  20 :       (direction de descente) stocke dans TAB2 (type 'TABLE').
  21 : 
  22 :     - Le mot clef 'AMOR' permet d'introduire le parametre de viscosite
  23 :       nu=FLOT1 (type 'FLOTTANT') de la methode L-M . Par defaut
  24 :       nu=0 (methode G-N).
  25 : 
  26 :     La fonction a minimiser est: 2F(X)={f(X)}.{f(X)}, and [J(X)]=df/dX
  27 :     
  28 :     La direction de descente est H solution de:
  29 :       [[transpose([J])[J]]+[nu*I]]{H}  = -{transpose([J]){f}}
  30 : 
  31 :     Un exemple d'utilisation est donne dans gane.dgibi.
  32 : 
  33 :     Commentaire :
  34 :     _____________
  35 : 
  36 :     TAB1   :  Table de type 'VECTEUR' contenant
  37 :               TAB1 . 0 = f(X)
  38 :               TAB1 . 1 = df/dX1
  39 :                    .
  40 :               TAB1 . n = df/dXn
  41 :               
  42 :     FLOT1  :  'FLOTTANT'=nu
  43 :               
  44 :     CHPO1  :  -{transpose([J]){f}}
  45 :     RIGI1  :  [transpose([J])[J]]
  46 :     
  47 :     TAB2   :  TABLE of type 'VECTEUR' contenant
  48 :               TAB2 . 1 = H1
  49 :                    .           
  50 :               TAB2 . n = Hn
  51 :