| |||
File: doug.out master thread Parsing DOUG control file: /home/olegus/work/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 grid_size = -1 coarse_method = -1 num_subdomains = 0 levels = 1 overlap = -1 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/work/doug-git/tests/testdir/rhs.txt info_file = /home/olegus/work/doug-git/tests/elemental/e4x4/doug_info.dat freedom_lists_file = /home/olegus/work/doug-git/tests/elemental/e4x4/doug_element.dat elemmat_rhs_file = /home/olegus/work/doug-git/tests/elemental/e4x4/doug_system.dat coords_file = /home/olegus/work/doug-git/tests/elemental/e4x4/doug_coord.dat freemap_file = /home/olegus/work/doug-git/tests/elemental/e4x4/doug_freemap.dat freedom_mask_file = /home/olegus/work/doug-git/tests/elemental/e4x4/./NOT.DEFINED.freedom_mask_file solution_file = /home/olegus/work/doug-git/tests/elemental/e4x4/solution.xdr start_vec_file = /home/olegus/work/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/work/doug-git/tests/testdir/rhs.txt info_file = /home/olegus/work/doug-git/tests/elemental/e4x4/doug_info.dat freedom_lists_file = /home/olegus/work/doug-git/tests/elemental/e4x4/doug_element.dat elemmat_rhs_file = /home/olegus/work/doug-git/tests/elemental/e4x4/doug_system.dat coords_file = /home/olegus/work/doug-git/tests/elemental/e4x4/doug_coord.dat freemap_file = /home/olegus/work/doug-git/tests/elemental/e4x4/doug_freemap.dat freedom_mask_file = /home/olegus/work/doug-git/tests/elemental/e4x4/./NOT.DEFINED.freedom_mask_file solution_file = /home/olegus/work/doug-git/tests/elemental/e4x4/solution.xdr start_vec_file = /home/olegus/work/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 : 32663 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(): 1.15768909454345703E-002 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 |
University of Tartu 2007-2008, info(REMOVEME)@dougdevel.org