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