|
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 = 2
overlap = -1
smoothers = 1
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 nodes coordinates ... 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: multilevel K-way (min. communication)
Number of edgecuts : 4
Global to local mapping: [M%gl_fmap] :size [ 9]:
[ 1]= 0
[ 2]= 0
[ 3]= 0
[ 4]= 1
[ 5]= 2
[ 6]= 3
[ 7]= 4
[ 8]= 5
[ 9]= 6
Local to global mapping: [M%lg_fmap] :size [ 6]:
[ 1]= 4
[ 2]= 5
[ 3]= 6
[ 4]= 7
[ 5]= 8
[ 6]= 9
RANK< 0> my neighbours [ 1]
Reading in element matrices and RHSs ...
Building freedoms permutation map ... done
Building freedoms permutation map ... done
number of packets in use: 2
done
System matrix:
-----------------------------
Sparse matrix:
nnz = 28
nrows = 6
ncols = 6
arrange type: not arranged
shape: square
symmetric structure: Yes
symmetric numerically: Yes
number of blocks = 1
inner freedoms bound = 22
mtx_bbs =
1, 15
8, 22
mtx_bbe =
7, 21
14, 28
-----------------------------
[ 1.823, 0.238, 0.000,-1.388,-0.456, 0.000]
[ 0.238, 1.823, 0.238,-0.456,-1.388,-0.456]
[ 0.000, 0.238, 1.823, 0.000,-0.456,-1.388]
[-1.388,-0.456, 0.000, 3.646, 0.477, 0.000]
[-0.456,-1.388,-0.456, 0.477, 3.646, 0.477]
[ 0.000,-0.456,-1.388, 0.000, 0.477, 3.646]
RHS assembled (local) :size [ 6]:
[ 1]= 0.85650069120523E+00
[ 2]= 0.63487983762098E+00
[ 3]=-0.40641738450114E+00
[ 4]= 0.15611135226129E+01
[ 5]= 0.15285419120539E+01
[ 6]=-0.58160286057994E+00
dot_intf_fmask :size [ 3]:
[ 1]= 1
[ 2]= 1
[ 3]= 1
dot_intf_fmap :size [ 3]:
[ 1]= 1
[ 2]= 2
[ 3]= 3
Building coarse grid
Initial coarse grid has dimensions 3 1 with a total of 3 elements
Removing 0 empty coarse elements
Global coarse mesh has:
Grid nodes: 8
Refinements: 0
Hanging nodes: 0
Sending parts of the coarse grid to other threads
Local coarse mesh has:
Grid nodes: 8
Refinements: 0
Hanging nodes: 0
Creating Restriction matrix
Restriction matrix has 24 elements
Cleaning Restriction matrix
Building coarse matrix
Building coarse matrix
Multiplying matrices with shapes: 6x6, 8x6, transposed:F,T
WARNING: SpMtx_arrange to 1 not needed...
WARNING: SpMtx_arrange to 2 not needed...
Multiplying matrices with shapes: 8x6, 6x8, transposed:
WARNING: SpMtx_arrange to 1 not needed...
Stripping the restriction matrix
Stripped restriction matrix has 24 elements
Transmitting local-to-global maps
Preconditioned conjugate gradient:
maxit = 300, tol = 0.10E-11
init_norm = 0.260785114693809E+01
Applying preconditioner
Preparing 2. level
Exchanging coarse matrix
Restricting into local coarse vector 8
Sending local coarse vector
Factorizing 1. level
symbolic analysis:
status: 0.
time: 0.00E+00 (sec)
estimates (upper bound) for numeric LU:
size of LU: 0.92E-03 (MB)
memory needed: 0.14E-01 (MB)
flop count: 0.85E+02
nnz (L): 17.
nnz (U): 21.
numeric factorization:
status: 0.
time: 0.00E+00
actual numeric LU statistics:
size of LU: 0.85E-03 (MB)
memory needed: 0.14E-01 (MB)
flop count: 0.65E+02
nnz (L): 17.
nnz (U): 17.
Solving 1. level
Solving 2. level
Receiving local coarse vector
factorising coarse matrix of size 8 and nnz: 56
symbolic analysis:
status: 0.
time: 0.00E+00 (sec)
estimates (upper bound) for numeric LU:
size of LU: 0.12E-02 (MB)
memory needed: 0.16E-01 (MB)
flop count: 0.25E+03
nnz (L): 32.
nnz (U): 36.
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.22E+03
nnz (L): 32.
nnz (U): 32.
1: res_norm= 0.207434776144955E+01
Applying preconditioner
Preparing 2. level
Restricting into local coarse vector 8
Sending local coarse vector
Solving 1. level
Solving 2. level
Receiving local coarse vector
2: res_norm= 0.170902945995823E+00
Applying preconditioner
Preparing 2. level
Restricting into local coarse vector 8
Sending local coarse vector
Solving 1. level
Solving 2. level
Receiving local coarse vector
3: res_norm= 0.954417210579805E-01
Applying preconditioner
Preparing 2. level
Restricting into local coarse vector 8
Sending local coarse vector
Solving 1. level
Solving 2. level
Receiving local coarse vector
4: res_norm= 0.212809244691061E-01
Applying preconditioner
Preparing 2. level
Restricting into local coarse vector 8
Sending local coarse vector
Solving 1. level
Solving 2. level
Receiving local coarse vector
5: res_norm= 0.547932431120933E-02
Applying preconditioner
Preparing 2. level
Restricting into local coarse vector 8
Sending local coarse vector
Solving 1. level
Solving 2. level
Receiving local coarse vector
6: res_norm= 0.539131608268427E-03
Applying preconditioner
Preparing 2. level
Restricting into local coarse vector 8
Sending local coarse vector
Solving 1. level
Solving 2. level
Receiving local coarse vector
7: res_norm= 0.333054705262478E-04
Applying preconditioner
Preparing 2. level
Restricting into local coarse vector 8
Sending local coarse vector
Solving 1. level
Solving 2. level
Receiving local coarse vector
8: res_norm= 0.152080848239328E-05
Applying preconditioner
Preparing 2. level
Restricting into local coarse vector 8
Sending local coarse vector
Solving 1. level
Solving 2. level
Receiving local coarse vector
9: res_norm= 0.164119672049414E-13
#it: 9 Cond#: 0.5248E+01
time spent in pcg(): 1.59406661987304688E-003
solution :size [ 9]:
[ 1]= 0.43329037790608E+00
[ 2]= 0.34041474735767E+00
[ 3]=-0.69611290397377E-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) = 2.31643185291875359E-014
DOUG: < 0> doug ended