Civil Engineering Reference
In-Depth Information
New dynamic real arrays:
col
column array
kc
conductivity matrix
x
contribution to conductivity matrix
kcx
y
contribution to conductivity matrix
kcy
contribution to conductivity matrix
row
row array
store pp
distributed penalty storage
val f
z
kcz
values of fixed freedoms
The closest serial equivalent to this program is Program 7.5. Geometrical restrictions
are to a quarter of a cuboidal box with zero potential on all its outer faces. A potential of
100.0 units can be fixed at the centre of the box or a flux of 10.0 units applied there, so
either
loaded_freedoms
or
fixed_freedoms
must be set to 1 and the other to 0.
Comparison of Programs 7.5 and 12.1 will show the familiar patterns of analysis type
and parallelisation process respectively. In the serial programs, loop
elements_1
con-
tains integration loop
gauss_pts_1
embedded within it whereas in the parallel version
these are separate as
elements_0
and
gauss_pts_1
. An extra loop
elements_1
is
necessary for formation of the preconditioner. For simplicity, in parallel a single
kc
matrix
is used. Loop
elements_2
is equivalent in both versions but in parallel, flow rates are
nels nxe nze nip
1000000 100 100 8
aa bb cc
0.01 0.01 0.01
kx ky kz
2.0 2.0 2.0
tol limit
1.0e-5 500
loaded_freedoms fixed_freedoms
1 0
Figure 12.17
Data for Program 12.3 example
This job ran on 8 processors
There are 1030301 nodes 30301 restrained and 1000000 equations
Time after setup is : 2.37000000000261934
The total load is 0.1000E+02
The number of iterations to convergence was 121
The potentials are :
Freedom Potential
9901 0.1748E+04
9902 0.1713E+03
9903 0.1586E+03
9904 0.9910E+02
This analysis took 382.870000000002619
Figure 12.18
Results from Program 12.3 example