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.9999998073
Solution in struct 6, block 33, dof_ent 1, posglob=33
373.9366408641
Solution in struct 38, block 11, dof_ent 4, posglob=11
357.6355205280
Solution in struct 40, block 57, dof_ent 7, posglob=57
357.6297768712
Solution in struct 42, block 19, dof_ent 5, posglob=19
399.9999998073
Solution in struct 48, block 15, dof_ent 2, posglob=15
300.0000000834
Solution in struct 52, block 51, dof_ent 9, posglob=51
300.0000000833
Solution in struct 62, block 1, dof_ent 3, posglob=1
300.0000000700
Solution in struct 64, block 64, dof_ent 8, posglob=64
300.0000000700
After linear solver in nonlinear iteration 1
Solution difference norm (u_k, prev. u_k): 25.703374570808 (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.9999997204
Solution in struct 6, block 33, dof_ent 1, posglob=33
367.9790512377
Solution in struct 38, block 11, dof_ent 4, posglob=11
352.1385292430
Solution in struct 40, block 57, dof_ent 7, posglob=57
352.1327516579
Solution in struct 42, block 19, dof_ent 5, posglob=19
399.9999997203
Solution in struct 48, block 15, dof_ent 2, posglob=15
300.0000001217
Solution in struct 52, block 51, dof_ent 9, posglob=51
300.0000001216
Solution in struct 62, block 1, dof_ent 3, posglob=1
300.0000001026
Solution in struct 64, block 64, dof_ent 8, posglob=64
300.0000001025
After linear solver in nonlinear iteration 2
Solution difference norm (u_k, prev. u_k): 6.006520162682 (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.9999997359
Solution in struct 6, block 33, dof_ent 1, posglob=33
368.4576909790
Solution in struct 38, block 11, dof_ent 4, posglob=11
352.2211510346
Solution in struct 40, block 57, dof_ent 7, posglob=57
352.2151210669
Solution in struct 42, block 19, dof_ent 5, posglob=19
399.9999997358
Solution in struct 48, block 15, dof_ent 2, posglob=15
300.0000001143
Solution in struct 52, block 51, dof_ent 9, posglob=51
300.0000001143
Solution in struct 62, block 1, dof_ent 3, posglob=1
300.0000000971
Solution in struct 64, block 64, dof_ent 8, posglob=64
300.0000000971
After linear solver in nonlinear iteration 3
Solution difference norm (u_k, prev. u_k): 0.549389026683 (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.9999997358
Solution in struct 6, block 33, dof_ent 1, posglob=33
368.5941877486
Solution in struct 38, block 11, dof_ent 4, posglob=11
352.4176767511
Solution in struct 40, block 57, dof_ent 7, posglob=57
352.4117489769
Solution in struct 42, block 19, dof_ent 5, posglob=19
399.9999997357
Solution in struct 48, block 15, dof_ent 2, posglob=15
300.0000001144
Solution in struct 52, block 51, dof_ent 9, posglob=51
300.0000001144
Solution in struct 62, block 1, dof_ent 3, posglob=1
300.0000000970
Solution in struct 64, block 64, dof_ent 8, posglob=64
300.0000000970
After linear solver in nonlinear iteration 4
Solution difference norm (u_k, prev. u_k): 0.203469624452 (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.9999997352
Solution in struct 6, block 33, dof_ent 1, posglob=33
368.5453338723
Solution in struct 38, block 11, dof_ent 4, posglob=11
352.3723603601
Solution in struct 40, block 57, dof_ent 7, posglob=57
352.3664194725
Solution in struct 42, block 19, dof_ent 5, posglob=19
399.9999997351
Solution in struct 48, block 15, dof_ent 2, posglob=15
300.0000001147
Solution in struct 52, block 51, dof_ent 9, posglob=51
300.0000001147
Solution in struct 62, block 1, dof_ent 3, posglob=1
300.0000000973
Solution in struct 64, block 64, dof_ent 8, posglob=64
300.0000000972
After linear solver in nonlinear iteration 5
Solution difference norm (u_k, prev. u_k): 0.049578949504 (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.9999997353
Solution in struct 6, block 33, dof_ent 1, posglob=33
368.5502793512
Solution in struct 38, block 11, dof_ent 4, posglob=11
352.3741456024
Solution in struct 40, block 57, dof_ent 7, posglob=57
352.3682023127
Solution in struct 42, block 19, dof_ent 5, posglob=19
399.9999997353
Solution in struct 48, block 15, dof_ent 2, posglob=15
300.0000001146
Solution in struct 52, block 51, dof_ent 9, posglob=51
300.0000001146
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 6
Solution difference norm (u_k, prev. u_k): 0.005236106696 (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.9999997353
Solution in struct 6, block 33, dof_ent 1, posglob=33
368.5511818222
Solution in struct 38, block 11, dof_ent 4, posglob=11
352.3755768568
Solution in struct 40, block 57, dof_ent 7, posglob=57
352.3696348619
Solution in struct 42, block 19, dof_ent 5, posglob=19
399.9999997353
Solution in struct 48, block 15, dof_ent 2, posglob=15
300.0000001146
Solution in struct 52, block 51, dof_ent 9, posglob=51
300.0000001146
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 7
Solution difference norm (u_k, prev. u_k): 0.001491723803 (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.9999997353
Solution in struct 6, block 33, dof_ent 1, posglob=33
368.5507921096
Solution in struct 38, block 11, dof_ent 4, posglob=11
352.3752011041
Solution in struct 40, block 57, dof_ent 7, posglob=57
352.3692589135
Solution in struct 42, block 19, dof_ent 5, posglob=19
399.9999997353
Solution in struct 48, block 15, dof_ent 2, posglob=15
300.0000001146
Solution in struct 52, block 51, dof_ent 9, posglob=51
300.0000001146
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 8
Solution difference norm (u_k, prev. u_k): 0.000399138027 (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.9999997353
Solution in struct 6, block 33, dof_ent 1, posglob=33
368.5508389023
Solution in struct 38, block 11, dof_ent 4, posglob=11
352.3752237929
Solution in struct 40, block 57, dof_ent 7, posglob=57
352.3692815874
Solution in struct 42, block 19, dof_ent 5, posglob=19
399.9999997353
Solution in struct 48, block 15, dof_ent 2, posglob=15
300.0000001146
Solution in struct 52, block 51, dof_ent 9, posglob=51
300.0000001146
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 9
Solution difference norm (u_k, prev. u_k): 0.000047263737 (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.9999997353
Solution in struct 6, block 33, dof_ent 1, posglob=33
368.5508445313
Solution in struct 38, block 11, dof_ent 4, posglob=11
352.3752339888
Solution in struct 40, block 57, dof_ent 7, posglob=57
352.3692917961
Solution in struct 42, block 19, dof_ent 5, posglob=19
399.9999997353
Solution in struct 48, block 15, dof_ent 2, posglob=15
300.0000001146
Solution in struct 52, block 51, dof_ent 9, posglob=51
300.0000001146
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 10
Solution difference norm (u_k, prev. u_k): 0.000010808544 (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.9999997353
Solution in struct 6, block 33, dof_ent 1, posglob=33
368.5508414636
Solution in struct 38, block 11, dof_ent 4, posglob=11
352.3752309135
Solution in struct 40, block 57, dof_ent 7, posglob=57
352.3692887186
Solution in struct 42, block 19, dof_ent 5, posglob=19
399.9999997353
Solution in struct 48, block 15, dof_ent 2, posglob=15
300.0000001146
Solution in struct 52, block 51, dof_ent 9, posglob=51
300.0000001146
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 11
Solution difference norm (u_k, prev. u_k): 0.000003172382 (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.9999997353
Solution in struct 6, block 33, dof_ent 1, posglob=33
368.5508418909
Solution in struct 38, block 11, dof_ent 4, posglob=11
352.3752311599
Solution in struct 40, block 57, dof_ent 7, posglob=57
352.3692889650
Solution in struct 42, block 19, dof_ent 5, posglob=19
399.9999997353
Solution in struct 48, block 15, dof_ent 2, posglob=15
300.0000001146
Solution in struct 52, block 51, dof_ent 9, posglob=51
300.0000001146
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 12
Solution difference norm (u_k, prev. u_k): 0.000000427268 (limit 0.000001000000)
Convergence in non-linear iterations: 0.000000427268 < 0.000001000000 - breaking.
---------------------------------------------------------------------------
After non-linear solver in time step 1
Solution difference norm (u_n, prev. u_n): 99.999999875509 (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.057822
