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: 35.000000
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
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.9999995096
Solution in struct 6, block 33, dof_ent 1, posglob=33
350.8440183856
Solution in struct 38, block 11, dof_ent 4, posglob=11
332.3578945700
Solution in struct 40, block 57, dof_ent 7, posglob=57
332.3516102934
Solution in struct 42, block 19, dof_ent 5, posglob=19
399.9999995096
Solution in struct 48, block 15, dof_ent 2, posglob=15
300.0000001911
Solution in struct 52, block 51, dof_ent 9, posglob=51
300.0000001910
Solution in struct 62, block 1, dof_ent 3, posglob=1
300.0000001686
Solution in struct 64, block 64, dof_ent 8, posglob=64
300.0000001686
After linear solver in nonlinear iteration 0
Solution difference norm (u_k, prev. u_k): 99.999999778416 (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.9999995096
Solution in struct 6, block 33, dof_ent 1, posglob=33
350.8440183856
Solution in struct 38, block 11, dof_ent 4, posglob=11
332.3578945700
Solution in struct 40, block 57, dof_ent 7, posglob=57
332.3516102934
Solution in struct 42, block 19, dof_ent 5, posglob=19
399.9999995096
Solution in struct 48, block 15, dof_ent 2, posglob=15
300.0000001911
Solution in struct 52, block 51, dof_ent 9, posglob=51
300.0000001910
Solution in struct 62, block 1, dof_ent 3, posglob=1
300.0000001686
Solution in struct 64, block 64, dof_ent 8, posglob=64
300.0000001686
After linear solver in nonlinear iteration 1
Solution difference norm (u_k, prev. u_k): 0.000000000000 (limit 0.000001000000)
Convergence in non-linear iterations: 0.000000000000 < 0.000001000000 - breaking.
---------------------------------------------------------------------------
After non-linear solver in time step 1
Solution difference norm (u_n, prev. u_n): 99.999999778416 (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.041004
