int18_law151_omp_accumulation.F File Reference

#include "implicit_f.inc"
#include "comlock.inc"
#include "task_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	int18_law151_omp_accumulation (multi_fvm)

Function/Subroutine Documentation

int18_law151_omp_accumulation()

subroutine int18_law151_omp_accumulation

(

type(multi_fvm_struct), intent(inout)

multi_fvm

)

Definition at line 32 of file int18_law151_omp_accumulation.F.

!$COMMENT
!       INT18_LAW151_OMP_ACCUMULATION description
!       accumulation of force for remote nodes on the 
!       main OpenMP thread
!       only useful when OMP_NUM_THREADS > 1
!       
!       INT18_LAW151_OMP_ACCUMULATION organization :
!           loop over the NTHREADS OpenMP task + NODFI phantom nodes 
!           & accumulation in the 1rst memory cell
!$ENDCOMMENT
C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE tri7box
      USE message_mod
      USE multi_fvm_mod
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
#include      "comlock.inc"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "task_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
        TYPE(MULTI_FVM_STRUCT), INTENT(INOUT) :: MULTI_FVM
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
        INTEGER :: J,N,NN
        INTEGER :: NODFI,SHIFT_FORCE_INT
!   -----------------------------------
        IF(nthread > 1) THEN
            DO nn=1,multi_fvm%NUMBER_INT18
                n = multi_fvm%INT18_LIST(nn)
                nodfi = multi_fvm%R_AFI(n)%NODFI
                DO j=2,nthread
                    shift_force_int = (j-1)*nodfi
                    multi_fvm%R_AFI(n)%R_FORCE_INT(1:3,1:6,1:nodfi) = 
     .               multi_fvm%R_AFI(n)%R_FORCE_INT(1:3,1:6,1:nodfi) 
     .               + multi_fvm%R_AFI(n)%R_FORCE_INT(1:3,1:6,1+shift_force_int:nodfi+shift_force_int)
                    multi_fvm%R_AFI(n)%R_FORCE_INT(1:3,1:6,1+shift_force_int:nodfi+shift_force_int) = 0.d+00                        
                ENDDO
            ENDDO
        ENDIF 
 
        RETURN