um1x.f

Go to the documentation of this file.

!                   ***************
                    SUBROUTINE um1x
!                   ***************
!
     &(x,d,s)
!
!***********************************************************************
! BIEF   V6P1                                   21/08/2010
!***********************************************************************
!
!brief    END OF THE OPERATIONS FOR DIAGONAL PRECONDITIONING
!+                WITH THE CONDENSED MATRIX  (OR ANY SIMILAR OPERATION).
!code
!+                   -1   PRIME
!+    OPERATION X = U    X
!+
!+    EXAMPLE OF A BLOCK OF 4:
!+
!+              (   I     D12  )
!+              (              )
!+         U =  (              )
!+              (              )
!+              (   0      I   )
!+                                   PRIME          PRIME         PRIME
!+              (   I    -D12  ) ( X1      )    ( X1    -   D12 X2 )
!+  -1  PRIME   (              ) (         )    (                  )
!+ U   X     =  (              ) (         )  = (                  )
!+              (              ) (   PRIME )    (                  )
!+              (   0      I   ) ( X2      )    ( X2               )
!
!history  J.M. HERVOUET (LNH)
!+        23/12/94
!+        V5P1
!+
!
!history  N.DURAND (HRW), S.E.BOURBAN (HRW)
!+        13/07/2010
!+        V6P0
!+   Translation of French comments within the FORTRAN sources into
!+   English comments
!
!history  N.DURAND (HRW), S.E.BOURBAN (HRW)
!+        21/08/2010
!+        V6P0
!+   Creation of DOXYGEN tags for automated documentation and
!+   cross-referencing of the FORTRAN sources
!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!| D              |<--| BLOCK OF DIAGONAL MATRICES
!| S              |-->| 2 : BLOCK WITH 4 MATRICES
!|                |   | 3 : BLOCK WITH 9 MATRICES
!| X              |<->| X AND X' (IN SITU TRANSFORMATION)
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
      USE bief, ex_um1x => um1x
!
      USE declarations_special
      IMPLICIT NONE
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
!     STRUCTURES OF VECTORS OR BLOCKS OF VECTORS
!
      INTEGER, INTENT(IN)           :: S
      TYPE(bief_obj), INTENT(INOUT) :: X
      TYPE(bief_obj), INTENT(IN)    :: D
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
!     BLOCKS OF 4 MATRICES:
!
      IF(s.EQ.2) THEN
!
!     BLOCKS OF 4 MATRICES:
!
        CALL um1x04(x%ADR(1)%P,x%ADR(2)%P,d%ADR(3)%P)
!
      ELSEIF(s.EQ.3) THEN
!
!     BLOCKS OF 9 MATRICES:
!
        CALL um1x09(x%ADR(1)%P,x%ADR(2)%P,x%ADR(3)%P,
     &              d%ADR(4)%P,d%ADR(5)%P,d%ADR(7)%P)
!
      ELSE
!
!-----------------------------------------------------------------------
!
!  ERROR
!
        WRITE(lu,200) s
200     FORMAT(1x,'UM1X (BIEF) : UNEXPECTED S :',1i6)
        CALL plante(1)
        stop
!
      ENDIF
!
!-----------------------------------------------------------------------
!                                                            -1
      RETURN
      END