GET_NEIGHBOUR_ELEM.f90 File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	get_neighbour_elem (omega_m, npt, n, id_mpi, n_qp)
	Computes neighbouring element.

Function/Subroutine Documentation

get_neighbour_elem()

subroutine get_neighbour_elem	(	integer*4, dimension(n_qp,0:3), intent(inout)	omega_m,
		integer*4, intent(in)	npt,
		integer*4, intent(out)	n,
		integer*4, intent(in)	id_mpi,
		integer*4, intent(in)	n_qp
	)

Computes neighbouring element.

Author

Ilario Mazzieri

Date

September, 2013

Version

1.0

Parameters

[in]	omega_m	omega_m(:,0) = block identification for the neighbouring element, omega_m(:,1) = neighbouring element number, omega_m(:,2) = neighbouring element face
[in]	npt	number of quadrature points lying on the neighbouring element
[in]	id_mpi	number of process mpi
[in]	n_qp	number of quadrature points
[out]	omega_m(,3)	label identifying the neighbournig element for each quadrature point
[out]	n	total neighbouring elements

Definition at line 32 of file GET_NEIGHBOUR_ELEM.f90.

          
     implicit none
     
     integer*4 :: i, lett, j     
     integer*4, intent(in) :: npt, id_mpi, n_qp
     integer*4, intent(out)  :: n
 
     integer*4, dimension(n_qp) :: omega_m_copy
     integer*4, dimension(n_qp,0:3), intent(inout) :: omega_m 
     
     
     n = 0
     i = 1
     omega_m_copy = omega_m(:,1)
     do while (i .le. npt)
        lett = omega_m_copy(i)
        
        if(lett .ne. 0) then
           n = n + 1
           omega_m(i,3) = n
        endif
        
        do j = i + 1, npt
           
           if (omega_m(j,1) .eq. lett) then
             omega_m(j,3) = n
             omega_m_copy(j) = 0
           endif
           
        enddo
        
        i = i + 1
        
     enddo
     
   
     
     
     

Referenced by make_dg_interface().

Here is the caller graph for this function: