|
DOUG development site www.dougdevel.org |
File: doug.out
master thread
Parsing DOUG control file: /home/olegus/work/doug-git/tests/testdir/DOUG-exec.ctl
Seperate aggregated RHS file is used. RHS from elemmat_rhs_file is discarded.
Control parameters:
solver = 2
method = 1
fine_method = 2
num_iters = -1
coarse_method = -1
num_subdomains = 1
levels = 1
overlap = 1
smoothers = 0
symmstruct = T
symmnumeric = T
input_type = 2
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 = 1
minasize1 = 2
minasize2 = 2
maxasize1 = -1
maxasize2 = 128
debug = 0
verbose = 3
plotting = 0
initial_guess = 2
interpolation_type = -1
solution_format = 2
start_vec_type = 2
assembled_rhs_format = 2
assembled_rhs_file = /home/olegus/work/doug-git/tests/assembled/lap16x16/rhs.xdr
assembled_mtx_file = /home/olegus/work/doug-git/tests/assembled/lap16x16/Lap16x16.txt
assembled_mtx_format = 0
Type of Graph set to default D_GRAPH_DUAL
##### Assembled input file: ##### /home/olegus/work/doug-git/tests/assembled/lap16x16/Lap16x16.txt
n,nnz: 225 1065
1 mat: 1 1 0.4000000000000000E+01
2 mat: 116 -0.1000000000000000E+01
1064 mat:**** -0.1000000000000000E+01
1065 mat:**** 0.4000000000000000E+01
##### Assembled RHS file: ##### /home/olegus/work/doug-git/tests/assembled/lap16x16/rhs.xdr
Cartesian Mesh!!!
Mesh object [info]:
# of elements (nell) = 0
# of global freedoms (ngf) = 225
# of local freedoms (nlf) = 225
# of spacial dimensions (nsd) = 2
max # of freedoms per elem (mfrelt) = 4
total # of nodes in a mesh (nnode) = 225
Building restriction matrix
Building coarse matrix
WARNING: SpMtx_arrange to 1 not needed...
WARNING: SpMtx_arrange to 1 not needed...
WARNING: SpMtx_arrange to 1 not needed...
WARNING: SpMtx_arrange to 1 not needed...
too small aggr 3 of size: 1 maxcw_sum: 0.35714285714285715
eater of 3 is: 2
too small aggr 5 of size: 1 maxcw_sum: 0.31250000000000000
eater of 5 is: 4
too small aggr 6 of size: 1 maxcw_sum: 0.50000000000000000
eater of 6 is: 4
too small aggr 7 of size: 1 maxcw_sum: 0.50000000000000000
eater of 7 is: 4
too small aggr 8 of size: 1 maxcw_sum: 0.50000000000000000
eater of 8 is: 2
too small aggr 9 of size: 1 maxcw_sum: 0.50000000000000000
eater of 9 is: 4
# too small aggrs: 6 #eaten: 6 #occupied: 0 # remaining: 0
# coarse aggregates: 3
WARNING: SpMtx_arrange to 1 not needed...
INFO: Sym. Gauss-Seidel preconditioner with 3 iterations
INFO: Initialize 3 submatrices
INFO: Submatrix: 364 nonzeros
INFO: Submatrix: 554 nonzeros
INFO: Submatrix: 321 nonzeros
calling pcg_weigs
Preconditioned conjugate gradient:
maxit = 300, tol = 0.10E-11
init_norm = 0.150000000000000E+02
1: res_norm= 0.304307179377777E+02
2: res_norm= 0.152885967082825E+02
3: res_norm= 0.118616097644779E+02
4: res_norm= 0.279430555019003E+01
5: res_norm= 0.115655770082006E+01
6: res_norm= 0.214200120718322E+00
7: res_norm= 0.590948320834870E-01
8: res_norm= 0.219523359877573E-01
9: res_norm= 0.168173452796694E-01
10: res_norm= 0.776575492138726E-02
11: res_norm= 0.211087224301419E-02
12: res_norm= 0.881823207361490E-03
13: res_norm= 0.311574211385333E-03
14: res_norm= 0.119363166675600E-03
15: res_norm= 0.175006803922876E-04
16: res_norm= 0.442952874368849E-05
17: res_norm= 0.249008131111249E-05
18: res_norm= 0.746944274586878E-06
19: res_norm= 0.303694362034986E-06
20: res_norm= 0.109968703689423E-06
21: res_norm= 0.233969609849927E-07
22: res_norm= 0.100948074073761E-07
23: res_norm= 0.106987246532645E-08
24: res_norm= 0.250228272135277E-09
25: res_norm= 0.543706591720671E-10
26: res_norm= 0.522204367598986E-11
#it: 26 Cond#: 0.1553E+02
time spent in pcg(): 6.23083114624023438E-003
...of which setup: 0.0000000000000000
...of which factorisation: 0.0000000000000000
solve time without setup: 6.23083114624023438E-003
DOUG: < 0> doug ended