PMD_Operator_1DD - A method to handel the user provided local operator matrix shape.


CALL PMD_Operator_1DD(Comm=comm, Dense=Ge,Operator=Op)

CALL PMD_Operator_1DD(Comm=comm, Low=L, Main=D, Upper=U,Operator=Op)


<IN>  comm     : PMD handel of type PMD_Comm_1D or PMD_Comm_2D.
<IN>  Ge       : REAL(Kind=R8 or R4) rank 2 array.
                 A general dense matrix corresponding to the local operator matrix.
<IN>  L, D, U  : REAL(Kind=R8 or R4) rank 1 or 2 array.
                 Usefull when the local operator matrix is tridiagonal. In this case, 
                 the Lower, the Main and the Upper diagonals have to be provided.
<OUT> Op       : PMD handel of type PMD_R4_Ge_Operator or PMD_R8_Ge_Operator
                 if a general dense local operator matrix has to be passed,
                 or of type PMD_R4_Tr_Operator or PMD_R8_Tr_Operator
                 if a tridiagonal local operator matrix has to be passed.
                 Returns a handel which point to the user provided local
                 operator matrix shape.


     USE PMD
     include 'mpif.h'
     INTEGER, PARAMETER              :: Nx=21, Ny=51
     TYPE(PMD_Comm_2D)               :: comm
     TYPE(PMD_R8_Ge_Operator)        :: Op
     REAL(kind=R8), DIMENSION(Nx,Ny) :: Ge
     CALL MPI_INIT( Info )
     CALL PMD_Init_1DD( MPI_COMM_WORLD, comm )

     Ge(:,:) = ...
     CALL PMD_Operator_1DD( comm, Dense=Ge, Operator=Op )


Back to the routines' list...
Last modified: Tue Nov 30 11:03:21 MET 1999