|
DOUG development site www.dougdevel.org |
File: doug.out
master thread
Parsing DOUG control file: /home/olegus/uni/projects/doug-git/tests/testdir/DOUG-exec.ctl
Seperate aggregated RHS file can not be read. RHS from elemmat_rhs_file is used instead.
Control parameters:
solver = 2
method = 1
fine_method = 1
num_iters = -1
coarse_method = -1
num_subdomains = 0
levels = 1
overlap = 2
smoothers = 0
symmstruct = T
symmnumeric = T
input_type = 1
matrix_type = 1
number_of_blocks = 1
strong1 = 6.700E-01
strong2 = 6.700E-01
solve_tolerance = 1.000E-12
solve_maxiters = 300
radius1 = 2
radius2 = 5
minasize1 = 2
minasize2 = 2
maxasize1 = -1
maxasize2 = 128
debug = 0
verbose = 10
plotting = 0
initial_guess = 2
interpolation_type = -1
solution_format = 2
start_vec_type = 2
assembled_rhs_format = 0
assembled_rhs_file = /home/olegus/uni/projects/doug-git/tests/testdir/rhs.txt
info_file = /home/olegus/uni/projects/doug-git/tests/elemental/e4x4/doug_info.dat
freedom_lists_file = /home/olegus/uni/projects/doug-git/tests/elemental/e4x4/doug_element.dat
elemmat_rhs_file = /home/olegus/uni/projects/doug-git/tests/elemental/e4x4/doug_system.dat
coords_file = /home/olegus/uni/projects/doug-git/tests/elemental/e4x4/doug_coord.dat
freemap_file = /home/olegus/uni/projects/doug-git/tests/elemental/e4x4/doug_freemap.dat
freedom_mask_file = /home/olegus/uni/projects/doug-git/tests/elemental/e4x4/./NOT.DEFINED.freedom_mask_file
solution_file = /home/olegus/uni/projects/doug-git/tests/elemental/e4x4/solution.xdr
start_vec_file = /home/olegus/uni/projects/doug-git/tests/elemental/e4x4/./NOT.DEFINED.start_vec_file
dump_matrix_file =
dump_matrix_only = F
maxcie = 4
maxnd = -1
cutbal = 3
Type of Graph set to default D_GRAPH_DUAL
master thread
Master control parameters:
solution_format = 2
start_vec_type = 2
assembled_rhs_format = 0
assembled_rhs_file = /home/olegus/uni/projects/doug-git/tests/testdir/rhs.txt
info_file = /home/olegus/uni/projects/doug-git/tests/elemental/e4x4/doug_info.dat
freedom_lists_file = /home/olegus/uni/projects/doug-git/tests/elemental/e4x4/doug_element.dat
elemmat_rhs_file = /home/olegus/uni/projects/doug-git/tests/elemental/e4x4/doug_system.dat
coords_file = /home/olegus/uni/projects/doug-git/tests/elemental/e4x4/doug_coord.dat
freemap_file = /home/olegus/uni/projects/doug-git/tests/elemental/e4x4/doug_freemap.dat
freedom_mask_file = /home/olegus/uni/projects/doug-git/tests/elemental/e4x4/./NOT.DEFINED.freedom_mask_file
solution_file = /home/olegus/uni/projects/doug-git/tests/elemental/e4x4/solution.xdr
start_vec_file = /home/olegus/uni/projects/doug-git/tests/elemental/e4x4/./NOT.DEFINED.start_vec_file
dump_matrix_file =
dump_matrix_only = F
maxcie = 4
maxnd = -1
cutbal = 3
Type of Graph set to default D_GRAPH_DUAL
Mesh object [info]:
# of elements (nell) = 16
# of global freedoms (ngf) = 9
# of local freedoms (nlf) = -1
# of spacial dimensions (nsd) = 2
max # of freedoms per elem (mfrelt) = 4
total # of nodes in a mesh (nnode) = 25
Reading in element nodes numbering ... done
Reading in freedoms' map ... done
Mesh object [element freedoms]:
| 1| 2| 3| 4| 5| 6| 7| 8| 9| 10|
1 1 2 3 1 1 2 3 4 4
0 2 3 0 4 2 3 6 7 5
0 0 0 0 0 4 5 0 0 7
0 0 0 0 0 5 6 0 0 8
| 11| 12| 13| 14| 15| 16|
5 6 7 7 8 9
6 9 0 8 9 0
8 0 0 0 0 0
9 0 0 0 0 0
Statistics of mesh assessment:
nactivef = 9
numbins = 2
hashsize = 202
hscale = 25.00000
Average connection : 4.00000
Standard deviation : 0.00000
Statistics on edge hash table:
Size of edge hash table : 99
Number of edges found : 24 ( 24% of hash table)
Graph partitioning: Number of edgecuts : 0
Global to local mapping: [M%gl_fmap] :size [ 9]:
[ 1]= 1
[ 2]= 2
[ 3]= 3
[ 4]= 4
[ 5]= 5
[ 6]= 6
[ 7]= 7
[ 8]= 8
[ 9]= 9
Local to global mapping: [M%lg_fmap] :size [ 9]:
[ 1]= 1
[ 2]= 2
[ 3]= 3
[ 4]= 4
[ 5]= 5
[ 6]= 6
[ 7]= 7
[ 8]= 8
[ 9]= 9
RANK< 0> my neighbours []
Reading in element matrices and RHSs ...
Building freedoms permutation map ... done
number of packets in use: 1
done
System matrix:
-----------------------------
Sparse matrix:
nnz = 49
nrows = 9
ncols = 9
arrange type: not arranged
shape: square
symmetric structure: Yes
symmetric numerically: Yes
number of blocks = 1
inner freedoms bound = 1
mtx_bbs =
1, 1
1, 1
mtx_bbe =
0, 0
0, 49
-----------------------------
[ 3.646, 0.477, 0.000,-1.388,-0.456, 0.000, 0.000, 0.000, 0.000]
[ 0.477, 3.646, 0.477,-0.456,-1.388,-0.456, 0.000, 0.000, 0.000]
[ 0.000, 0.477, 3.646, 0.000,-0.456,-1.388, 0.000, 0.000, 0.000]
[-1.388,-0.456, 0.000, 3.646, 0.477, 0.000,-1.388,-0.456, 0.000]
[-0.456,-1.388,-0.456, 0.477, 3.646, 0.477,-0.456,-1.388,-0.456]
[ 0.000,-0.456,-1.388, 0.000, 0.477, 3.646, 0.000,-0.456,-1.388]
[ 0.000, 0.000, 0.000,-1.388,-0.456, 0.000, 3.646, 0.477, 0.000]
[ 0.000, 0.000, 0.000,-0.456,-1.388,-0.456, 0.477, 3.646, 0.477]
[ 0.000, 0.000, 0.000, 0.000,-0.456,-1.388, 0.000, 0.477, 3.646]
RHS assembled (local) :size [ 9]:
[ 1]= 0.49288054439675E+00
[ 2]= 0.36534695559068E+00
[ 3]=-0.10420819583682E+00
[ 4]= 0.85650069120523E+00
[ 5]= 0.63487983762098E+00
[ 6]=-0.40641738450114E+00
[ 7]= 0.15611135226129E+01
[ 8]= 0.15285419120539E+01
[ 9]=-0.58160286057994E+00
dot_intf_fmask :size [ 0]:
dot_intf_fmap :size [ 0]:
Preconditioned conjugate gradient:
maxit = 300, tol = 0.10E-11
init_norm = 0.260785114693809E+01
Applying preconditioner
Factorizing 1. level
symbolic analysis:
status: 0.
time: 0.00E+00 (sec)
estimates (upper bound) for numeric LU:
size of LU: 0.14E-02 (MB)
memory needed: 0.17E-01 (MB)
flop count: 0.30E+03
nnz (L): 35.
nnz (U): 45.
numeric factorization:
status: 0.
time: 0.00E+00
actual numeric LU statistics:
size of LU: 0.12E-02 (MB)
memory needed: 0.16E-01 (MB)
flop count: 0.14E+03
nnz (L): 30.
nnz (U): 30.
Solving 1. level
1: res_norm= 0.538915796403849E-15
#it: 1 Cond#: 0.1000E+01
time spent in pcg(): 8.76903533935546875E-004
solution :size [ 9]:
[ 1]= 0.43329037790609E+00
[ 2]= 0.34041474735768E+00
[ 3]=-0.69611290397381E-01
[ 4]= 0.70845466779722E+00
[ 5]= 0.58304577204737E+00
[ 6]=-0.18226492164993E+00
[ 7]= 0.68592599293014E+00
[ 8]= 0.64877203480026E+00
[ 9]=-0.24085607030253E+00
dsqrt(res_norm) = 5.38915796403849330E-016
DOUG: < 0> doug ended