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: -1.000000 ambient temperature: 300.000000 no file with materials configuration (data in problem_heat.dat file) single material data (not temperature dependent) thermal_conductivity: 35.000000 density: 7000.000000 specific_heat: 700.000000 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.040047