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

Elements with Neuman BC´s:
Element 600 Prescribed values:
-10.0000000000000 0.000000000000000E+000
0.000000000000000E+000
-10.0000000000000 0.000000000000000E+000
0.000000000000000E+000

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

These results are the same as those produced by Program 7.1 to 5 decimal places.

8.3 PROGRAM 8.2 : REPLACING ASSEMBLY BY THE

ELEMENT-BY-ELEMENT PROCEDURE

PROGRAM EBE_BEM
!------------------------------------------------------
! General purpose BEM program for solving elasticity problems
! This version EBE with bicgstab(l)
!------------------------------------------------------
USE bem_lib ; IMPLICIT NONE ! N_dof replaces Ndof
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(:)
REAL(iwp), ALLOCATABLE :: store_dUe(:,:,:),store_dTe(:,:,:)
REAL(iwp), ALLOCATABLE :: F(:) ! global RHS
REAL(iwp), ALLOCATABLE :: u1(:) ! global vector of unknowns