Exemples
========


Résolution par la méthode directe MULT_FRONT
--------------------------------------------


*   Constitution de la matrice assemblée et du second membre :

On a calculé auparavant les termes élémentaires KEL, FEL.


.. code-block:: text
NU    =NUME_DDL( MATR_RIGI=KEL )
K    =ASSE_MATRICE( MATR_ELEM=KEL,    NUME_DDL=NU, )
F    =ASSE_VECTEUR( MATR_ELEM=FEL,    NUME_DDL=NU, )


*   Factorisation:


.. code-block:: text
K    =FACTORISER( reuse=K,    MATR_ASSE=K, METHODE='MULT_FRONT',)


*   Résolution:


.. code-block:: text
U    =RESOUDRE( MATR=K,    CHAM_NO=F,)


*   pour l’utilisation des charges cinématiques (avec élimination des degrés de liberté imposés), voir l’exemple donné dans la commande AFFE_CHAR_CINE [:external:ref:`U4.44.03 <U4.44.03>`].


Résolution par la méthode MUMPS
-------------------------------


.. code-block:: text

    NU        =    NUME_DDL( MATR_RIGI= KEL)


    K        =    ASSE_MATRICE    (    MATR_ELEM= KEL, NUME_DDL= NU)

    F        =    ASSE_VECTEUR    (    VECT_ELEM= FEL, NUME_DDL= NU )

    K        =    FACTORISER    (    reuse= K, MATR_ASSE= K, METHODE= 'MUMPS' )

    dep        =    RESOUDRE        (    CHAM_NO = F    ,    MATR= K)


Résolution par la méthode du gradient conjugué pré conditionné
--------------------------------------------------------------


.. code-block:: text

    NU        =    NUME_DDL( MATR_RIGI= KEL)


    K        =    ASSE_MATRICE    (    MATR_ELEM= KEL, NUME_DDL= NU )

    F        =    ASSE_VECTEUR    (    VECT_ELEM= FEL, NUME_DDL= NU )

    KPREC    =    FACTORISER    (    MATR_ASSE= K, METHODE= 'GCPC',
PRE_COND='LDLT_INC' )

    dep        =    RESOUDRE        (    CHAM_NO = F    ,    MATR= K, MATR_PREC= KPREC,
NMAX_ITER= 1000    ,    RESI_RELA= 1e-07
)


Résolution par la méthode PETSC
-------------------------------


.. code-block:: text

    NU        =    NUME_DDL( MATR_RIGI= KEL)


    K        =    ASSE_MATRICE    (    MATR_ELEM= KEL, NUME_DDL= NU )

    F        =    ASSE_VECTEUR    (    VECT_ELEM= FEL, NUME_DDL= NU )

    K        =    FACTORISER    (    reuse=K, MATR_ASSE= K, METHODE= 'PETSC')

    dep        =    RESOUDRE        (    CHAM_NO = F    ,    MATR= K, MATR_PREC= K,
ALGORITHME='GMRES',
NMAX_ITER= 1000    ,    RESI_RELA= 1e-07)