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.9999998819 Solution in struct 6, block 33, dof_ent 1, posglob=33 387.0342805637 Solution in struct 38, block 11, dof_ent 4, posglob=11 373.3918288190 Solution in struct 40, block 57, dof_ent 7, posglob=57 373.3843029598 Solution in struct 42, block 19, dof_ent 5, posglob=19 399.9999998819 Solution in struct 48, block 15, dof_ent 2, posglob=15 300.0000000470 Solution in struct 52, block 51, dof_ent 9, posglob=51 300.0000000470 Solution in struct 62, block 1, dof_ent 3, posglob=1 300.0000000399 Solution in struct 64, block 64, dof_ent 8, posglob=64 300.0000000399 After linear solver in nonlinear iteration 1 Solution difference norm (u_k, prev. u_k): 41.668826421575 (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.9999997670 Solution in struct 6, block 33, dof_ent 1, posglob=33 376.5720549669 Solution in struct 38, block 11, dof_ent 4, posglob=11 368.0184427079 Solution in struct 40, block 57, dof_ent 7, posglob=57 368.0168453665 Solution in struct 42, block 19, dof_ent 5, posglob=19 399.9999997669 Solution in struct 48, block 15, dof_ent 2, posglob=15 300.0000001235 Solution in struct 52, block 51, dof_ent 9, posglob=51 300.0000001235 Solution in struct 62, block 1, dof_ent 3, posglob=1 300.0000000778 Solution in struct 64, block 64, dof_ent 8, posglob=64 300.0000000778 After linear solver in nonlinear iteration 2 Solution difference norm (u_k, prev. u_k): 22.051228235256 (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.9999996749 Solution in struct 6, block 33, dof_ent 1, posglob=33 367.4632974447 Solution in struct 38, block 11, dof_ent 4, posglob=11 354.9155597059 Solution in struct 40, block 57, dof_ent 7, posglob=57 354.9098631058 Solution in struct 42, block 19, dof_ent 5, posglob=19 399.9999996749 Solution in struct 48, block 15, dof_ent 2, posglob=15 300.0000001481 Solution in struct 52, block 51, dof_ent 9, posglob=51 300.0000001481 Solution in struct 62, block 1, dof_ent 3, posglob=1 300.0000001253 Solution in struct 64, block 64, dof_ent 8, posglob=64 300.0000001253 After linear solver in nonlinear iteration 3 Solution difference norm (u_k, prev. u_k): 15.833593412375 (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.9999997763 Solution in struct 6, block 33, dof_ent 1, posglob=33 377.5864737367 Solution in struct 38, block 11, dof_ent 4, posglob=11 366.6406336459 Solution in struct 40, block 57, dof_ent 7, posglob=57 366.6376884087 Solution in struct 42, block 19, dof_ent 5, posglob=19 399.9999997763 Solution in struct 48, block 15, dof_ent 2, posglob=15 300.0000001021 Solution in struct 52, block 51, dof_ent 9, posglob=51 300.0000001021 Solution in struct 62, block 1, dof_ent 3, posglob=1 300.0000000820 Solution in struct 64, block 64, dof_ent 8, posglob=64 300.0000000820 After linear solver in nonlinear iteration 4 Solution difference norm (u_k, prev. u_k): 11.927967023985 (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.9999997005 Solution in struct 6, block 33, dof_ent 1, posglob=33 370.0080694900 Solution in struct 38, block 11, dof_ent 4, posglob=11 358.4316293480 Solution in struct 40, block 57, dof_ent 7, posglob=57 358.4273122719 Solution in struct 42, block 19, dof_ent 5, posglob=19 399.9999997005 Solution in struct 48, block 15, dof_ent 2, posglob=15 300.0000001421 Solution in struct 52, block 51, dof_ent 9, posglob=51 300.0000001421 Solution in struct 62, block 1, dof_ent 3, posglob=1 300.0000001124 Solution in struct 64, block 64, dof_ent 8, posglob=64 300.0000001124 After linear solver in nonlinear iteration 5 Solution difference norm (u_k, prev. u_k): 8.345327108231 (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.9999997506 Solution in struct 6, block 33, dof_ent 1, posglob=33 375.0576010559 Solution in struct 38, block 11, dof_ent 4, posglob=11 363.8531160391 Solution in struct 40, block 57, dof_ent 7, posglob=57 363.8491369104 Solution in struct 42, block 19, dof_ent 5, posglob=19 399.9999997506 Solution in struct 48, block 15, dof_ent 2, posglob=15 300.0000001153 Solution in struct 52, block 51, dof_ent 9, posglob=51 300.0000001153 Solution in struct 62, block 1, dof_ent 3, posglob=1 300.0000000942 Solution in struct 64, block 64, dof_ent 8, posglob=64 300.0000000942 After linear solver in nonlinear iteration 6 Solution difference norm (u_k, prev. u_k): 5.545989739996 (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.9999997163 Solution in struct 6, block 33, dof_ent 1, posglob=33 371.5907722327 Solution in struct 38, block 11, dof_ent 4, posglob=11 360.2526414478 Solution in struct 40, block 57, dof_ent 7, posglob=57 360.2483215230 Solution in struct 42, block 19, dof_ent 5, posglob=19 399.9999997163 Solution in struct 48, block 15, dof_ent 2, posglob=15 300.0000001343 Solution in struct 52, block 51, dof_ent 9, posglob=51 300.0000001342 Solution in struct 62, block 1, dof_ent 3, posglob=1 300.0000001061 Solution in struct 64, block 64, dof_ent 8, posglob=64 300.0000001061 After linear solver in nonlinear iteration 7 Solution difference norm (u_k, prev. u_k): 3.695900724273 (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.9999997383 Solution in struct 6, block 33, dof_ent 1, posglob=33 373.8270427609 Solution in struct 38, block 11, dof_ent 4, posglob=11 362.5752739362 Solution in struct 40, block 57, dof_ent 7, posglob=57 362.5714953318 Solution in struct 42, block 19, dof_ent 5, posglob=19 399.9999997383 Solution in struct 48, block 15, dof_ent 2, posglob=15 300.0000001219 Solution in struct 52, block 51, dof_ent 9, posglob=51 300.0000001218 Solution in struct 62, block 1, dof_ent 3, posglob=1 300.0000000987 Solution in struct 64, block 64, dof_ent 8, posglob=64 300.0000000987 After linear solver in nonlinear iteration 8 Solution difference norm (u_k, prev. u_k): 2.399732208651 (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.9999997240 Solution in struct 6, block 33, dof_ent 1, posglob=33 372.3727688657 Solution in struct 38, block 11, dof_ent 4, posglob=11 361.0701991657 Solution in struct 40, block 57, dof_ent 7, posglob=57 361.0659182065 Solution in struct 42, block 19, dof_ent 5, posglob=19 399.9999997239 Solution in struct 48, block 15, dof_ent 2, posglob=15 300.0000001299 Solution in struct 52, block 51, dof_ent 9, posglob=51 300.0000001299 Solution in struct 62, block 1, dof_ent 3, posglob=1 300.0000001036 Solution in struct 64, block 64, dof_ent 8, posglob=64 300.0000001035 After linear solver in nonlinear iteration 9 Solution difference norm (u_k, prev. u_k): 1.552832276081 (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.9999997332 Solution in struct 6, block 33, dof_ent 1, posglob=33 373.3064736863 Solution in struct 38, block 11, dof_ent 4, posglob=11 362.0351198383 Solution in struct 40, block 57, dof_ent 7, posglob=57 362.0312901462 Solution in struct 42, block 19, dof_ent 5, posglob=19 399.9999997331 Solution in struct 48, block 15, dof_ent 2, posglob=15 300.0000001247 Solution in struct 52, block 51, dof_ent 9, posglob=51 300.0000001247 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.0000001004 After linear solver in nonlinear iteration 10 Solution difference norm (u_k, prev. u_k): 0.996379051914 (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.9999997272 Solution in struct 6, block 33, dof_ent 1, posglob=33 372.7050093747 Solution in struct 38, block 11, dof_ent 4, posglob=11 361.4151189397 Solution in struct 40, block 57, dof_ent 7, posglob=57 361.4108894514 Solution in struct 42, block 19, dof_ent 5, posglob=19 399.9999997272 Solution in struct 48, block 15, dof_ent 2, posglob=15 300.0000001281 Solution in struct 52, block 51, dof_ent 9, posglob=51 300.0000001280 Solution in struct 62, block 1, dof_ent 3, posglob=1 300.0000001025 Solution in struct 64, block 64, dof_ent 8, posglob=64 300.0000001024 After linear solver in nonlinear iteration 11 Solution difference norm (u_k, prev. u_k): 0.639576341159 (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.9999997310 Solution in struct 6, block 33, dof_ent 1, posglob=33 373.0902664584 Solution in struct 38, block 11, dof_ent 4, posglob=11 361.8114995921 Solution in struct 40, block 57, dof_ent 7, posglob=57 361.8075818085 Solution in struct 42, block 19, dof_ent 5, posglob=19 399.9999997310 Solution in struct 48, block 15, dof_ent 2, posglob=15 300.0000001259 Solution in struct 52, block 51, dof_ent 9, posglob=51 300.0000001259 Solution in struct 62, block 1, dof_ent 3, posglob=1 300.0000001012 Solution in struct 64, block 64, dof_ent 8, posglob=64 300.0000001012 After linear solver in nonlinear iteration 12 Solution difference norm (u_k, prev. u_k): 0.408844941220 (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.9999997286 Solution in struct 6, block 33, dof_ent 1, posglob=33 372.8434763080 Solution in struct 38, block 11, dof_ent 4, posglob=11 361.5574177852 Solution in struct 40, block 57, dof_ent 7, posglob=57 361.5532727710 Solution in struct 42, block 19, dof_ent 5, posglob=19 399.9999997286 Solution in struct 48, block 15, dof_ent 2, posglob=15 300.0000001273 Solution in struct 52, block 51, dof_ent 9, posglob=51 300.0000001273 Solution in struct 62, block 1, dof_ent 3, posglob=1 300.0000001020 Solution in struct 64, block 64, dof_ent 8, posglob=64 300.0000001020 After linear solver in nonlinear iteration 13 Solution difference norm (u_k, prev. u_k): 0.262146805633 (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.9999997302 Solution in struct 6, block 33, dof_ent 1, posglob=33 373.0014819334 Solution in struct 38, block 11, dof_ent 4, posglob=11 361.7200123031 Solution in struct 40, block 57, dof_ent 7, posglob=57 361.7160258389 Solution in struct 42, block 19, dof_ent 5, posglob=19 399.9999997301 Solution in struct 48, block 15, dof_ent 2, posglob=15 300.0000001264 Solution in struct 52, block 51, dof_ent 9, posglob=51 300.0000001264 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.0000001014 After linear solver in nonlinear iteration 14 Solution difference norm (u_k, prev. u_k): 0.167749846696 (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.9999997292 Solution in struct 6, block 33, dof_ent 1, posglob=33 372.9002905695 Solution in struct 38, block 11, dof_ent 4, posglob=11 361.6158510200 Solution in struct 40, block 57, dof_ent 7, posglob=57 361.6117565058 Solution in struct 42, block 19, dof_ent 5, posglob=19 399.9999997291 Solution in struct 48, block 15, dof_ent 2, posglob=15 300.0000001270 Solution in struct 52, block 51, dof_ent 9, posglob=51 300.0000001269 Solution in struct 62, block 1, dof_ent 3, posglob=1 300.0000001018 Solution in struct 64, block 64, dof_ent 8, posglob=64 300.0000001018 After linear solver in nonlinear iteration 15 Solution difference norm (u_k, prev. u_k): 0.107478373910 (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.9999997298 Solution in struct 6, block 33, dof_ent 1, posglob=33 372.9650844206 Solution in struct 38, block 11, dof_ent 4, posglob=11 361.6825332747 Solution in struct 40, block 57, dof_ent 7, posglob=57 361.6785110489 Solution in struct 42, block 19, dof_ent 5, posglob=19 399.9999997298 Solution in struct 48, block 15, dof_ent 2, posglob=15 300.0000001266 Solution in struct 52, block 51, dof_ent 9, posglob=51 300.0000001266 Solution in struct 62, block 1, dof_ent 3, posglob=1 300.0000001016 Solution in struct 64, block 64, dof_ent 8, posglob=64 300.0000001016 After linear solver in nonlinear iteration 16 Solution difference norm (u_k, prev. u_k): 0.068805813708 (limit 0.000001000000) --------------------------------------------------------------------------- Solving nonlinear iteration 17 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.9235910842 Solution in struct 38, block 11, dof_ent 4, posglob=11 361.6398254161 Solution in struct 40, block 57, dof_ent 7, posglob=57 361.6357553568 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.0000001269 Solution in struct 52, block 51, dof_ent 9, posglob=51 300.0000001268 Solution in struct 62, block 1, dof_ent 3, posglob=1 300.0000001017 Solution in struct 64, block 64, dof_ent 8, posglob=64 300.0000001017 After linear solver in nonlinear iteration 17 Solution difference norm (u_k, prev. u_k): 0.044070833229 (limit 0.000001000000) --------------------------------------------------------------------------- Solving nonlinear iteration 18 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.9501610331 Solution in struct 38, block 11, dof_ent 4, posglob=11 361.6671708332 Solution in struct 40, block 57, dof_ent 7, posglob=57 361.6631321419 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.0000001267 Solution in struct 52, block 51, dof_ent 9, posglob=51 300.0000001267 Solution in struct 62, block 1, dof_ent 3, posglob=1 300.0000001016 Solution in struct 64, block 64, dof_ent 8, posglob=64 300.0000001016 After linear solver in nonlinear iteration 18 Solution difference norm (u_k, prev. u_k): 0.028218343437 (limit 0.000001000000) --------------------------------------------------------------------------- Solving nonlinear iteration 19 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.9999997295 Solution in struct 6, block 33, dof_ent 1, posglob=33 372.9331462750 Solution in struct 38, block 11, dof_ent 4, posglob=11 361.6496586148 Solution in struct 40, block 57, dof_ent 7, posglob=57 361.6455994710 Solution in struct 42, block 19, dof_ent 5, posglob=19 399.9999997294 Solution in struct 48, block 15, dof_ent 2, posglob=15 300.0000001268 Solution in struct 52, block 51, dof_ent 9, posglob=51 300.0000001268 Solution in struct 62, block 1, dof_ent 3, posglob=1 300.0000001017 Solution in struct 64, block 64, dof_ent 8, posglob=64 300.0000001017 After linear solver in nonlinear iteration 19 Solution difference norm (u_k, prev. u_k): 0.018071845284 (limit 0.000001000000) --------------------------------------------------------------------------- Solving nonlinear iteration 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.9999997296 Solution in struct 6, block 33, dof_ent 1, posglob=33 372.9440417946 Solution in struct 38, block 11, dof_ent 4, posglob=11 361.6608723060 Solution in struct 40, block 57, dof_ent 7, posglob=57 361.6568264351 Solution in struct 42, block 19, dof_ent 5, posglob=19 399.9999997295 Solution in struct 48, block 15, dof_ent 2, posglob=15 300.0000001267 Solution in struct 52, block 51, dof_ent 9, posglob=51 300.0000001267 Solution in struct 62, block 1, dof_ent 3, posglob=1 300.0000001017 Solution in struct 64, block 64, dof_ent 8, posglob=64 300.0000001016 After linear solver in nonlinear iteration 20 Solution difference norm (u_k, prev. u_k): 0.011572132054 (limit 0.000001000000) --------------------------------------------------------------------------- Solving nonlinear iteration 21 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.9999997295 Solution in struct 6, block 33, dof_ent 1, posglob=33 372.9370646191 Solution in struct 38, block 11, dof_ent 4, posglob=11 361.6536912384 Solution in struct 40, block 57, dof_ent 7, posglob=57 361.6496367815 Solution in struct 42, block 19, dof_ent 5, posglob=19 399.9999997295 Solution in struct 48, block 15, dof_ent 2, posglob=15 300.0000001268 Solution in struct 52, block 51, dof_ent 9, posglob=51 300.0000001267 Solution in struct 62, block 1, dof_ent 3, posglob=1 300.0000001017 Solution in struct 64, block 64, dof_ent 8, posglob=64 300.0000001017 After linear solver in nonlinear iteration 21 Solution difference norm (u_k, prev. u_k): 0.007410733933 (limit 0.000001000000) --------------------------------------------------------------------------- Solving nonlinear iteration 22 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.9415325473 Solution in struct 38, block 11, dof_ent 4, posglob=11 361.6582896669 Solution in struct 40, block 57, dof_ent 7, posglob=57 361.6542407498 Solution in struct 42, block 19, dof_ent 5, posglob=19 399.9999997295 Solution in struct 48, block 15, dof_ent 2, posglob=15 300.0000001268 Solution in struct 52, block 51, dof_ent 9, posglob=51 300.0000001267 Solution in struct 62, block 1, dof_ent 3, posglob=1 300.0000001017 Solution in struct 64, block 64, dof_ent 8, posglob=64 300.0000001016 After linear solver in nonlinear iteration 22 Solution difference norm (u_k, prev. u_k): 0.004745524317 (limit 0.000001000000) --------------------------------------------------------------------------- Solving nonlinear iteration 23 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.9999997295 Solution in struct 6, block 33, dof_ent 1, posglob=33 372.9386714230 Solution in struct 38, block 11, dof_ent 4, posglob=11 361.6553449502 Solution in struct 40, block 57, dof_ent 7, posglob=57 361.6512924650 Solution in struct 42, block 19, dof_ent 5, posglob=19 399.9999997295 Solution in struct 48, block 15, dof_ent 2, posglob=15 300.0000001268 Solution in struct 52, block 51, dof_ent 9, posglob=51 300.0000001267 Solution in struct 62, block 1, dof_ent 3, posglob=1 300.0000001017 Solution in struct 64, block 64, dof_ent 8, posglob=64 300.0000001017 After linear solver in nonlinear iteration 23 Solution difference norm (u_k, prev. u_k): 0.003038940670 (limit 0.000001000000) --------------------------------------------------------------------------- Solving nonlinear iteration 24 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.9999997295 Solution in struct 6, block 33, dof_ent 1, posglob=33 372.9405035907 Solution in struct 38, block 11, dof_ent 4, posglob=11 361.6572306365 Solution in struct 40, block 57, dof_ent 7, posglob=57 361.6531804463 Solution in struct 42, block 19, dof_ent 5, posglob=19 399.9999997295 Solution in struct 48, block 15, dof_ent 2, posglob=15 300.0000001268 Solution in struct 52, block 51, dof_ent 9, posglob=51 300.0000001267 Solution in struct 62, block 1, dof_ent 3, posglob=1 300.0000001017 Solution in struct 64, block 64, dof_ent 8, posglob=64 300.0000001017 After linear solver in nonlinear iteration 24 Solution difference norm (u_k, prev. u_k): 0.001946031693 (limit 0.000001000000) --------------------------------------------------------------------------- Solving nonlinear iteration 25 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.9999997295 Solution in struct 6, block 33, dof_ent 1, posglob=33 372.9393303275 Solution in struct 38, block 11, dof_ent 4, posglob=11 361.6560230979 Solution in struct 40, block 57, dof_ent 7, posglob=57 361.6519714331 Solution in struct 42, block 19, dof_ent 5, posglob=19 399.9999997295 Solution in struct 48, block 15, dof_ent 2, posglob=15 300.0000001268 Solution in struct 52, block 51, dof_ent 9, posglob=51 300.0000001267 Solution in struct 62, block 1, dof_ent 3, posglob=1 300.0000001017 Solution in struct 64, block 64, dof_ent 8, posglob=64 300.0000001017 After linear solver in nonlinear iteration 25 Solution difference norm (u_k, prev. u_k): 0.001246188164 (limit 0.000001000000) --------------------------------------------------------------------------- Solving nonlinear iteration 26 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.9999997295 Solution in struct 6, block 33, dof_ent 1, posglob=33 372.9400816475 Solution in struct 38, block 11, dof_ent 4, posglob=11 361.6567963649 Solution in struct 40, block 57, dof_ent 7, posglob=57 361.6527456466 Solution in struct 42, block 19, dof_ent 5, posglob=19 399.9999997295 Solution in struct 48, block 15, dof_ent 2, posglob=15 300.0000001268 Solution in struct 52, block 51, dof_ent 9, posglob=51 300.0000001267 Solution in struct 62, block 1, dof_ent 3, posglob=1 300.0000001017 Solution in struct 64, block 64, dof_ent 8, posglob=64 300.0000001017 After linear solver in nonlinear iteration 26 Solution difference norm (u_k, prev. u_k): 0.000798018678 (limit 0.000001000000) --------------------------------------------------------------------------- Solving nonlinear iteration 27 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.9999997295 Solution in struct 6, block 33, dof_ent 1, posglob=33 372.9396005257 Solution in struct 38, block 11, dof_ent 4, posglob=11 361.6563011882 Solution in struct 40, block 57, dof_ent 7, posglob=57 361.6522498627 Solution in struct 42, block 19, dof_ent 5, posglob=19 399.9999997295 Solution in struct 48, block 15, dof_ent 2, posglob=15 300.0000001268 Solution in struct 52, block 51, dof_ent 9, posglob=51 300.0000001267 Solution in struct 62, block 1, dof_ent 3, posglob=1 300.0000001017 Solution in struct 64, block 64, dof_ent 8, posglob=64 300.0000001017 After linear solver in nonlinear iteration 27 Solution difference norm (u_k, prev. u_k): 0.000511028243 (limit 0.000001000000) --------------------------------------------------------------------------- Solving nonlinear iteration 28 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.9999997295 Solution in struct 6, block 33, dof_ent 1, posglob=33 372.9399086208 Solution in struct 38, block 11, dof_ent 4, posglob=11 361.6566182832 Solution in struct 40, block 57, dof_ent 7, posglob=57 361.6525673473 Solution in struct 42, block 19, dof_ent 5, posglob=19 399.9999997295 Solution in struct 48, block 15, dof_ent 2, posglob=15 300.0000001268 Solution in struct 52, block 51, dof_ent 9, posglob=51 300.0000001267 Solution in struct 62, block 1, dof_ent 3, posglob=1 300.0000001017 Solution in struct 64, block 64, dof_ent 8, posglob=64 300.0000001017 After linear solver in nonlinear iteration 28 Solution difference norm (u_k, prev. u_k): 0.000327246652 (limit 0.000001000000) --------------------------------------------------------------------------- Solving nonlinear iteration 29 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.9999997295 Solution in struct 6, block 33, dof_ent 1, posglob=33 372.9397113262 Solution in struct 38, block 11, dof_ent 4, posglob=11 361.6564152252 Solution in struct 40, block 57, dof_ent 7, posglob=57 361.6523640397 Solution in struct 42, block 19, dof_ent 5, posglob=19 399.9999997295 Solution in struct 48, block 15, dof_ent 2, posglob=15 300.0000001268 Solution in struct 52, block 51, dof_ent 9, posglob=51 300.0000001267 Solution in struct 62, block 1, dof_ent 3, posglob=1 300.0000001017 Solution in struct 64, block 64, dof_ent 8, posglob=64 300.0000001017 After linear solver in nonlinear iteration 29 Solution difference norm (u_k, prev. u_k): 0.000209558968 (limit 0.000001000000) --------------------------------------------------------------------------- Solving nonlinear iteration 30 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.9999997295 Solution in struct 6, block 33, dof_ent 1, posglob=33 372.9398376677 Solution in struct 38, block 11, dof_ent 4, posglob=11 361.6565452574 Solution in struct 40, block 57, dof_ent 7, posglob=57 361.6524942317 Solution in struct 42, block 19, dof_ent 5, posglob=19 399.9999997295 Solution in struct 48, block 15, dof_ent 2, posglob=15 300.0000001268 Solution in struct 52, block 51, dof_ent 9, posglob=51 300.0000001267 Solution in struct 62, block 1, dof_ent 3, posglob=1 300.0000001017 Solution in struct 64, block 64, dof_ent 8, posglob=64 300.0000001017 After linear solver in nonlinear iteration 30 Solution difference norm (u_k, prev. u_k): 0.000134195149 (limit 0.000001000000) --------------------------------------------------------------------------- Solving nonlinear iteration 31 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.9999997295 Solution in struct 6, block 33, dof_ent 1, posglob=33 372.9397567626 Solution in struct 38, block 11, dof_ent 4, posglob=11 361.6564619890 Solution in struct 40, block 57, dof_ent 7, posglob=57 361.6524108604 Solution in struct 42, block 19, dof_ent 5, posglob=19 399.9999997295 Solution in struct 48, block 15, dof_ent 2, posglob=15 300.0000001268 Solution in struct 52, block 51, dof_ent 9, posglob=51 300.0000001267 Solution in struct 62, block 1, dof_ent 3, posglob=1 300.0000001017 Solution in struct 64, block 64, dof_ent 8, posglob=64 300.0000001017 After linear solver in nonlinear iteration 31 Solution difference norm (u_k, prev. u_k): 0.000085934775 (limit 0.000001000000) --------------------------------------------------------------------------- Solving nonlinear iteration 32 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.9999997295 Solution in struct 6, block 33, dof_ent 1, posglob=33 372.9398085716 Solution in struct 38, block 11, dof_ent 4, posglob=11 361.6565153114 Solution in struct 40, block 57, dof_ent 7, posglob=57 361.6524642488 Solution in struct 42, block 19, dof_ent 5, posglob=19 399.9999997295 Solution in struct 48, block 15, dof_ent 2, posglob=15 300.0000001268 Solution in struct 52, block 51, dof_ent 9, posglob=51 300.0000001267 Solution in struct 62, block 1, dof_ent 3, posglob=1 300.0000001017 Solution in struct 64, block 64, dof_ent 8, posglob=64 300.0000001017 After linear solver in nonlinear iteration 32 Solution difference norm (u_k, prev. u_k): 0.000055029841 (limit 0.000001000000) --------------------------------------------------------------------------- Solving nonlinear iteration 33 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.9999997295 Solution in struct 6, block 33, dof_ent 1, posglob=33 372.9397753949 Solution in struct 38, block 11, dof_ent 4, posglob=11 361.6564811654 Solution in struct 40, block 57, dof_ent 7, posglob=57 361.6524300607 Solution in struct 42, block 19, dof_ent 5, posglob=19 399.9999997295 Solution in struct 48, block 15, dof_ent 2, posglob=15 300.0000001268 Solution in struct 52, block 51, dof_ent 9, posglob=51 300.0000001267 Solution in struct 62, block 1, dof_ent 3, posglob=1 300.0000001017 Solution in struct 64, block 64, dof_ent 8, posglob=64 300.0000001017 After linear solver in nonlinear iteration 33 Solution difference norm (u_k, prev. u_k): 0.000035239192 (limit 0.000001000000) --------------------------------------------------------------------------- Solving nonlinear iteration 34 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.9999997295 Solution in struct 6, block 33, dof_ent 1, posglob=33 372.9397966401 Solution in struct 38, block 11, dof_ent 4, posglob=11 361.6565030313 Solution in struct 40, block 57, dof_ent 7, posglob=57 361.6524519539 Solution in struct 42, block 19, dof_ent 5, posglob=19 399.9999997295 Solution in struct 48, block 15, dof_ent 2, posglob=15 300.0000001268 Solution in struct 52, block 51, dof_ent 9, posglob=51 300.0000001267 Solution in struct 62, block 1, dof_ent 3, posglob=1 300.0000001017 Solution in struct 64, block 64, dof_ent 8, posglob=64 300.0000001017 After linear solver in nonlinear iteration 34 Solution difference norm (u_k, prev. u_k): 0.000022566366 (limit 0.000001000000) --------------------------------------------------------------------------- Solving nonlinear iteration 35 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.9999997295 Solution in struct 6, block 33, dof_ent 1, posglob=33 372.9397830352 Solution in struct 38, block 11, dof_ent 4, posglob=11 361.6564890289 Solution in struct 40, block 57, dof_ent 7, posglob=57 361.6524379338 Solution in struct 42, block 19, dof_ent 5, posglob=19 399.9999997295 Solution in struct 48, block 15, dof_ent 2, posglob=15 300.0000001268 Solution in struct 52, block 51, dof_ent 9, posglob=51 300.0000001267 Solution in struct 62, block 1, dof_ent 3, posglob=1 300.0000001017 Solution in struct 64, block 64, dof_ent 8, posglob=64 300.0000001017 After linear solver in nonlinear iteration 35 Solution difference norm (u_k, prev. u_k): 0.000014451187 (limit 0.000001000000) --------------------------------------------------------------------------- Solving nonlinear iteration 36 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.9999997295 Solution in struct 6, block 33, dof_ent 1, posglob=33 372.9397917473 Solution in struct 38, block 11, dof_ent 4, posglob=11 361.6564979955 Solution in struct 40, block 57, dof_ent 7, posglob=57 361.6524469122 Solution in struct 42, block 19, dof_ent 5, posglob=19 399.9999997295 Solution in struct 48, block 15, dof_ent 2, posglob=15 300.0000001268 Solution in struct 52, block 51, dof_ent 9, posglob=51 300.0000001267 Solution in struct 62, block 1, dof_ent 3, posglob=1 300.0000001017 Solution in struct 64, block 64, dof_ent 8, posglob=64 300.0000001017 After linear solver in nonlinear iteration 36 Solution difference norm (u_k, prev. u_k): 0.000009254503 (limit 0.000001000000) --------------------------------------------------------------------------- Solving nonlinear iteration 37 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.9999997295 Solution in struct 6, block 33, dof_ent 1, posglob=33 372.9397861684 Solution in struct 38, block 11, dof_ent 4, posglob=11 361.6564922538 Solution in struct 40, block 57, dof_ent 7, posglob=57 361.6524411625 Solution in struct 42, block 19, dof_ent 5, posglob=19 399.9999997295 Solution in struct 48, block 15, dof_ent 2, posglob=15 300.0000001268 Solution in struct 52, block 51, dof_ent 9, posglob=51 300.0000001267 Solution in struct 62, block 1, dof_ent 3, posglob=1 300.0000001017 Solution in struct 64, block 64, dof_ent 8, posglob=64 300.0000001017 After linear solver in nonlinear iteration 37 Solution difference norm (u_k, prev. u_k): 0.000005926521 (limit 0.000001000000) --------------------------------------------------------------------------- Solving nonlinear iteration 38 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.9999997295 Solution in struct 6, block 33, dof_ent 1, posglob=33 372.9397897410 Solution in struct 38, block 11, dof_ent 4, posglob=11 361.6564959307 Solution in struct 40, block 57, dof_ent 7, posglob=57 361.6524448445 Solution in struct 42, block 19, dof_ent 5, posglob=19 399.9999997295 Solution in struct 48, block 15, dof_ent 2, posglob=15 300.0000001268 Solution in struct 52, block 51, dof_ent 9, posglob=51 300.0000001267 Solution in struct 62, block 1, dof_ent 3, posglob=1 300.0000001017 Solution in struct 64, block 64, dof_ent 8, posglob=64 300.0000001017 After linear solver in nonlinear iteration 38 Solution difference norm (u_k, prev. u_k): 0.000003795140 (limit 0.000001000000) --------------------------------------------------------------------------- Solving nonlinear iteration 39 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.9999997295 Solution in struct 6, block 33, dof_ent 1, posglob=33 372.9397874530 Solution in struct 38, block 11, dof_ent 4, posglob=11 361.6564935758 Solution in struct 40, block 57, dof_ent 7, posglob=57 361.6524424863 Solution in struct 42, block 19, dof_ent 5, posglob=19 399.9999997295 Solution in struct 48, block 15, dof_ent 2, posglob=15 300.0000001268 Solution in struct 52, block 51, dof_ent 9, posglob=51 300.0000001267 Solution in struct 62, block 1, dof_ent 3, posglob=1 300.0000001017 Solution in struct 64, block 64, dof_ent 8, posglob=64 300.0000001017 After linear solver in nonlinear iteration 39 Solution difference norm (u_k, prev. u_k): 0.000002430644 (limit 0.000001000000) --------------------------------------------------------------------------- Solving nonlinear iteration 40 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.9999997295 Solution in struct 6, block 33, dof_ent 1, posglob=33 372.9397889185 Solution in struct 38, block 11, dof_ent 4, posglob=11 361.6564950840 Solution in struct 40, block 57, dof_ent 7, posglob=57 361.6524439966 Solution in struct 42, block 19, dof_ent 5, posglob=19 399.9999997295 Solution in struct 48, block 15, dof_ent 2, posglob=15 300.0000001268 Solution in struct 52, block 51, dof_ent 9, posglob=51 300.0000001267 Solution in struct 62, block 1, dof_ent 3, posglob=1 300.0000001017 Solution in struct 64, block 64, dof_ent 8, posglob=64 300.0000001017 After linear solver in nonlinear iteration 40 Solution difference norm (u_k, prev. u_k): 0.000001556668 (limit 0.000001000000) --------------------------------------------------------------------------- Solving nonlinear iteration 41 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.9999997295 Solution in struct 6, block 33, dof_ent 1, posglob=33 372.9397879802 Solution in struct 38, block 11, dof_ent 4, posglob=11 361.6564941184 Solution in struct 40, block 57, dof_ent 7, posglob=57 361.6524430295 Solution in struct 42, block 19, dof_ent 5, posglob=19 399.9999997295 Solution in struct 48, block 15, dof_ent 2, posglob=15 300.0000001268 Solution in struct 52, block 51, dof_ent 9, posglob=51 300.0000001267 Solution in struct 62, block 1, dof_ent 3, posglob=1 300.0000001017 Solution in struct 64, block 64, dof_ent 8, posglob=64 300.0000001017 After linear solver in nonlinear iteration 41 Solution difference norm (u_k, prev. u_k): 0.000000996845 (limit 0.000001000000) Convergence in non-linear iterations: 0.000000996845 < 0.000001000000 - breaking. --------------------------------------------------------------------------- After non-linear solver in time step 1 Solution difference norm (u_n, prev. u_n): 99.999999877864 (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.098262