| |||
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 = 0 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: 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 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.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 1: res_norm= 0.164855640588882E+01 Applying preconditioner Solving 1. level 2: res_norm= 0.260700127667448E+00 Applying preconditioner Solving 1. level 3: res_norm= 0.473176977323666E-01 Applying preconditioner Solving 1. level 4: res_norm= 0.664602591981273E-02 Applying preconditioner Solving 1. level 5: res_norm= 0.595796303097263E-03 Applying preconditioner Solving 1. level 6: res_norm= 0.414437390960596E-04 Applying preconditioner Solving 1. level 7: res_norm= 0.114644576216608E-15 #it: 7 Cond#: 0.2795E+01 time spent in pcg(): 1.17301940917968750E-003 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) = 1.22690946876004807E-015 DOUG: < 0> doug ended |
University of Tartu 2007-2008, info(REMOVEME)@dougdevel.org