Accueil | Cast3M

   1 : $$$$ KPRE     NOTICE  GOUNAND   25/03/12    21:15:05     12194          
   2 :                                              DATE     25/03/12
   3 : 
   4 :   Operateur KPRE                          Voir aussi : CMCT
   5 :     --------------  
   6 : 
   7 :     Objet :
   8 :     _______
   9 : 
  10 :     L'operateur KPRE construit une matrice de preconditionnement.
  11 : 
  12 :     Etant donnee la matrice 2x2 par blocs  ( D  B^t )
  13 :                                            ( C  E   )
  14 :     ou les sous-matrices B, C, D et E sont des matrices definies sur les
  15 :     elements de MAIL1 (D et E carrees), l'operateur KPRE construit la
  16 :     matrice :
  17 :                      S* = \sum_e  E_e - C_e D_e^-1 B_e^t
  18 :     ou e sont les elements de MAIL1.
  19 : 
  20 :     La matrice S* peut etre vue comme une approximation du complement
  21 :     de Schur :
  22 :                      S = E - C D^-1 B^t
  23 :     definie element par element.
  24 : 
  25 :     Commentaire :
  26 :     _____________
  27 : 
  28 :       RIGC, RIGD   : Matrices (type RIGIDITE)
  29 :       RIGB, RIGE
  30 : 
  31 :       MAIL1  : Maillage sous-tendant les matrices (type MAILLAGE)
  32 : 
  33 :       RIG2   : Matrice resultat S* (type RIGIDITE)
  34 : 
  35 :     Remarques :
  36 :     ___________
  37 : 
  38 :      1) Si RIGE n'est pas precisee, le bloc correspondant est nul.
  39 : 
  40 :      2) Si MAIL1 n'est pas precise, on prend le maillage sous-tendant
  41 :         la matrice D.
  42 :