diffusion_mod Module Reference

Data Types
type	t_diffusion

Functions/Subroutines
subroutine	init_diffusion (this, linsol, nrhs, iparg, ale_connect)
subroutine	solve_diffusion (this)
subroutine	get_solution (this, sol, glob_dim)
subroutine	terminate_diffusion (this)

Function/Subroutine Documentation

get_solution()

subroutine diffusion_mod::get_solution	(	class(t_diffusion), intent(inout), target	this,
		double precision, dimension(:), intent(inout), pointer	sol,
		integer, intent(out)	glob_dim )

Definition at line 208 of file diffusion_mod.F.

C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
        USE spmd_comm_world_mod, ONLY : spmd_comm_world
#include      "implicit_f.inc" 
      CLASS(T_DIFFUSION), INTENT(INOUT), TARGET :: this
      INTEGER, INTENT(OUT) :: GLOB_DIM
      
      DOUBLE PRECISION, DIMENSION(:), INTENT(INOUT), POINTER :: SOL
      sol => this%SOL
      glob_dim = this%LINEAR_SOLVER%GET_GLOBAL_DIM()

init_diffusion()

subroutine diffusion_mod::init_diffusion	(	class(t_diffusion)	this,
		class(t_linear_solver), intent(in), target	linsol,
		integer, intent(in)	nrhs,
		integer, dimension(nparg, ngroup), intent(in)	iparg,
		type(t_ale_connectivity), intent(in)	ale_connect )

Definition at line 81 of file diffusion_mod.F.

      USE linear_solver_mod
      USE ale_connectivity_mod
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
        USE spmd_comm_world_mod, ONLY : spmd_comm_world
#include      "implicit_f.inc" 
C-----------------------------------------------
C     C o m m o n   B l o c k s
C-----------------------------------------------
!     NSPMD
#include      "com01_c.inc"
!     NPARG
#include      "param_c.inc"
!     NUMELS
#include      "com04_c.inc"
!     NSVOIS
#include      "spmd_c.inc"
C-----------------------------------------------
C   M e s s a g e   P a s s i n g
C-----------------------------------------------
#include "spmd.inc"
C-----------------------------------------------
C     D u m m y   A r g u m e n t s
C-----------------------------------------------
      CLASS(T_DIFFUSION) :: this
      CLASS(T_LINEAR_SOLVER), INTENT(IN), TARGET :: LINSOL
      INTEGER, DIMENSION(NPARG, NGROUP), INTENT(IN) :: IPARG
      INTEGER, INTENT(IN) :: NRHS
      TYPE(T_ALE_CONNECTIVITY), INTENT(IN) :: ALE_CONNECT
C-----------------------------------------------
C     L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER :: MAT_NZ, NG, II, NEL, ITY, MATLAW, NFT, I, JJ, KK
      INTEGER :: GLOB_DIM, MAX_ID, IAD, LGTH
#ifdef MPI
      INTEGER IERR
#endif
 
      this%NRHS = nrhs
!     number of non zero values for the diffusion matrix elem <-> elem matrix
      mat_nz = 0
      DO ng = 1, ngroup
         matlaw = iparg(1, ng)
         IF (matlaw == 151) THEN
            nel = iparg(2, ng)
            nft = iparg(3, ng)
            ity = iparg(5, ng)
            DO ii = 1, nel
               i = ii + nft
               mat_nz = mat_nz + 1
               iad = ale_connect%EE_CONNECT%IAD_CONNECT(i)
               lgth = ale_connect%EE_CONNECT%IAD_CONNECT(i+1)-ale_connect%EE_CONNECT%IAD_CONNECT(i)
               DO jj = 1, lgth
                  kk = ale_connect%EE_CONNECT%CONNECTED(iad + jj - 1)
                  IF (kk > 0) THEN
                     mat_nz = mat_nz + 1
                  ENDIF
               ENDDO
            ENDDO
         ENDIF
      ENDDO
!     create matrix
      CALL this%MAT%MATRIX_CREATE(mat_nz)
!     create rhs
      CALL this%RHS%CREATE(nrhs * numels)
!     global system size
      glob_dim = numels
      max_id = maxval(ale_connect%IDGLOB%ID(1:numels + nsvois))
#ifdef MPI
      IF (nspmd > 1) THEN
         CALL mpi_allreduce(max_id, glob_dim, 1, mpi_int, mpi_max, spmd_comm_world, ierr)
      ENDIF
#endif
!     solution
      ALLOCATE(this%SOL(3 * glob_dim))
!     init solver
      this%LINEAR_SOLVER => linsol
      CALL this%LINEAR_SOLVER%INIT_SOLVER(glob_dim)
!     set rhs and matrix
      CALL this%LINEAR_SOLVER%SET_RHS(3, this%RHS)
      CALL this%LINEAR_SOLVER%SET_MATRIX(this%MAT)
!     diffusion coefficient, outlet flag
      this%OUTLET_FLAGGED = .false.
      IF (n2d == 0) THEN
         ALLOCATE(this%NU(numels + nsvois))
         this%NU(1:numels + nsvois) = zero
         ALLOCATE(this%FLAG_OUTLET(6 * numels))
         this%FLAG_OUTLET(1:6 * numels) = 0
      ELSE
         ALLOCATE(this%NU(numelq + numeltg + nqvois + ntgvois))
         this%NU(1:numelq + numeltg + nqvois + ntgvois) = zero
         ALLOCATE(this%FLAG_OUTLET(4 * numelq + 3 * numeltg))
         this%FLAG_OUTLET(1:4 * numelq + 3 * numeltg) = 0
      ENDIF
 

solve_diffusion()

subroutine diffusion_mod::solve_diffusion

(

class(t_diffusion)

this

)

Definition at line 187 of file diffusion_mod.F.

C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
        USE spmd_comm_world_mod, ONLY : spmd_comm_world
#include      "implicit_f.inc" 
      CLASS(T_DIFFUSION) :: this
      INTEGER :: GLOB_DIM
      
      glob_dim = this%LINEAR_SOLVER%GET_GLOBAL_DIM()
      CALL this%LINEAR_SOLVER%SOLVE(this%SOL, this%NRHS * glob_dim)
      

terminate_diffusion()

subroutine diffusion_mod::terminate_diffusion

(

class(t_diffusion), intent(inout)

this

)

Definition at line 227 of file diffusion_mod.F.

C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
        USE spmd_comm_world_mod, ONLY : spmd_comm_world
#include      "implicit_f.inc" 
      CLASS(T_DIFFUSION), INTENT(INOUT) :: this
      
    
      CALL this%LINEAR_SOLVER%TERMINATE()
      CALL this%MAT%MATRIX_DESTROY()
      CALL this%RHS%DESTROY()
      IF (ALLOCATED(this%SOL)) DEALLOCATE(this%SOL)
      IF (ALLOCATED(this%NU)) DEALLOCATE(this%NU)
      IF (ALLOCATED(this%FLAG_OUTLET)) DEALLOCATE(this%FLAG_OUTLET)