Opening HEAT problem file ./problem_heat.dat
Warning: time_step_length_nonl_control not found in problem file. Assuming defaults.
HEAT problem "problem_nieliniowy" settings :
CONTROL PARAMETERS:
mesh_type: j
mesh_file_in: A.jk
mesh_file_out: mesh
field_file_in: i
field_file_out: field_heat
field_function:
penalty for Dirichlet BCs: 10000000.000000
reference_temperature: 300.000000
ambient temperature: 300.000000
materials configuration file: materials.dat
HEAT material database or file read
material data at reference temperature
material number: 0 (data at reference temperature)
reference temperature: 300.000000
thermal_conductivity: 3.500000
density: 7500.000000
specific_heat: 650.000000
material number: 1 (data at reference temperature)
reference temperature: 300.000000
thermal_conductivity: 34.728838
density: 7847.175017
specific_heat: 574.520386
TIME INTEGRATION PARAMETERS:
type: 1
implicitness parameter: 1.000000
current timestep: 0
current time: 0.000000
current timestep_length: 1000000000.000000
previous timestep_length: 1000000000.000000
final time: 50000.000000
final timestep: 2
convergence criterion type: 0
error tolerance (n-epsilon): 0.000001
monitoring level: 0
graph_dump_intv: 0
full_dump_intv: 0
NONLINEAR SOLVER PARAMETERS:
type: 0
max_nonl_iter: 100
convergence criterion type: 0
error tolerance (k-epsilon): 0.000001
monitoring level: 3
LINEAR SOLVER PARAMETERS:
solver type: 0
monitoring level: 0
no special file with solver parameters
ADAPTATION PARAMETERS:
adapt_type: 2
adapt_interval: 0
adapt_eps: 0.110000
adapt_ratio: 0.200000
monitoring level: 2
BOUNDARY CONDITIONS
boundary conditions configuration file: bc_heat.dat
number of BCs: 12
BC ID: 21
BC type: 4 - BC_HEAT_NORMAL_HEAT_FLUX
flux: 0.000000
BC ID: 22
BC type: 4 - BC_HEAT_NORMAL_HEAT_FLUX
flux: 0.000000
BC ID: 2
BC type: 3 - BC_HEAT_ISOTHERMAL
temperature: 300.000000
BC ID: 3
BC type: 4 - BC_HEAT_NORMAL_HEAT_FLUX
flux: 0.000000
BC ID: 4
BC type: 4 - BC_HEAT_NORMAL_HEAT_FLUX
flux: 0.000000
BC ID: 5
BC type: 3 - BC_HEAT_ISOTHERMAL
temperature: 400.000000
BC ID: 6
BC type: 4 - BC_HEAT_NORMAL_HEAT_FLUX
flux: 0.000000
BC ID: 7
BC type: 4 - BC_HEAT_NORMAL_HEAT_FLUX
flux: 0.000000
BC ID: 8
BC type: 3 - BC_HEAT_ISOTHERMAL
temperature: 300.000000
BC ID: 9
BC type: 4 - BC_HEAT_NORMAL_HEAT_FLUX
flux: 0.000000
BC ID: 10
BC type: 4 - BC_HEAT_NORMAL_HEAT_FLUX
flux: 0.000000
BC ID: 1
BC type: 4 - BC_HEAT_NORMAL_HEAT_FLUX
flux: 0.000000
MESH
Read mesh file A.jk
Material data specified for all elements
All elements assigned to single material data
Boundary conditions specified for all boundary faces
After reading initial mesh data.
Mesh entities (number of active, maximal index):
Elements: nrel 42, nmel 42
Faces: nrfa 158, nmfa 158
Edges: nred 181, nmed 181
Nodes: nrno 66, nmno 66
Maximal generation level set to 9, maximal generation difference set to 1
Initializing heat field with initial_condition function
Approximation type: LINEAR
NR DOF ENTS/VERTs: 66 ; nreq = 1 ; nr_sol = 3
After create CONSTRAINTS DATA:
NR DOF ENTS/VERTs: 66 ; nreq = 1 ; nr_sol = 3
Choose a command from the menu:
t - solve the problem (time integration)
s - solve single heat problem
e - compute error for test problems with known exact solution
z - compute error (ZZ aproximation)
p - postprocessing
g - launch graphics module
m - perform uniform mesh refinement
a - perform automatic mesh adaptation
r - print profile
d - dump out data
v - write ParaView graphics data
c - change control data
h - perform random test of h-adaptivity
q - exit the program
t
Beginning time integration of heat problem
**************** BEGINNING TIME INTEGRATION ***************
Time integration will stop whichever comes first:
final_time: 50000.000000
final_timestep: 2
error less than time_integration_tolerance: 0.000001
Type [Ctrl-C] to manually break time integration.
Entering linear solver module for solver 0
lsr_mkb_init input parameters: Solver_type 0, Parallel 0, Monitor 0
No special file with solver parameters
Assigned solver IDs for heat 0
***************************************************************************
***************************************************************************
Solving time step 1 (step_length: 1000000000.000000, initial time: 0.000000)
---------------------------------------------------------------------------
Solving nonlinear iteration 0
Entering RCM renumbering !!!!
Entering lsr_mkb_create_matrix: solver 0, level 0, storage_type -1
nrblocks 66, max_sm_size 6, nrdof_glob 66, block_size 1
Entering lar_allocate_SM_and_LV: storage_type: requested -1, final 2
nrblocks 66, max_sm_size 6, nrdof_glob 66, block_size 1
Allocated CRS matrix 0: nglob = 66, nloc=66 nnz = 724, half bandwidth = 20
Entering sir_solve (with initial guess): Solver 0, Monitor 0, nrdofs_glob 66
Entering lsr_mkb_clear_matrix: solver 0, level 0, comp_type 0
Entering lar_initialize_SM_and_LV: storage_type: 2, scope 0
matrix id: in interface 0, internally in implementation 0
Initialized CRS matrix 0: n = 66, nnz = 724, half bandwidth = 20
The average number of non-zeros in a single row: 10
Entering lsr_mkb_solve: Solver_type 0, Solver id 0, Comp_type 0
Solution in struct 0, block 35, dof_ent 6, posglob=35
399.9999999510
Solution in struct 6, block 33, dof_ent 1, posglob=33
350.8440184903
Solution in struct 38, block 11, dof_ent 4, posglob=11
332.3578945818
Solution in struct 40, block 57, dof_ent 7, posglob=57
332.3516103050
Solution in struct 42, block 19, dof_ent 5, posglob=19
399.9999999510
Solution in struct 48, block 15, dof_ent 2, posglob=15
300.0000000191
Solution in struct 52, block 51, dof_ent 9, posglob=51
300.0000000191
Solution in struct 62, block 1, dof_ent 3, posglob=1
300.0000000169
Solution in struct 64, block 64, dof_ent 8, posglob=64
300.0000000169
After linear solver in nonlinear iteration 0
Solution difference norm (u_k, prev. u_k): 99.999999977842 (limit 0.000001000000)
---------------------------------------------------------------------------
Solving nonlinear iteration 1
Entering sir_solve (with initial guess): Solver 0, Monitor 0, nrdofs_glob 66
Entering lsr_mkb_clear_matrix: solver 0, level 0, comp_type 0
Entering lar_initialize_SM_and_LV: storage_type: 2, scope 0
matrix id: in interface 0, internally in implementation 0
Initialized CRS matrix 0: n = 66, nnz = 724, half bandwidth = 20
The average number of non-zeros in a single row: 10
Entering lsr_mkb_solve: Solver_type 0, Solver id 0, Comp_type 0
Solution in struct 0, block 35, dof_ent 6, posglob=35
399.9999998659
Solution in struct 6, block 33, dof_ent 1, posglob=33
384.1005805230
Solution in struct 38, block 11, dof_ent 4, posglob=11
370.9337553779
Solution in struct 40, block 57, dof_ent 7, posglob=57
370.9295046777
Solution in struct 42, block 19, dof_ent 5, posglob=19
399.9999998659
Solution in struct 48, block 15, dof_ent 2, posglob=15
300.0000000513
Solution in struct 52, block 51, dof_ent 9, posglob=51
300.0000000513
Solution in struct 62, block 1, dof_ent 3, posglob=1
300.0000000470
Solution in struct 64, block 64, dof_ent 8, posglob=64
300.0000000470
After linear solver in nonlinear iteration 1
Solution difference norm (u_k, prev. u_k): 39.761831371248 (limit 0.000001000000)
---------------------------------------------------------------------------
Solving nonlinear iteration 2
Entering sir_solve (with initial guess): Solver 0, Monitor 0, nrdofs_glob 66
Entering lsr_mkb_clear_matrix: solver 0, level 0, comp_type 0
Entering lar_initialize_SM_and_LV: storage_type: 2, scope 0
matrix id: in interface 0, internally in implementation 0
Initialized CRS matrix 0: n = 66, nnz = 724, half bandwidth = 20
The average number of non-zeros in a single row: 10
Entering lsr_mkb_solve: Solver_type 0, Solver id 0, Comp_type 0
Solution in struct 0, block 35, dof_ent 6, posglob=35
399.9999997122
Solution in struct 6, block 33, dof_ent 1, posglob=33
371.1526939818
Solution in struct 38, block 11, dof_ent 4, posglob=11
359.6645755124
Solution in struct 40, block 57, dof_ent 7, posglob=57
359.6623772464
Solution in struct 42, block 19, dof_ent 5, posglob=19
399.9999997121
Solution in struct 48, block 15, dof_ent 2, posglob=15
300.0000001397
Solution in struct 52, block 51, dof_ent 9, posglob=51
300.0000001397
Solution in struct 62, block 1, dof_ent 3, posglob=1
300.0000001046
Solution in struct 64, block 64, dof_ent 8, posglob=64
300.0000001046
After linear solver in nonlinear iteration 2
Solution difference norm (u_k, prev. u_k): 12.947976052720 (limit 0.000001000000)
---------------------------------------------------------------------------
Solving nonlinear iteration 3
Entering sir_solve (with initial guess): Solver 0, Monitor 0, nrdofs_glob 66
Entering lsr_mkb_clear_matrix: solver 0, level 0, comp_type 0
Entering lar_initialize_SM_and_LV: storage_type: 2, scope 0
matrix id: in interface 0, internally in implementation 0
Initialized CRS matrix 0: n = 66, nnz = 724, half bandwidth = 20
The average number of non-zeros in a single row: 10
Entering lsr_mkb_solve: Solver_type 0, Solver id 0, Comp_type 0
Solution in struct 0, block 35, dof_ent 6, posglob=35
399.9999997220
Solution in struct 6, block 33, dof_ent 1, posglob=33
371.6179253984
Solution in struct 38, block 11, dof_ent 4, posglob=11
357.7290961327
Solution in struct 40, block 57, dof_ent 7, posglob=57
357.7238299576
Solution in struct 42, block 19, dof_ent 5, posglob=19
399.9999997219
Solution in struct 48, block 15, dof_ent 2, posglob=15
300.0000001222
Solution in struct 52, block 51, dof_ent 9, posglob=51
300.0000001222
Solution in struct 62, block 1, dof_ent 3, posglob=1
300.0000001043
Solution in struct 64, block 64, dof_ent 8, posglob=64
300.0000001043
After linear solver in nonlinear iteration 3
Solution difference norm (u_k, prev. u_k): 3.786491744812 (limit 0.000001000000)
---------------------------------------------------------------------------
Solving nonlinear iteration 4
Entering sir_solve (with initial guess): Solver 0, Monitor 0, nrdofs_glob 66
Entering lsr_mkb_clear_matrix: solver 0, level 0, comp_type 0
Entering lar_initialize_SM_and_LV: storage_type: 2, scope 0
matrix id: in interface 0, internally in implementation 0
Initialized CRS matrix 0: n = 66, nnz = 724, half bandwidth = 20
The average number of non-zeros in a single row: 10
Entering lsr_mkb_solve: Solver_type 0, Solver id 0, Comp_type 0
Solution in struct 0, block 35, dof_ent 6, posglob=35
399.9999997369
Solution in struct 6, block 33, dof_ent 1, posglob=33
373.1974678572
Solution in struct 38, block 11, dof_ent 4, posglob=11
359.8569255609
Solution in struct 40, block 57, dof_ent 7, posglob=57
359.8532924599
Solution in struct 42, block 19, dof_ent 5, posglob=19
399.9999997368
Solution in struct 48, block 15, dof_ent 2, posglob=15
300.0000001175
Solution in struct 52, block 51, dof_ent 9, posglob=51
300.0000001174
Solution in struct 62, block 1, dof_ent 3, posglob=1
300.0000000972
Solution in struct 64, block 64, dof_ent 8, posglob=64
300.0000000972
After linear solver in nonlinear iteration 4
Solution difference norm (u_k, prev. u_k): 2.180315184417 (limit 0.000001000000)
---------------------------------------------------------------------------
Solving nonlinear iteration 5
Entering sir_solve (with initial guess): Solver 0, Monitor 0, nrdofs_glob 66
Entering lsr_mkb_clear_matrix: solver 0, level 0, comp_type 0
Entering lar_initialize_SM_and_LV: storage_type: 2, scope 0
matrix id: in interface 0, internally in implementation 0
Initialized CRS matrix 0: n = 66, nnz = 724, half bandwidth = 20
The average number of non-zeros in a single row: 10
Entering lsr_mkb_solve: Solver_type 0, Solver id 0, Comp_type 0
Solution in struct 0, block 35, dof_ent 6, posglob=35
399.9999997265
Solution in struct 6, block 33, dof_ent 1, posglob=33
372.2969741692
Solution in struct 38, block 11, dof_ent 4, posglob=11
358.9286370807
Solution in struct 40, block 57, dof_ent 7, posglob=57
358.9233482084
Solution in struct 42, block 19, dof_ent 5, posglob=19
399.9999997264
Solution in struct 48, block 15, dof_ent 2, posglob=15
300.0000001220
Solution in struct 52, block 51, dof_ent 9, posglob=51
300.0000001219
Solution in struct 62, block 1, dof_ent 3, posglob=1
300.0000001015
Solution in struct 64, block 64, dof_ent 8, posglob=64
300.0000001015
After linear solver in nonlinear iteration 5
Solution difference norm (u_k, prev. u_k): 0.934993227910 (limit 0.000001000000)
---------------------------------------------------------------------------
Solving nonlinear iteration 6
Entering sir_solve (with initial guess): Solver 0, Monitor 0, nrdofs_glob 66
Entering lsr_mkb_clear_matrix: solver 0, level 0, comp_type 0
Entering lar_initialize_SM_and_LV: storage_type: 2, scope 0
matrix id: in interface 0, internally in implementation 0
Initialized CRS matrix 0: n = 66, nnz = 724, half bandwidth = 20
The average number of non-zeros in a single row: 10
Entering lsr_mkb_solve: Solver_type 0, Solver id 0, Comp_type 0
Solution in struct 0, block 35, dof_ent 6, posglob=35
399.9999997306
Solution in struct 6, block 33, dof_ent 1, posglob=33
372.6251893738
Solution in struct 38, block 11, dof_ent 4, posglob=11
359.2205244219
Solution in struct 40, block 57, dof_ent 7, posglob=57
359.2159137218
Solution in struct 42, block 19, dof_ent 5, posglob=19
399.9999997306
Solution in struct 48, block 15, dof_ent 2, posglob=15
300.0000001201
Solution in struct 52, block 51, dof_ent 9, posglob=51
300.0000001200
Solution in struct 62, block 1, dof_ent 3, posglob=1
300.0000000999
Solution in struct 64, block 64, dof_ent 8, posglob=64
300.0000000999
After linear solver in nonlinear iteration 6
Solution difference norm (u_k, prev. u_k): 0.328853724129 (limit 0.000001000000)
---------------------------------------------------------------------------
Solving nonlinear iteration 7
Entering sir_solve (with initial guess): Solver 0, Monitor 0, nrdofs_glob 66
Entering lsr_mkb_clear_matrix: solver 0, level 0, comp_type 0
Entering lar_initialize_SM_and_LV: storage_type: 2, scope 0
matrix id: in interface 0, internally in implementation 0
Initialized CRS matrix 0: n = 66, nnz = 724, half bandwidth = 20
The average number of non-zeros in a single row: 10
Entering lsr_mkb_solve: Solver_type 0, Solver id 0, Comp_type 0
Solution in struct 0, block 35, dof_ent 6, posglob=35
399.9999997294
Solution in struct 6, block 33, dof_ent 1, posglob=33
372.5334194324
Solution in struct 38, block 11, dof_ent 4, posglob=11
359.1486514265
Solution in struct 40, block 57, dof_ent 7, posglob=57
359.1440206951
Solution in struct 42, block 19, dof_ent 5, posglob=19
399.9999997293
Solution in struct 48, block 15, dof_ent 2, posglob=15
300.0000001207
Solution in struct 52, block 51, dof_ent 9, posglob=51
300.0000001206
Solution in struct 62, block 1, dof_ent 3, posglob=1
300.0000001004
Solution in struct 64, block 64, dof_ent 8, posglob=64
300.0000001003
After linear solver in nonlinear iteration 7
Solution difference norm (u_k, prev. u_k): 0.091983489216 (limit 0.000001000000)
---------------------------------------------------------------------------
Solving nonlinear iteration 8
Entering sir_solve (with initial guess): Solver 0, Monitor 0, nrdofs_glob 66
Entering lsr_mkb_clear_matrix: solver 0, level 0, comp_type 0
Entering lar_initialize_SM_and_LV: storage_type: 2, scope 0
matrix id: in interface 0, internally in implementation 0
Initialized CRS matrix 0: n = 66, nnz = 724, half bandwidth = 20
The average number of non-zeros in a single row: 10
Entering lsr_mkb_solve: Solver_type 0, Solver id 0, Comp_type 0
Solution in struct 0, block 35, dof_ent 6, posglob=35
399.9999997297
Solution in struct 6, block 33, dof_ent 1, posglob=33
372.5530005597
Solution in struct 38, block 11, dof_ent 4, posglob=11
359.1612305597
Solution in struct 40, block 57, dof_ent 7, posglob=57
359.1565558073
Solution in struct 42, block 19, dof_ent 5, posglob=19
399.9999997296
Solution in struct 48, block 15, dof_ent 2, posglob=15
300.0000001205
Solution in struct 52, block 51, dof_ent 9, posglob=51
300.0000001205
Solution in struct 62, block 1, dof_ent 3, posglob=1
300.0000001003
Solution in struct 64, block 64, dof_ent 8, posglob=64
300.0000001002
After linear solver in nonlinear iteration 8
Solution difference norm (u_k, prev. u_k): 0.020078299870 (limit 0.000001000000)
---------------------------------------------------------------------------
Solving nonlinear iteration 9
Entering sir_solve (with initial guess): Solver 0, Monitor 0, nrdofs_glob 66
Entering lsr_mkb_clear_matrix: solver 0, level 0, comp_type 0
Entering lar_initialize_SM_and_LV: storage_type: 2, scope 0
matrix id: in interface 0, internally in implementation 0
Initialized CRS matrix 0: n = 66, nnz = 724, half bandwidth = 20
The average number of non-zeros in a single row: 10
Entering lsr_mkb_solve: Solver_type 0, Solver id 0, Comp_type 0
Solution in struct 0, block 35, dof_ent 6, posglob=35
399.9999997296
Solution in struct 6, block 33, dof_ent 1, posglob=33
372.5505201969
Solution in struct 38, block 11, dof_ent 4, posglob=11
359.1606487152
Solution in struct 40, block 57, dof_ent 7, posglob=57
359.1559941897
Solution in struct 42, block 19, dof_ent 5, posglob=19
399.9999997296
Solution in struct 48, block 15, dof_ent 2, posglob=15
300.0000001205
Solution in struct 52, block 51, dof_ent 9, posglob=51
300.0000001205
Solution in struct 62, block 1, dof_ent 3, posglob=1
300.0000001003
Solution in struct 64, block 64, dof_ent 8, posglob=64
300.0000001002
After linear solver in nonlinear iteration 9
Solution difference norm (u_k, prev. u_k): 0.003052243251 (limit 0.000001000000)
---------------------------------------------------------------------------
Solving nonlinear iteration 10
Entering sir_solve (with initial guess): Solver 0, Monitor 0, nrdofs_glob 66
Entering lsr_mkb_clear_matrix: solver 0, level 0, comp_type 0
Entering lar_initialize_SM_and_LV: storage_type: 2, scope 0
matrix id: in interface 0, internally in implementation 0
Initialized CRS matrix 0: n = 66, nnz = 724, half bandwidth = 20
The average number of non-zeros in a single row: 10
Entering lsr_mkb_solve: Solver_type 0, Solver id 0, Comp_type 0
Solution in struct 0, block 35, dof_ent 6, posglob=35
399.9999997296
Solution in struct 6, block 33, dof_ent 1, posglob=33
372.5502113067
Solution in struct 38, block 11, dof_ent 4, posglob=11
359.1599514492
Solution in struct 40, block 57, dof_ent 7, posglob=57
359.1552914200
Solution in struct 42, block 19, dof_ent 5, posglob=19
399.9999997296
Solution in struct 48, block 15, dof_ent 2, posglob=15
300.0000001205
Solution in struct 52, block 51, dof_ent 9, posglob=51
300.0000001205
Solution in struct 62, block 1, dof_ent 3, posglob=1
300.0000001003
Solution in struct 64, block 64, dof_ent 8, posglob=64
300.0000001002
After linear solver in nonlinear iteration 10
Solution difference norm (u_k, prev. u_k): 0.000777424791 (limit 0.000001000000)
---------------------------------------------------------------------------
Solving nonlinear iteration 11
Entering sir_solve (with initial guess): Solver 0, Monitor 0, nrdofs_glob 66
Entering lsr_mkb_clear_matrix: solver 0, level 0, comp_type 0
Entering lar_initialize_SM_and_LV: storage_type: 2, scope 0
matrix id: in interface 0, internally in implementation 0
Initialized CRS matrix 0: n = 66, nnz = 724, half bandwidth = 20
The average number of non-zeros in a single row: 10
Entering lsr_mkb_solve: Solver_type 0, Solver id 0, Comp_type 0
Solution in struct 0, block 35, dof_ent 6, posglob=35
399.9999997296
Solution in struct 6, block 33, dof_ent 1, posglob=33
372.5505575830
Solution in struct 38, block 11, dof_ent 4, posglob=11
359.1603407188
Solution in struct 40, block 57, dof_ent 7, posglob=57
359.1556815643
Solution in struct 42, block 19, dof_ent 5, posglob=19
399.9999997296
Solution in struct 48, block 15, dof_ent 2, posglob=15
300.0000001205
Solution in struct 52, block 51, dof_ent 9, posglob=51
300.0000001205
Solution in struct 62, block 1, dof_ent 3, posglob=1
300.0000001003
Solution in struct 64, block 64, dof_ent 8, posglob=64
300.0000001002
After linear solver in nonlinear iteration 11
Solution difference norm (u_k, prev. u_k): 0.000396094464 (limit 0.000001000000)
---------------------------------------------------------------------------
Solving nonlinear iteration 12
Entering sir_solve (with initial guess): Solver 0, Monitor 0, nrdofs_glob 66
Entering lsr_mkb_clear_matrix: solver 0, level 0, comp_type 0
Entering lar_initialize_SM_and_LV: storage_type: 2, scope 0
matrix id: in interface 0, internally in implementation 0
Initialized CRS matrix 0: n = 66, nnz = 724, half bandwidth = 20
The average number of non-zeros in a single row: 10
Entering lsr_mkb_solve: Solver_type 0, Solver id 0, Comp_type 0
Solution in struct 0, block 35, dof_ent 6, posglob=35
399.9999997296
Solution in struct 6, block 33, dof_ent 1, posglob=33
372.5504110459
Solution in struct 38, block 11, dof_ent 4, posglob=11
359.1602033617
Solution in struct 40, block 57, dof_ent 7, posglob=57
359.1555442528
Solution in struct 42, block 19, dof_ent 5, posglob=19
399.9999997296
Solution in struct 48, block 15, dof_ent 2, posglob=15
300.0000001205
Solution in struct 52, block 51, dof_ent 9, posglob=51
300.0000001205
Solution in struct 62, block 1, dof_ent 3, posglob=1
300.0000001003
Solution in struct 64, block 64, dof_ent 8, posglob=64
300.0000001002
After linear solver in nonlinear iteration 12
Solution difference norm (u_k, prev. u_k): 0.000146537146 (limit 0.000001000000)
---------------------------------------------------------------------------
Solving nonlinear iteration 13
Entering sir_solve (with initial guess): Solver 0, Monitor 0, nrdofs_glob 66
Entering lsr_mkb_clear_matrix: solver 0, level 0, comp_type 0
Entering lar_initialize_SM_and_LV: storage_type: 2, scope 0
matrix id: in interface 0, internally in implementation 0
Initialized CRS matrix 0: n = 66, nnz = 724, half bandwidth = 20
The average number of non-zeros in a single row: 10
Entering lsr_mkb_solve: Solver_type 0, Solver id 0, Comp_type 0
Solution in struct 0, block 35, dof_ent 6, posglob=35
399.9999997296
Solution in struct 6, block 33, dof_ent 1, posglob=33
372.5504561951
Solution in struct 38, block 11, dof_ent 4, posglob=11
359.1602406229
Solution in struct 40, block 57, dof_ent 7, posglob=57
359.1555814195
Solution in struct 42, block 19, dof_ent 5, posglob=19
399.9999997296
Solution in struct 48, block 15, dof_ent 2, posglob=15
300.0000001205
Solution in struct 52, block 51, dof_ent 9, posglob=51
300.0000001205
Solution in struct 62, block 1, dof_ent 3, posglob=1
300.0000001003
Solution in struct 64, block 64, dof_ent 8, posglob=64
300.0000001002
After linear solver in nonlinear iteration 13
Solution difference norm (u_k, prev. u_k): 0.000045149188 (limit 0.000001000000)
---------------------------------------------------------------------------
Solving nonlinear iteration 14
Entering sir_solve (with initial guess): Solver 0, Monitor 0, nrdofs_glob 66
Entering lsr_mkb_clear_matrix: solver 0, level 0, comp_type 0
Entering lar_initialize_SM_and_LV: storage_type: 2, scope 0
matrix id: in interface 0, internally in implementation 0
Initialized CRS matrix 0: n = 66, nnz = 724, half bandwidth = 20
The average number of non-zeros in a single row: 10
Entering lsr_mkb_solve: Solver_type 0, Solver id 0, Comp_type 0
Solution in struct 0, block 35, dof_ent 6, posglob=35
399.9999997296
Solution in struct 6, block 33, dof_ent 1, posglob=33
372.5504454015
Solution in struct 38, block 11, dof_ent 4, posglob=11
359.1602329998
Solution in struct 40, block 57, dof_ent 7, posglob=57
359.1555738379
Solution in struct 42, block 19, dof_ent 5, posglob=19
399.9999997296
Solution in struct 48, block 15, dof_ent 2, posglob=15
300.0000001205
Solution in struct 52, block 51, dof_ent 9, posglob=51
300.0000001205
Solution in struct 62, block 1, dof_ent 3, posglob=1
300.0000001003
Solution in struct 64, block 64, dof_ent 8, posglob=64
300.0000001002
After linear solver in nonlinear iteration 14
Solution difference norm (u_k, prev. u_k): 0.000010949884 (limit 0.000001000000)
---------------------------------------------------------------------------
Solving nonlinear iteration 15
Entering sir_solve (with initial guess): Solver 0, Monitor 0, nrdofs_glob 66
Entering lsr_mkb_clear_matrix: solver 0, level 0, comp_type 0
Entering lar_initialize_SM_and_LV: storage_type: 2, scope 0
matrix id: in interface 0, internally in implementation 0
Initialized CRS matrix 0: n = 66, nnz = 724, half bandwidth = 20
The average number of non-zeros in a single row: 10
Entering lsr_mkb_solve: Solver_type 0, Solver id 0, Comp_type 0
Solution in struct 0, block 35, dof_ent 6, posglob=35
399.9999997296
Solution in struct 6, block 33, dof_ent 1, posglob=33
372.5504471904
Solution in struct 38, block 11, dof_ent 4, posglob=11
359.1602338426
Solution in struct 40, block 57, dof_ent 7, posglob=57
359.1555746689
Solution in struct 42, block 19, dof_ent 5, posglob=19
399.9999997296
Solution in struct 48, block 15, dof_ent 2, posglob=15
300.0000001205
Solution in struct 52, block 51, dof_ent 9, posglob=51
300.0000001205
Solution in struct 62, block 1, dof_ent 3, posglob=1
300.0000001003
Solution in struct 64, block 64, dof_ent 8, posglob=64
300.0000001002
After linear solver in nonlinear iteration 15
Solution difference norm (u_k, prev. u_k): 0.000001976939 (limit 0.000001000000)
---------------------------------------------------------------------------
Solving nonlinear iteration 16
Entering sir_solve (with initial guess): Solver 0, Monitor 0, nrdofs_glob 66
Entering lsr_mkb_clear_matrix: solver 0, level 0, comp_type 0
Entering lar_initialize_SM_and_LV: storage_type: 2, scope 0
matrix id: in interface 0, internally in implementation 0
Initialized CRS matrix 0: n = 66, nnz = 724, half bandwidth = 20
The average number of non-zeros in a single row: 10
Entering lsr_mkb_solve: Solver_type 0, Solver id 0, Comp_type 0
Solution in struct 0, block 35, dof_ent 6, posglob=35
399.9999997296
Solution in struct 6, block 33, dof_ent 1, posglob=33
372.5504471529
Solution in struct 38, block 11, dof_ent 4, posglob=11
359.1602340235
Solution in struct 40, block 57, dof_ent 7, posglob=57
359.1555748521
Solution in struct 42, block 19, dof_ent 5, posglob=19
399.9999997296
Solution in struct 48, block 15, dof_ent 2, posglob=15
300.0000001205
Solution in struct 52, block 51, dof_ent 9, posglob=51
300.0000001205
Solution in struct 62, block 1, dof_ent 3, posglob=1
300.0000001003
Solution in struct 64, block 64, dof_ent 8, posglob=64
300.0000001002
After linear solver in nonlinear iteration 16
Solution difference norm (u_k, prev. u_k): 0.000000263265 (limit 0.000001000000)
Convergence in non-linear iterations: 0.000000263265 < 0.000001000000 - breaking.
---------------------------------------------------------------------------
After non-linear solver in time step 1
Solution difference norm (u_n, prev. u_n): 99.999999877692 (limit 0.000001000000)
Final time reached. Stopping.
Entering lsr_mkb_free_matrix: Solver_type 0, Solver id 0
Entering lar_free_SM_and_LV: storage_type: 2
matrix id: in interface 0, internally in implementation 0
Entering lsr_mkb_destroy: Solver_type 0, Solver id 0
Leaving linear solver 0 (sir_destroy)
**************** LEAVING TIME INTEGRATION ***************
Execution time total: 0.064610
