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 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 =     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/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:  multilevel K-way (min. communication)
 Number of edgecuts :     8

Global to local mapping: [M%gl_fmap] :size [     9]:
 [     1]=               4
 [     2]=               1
 [     3]=               0
 [     4]=               2
 [     5]=               3
 [     6]=               0
 [     7]=               0
 [     8]=               0
 [     9]=               0

Local to global mapping: [M%lg_fmap] :size [     4]:
 [     1]=               2
 [     2]=               4
 [     3]=               5
 [     4]=               1

RANK<  0> my neighbours [  1,   2,   3]

Reading in element matrices and RHSs ... 
Building freedoms permutation map ...  done

Building freedoms permutation map ...  done
 number of packets in use:            4
 done

System matrix:
 -----------------------------
 Sparse matrix:
   nnz =           16
 nrows =            4
 ncols =            4
 arrange type:  not arranged
 shape:  square
 symmetric structure:   Yes
 symmetric numerically: Yes
 number of blocks     =            1
 inner freedoms bound =           16
 mtx_bbs = 
       1,       13
      10,       16
 mtx_bbe = 
       9,       15
      12,       16
 -----------------------------
[ 1.823,-0.456,-0.694, 0.477]
[-0.456, 1.823, 0.238,-1.388]
[-0.694, 0.238, 0.912,-0.456]
[ 0.477,-1.388,-0.456, 3.646]

RHS assembled (local)  :size [     4]:
 [     1]= 0.36534695559068E+00
 [     2]= 0.85650069120523E+00
 [     3]= 0.63487983762098E+00
 [     4]= 0.49288054439675E+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.69E-03 (MB)
   memory needed:   0.12E-01 (MB)
   flop count:      0.34E+02
   nnz (L):              10.
   nnz (U):              10.
numeric factorization:
   status:     0.
   time:      0.00E+00
   actual numeric LU statistics:
   size of LU:      0.68E-03 (MB)
   memory needed:   0.12E-01 (MB)
   flop count:      0.34E+02
   nnz (L):              10.
   nnz (U):              10.
 Solving 1. level
    1: res_norm= 0.172622962226937E+01
 Applying preconditioner
 Solving 1. level
    2: res_norm= 0.856022663717325E+00
 Applying preconditioner
 Solving 1. level
    3: res_norm= 0.227482145797073E+00
 Applying preconditioner
 Solving 1. level
    4: res_norm= 0.437132065572812E-01
 Applying preconditioner
 Solving 1. level
    5: res_norm= 0.102801226712442E-01
 Applying preconditioner
 Solving 1. level
    6: res_norm= 0.241456376167215E-02
 Applying preconditioner
 Solving 1. level
    7: res_norm= 0.651077421166873E-04
 Applying preconditioner
 Solving 1. level
    8: res_norm= 0.917745683664503E-05
 Applying preconditioner
 Solving 1. level
    9: res_norm= 0.653604421191989E-14
#it:  9 Cond#: 0.5017E+01
 time spent in pcg():  2.92739868164062500E-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) =  1.41661151596772975E-014

DOUG: <  0> doug ended

University of Tartu 2007-2008, info(REMOVEME)@dougdevel.org