Element-by-element techniques and Parallel Programming - The Boundary Element Method with Programming - page 215

It took BiCGSTAB_L 4 iterations to converge

-0.00508 0.00508 0.00508 -0.00787 0.00123
0.00787 0.00079 -0.00079 0.02236 -0.00123
0.00787 0.00787
0.00000 0.00000 10.00000 0.00000 0.00000
10.00000 0.00000 0.00000 10.00000 0.00000
0.00000 10.00000
-0.02236 0.00079 0.00079 -0.00787 -0.00123 -
0.00787 -0.00508 -0.00508 -0.00508 -0.00787 -
0.00787 -0.00123
-10.00000 0.00000 0.00000 -10.00000
0.00000 0.00000 -10.00000 0.00000 0.00000 -
10.00000 0.00000 0.00000

The output can be seen to be the same as for Program 8.1 to 5 significant figures. The
iterative algorithm converged in 4 iterations in this case.

8.4 PROGRAM 8.3 : PARALLELISING THE

ELEMENT_BY_ELEMENT PROCEDURE

PROGRAM PARALLEL_BEM
!------------------------------------------------------
! General purpose BEM program for solving elasticity problems
! This version parallel with bicgstab(l)
!------------------------------------------------------
USE bem_lib_p; USE precision; USE timing; USE utility; USE
mp_module
USE global_variables1; USE gather_scatter6
IMPLICIT NONE ! Ndof changed to N_dof
INTEGER, ALLOCATABLE :: Inci(:,:) ! Element Incidences
INTEGER, ALLOCATABLE :: BCode(:,:), NCode(:) ! Element BC´s
INTEGER, ALLOCATABLE :: Ldest(:,:) ! Element destination vector
INTEGER, ALLOCATABLE :: Ndest(:,:) ! Node destination vector
REAL(iwp), ALLOCATABLE :: Elres_u(:,:) ! Element results , u
REAL(iwp), ALLOCATABLE :: Elres_t(:,:) ! Element results , t
REAL(iwp), ALLOCATABLE :: Elcor(:,:) ! Element coordinates
REAL(iwp), ALLOCATABLE :: xP(:,:) ! Node co-ordinates
REAL(iwp),ALLOCATABLE :: &
dUe(:,:),dTe(:,:),lhs(:,:),Diag(:,:),pmul(:)
REAL(iwp), ALLOCATABLE :: km(:,:),qmul(:), Diag1(:,:)
REAL(iwp), ALLOCATABLE :: store_dUe_pp(:,:,:),store_dTe_pp(:,:,:)
REAL(iwp), ALLOCATABLE :: F(:), F1(:) ! global RHS
REAL(iwp), ALLOCATABLE :: u1(:), y_cop(:) ! vector of unknowns
CHARACTER (LEN=80) :: Title
INTEGER :: Cdim,m,n,Nodel,Nel,N_dof,Toa, N_tot