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