<<< directory: c:\Python27\lib\site-packages\sfepy\tests, test files: 90 <<< c:\Python27\lib\site-packages\sfepy\tests\test_assembling.py sfepy: left over: ['nm', '__builtins__', '_filename', '__file__', 'sps', '__package__', '__doc__', '__name__', 'TestCommon', 'verbose'] >>> test instance prepared (4 test(s)) ... assembled: [[ 1. 1. 1. 0. 0.] [ 1. 1. 1. 0. 0.] [ 1. 1. 3. 2. 2.] [ 0. 0. 2. 2. 2.] [ 0. 0. 2. 2. 2.]] ... expected: [[ 1. 1. 1. 0. 0.] [ 1. 1. 1. 0. 0.] [ 1. 1. 3. 2. 2.] [ 0. 0. 2. 2. 2.] [ 0. 0. 2. 2. 2.]] ... ||assembled - expected||: 0.000000e+00 +++ test_assemble_matrix: ok ... assembled: [[ 2.-3.j 2.-3.j 2.-3.j 0.+0.j 0.+0.j] [ 2.-3.j 2.-3.j 2.-3.j 0.+0.j 0.+0.j] [ 2.-3.j 2.-3.j 6.-9.j 4.-6.j 4.-6.j] [ 0.+0.j 0.+0.j 4.-6.j 4.-6.j 4.-6.j] [ 0.+0.j 0.+0.j 4.-6.j 4.-6.j 4.-6.j]] ... expected: [[ 2.-3.j 2.-3.j 2.-3.j 0.+0.j 0.+0.j] [ 2.-3.j 2.-3.j 2.-3.j 0.+0.j 0.+0.j] [ 2.-3.j 2.-3.j 6.-9.j 4.-6.j 4.-6.j] [ 0.+0.j 0.+0.j 4.-6.j 4.-6.j 4.-6.j] [ 0.+0.j 0.+0.j 4.-6.j 4.-6.j 4.-6.j]] ... ||assembled - expected||: 0.000000e+00 +++ test_assemble_matrix_complex: ok ... assembled: [ 1. 1. 3. 2. 2.] ... expected: [ 1. 1. 3. 2. 2.] ... ||assembled - expected||: 0.000000e+00 +++ test_assemble_vector: ok ... assembled: [ 2.-3.j 2.-3.j 6.-9.j 4.-6.j 4.-6.j] ... expected: [ 2.-3.j 2.-3.j 6.-9.j 4.-6.j 4.-6.j] ... ||assembled - expected||: 0.000000e+00 +++ test_assemble_vector_complex: ok >>> all passed in 0.20 s <<< c:\Python27\lib\site-packages\sfepy\tests\test_base.py sfepy: left over: ['verbose', '__builtins__', '_filename', '__file__', '__package__', '__doc__', '__name__', 'assert_', 'TestCommon'] >>> test instance prepared (4 test(s)) +++ test_resolve_deps: ok +++ test_struct_add: ok +++ test_struct_i_add: ok +++ test_verbose_output: ok >>> all passed in 0.00 s <<< c:\Python27\lib\site-packages\sfepy\tests\test_cmesh.py sfepy: left over: ['data_dir', 'nm', '__builtins__', '_filename', '__file__', '__package__', 'expected', '__doc__', '__name__', 'os', 'TestCommon', 'verbose'] >>> test instance prepared (1 test(s)) ... mesh: 1_2_2.mesh sfepy: reading mesh [line2, tri3, quad4, tetra4, hexa8] (c:\Python27\lib\site-packages\sfepy/meshes/elements/1_2_2.mesh)... sfepy: ...done in 0.00 s ... dim: 1 ... n_vertex: 3, n_edge: 2, n_face: 0, n_cell: 0 ... (0, 0) : (3, 4) ... (0, 1) : (3, 4) ... (1, 1) : (2, 2) ... mesh: 2_3_2.mesh sfepy: reading mesh [line2, tri3, quad4, tetra4, hexa8] (c:\Python27\lib\site-packages\sfepy/meshes/elements/2_3_2.mesh)... sfepy: ...done in 0.00 s ... dim: 2 ... n_vertex: 4, n_edge: 5, n_face: 2, n_cell: 0 ... (0, 0) : (4, 10) ... (0, 1) : (4, 10) ... (0, 2) : (4, 6) ... (1, 0) : (5, 10) ... (1, 1) : (5, 16) ... (1, 2) : (5, 6) ... (2, 1) : (2, 6) ... (2, 2) : (2, 2) ... mesh: 2_4_2.mesh sfepy: reading mesh [line2, tri3, quad4, tetra4, hexa8] (c:\Python27\lib\site-packages\sfepy/meshes/elements/2_4_2.mesh)... sfepy: ...done in 0.00 s ... dim: 2 ... n_vertex: 6, n_edge: 7, n_face: 2, n_cell: 0 ... (0, 0) : (6, 22) ... (0, 1) : (6, 14) ... (0, 2) : (6, 8) ... (1, 0) : (7, 14) ... (1, 1) : (7, 20) ... (1, 2) : (7, 8) ... (2, 1) : (2, 8) ... (2, 2) : (2, 2) ... mesh: 3_4_2.mesh sfepy: reading mesh [line2, tri3, quad4, tetra4, hexa8] (c:\Python27\lib\site-packages\sfepy/meshes/elements/3_4_2.mesh)... sfepy: ...done in 0.00 s ... dim: 3 ... n_vertex: 5, n_edge: 9, n_face: 7, n_cell: 2 ... (0, 0) : (5, 18) ... (0, 1) : (5, 18) ... (0, 2) : (5, 21) ... (0, 3) : (5, 8) ... (1, 0) : (9, 18) ... (1, 1) : (9, 48) ... (1, 2) : (9, 21) ... (1, 3) : (9, 12) ... (2, 0) : (7, 21) ... (2, 1) : (7, 21) ... (2, 2) : (7, 42) ... (2, 3) : (7, 8) ... (3, 1) : (2, 12) ... (3, 2) : (2, 8) ... (3, 3) : (2, 2) ... mesh: 3_8_2.mesh sfepy: reading mesh [line2, tri3, quad4, tetra4, hexa8] (c:\Python27\lib\site-packages\sfepy/meshes/elements/3_8_2.mesh)... sfepy: ...done in 0.00 s ... dim: 3 ... n_vertex: 12, n_edge: 20, n_face: 11, n_cell: 2 ... (0, 0) : (12, 100) ... (0, 1) : (12, 40) ... (0, 2) : (12, 44) ... (0, 3) : (12, 16) ... (1, 0) : (20, 40) ... (1, 1) : (20, 96) ... (1, 2) : (20, 44) ... (1, 3) : (20, 24) ... (2, 0) : (11, 44) ... (2, 1) : (11, 44) ... (2, 2) : (11, 72) ... (2, 3) : (11, 12) ... (3, 1) : (2, 24) ... (3, 2) : (2, 12) ... (3, 3) : (2, 2) ... mesh: square_triquad.mesh sfepy: reading mesh [line2, tri3, quad4, tetra4, hexa8] (c:\Python27\lib\site-packages\sfepy/meshes/2d/special/square_triquad.mesh)... sfepy: ...done in 0.00 s ... dim: 2 ... n_vertex: 470, n_edge: 1127, n_face: 658, n_cell: 0 ... (0, 0) : (470, 3054) ... (0, 1) : (470, 2254) ... (0, 2) : (470, 2174) ... (1, 0) : (1127, 2254) ... (1, 1) : (1127, 9174) ... (1, 2) : (1127, 2174) ... (2, 1) : (658, 2174) ... (2, 2) : (658, 6686) +++ test_cmesh_counts: ok >>> all passed in 0.01 s <<< c:\Python27\lib\site-packages\sfepy\tests\test_domain.py sfepy: left over: ['expected_conn', 'expected_coors', 'verbose', '__builtins__', '__file__', 'compare_mesh', 'FEDomain', 'refine', '__name__', 'TestCommon', 'data_dir', 'nm', '_filename', '__package__', 'Mesh', '__doc__', 'op'] sfepy: reading mesh [line2, tri3, quad4, tetra4, hexa8] (c:\Python27\lib\site-packages\sfepy/meshes/various_formats/small3d.mesh)... sfepy: ...done in 0.00 s >>> test instance prepared (7 test(s)) ... unique edges: True ... unique faces: True +++ test_facets: ok sfepy: reading mesh [line2, tri3, quad4, tetra4, hexa8] (c:\Python27\lib\site-packages\sfepy/meshes/elements/2_3_1.mesh)... sfepy: ...done in 0.00 s array([[ 0. , 0. ], [ 1. , 0. ], [ 0. , 1. ], [ 0.5, 0. ], [ 0.5, 0.5], [ 0. , 0.5]]) array([[0, 3, 5], [3, 4, 5], [1, 4, 3], [2, 5, 4]]) +++ test_refine_2_3: ok sfepy: reading mesh [line2, tri3, quad4, tetra4, hexa8] (c:\Python27\lib\site-packages\sfepy/meshes/elements/2_4_1.mesh)... sfepy: ...done in 0.00 s array([[ 0. , 0. ], [ 1. , 0. ], [ 1. , 1. ], [ 0. , 1. ], [ 0.5, 0. ], [ 1. , 0.5], [ 0.5, 1. ], [ 0. , 0.5], [ 0.5, 0.5]]) array([[0, 4, 8, 7], [1, 5, 8, 4], [2, 6, 8, 5], [3, 7, 8, 6]]) +++ test_refine_2_4: ok sfepy: reading mesh [line2, tri3, quad4, tetra4, hexa8] (c:\Python27\lib\site-packages\sfepy/meshes/elements/3_4_1.mesh)... sfepy: ...done in 0.00 s array([[ 0. , 0. , 0. ], [ 1. , 0. , 0. ], [ 0. , 1. , 0. ], [ 0. , 0. , 1. ], [ 0.5, 0. , 0. ], [ 0.5, 0.5, 0. ], [ 0. , 0.5, 0. ], [ 0. , 0. , 0.5], [ 0.5, 0. , 0.5], [ 0. , 0.5, 0.5]]) array([[0, 4, 6, 7], [4, 1, 5, 8], [6, 5, 2, 9], [7, 8, 9, 3], [4, 6, 7, 8], [4, 6, 8, 5], [6, 7, 8, 9], [6, 5, 9, 8]]) +++ test_refine_3_4: ok sfepy: reading mesh [line2, tri3, quad4, tetra4, hexa8] (c:\Python27\lib\site-packages\sfepy/meshes/elements/3_8_1.mesh)... sfepy: ...done in 0.00 s array([[ 0. , 0. , 0. ], [ 1. , 0. , 0. ], [ 1. , 1. , 0. ], [ 0. , 1. , 0. ], [ 0. , 0. , 1. ], [ 1. , 0. , 1. ], [ 1. , 1. , 1. ], [ 0. , 1. , 1. ], [ 0.5, 0. , 0. ], [ 1. , 0.5, 0. ], [ 0.5, 1. , 0. ], [ 0. , 0.5, 0. ], [ 0.5, 0. , 1. ], [ 1. , 0.5, 1. ], [ 0.5, 1. , 1. ], [ 0. , 0.5, 1. ], [ 0. , 0. , 0.5], [ 1. , 0. , 0.5], [ 1. , 1. , 0.5], [ 0. , 1. , 0.5], [ 0.5, 0.5, 0. ], [ 0. , 0.5, 0.5], [ 0.5, 0. , 0.5], [ 0.5, 0.5, 1. ], [ 1. , 0.5, 0.5], [ 0.5, 1. , 0.5], [ 0.5, 0.5, 0.5]]) array([[ 0, 8, 20, 11, 16, 22, 26, 21], [ 1, 9, 20, 8, 17, 24, 26, 22], [ 2, 10, 20, 9, 18, 25, 26, 24], [ 3, 11, 20, 10, 19, 21, 26, 25], [ 4, 15, 23, 12, 16, 21, 26, 22], [ 5, 12, 23, 13, 17, 22, 26, 24], [ 6, 13, 23, 14, 18, 24, 26, 25], [ 7, 14, 23, 15, 19, 25, 26, 21]]) +++ test_refine_3_8: ok sfepy: guessing abaqus sfepy: reading mesh [tri3, quad4, tetra4, hexa8] (c:\Python27\lib\site-packages\sfepy/meshes/various_formats/abaqus_hex.inp)... sfepy: ...done in 0.00 s +++ test_refine_hexa: ok +++ test_refine_tetra: ok >>> all passed in 0.38 s <<< c:\Python27\lib\site-packages\sfepy\tests\test_ebcs.py sfepy: left over: ['verbose', '__builtins__', '__doc__', '__name__', 'TestCommon', 'data_dir', 'match_y_line', '__package__', '_filename', '__file__', 'op'] sfepy: reading mesh [line2, tri3, quad4, tetra4, hexa8] (c:\Python27\lib\site-packages\sfepy\meshes\2d\square_unit_tri.mesh)... sfepy: ...done in 0.00 s sfepy: creating regions... sfepy: Right sfepy: LeftStrip sfepy: Omega sfepy: RightStrip sfepy: Left sfepy: RightFix sfepy: LeftFix sfepy: ...done in 0.02 s sfepy: equation "eq": sfepy: dw_diffusion.2.Omega( m.K, q, p ) = 0 sfepy: using solvers: ts: no ts nls: newton ls: ls sfepy: updating variables... sfepy: ...done sfepy: setting up dof connectivities... sfepy: ...done in 0.00 s sfepy: matrix shape: (460L, 460L) sfepy: assembling matrix graph... sfepy: ...done in 0.00 s sfepy: matrix structural nonzeros: 4938 (2.33e-02% fill) sfepy: updating materials... sfepy: m sfepy: ...done in 0.00 s sfepy: nls: iter: 0, residual: 8.615423e+01 (rel: 1.000000e+00) sfepy: rezidual: 0.00 [s] sfepy: solve: 0.00 [s] sfepy: matrix: 0.00 [s] sfepy: nls: iter: 1, residual: 1.031171e-13 (rel: 1.196890e-15) >>> test instance prepared (1 test(s)) sfepy: saving ebc... sfepy: ...done sfepy: saving ebc... sfepy: ...done +++ test_save_ebc: ok >>> all passed in 0.03 s <<< c:\Python27\lib\site-packages\sfepy\tests\test_eigenvalue_solvers.py sfepy: left over: ['verbose', 'Solver', '__builtins__', 'dims', 'shape', '__file__', '__name__', 'TestCommon', 'nm', '_filename', 'UserMeshIO', 'gen_block_mesh', '__package__', 'eigs_expected', 'time', '__doc__', 'mesh_hook'] sfepy: reading mesh [user] (function:mesh_hook)... sfepy: ...done in 0.00 s sfepy: creating regions... sfepy: Omega sfepy: Gamma sfepy: ...done in 0.00 s sfepy: equation "eq": sfepy: dw_laplace.i.Omega(s, t) = 0 sfepy: updating variables... sfepy: ...done sfepy: setting up dof connectivities... sfepy: ...done in 0.00 s sfepy: matrix shape: (361L, 361L) sfepy: assembling matrix graph... sfepy: ...done in 0.00 s sfepy: matrix structural nonzeros: 3025 (2.32e-02% fill) >>> test instance prepared (1 test(s)) ... evp0 ... [ 0.04904454+0.j 0.12170685+0.j 0.12170685+0.j 0.19257998+0.j 0.24082108+0.j] ... evp0h ... [ 0.04904454 0.12170685 0.12170685 0.19257998 0.24082108] ... evp1 ... [ 0.04904454 0.12170685 0.12170685 0.19257998 0.24082108] ... evp2 ... [ 0.04904454 0.12170685 0.12170685 0.19257998 0.24082108] ... evp3 ... [ 0.04904454 0.12170685 0.12170685 0.19257998 0.24082108] ... solution times: ... 0.01 [s] : evp3 eig.pysparse ... 0.01 [s] : evp0h eig.scipy ... 0.02 [s] : evp0 eig.scipy ... 0.04 [s] : evp1 eig.sgscipy ... 0.24 [s] : evp2 eig.scipy_lobpcg +++ test_eigenvalue_solvers: ok >>> all passed in 0.34 s <<< c:\Python27\lib\site-packages\sfepy\tests\test_elasticity_small_strain.py sfepy: left over: ['verbose', '__builtins__', '__file__', 'all_your_bases', 'equations_general', '__name__', 'equations_iso', 'get_pars', 'TestCommon', 'data_dir', 'name', '_filename', '__package__', '__doc__'] >>> test instance prepared (2 test(s)) [('test_get_solution', ), ('test_linear_terms', )] ... mesh: c:\Python27\lib\site-packages\sfepy/meshes/3d/cube_medium_tetra.mesh, base: 1 ... isotropic sfepy: reading mesh [line2, tri3, quad4, tetra4, hexa8] (c:\Python27\lib\site-packages\sfepy\meshes\3d\cube_medium_tetra.mesh)... sfepy: ...done in 0.01 s sfepy: creating regions... sfepy: Top sfepy: Omega sfepy: Bottom sfepy: ...done in 0.01 s sfepy: equation "balance_of_forces": sfepy: dw_lin_elastic_iso.i.Omega( solid.lam, solid.mu, v, u ) = dw_point_lspring.i.Bottom( spring.pars, v, u ) sfepy: using solvers: ts: no ts nls: newton ls: ls sfepy: updating variables... sfepy: ...done sfepy: setting up dof connectivities... sfepy: ...done in 0.00 s sfepy: matrix shape: (1287L, 1287L) sfepy: assembling matrix graph... sfepy: ...done in 0.02 s sfepy: matrix structural nonzeros: 45993 (2.78e-02% fill) sfepy: updating materials... sfepy: spring sfepy: solid sfepy: ...done in 0.00 s sfepy: nls: iter: 0, residual: 1.389999e+00 (rel: 1.000000e+00) sfepy: rezidual: 0.00 [s] sfepy: solve: 0.04 [s] sfepy: matrix: 0.02 [s] sfepy: nls: iter: 1, residual: 1.832385e-15 (rel: 1.318264e-15) ... converged: True ... general sfepy: reading mesh [line2, tri3, quad4, tetra4, hexa8] (c:\Python27\lib\site-packages\sfepy\meshes\3d\cube_medium_tetra.mesh)... sfepy: ...done in 0.01 s sfepy: creating regions... sfepy: Top sfepy: Omega sfepy: Bottom sfepy: ...done in 0.02 s sfepy: equation "balance_of_forces": sfepy: dw_lin_elastic.i.Omega( solid.Dijkl, v, u ) = dw_point_lspring.i.Bottom( spring.pars, v, u ) sfepy: using solvers: ts: no ts nls: newton ls: ls sfepy: updating variables... sfepy: ...done sfepy: setting up dof connectivities... sfepy: ...done in 0.00 s sfepy: matrix shape: (1287L, 1287L) sfepy: assembling matrix graph... sfepy: ...done in 0.02 s sfepy: matrix structural nonzeros: 45993 (2.78e-02% fill) sfepy: updating materials... sfepy: spring sfepy: solid sfepy: ...done in 0.00 s sfepy: nls: iter: 0, residual: 1.389999e+00 (rel: 1.000000e+00) sfepy: rezidual: 0.01 [s] sfepy: solve: 0.05 [s] sfepy: matrix: 0.02 [s] sfepy: nls: iter: 1, residual: 1.831943e-15 (rel: 1.317946e-15) ... converged: True ... mesh: c:\Python27\lib\site-packages\sfepy/meshes/3d/cube_medium_tetra.mesh, base: 2 ... isotropic sfepy: reading mesh [line2, tri3, quad4, tetra4, hexa8] (c:\Python27\lib\site-packages\sfepy\meshes\3d\cube_medium_tetra.mesh)... sfepy: ...done in 0.01 s sfepy: creating regions... sfepy: Top sfepy: Omega sfepy: Bottom sfepy: ...done in 0.01 s sfepy: equation "balance_of_forces": sfepy: dw_lin_elastic_iso.i.Omega( solid.lam, solid.mu, v, u ) = dw_point_lspring.i.Bottom( spring.pars, v, u ) sfepy: using solvers: ts: no ts nls: newton ls: ls sfepy: updating variables... sfepy: ...done sfepy: setting up dof connectivities... sfepy: ...done in 0.00 s sfepy: matrix shape: (8622L, 8622L) sfepy: assembling matrix graph... sfepy: ...done in 0.17 s sfepy: matrix structural nonzeros: 647448 (8.71e-03% fill) sfepy: updating materials... sfepy: solid sfepy: spring sfepy: ...done in 0.00 s sfepy: nls: iter: 0, residual: 1.924264e+00 (rel: 1.000000e+00) sfepy: rezidual: 0.01 [s] sfepy: solve: 3.10 [s] sfepy: matrix: 0.16 [s] sfepy: nls: iter: 1, residual: 4.714560e-15 (rel: 2.450058e-15) ... converged: True ... general sfepy: reading mesh [line2, tri3, quad4, tetra4, hexa8] (c:\Python27\lib\site-packages\sfepy\meshes\3d\cube_medium_tetra.mesh)... sfepy: ...done in 0.01 s sfepy: creating regions... sfepy: Top sfepy: Omega sfepy: Bottom sfepy: ...done in 0.01 s sfepy: equation "balance_of_forces": sfepy: dw_lin_elastic.i.Omega( solid.Dijkl, v, u ) = dw_point_lspring.i.Bottom( spring.pars, v, u ) sfepy: using solvers: ts: no ts nls: newton ls: ls sfepy: updating variables... sfepy: ...done sfepy: setting up dof connectivities... sfepy: ...done in 0.00 s sfepy: matrix shape: (8622L, 8622L) sfepy: assembling matrix graph... sfepy: ...done in 0.17 s sfepy: matrix structural nonzeros: 647448 (8.71e-03% fill) sfepy: updating materials... sfepy: spring sfepy: solid sfepy: ...done in 0.00 s sfepy: nls: iter: 0, residual: 1.924264e+00 (rel: 1.000000e+00) sfepy: rezidual: 0.01 [s] sfepy: solve: 3.11 [s] sfepy: matrix: 0.18 [s] sfepy: nls: iter: 1, residual: 4.715117e-15 (rel: 2.450347e-15) ... converged: True ... mesh: c:\Python27\lib\site-packages\sfepy/meshes/3d/cube_medium_hexa.mesh, base: 1 ... isotropic sfepy: reading mesh [line2, tri3, quad4, tetra4, hexa8] (c:\Python27\lib\site-packages\sfepy\meshes\3d\cube_medium_hexa.mesh)... sfepy: ...done in 0.02 s sfepy: creating regions... sfepy: Top sfepy: Omega sfepy: Bottom sfepy: ...done in 0.01 s sfepy: equation "balance_of_forces": sfepy: dw_lin_elastic_iso.i.Omega( solid.lam, solid.mu, v, u ) = dw_point_lspring.i.Bottom( spring.pars, v, u ) sfepy: using solvers: ts: no ts nls: newton ls: ls sfepy: updating variables... sfepy: ...done sfepy: setting up dof connectivities... sfepy: ...done in 0.00 s sfepy: matrix shape: (3872L, 3872L) sfepy: assembling matrix graph... sfepy: ...done in 0.05 s sfepy: matrix structural nonzeros: 257548 (1.72e-02% fill) sfepy: updating materials... sfepy: solid sfepy: spring sfepy: ...done in 0.01 s sfepy: nls: iter: 0, residual: 1.220492e+00 (rel: 1.000000e+00) sfepy: rezidual: 0.01 [s] sfepy: solve: 0.69 [s] sfepy: matrix: 0.06 [s] sfepy: nls: iter: 1, residual: 1.810122e-15 (rel: 1.483109e-15) ... converged: True ... general sfepy: reading mesh [line2, tri3, quad4, tetra4, hexa8] (c:\Python27\lib\site-packages\sfepy\meshes\3d\cube_medium_hexa.mesh)... sfepy: ...done in 0.02 s sfepy: creating regions... sfepy: Top sfepy: Omega sfepy: Bottom sfepy: ...done in 0.01 s sfepy: equation "balance_of_forces": sfepy: dw_lin_elastic.i.Omega( solid.Dijkl, v, u ) = dw_point_lspring.i.Bottom( spring.pars, v, u ) sfepy: using solvers: ts: no ts nls: newton ls: ls sfepy: updating variables... sfepy: ...done sfepy: setting up dof connectivities... sfepy: ...done in 0.00 s sfepy: matrix shape: (3872L, 3872L) sfepy: assembling matrix graph... sfepy: ...done in 0.04 s sfepy: matrix structural nonzeros: 257548 (1.72e-02% fill) sfepy: updating materials... sfepy: solid sfepy: spring sfepy: ...done in 0.00 s sfepy: nls: iter: 0, residual: 1.220492e+00 (rel: 1.000000e+00) sfepy: rezidual: 0.01 [s] sfepy: solve: 0.65 [s] sfepy: matrix: 0.08 [s] sfepy: nls: iter: 1, residual: 1.809848e-15 (rel: 1.482885e-15) ... converged: True +++ test_get_solution: ok ... ||isotropic - general||: 0.000000e+00 ... ||isotropic - general||: 0.000000e+00 ... ||isotropic - general||: 0.000000e+00 +++ test_linear_terms: ok >>> all passed in 9.41 s <<< c:\Python27\lib\site-packages\sfepy\tests\test_fem.py sfepy: left over: ['verbose', 'nm', '__builtins__', '_filename', '__file__', '__doc__', '__package__', 'sfepy', 'test_bases', '__name__', 'TestCommon', 'op'] >>> test instance prepared (2 test(s)) ... 2_3 order 0 (n_nod: 1): True ... 2_3 order 1 (n_nod: 3): True ... 2_3 order 2 (n_nod: 6): True ... 2_3 order 3 (n_nod: 10): True ... 2_3 order 4 (n_nod: 15): True ... 2_3 order 5 (n_nod: 21): True ... 2_3 order 6 (n_nod: 28): True ... 2_3 order 7 (n_nod: 36): True ... 2_3 order 8 (n_nod: 45): True ... 2_3 order 9 (n_nod: 55): True ... 2_3 order 10 (n_nod: 66): True ... 2_4 order 0 (n_nod: 1): True ... 2_4 order 1 (n_nod: 4): True ... 2_4 order 2 (n_nod: 9): True ... 2_4 order 3 (n_nod: 16): True ... 2_4 order 4 (n_nod: 25): True ... 2_4 order 5 (n_nod: 36): True ... 2_4 order 6 (n_nod: 49): True ... 2_4 order 7 (n_nod: 64): True ... 2_4 order 8 (n_nod: 81): True ... 2_4 order 9 (n_nod: 100): True ... 2_4 order 10 (n_nod: 121): True ... 3_4 order 0 (n_nod: 1): True ... 3_4 order 1 (n_nod: 4): True ... 3_4 order 2 (n_nod: 10): True ... 3_4 order 3 (n_nod: 20): True ... 3_4 order 4 (n_nod: 35): True ... 3_4 order 5 (n_nod: 56): True ... 3_4 order 6 (n_nod: 84): True ... 3_4 order 7 (n_nod: 120): True ... 3_4 order 8 (n_nod: 165): True ... 3_4 order 9 (n_nod: 220): True ... 3_4 order 10 (n_nod: 286): True ... 3_8 order 0 (n_nod: 1): True ... 3_8 order 1 (n_nod: 8): True ... 3_8 order 2 (n_nod: 27): True ... 3_8 order 3 (n_nod: 64): True ... 3_8 order 4 (n_nod: 125): True ... 3_8 order 5 (n_nod: 216): True ... 3_8 order 6 (n_nod: 343): True ... 3_8 order 7 (n_nod: 512): True ... 3_8 order 8 (n_nod: 729): True ... 3_8 order 9 (n_nod: 1000): True ... 3_8 order 10 (n_nod: 1331): True +++ test_base_functions_delta: ok ... 2_3_P1: True ... 2_3_P1_grad: True ... 2_4_Q1: True ... 2_4_Q1_grad: True ... 3_4_P0: True ... 3_4_P0_grad: True ... 3_4_P2: True ... 3_4_P2B: True ... 3_4_P2B_grad: True ... 3_4_P2_grad: True ... 3_8_Q0: True ... 3_8_Q0_grad: True ... 3_8_Q1: True ... 3_8_Q1_grad: True +++ test_base_functions_values: ok >>> all passed in 0.86 s <<< c:\Python27\lib\site-packages\sfepy\tests\test_functions.py sfepy: left over: ['get_p_edge', 'get_circle', '__builtins__', '__file__', '__name__', 'assert_', 'get_pars', 'TestCommon', 'wx', 'data_dir', 'nm', '__package__', 'verbose', '_filename', '__doc__'] sfepy: reading mesh [line2, tri3, quad4, tetra4, hexa8] (c:\Python27\lib\site-packages\sfepy\meshes\2d\square_unit_tri.mesh)... sfepy: ...done in 0.00 s sfepy: creating regions... sfepy: Circle sfepy: Right sfepy: Omega sfepy: Left sfepy: ...done in 0.01 s sfepy: equation "e1": sfepy: dw_laplace.2.Omega( mf3.a, q, p ) = 0 sfepy: using solvers: ts: no ts nls: newton ls: ls >>> test instance prepared (3 test(s)) sfepy: equation "e1": sfepy: dw_laplace.2.Omega( mf3.a, q, p ) = 0 sfepy: updating variables... sfepy: ...done sfepy: setting up dof connectivities... sfepy: ...done in 0.00 s sfepy: matrix shape: (451L, 451L) sfepy: assembling matrix graph... sfepy: ...done in 0.00 s sfepy: matrix structural nonzeros: 4777 (2.35e-02% fill) sfepy: updating materials... sfepy: mf3 sfepy: ...done in 0.00 s sfepy: nls: iter: 0, residual: 7.391970e+01 (rel: 1.000000e+00) sfepy: rezidual: 0.00 [s] sfepy: solve: 0.00 [s] sfepy: matrix: 0.00 [s] sfepy: nls: iter: 1, residual: 8.203528e-14 (rel: 1.109789e-15) ... ||state_left - bc_left||: 0.000000e+00 ... ||state_right - bc_right||: 0.000000e+00 +++ test_ebc_functions: ok sfepy: updating materials... sfepy: mf3 sfepy: ...done in 0.00 s +++ test_material_functions: ok sfepy: saving regions... sfepy: Circle sfepy: ...done +++ test_region_functions: ok >>> all passed in 0.04 s <<< c:\Python27\lib\site-packages\sfepy\tests\test_high_level.py sfepy: left over: ['verbose', 'nm', '__builtins__', '_filename', '__file__', '__package__', '__doc__', 'fix_u_fun', 'time', '__name__', 'TestCommon', 'op'] sfepy: reading mesh [line2, tri3, quad4, tetra4, hexa8] (c:\Python27\lib\site-packages\sfepy\meshes\2d\rectangle_tri.mesh)... sfepy: ...done in 0.00 s >>> test instance prepared (4 test(s)) sfepy: updating variables... sfepy: ...done sfepy: setting up dof connectivities... sfepy: ...done in 0.00 s sfepy: matrix shape: (1815L, 1815L) sfepy: assembling matrix graph... sfepy: ...done in 0.00 s sfepy: matrix structural nonzeros: 39145 (1.19e-02% fill) sfepy: updating materials... sfepy: m sfepy: f sfepy: ...done in 0.00 s sfepy: nls: iter: 0, residual: 2.979640e+00 (rel: 1.000000e+00) sfepy: rezidual: 0.00 [s] sfepy: solve: 0.01 [s] sfepy: matrix: 0.00 [s] sfepy: nls: iter: 1, residual: 2.472617e-14 (rel: 8.298376e-15) +++ test_solving: ok +++ test_term_arithmetics: ok sfepy: setting up dof connectivities... sfepy: ...done in 0.00 s ... volume: 2000.00000000 == 2000.0 +++ test_term_evaluation: ok +++ test_variables: ok >>> all passed in 0.05 s <<< c:\Python27\lib\site-packages\sfepy\tests\test_homogenization_perfusion.py sfepy: left over: ['verbose', '__builtins__', '_filename', '__file__', '__package__', 'input_name', '__doc__', '__name__', 'TestCommon'] >>> test instance prepared (1 test(s)) ../examples/homogenization/perfusion_micro.py sfepy: left over: ['ch', 'verbose', 'cb', 'ii', 'io2', 'ebcs_gamma', 'ipm2', 'all_periodicY', 'io', 'val2', 'all_periodicYM', 'requirements', 'nm', 'val', '_filename', 'data_dir', 'ch2', '__package__', 'set_corrpis', 'param_h', 'aux', 'i_io', 'match_x_plane', '__doc__', '__builtins__', '__file__', 'get_channel', 'aux_bY', 'io12', 'get_mats', 'recovery_perf', 'fe', 'i_io2', 'ireg', 'ipm', 'matk2', 'ebcs_eta', 'matk1', 'set_corr_cc', 'geoms', 'match_y_plane', 'pb_def', '__name__', 'ion', 'set_corr_S', 'reg_io', 'eps0', 'coefs'] homogen: reading mesh [line2, tri3, quad4, tetra4, hexa8] (c:\Python27\lib\site-packages\sfepy\meshes\3d\perfusion_micro3d.mesh)... homogen: ...done in 0.08 s homogen: creating regions... homogen: Pr homogen: YM homogen: YB homogen: fixedYA homogen: Pn homogen: Pl homogen: YA homogen: Y homogen: Pf homogen: bYp homogen: bYm homogen: fixedYB homogen: bYB_3 homogen: bYB_2 homogen: PrYM homogen: bYMB homogen: PlYB homogen: PlYM homogen: PrYB homogen: PfYB homogen: PfYM homogen: PfYA homogen: bYA_1 homogen: bYA_2 homogen: bYA_3 homogen: PnYM homogen: PnYA homogen: PnYB homogen: bYB_1 homogen: bYMA homogen: PrYA homogen: PlYA homogen: bYMp homogen: bYMm homogen: bYA homogen: bYMchs homogen: bYB homogen: bYMpm homogen: bYA_n3 homogen: bYA_n2 homogen: bYA_n1 homogen: bYB_n2 homogen: YMmchs homogen: bYB_n3 homogen: bYB_n1 homogen: ...done in 0.16 s homogen: using solvers: ts: no ts nls: newton ls: ls homogen: equation "tmp": homogen: d_surface.iS.bYMp(pM) homogen: updating materials... homogen: ...done in 0.00 s homogen: equation "tmp": homogen: d_surface.iS.bYp(vol_all) homogen: updating materials... homogen: ...done in 0.00 s homogen: equation "tmp": homogen: d_surface.iS.bYB_3(pB) homogen: updating materials... homogen: ...done in 0.00 s homogen: equation "tmp": homogen: d_surface.iS.bYB_2(pB) homogen: updating materials... homogen: ...done in 0.00 s homogen: equation "tmp": homogen: d_surface.iS.bYB_1(pB) homogen: updating materials... homogen: ...done in 0.00 s homogen: equation "tmp": homogen: d_surface.iS.bYm(vol_all) homogen: updating materials... homogen: ...done in 0.00 s homogen: equation "tmp": homogen: d_surface.iS.bYA_1(pA) homogen: updating materials... homogen: ...done in 0.00 s homogen: equation "tmp": homogen: d_surface.iS.bYA_2(pA) homogen: updating materials... homogen: ...done in 0.00 s homogen: equation "tmp": homogen: d_surface.iS.bYA_3(pA) homogen: updating materials... homogen: ...done in 0.00 s homogen: equation "tmp": homogen: d_surface.iS.bYMm(pM) homogen: updating materials... homogen: ...done in 0.00 s homogen: equation "tmp": homogen: d_volume.iV.Y(vol_all) homogen: updating materials... homogen: ...done in 0.00 s homogen: volume: bYMp = 0.96 homogen: volume: bYp = 1.00 homogen: volume: one = 1.00 homogen: volume: bYB_3 = 0.01 homogen: volume: bYB_2 = 0.01 homogen: volume: bYB_1 = 0.03 homogen: volume: bYm = 1.00 homogen: volume: bYA_1 = 0.03 homogen: volume: bYA_2 = 0.01 homogen: volume: bYA_3 = 0.01 homogen: volume: bYMm = 0.96 homogen: volume: total = 1.00 he: computing S_testB_1... he: computing dependency pis_B... he: ...done he: computing dependency corrs_piB... he: equation "eq_pi": he: dw_diffusion.iV.YB(mat2B.k, qB, pB) + dw_volume_dot.iV.YB(qB, lsB) = - dw_diffusion.iV.YB(mat2B.k, qB, Pi_B) he: equation "eq_imv": he: dw_volume_dot.iV.YB(lvB, pB) = 0 he: updating variables... he: ...done he: setting up dof connectivities... he: ...done in 0.00 s he: lcbc: imvB he: dofs: total 1406, free 694, constrained 712, new 1 he: -> active 695 he: matrix shape: (1406L, 1406L) he: assembling matrix graph... he: ...done in 0.01 s he: matrix structural nonzeros: 22052 (1.12e-02% fill) he: updating materials... he: mat2B he: ...done in 0.00 s he: using solvers: ts: no ts nls: newton ls: ls he: updating materials... he: mat2B he: ...done in 0.00 s he: nls: iter: 0, residual: 6.927881e-02 (rel: 1.000000e+00) he: rezidual: 0.01 [s] he: solve: 0.00 [s] he: matrix: 0.01 [s] he: nls: iter: 1, residual: 1.681126e-16 (rel: 2.426609e-15) he: updating materials... he: mat2B he: ...done in 0.00 s he: nls: iter: 0, residual: 6.194610e-02 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.00 [s] he: matrix: 0.01 [s] he: nls: iter: 1, residual: 3.941807e-16 (rel: 6.363286e-15) he: ...done he: equation "tmp": he: 1.000000e+00 * ev_surface_integrate.iS.bYB_1(corr1_B) he: updating materials... he: ...done in 0.00 s he: ...done he: computing S_testB_2... he: equation "tmp": he: 1.000000e+00 * ev_surface_integrate.iS.bYB_2(corr1_B) he: updating materials... he: ...done in 0.00 s he: ...done he: computing S_testB_3... he: equation "tmp": he: 1.000000e+00 * ev_surface_integrate.iS.bYB_3(corr1_B) he: updating materials... he: ...done in 0.00 s he: ...done he: computing Hmp... he: computing dependency corrs_gamma_m... he: equation "eq_gamma_pm": he: dw_diffusion.iV.YM(mat2M.k, qM, pM) = 1.000000e+00 * dw_surface_integrate.iS.bYMm(qM) he: updating variables... he: ...done he: setting up dof connectivities... he: ...done in 0.00 s he: matrix shape: (1992L, 1992L) he: assembling matrix graph... he: ...done in 0.01 s he: matrix structural nonzeros: 25154 (6.34e-03% fill) he: updating materials... he: mat2M he: ...done in 0.01 s he: using solvers: ts: no ts nls: newton ls: ls he: updating materials... he: mat2M he: ...done in 0.00 s he: nls: iter: 0, residual: 6.480876e-02 (rel: 1.000000e+00) he: rezidual: 0.02 [s] he: solve: 0.06 [s] he: matrix: 0.02 [s] he: nls: iter: 1, residual: 3.815105e-16 (rel: 5.886712e-15) he: ...done he: equation "tmp": he: ev_surface_integrate.iS.bYMp(corr_M) he: updating materials... he: ...done in 0.00 s he: ...done he: computing Hpp... he: computing dependency corrs_gamma_p... he: equation "eq_gamma_pm": he: dw_diffusion.iV.YM(mat2M.k, qM, pM) = 1.000000e+00 * dw_surface_integrate.iS.bYMp(qM) he: updating variables... he: ...done he: setting up dof connectivities... he: ...done in 0.00 s he: matrix shape: (1992L, 1992L) he: assembling matrix graph... he: ...done in 0.01 s he: matrix structural nonzeros: 25154 (6.34e-03% fill) he: updating materials... he: mat2M he: ...done in 0.01 s he: using solvers: ts: no ts nls: newton ls: ls he: updating materials... he: mat2M he: ...done in 0.00 s he: nls: iter: 0, residual: 7.140562e-02 (rel: 1.000000e+00) he: rezidual: 0.01 [s] he: solve: 0.06 [s] he: matrix: 0.02 [s] he: nls: iter: 1, residual: 4.097167e-16 (rel: 5.737878e-15) he: ...done he: equation "tmp": he: ev_surface_integrate.iS.bYMp(corr_M) he: updating materials... he: ...done in 0.00 s he: ...done he: computing RB_1_3... he: computing dependency corrs_gamma_B_3... he: equation "eq_gamma": he: dw_diffusion.iV.YB(mat2B.k, qB, pB) + dw_volume_dot.iV.YB(qB, lsB) = 1.000000e+00 * dw_surface_integrate.iS.bYB_3(qB) he: equation "eq_imv": he: dw_volume_dot.iV.YB(lvB, pB) = 0 he: updating variables... he: ...done he: setting up dof connectivities... he: ...done in 0.00 s he: lcbc: imvB he: dofs: total 1406, free 694, constrained 712, new 1 he: -> active 695 he: matrix shape: (1406L, 1406L) he: assembling matrix graph... he: ...done in 0.01 s he: matrix structural nonzeros: 22052 (1.12e-02% fill) he: updating materials... he: mat2B he: ...done in 0.00 s he: using solvers: ts: no ts nls: newton ls: ls he: updating materials... he: mat2B he: ...done in 0.00 s he: nls: iter: 0, residual: 2.886740e-03 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.00 [s] he: matrix: 0.01 [s] he: nls: iter: 1, residual: 1.933140e-16 (rel: 6.696620e-14) he: ...done he: equation "tmp": he: ev_surface_integrate.iS.bYB_1(corr1_B) he: updating materials... he: ...done in 0.00 s he: ...done he: computing RB_1_2... he: computing dependency corrs_gamma_B_2... he: equation "eq_gamma": he: dw_diffusion.iV.YB(mat2B.k, qB, pB) + dw_volume_dot.iV.YB(qB, lsB) = 1.000000e+00 * dw_surface_integrate.iS.bYB_2(qB) he: equation "eq_imv": he: dw_volume_dot.iV.YB(lvB, pB) = 0 he: updating variables... he: ...done he: setting up dof connectivities... he: ...done in 0.00 s he: lcbc: imvB he: dofs: total 1406, free 694, constrained 712, new 1 he: -> active 695 he: matrix shape: (1406L, 1406L) he: assembling matrix graph... he: ...done in 0.01 s he: matrix structural nonzeros: 22052 (1.12e-02% fill) he: updating materials... he: mat2B he: ...done in 0.00 s he: using solvers: ts: no ts nls: newton ls: ls he: updating materials... he: mat2B he: ...done in 0.00 s he: nls: iter: 0, residual: 3.051590e-03 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.00 [s] he: matrix: 0.01 [s] he: nls: iter: 1, residual: 1.106256e-16 (rel: 3.625177e-14) he: ...done he: equation "tmp": he: ev_surface_integrate.iS.bYB_1(corr1_B) he: updating materials... he: ...done in 0.00 s he: ...done he: computing PA_2... he: computing dependency pis_A... he: ...done he: computing dependency corrs_piA... he: equation "eq_pi": he: dw_diffusion.iV.YA(mat2A.k, qA, pA) + dw_volume_dot.iV.YA(qA, lsA) = - dw_diffusion.iV.YA(mat2A.k, qA, Pi_A) he: equation "eq_imv": he: dw_volume_dot.iV.YA(lvA, pA) = 0 he: updating variables... he: ...done he: setting up dof connectivities... he: ...done in 0.00 s he: lcbc: imvA he: dofs: total 1230, free 606, constrained 624, new 1 he: -> active 607 he: matrix shape: (1230L, 1230L) he: assembling matrix graph... he: ...done in 0.01 s he: matrix structural nonzeros: 19076 (1.26e-02% fill) he: updating materials... he: mat2A he: ...done in 0.00 s he: using solvers: ts: no ts nls: newton ls: ls he: updating materials... he: mat2A he: ...done in 0.00 s he: nls: iter: 0, residual: 3.209525e-02 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.00 [s] he: matrix: 0.01 [s] he: nls: iter: 1, residual: 6.987242e-17 (rel: 2.177033e-15) he: updating materials... he: mat2A he: ...done in 0.00 s he: nls: iter: 0, residual: 2.840559e-02 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.00 [s] he: matrix: 0.01 [s] he: nls: iter: 1, residual: 1.195766e-16 (rel: 4.209615e-15) he: ...done he: equation "tmp": he: ev_surface_integrate.iS.bYA_2(corr1_A) he: updating materials... he: ...done in 0.00 s he: ...done he: computing PA_3... he: equation "tmp": he: ev_surface_integrate.iS.bYA_3(corr1_A) he: updating materials... he: ...done in 0.00 s he: ...done he: computing Hpm... he: equation "tmp": he: ev_surface_integrate.iS.bYMm(corr_M) he: updating materials... he: ...done in 0.00 s he: ...done he: computing S_testA_3... he: equation "tmp": he: 1.000000e+00 * ev_surface_integrate.iS.bYA_3(corr1_A) he: updating materials... he: ...done in 0.00 s he: ...done he: computing S_testA_2... he: equation "tmp": he: 1.000000e+00 * ev_surface_integrate.iS.bYA_2(corr1_A) he: updating materials... he: ...done in 0.00 s he: ...done he: computing S_testA_1... he: equation "tmp": he: 1.000000e+00 * ev_surface_integrate.iS.bYA_1(corr1_A) he: updating materials... he: ...done in 0.00 s he: ...done he: computing FpA... he: computing dependency corrs_oneA... he: ...done he: equation "tmp": he: dw_diffusion.iV.YM(mat2M.k, corr1_M, corr2_M) - 1.000000e+00 * ev_surface_integrate.iS.bYMp(corr1_M) he: updating materials... he: mat2M he: ...done in 0.01 s he: ...done he: computing RB_3_3... he: equation "tmp": he: ev_surface_integrate.iS.bYB_3(corr1_B) he: updating materials... he: ...done in 0.00 s he: ...done he: computing RB_3_2... he: equation "tmp": he: ev_surface_integrate.iS.bYB_3(corr1_B) he: updating materials... he: ...done in 0.00 s he: ...done he: computing PB_3... he: equation "tmp": he: ev_surface_integrate.iS.bYB_3(corr1_B) he: updating materials... he: ...done in 0.00 s he: ...done he: computing PB_2... he: equation "tmp": he: ev_surface_integrate.iS.bYB_2(corr1_B) he: updating materials... he: ...done in 0.00 s he: ...done he: computing PB_1... he: equation "tmp": he: ev_surface_integrate.iS.bYB_1(corr1_B) he: updating materials... he: ...done in 0.00 s he: ...done he: computing Hmm... he: equation "tmp": he: ev_surface_integrate.iS.bYMm(corr_M) he: updating materials... he: ...done in 0.00 s he: ...done he: computing SB_1... he: computing dependency corrs_gamma_B_1... he: equation "eq_gamma": he: dw_diffusion.iV.YB(mat2B.k, qB, pB) + dw_volume_dot.iV.YB(qB, lsB) = 1.000000e+00 * dw_surface_integrate.iS.bYB_1(qB) he: equation "eq_imv": he: dw_volume_dot.iV.YB(lvB, pB) = 0 he: updating variables... he: ...done he: setting up dof connectivities... he: ...done in 0.00 s he: lcbc: imvB he: dofs: total 1406, free 694, constrained 712, new 1 he: -> active 695 he: matrix shape: (1406L, 1406L) he: assembling matrix graph... he: ...done in 0.01 s he: matrix structural nonzeros: 22052 (1.12e-02% fill) he: updating materials... he: mat2B he: ...done in 0.00 s he: using solvers: ts: no ts nls: newton ls: ls he: updating materials... he: mat2B he: ...done in 0.00 s he: nls: iter: 0, residual: 9.945233e-03 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.00 [s] he: matrix: 0.01 [s] he: nls: iter: 1, residual: 2.087191e-16 (rel: 2.098685e-14) he: ...done he: equation "tmp": he: dw_diffusion.iV.YB(mat2B.k,corr1_B,corr2_B) he: updating materials... he: mat2B he: ...done in 0.00 s he: ...done he: computing FmB... he: computing dependency corrs_oneB... he: ...done he: equation "tmp": he: dw_diffusion.iV.YM(mat2M.k, corr1_M, corr2_M) - 1.000000e+00 * ev_surface_integrate.iS.bYMm(corr1_M) he: updating materials... he: mat2M he: ...done in 0.01 s he: ...done he: computing RB_2_1... he: equation "tmp": he: ev_surface_integrate.iS.bYB_2(corr1_B) he: updating materials... he: ...done in 0.00 s he: ...done he: computing RB_2_2... he: equation "tmp": he: ev_surface_integrate.iS.bYB_2(corr1_B) he: updating materials... he: ...done in 0.00 s he: ...done he: computing FmA... he: equation "tmp": he: dw_diffusion.iV.YM(mat2M.k, corr1_M, corr2_M) - 1.000000e+00 * ev_surface_integrate.iS.bYMm(corr1_M) he: updating materials... he: mat2M he: ...done in 0.01 s he: ...done he: computing EmB... he: computing dependency corrs_etaB... he: equation "eq_eta": he: dw_diffusion.iV.YM(mat2M.k, qM, pM) = 0 he: updating variables... he: ...done he: setting up dof connectivities... he: ...done in 0.00 s he: matrix shape: (1992L, 1992L) he: assembling matrix graph... he: ...done in 0.01 s he: matrix structural nonzeros: 25154 (6.34e-03% fill) he: updating materials... he: mat2M he: ...done in 0.01 s he: using solvers: ts: no ts nls: newton ls: ls he: updating materials... he: mat2M he: ...done in 0.00 s he: nls: iter: 0, residual: 1.923779e-01 (rel: 1.000000e+00) he: rezidual: 0.01 [s] he: solve: 0.07 [s] he: matrix: 0.02 [s] he: nls: iter: 1, residual: 4.688755e-16 (rel: 2.437263e-15) he: ...done he: equation "tmp": he: ev_surface_integrate.iS.bYMm(corr_M) he: updating materials... he: ...done in 0.00 s he: ...done he: computing SA_1... he: computing dependency corrs_gamma_A_1... he: equation "eq_gamma": he: dw_diffusion.iV.YA(mat2A.k, qA, pA) + dw_volume_dot.iV.YA(qA, lsA) = 1.000000e+00 * dw_surface_integrate.iS.bYA_1(qA) he: equation "eq_imv": he: dw_volume_dot.iV.YA(lvA, pA) = 0 he: updating variables... he: ...done he: setting up dof connectivities... he: ...done in 0.00 s he: lcbc: imvA he: dofs: total 1230, free 606, constrained 624, new 1 he: -> active 607 he: matrix shape: (1230L, 1230L) he: assembling matrix graph... he: ...done in 0.01 s he: matrix structural nonzeros: 19076 (1.26e-02% fill) he: updating materials... he: mat2A he: ...done in 0.00 s he: using solvers: ts: no ts nls: newton ls: ls he: updating materials... he: mat2A he: ...done in 0.00 s he: nls: iter: 0, residual: 9.295401e-03 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.00 [s] he: matrix: 0.01 [s] he: nls: iter: 1, residual: 2.358114e-16 (rel: 2.536861e-14) he: ...done he: equation "tmp": he: dw_diffusion.iV.YA(mat2A.k,corr1_A,corr2_A) he: updating materials... he: mat2A he: ...done in 0.00 s he: ...done he: computing SA_3... he: computing dependency corrs_gamma_A_3... he: equation "eq_gamma": he: dw_diffusion.iV.YA(mat2A.k, qA, pA) + dw_volume_dot.iV.YA(qA, lsA) = 1.000000e+00 * dw_surface_integrate.iS.bYA_3(qA) he: equation "eq_imv": he: dw_volume_dot.iV.YA(lvA, pA) = 0 he: updating variables... he: ...done he: setting up dof connectivities... he: ...done in 0.00 s he: lcbc: imvA he: dofs: total 1230, free 606, constrained 624, new 1 he: -> active 607 he: matrix shape: (1230L, 1230L) he: assembling matrix graph... he: ...done in 0.01 s he: matrix structural nonzeros: 19076 (1.26e-02% fill) he: updating materials... he: mat2A he: ...done in 0.00 s he: using solvers: ts: no ts nls: newton ls: ls he: updating materials... he: mat2A he: ...done in 0.00 s he: nls: iter: 0, residual: 2.886740e-03 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.00 [s] he: matrix: 0.01 [s] he: nls: iter: 1, residual: 1.172788e-16 (rel: 4.062671e-14) he: ...done he: equation "tmp": he: dw_diffusion.iV.YA(mat2A.k,corr1_A,corr2_A) he: updating materials... he: mat2A he: ...done in 0.00 s he: ...done he: computing SA_2... he: computing dependency corrs_gamma_A_2... he: equation "eq_gamma": he: dw_diffusion.iV.YA(mat2A.k, qA, pA) + dw_volume_dot.iV.YA(qA, lsA) = 1.000000e+00 * dw_surface_integrate.iS.bYA_2(qA) he: equation "eq_imv": he: dw_volume_dot.iV.YA(lvA, pA) = 0 he: updating variables... he: ...done he: setting up dof connectivities... he: ...done in 0.00 s he: lcbc: imvA he: dofs: total 1230, free 606, constrained 624, new 1 he: -> active 607 he: matrix shape: (1230L, 1230L) he: assembling matrix graph... he: ...done in 0.01 s he: matrix structural nonzeros: 19076 (1.26e-02% fill) he: updating materials... he: mat2A he: ...done in 0.00 s he: using solvers: ts: no ts nls: newton ls: ls he: updating materials... he: mat2A he: ...done in 0.00 s he: nls: iter: 0, residual: 3.258425e-03 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.00 [s] he: matrix: 0.01 [s] he: nls: iter: 1, residual: 6.348575e-17 (rel: 1.948357e-14) he: ...done he: equation "tmp": he: dw_diffusion.iV.YA(mat2A.k,corr1_A,corr2_A) he: updating materials... he: mat2A he: ...done in 0.00 s he: ...done he: computing RA_2_3... he: equation "tmp": he: ev_surface_integrate.iS.bYA_2(corr1_A) he: updating materials... he: ...done in 0.00 s he: ...done he: computing RA_2_2... he: equation "tmp": he: ev_surface_integrate.iS.bYA_2(corr1_A) he: updating materials... he: ...done in 0.00 s he: ...done he: computing GA... he: computing dependency corrs_etaA... he: equation "eq_eta": he: dw_diffusion.iV.YM(mat2M.k, qM, pM) = 0 he: updating variables... he: ...done he: setting up dof connectivities... he: ...done in 0.00 s he: matrix shape: (1992L, 1992L) he: assembling matrix graph... he: ...done in 0.01 s he: matrix structural nonzeros: 25154 (6.34e-03% fill) he: updating materials... he: mat2M he: ...done in 0.01 s he: using solvers: ts: no ts nls: newton ls: ls he: updating materials... he: mat2M he: ...done in 0.00 s he: nls: iter: 0, residual: 1.862908e-01 (rel: 1.000000e+00) he: rezidual: 0.01 [s] he: solve: 0.07 [s] he: matrix: 0.02 [s] he: nls: iter: 1, residual: 4.329536e-16 (rel: 2.324074e-15) he: ...done he: equation "tmp": he: dw_diffusion.iV.YM(mat2M.k, corr1_M, corr2_M) he: updating materials... he: mat2M he: ...done in 0.01 s he: ...done he: computing RB_3_1... he: equation "tmp": he: ev_surface_integrate.iS.bYB_3(corr1_B) he: updating materials... he: ...done in 0.00 s he: ...done he: computing KB... he: equation "tmp": he: dw_diffusion.iV.YB(mat2B.k, corr1_B, corr2_B) he: updating materials... he: mat2B he: ...done in 0.00 s he: ...done he: computing KA... he: equation "tmp": he: dw_diffusion.iV.YA(mat2A.k, corr1_A, corr2_A) he: updating materials... he: mat2A he: ...done in 0.00 s he: ...done he: computing EmA... he: equation "tmp": he: ev_surface_integrate.iS.bYMm(corr_M) he: updating materials... he: ...done in 0.00 s he: ...done he: computing PA_1... he: equation "tmp": he: ev_surface_integrate.iS.bYA_1(corr1_A) he: updating materials... he: ...done in 0.00 s he: ...done he: computing EpB... he: equation "tmp": he: ev_surface_integrate.iS.bYMp(corr_M) he: updating materials... he: ...done in 0.00 s he: ...done he: computing FpB... he: equation "tmp": he: dw_diffusion.iV.YM(mat2M.k, corr1_M, corr2_M) - 1.000000e+00 * ev_surface_integrate.iS.bYMp(corr1_M) he: updating materials... he: mat2M he: ...done in 0.01 s he: ...done he: computing RA_3_2... he: equation "tmp": he: ev_surface_integrate.iS.bYA_3(corr1_A) he: updating materials... he: ...done in 0.00 s he: ...done he: computing RA_3_3... he: equation "tmp": he: ev_surface_integrate.iS.bYA_3(corr1_A) he: updating materials... he: ...done in 0.00 s he: ...done he: computing RA_3_1... he: equation "tmp": he: ev_surface_integrate.iS.bYA_3(corr1_A) he: updating materials... he: ...done in 0.00 s he: ...done he: computing RA_1_1... he: equation "tmp": he: ev_surface_integrate.iS.bYA_1(corr1_A) he: updating materials... he: ...done in 0.00 s he: ...done he: computing RA_1_2... he: equation "tmp": he: ev_surface_integrate.iS.bYA_1(corr1_A) he: updating materials... he: ...done in 0.00 s he: ...done he: computing RA_1_3... he: equation "tmp": he: ev_surface_integrate.iS.bYA_1(corr1_A) he: updating materials... he: ...done in 0.00 s he: ...done he: computing RA_2_1... he: equation "tmp": he: ev_surface_integrate.iS.bYA_2(corr1_A) he: updating materials... he: ...done in 0.00 s he: ...done he: computing RB_1_1... he: equation "tmp": he: ev_surface_integrate.iS.bYB_1(corr1_B) he: updating materials... he: ...done in 0.00 s he: ...done he: computing vol_bYMpm... he: equation "tmp": he: d_surface.iS.bYMp(pM) he: updating materials... he: ...done in 0.00 s he: equation "tmp": he: d_surface.iS.bYMm(pM) he: updating materials... he: ...done in 0.00 s he: ...done he: computing RB_2_3... he: equation "tmp": he: ev_surface_integrate.iS.bYB_2(corr1_B) he: updating materials... he: ...done in 0.00 s he: ...done he: computing SB_2... he: equation "tmp": he: dw_diffusion.iV.YB(mat2B.k,corr1_B,corr2_B) he: updating materials... he: mat2B he: ...done in 0.00 s he: ...done he: computing SB_3... he: equation "tmp": he: dw_diffusion.iV.YB(mat2B.k,corr1_B,corr2_B) he: updating materials... he: mat2B he: ...done in 0.00 s he: ...done he: computing GB... he: equation "tmp": he: dw_diffusion.iV.YM(mat2M.k, corr1_M, corr2_M) he: updating materials... he: mat2M he: ...done in 0.01 s he: ...done he: computing EpA... he: equation "tmp": he: ev_surface_integrate.iS.bYMp(corr_M) he: updating materials... he: ...done in 0.00 s he: ...done ... allowed error: abs = 1.000000e-09 ... |G^A - G^B|: 2.220446e-16 ... |F^{A+} + F^{B+} - -bYM^+|: 2.220446e-16 ... |F^{A-} + F^{B-} - -bYM^-|: 3.330669e-16 ... |H^{+-} - H^{-+}|: 2.602085e-18 ... |E^{A-} - -F^{A-}|: 1.110223e-16 ... |E^{A+} - -F^{A+}|: 1.665335e-16 ... |E^{B-} - -F^{B-}|: 1.110223e-16 ... |E^{B+} - -F^{B+}|: 1.110223e-16 ... ||S_test - S||: 4.165369e-18 ... ||S_test - S||: 4.661124e-18 ... ||S_test - S||: 2.174494e-18 ... ||S_test - S||: 2.607162e-18 ... ||S_test - S||: 4.336809e-18 ... ||S_test - S||: 3.576224e-18 +++ test_solution: ok >>> all passed in 7.36 s <<< c:\Python27\lib\site-packages\sfepy\tests\test_hyperelastic_tlul.py he: left over: ['input_names', 'verbose', '__builtins__', '_filename', '__file__', '__package__', '__doc__', '__name__', 'TestCommon', 'output_name_trunk', 'NLSStatus'] >>> test instance prepared (1 test(s)) he: left over: ['verbose', 'centre', '__builtins__', '__doc__', '__name__', 'data_dir', 'nm', 'stress_strain', '__package__', 'rotate_yz', '_filename', '__file__'] ... hyperelastic formulation: UL he: reading mesh [line2, tri3, quad4, tetra4, hexa8] (c:\Python27\lib\site-packages\sfepy\meshes\3d\cylinder.mesh)... he: ...done in 0.02 s he: creating regions... he: Right he: Omega he: Left he: ...done in 0.01 s he: equation "balance": he: dw_ul_he_neohook.3.Omega( solid.mu, v, u ) + dw_ul_he_mooney_rivlin.3.Omega(solid.kappa, v, u) + dw_ul_bulk_penalty.3.Omega( solid.K, v, u ) = 0 he: using solvers: ts: ts nls: newton ls: ls he: updating variables... he: ...done angle: 0.0 he: setting up dof connectivities... he: ...done in 0.00 s he: matrix shape: (900L, 900L) he: assembling matrix graph... he: ...done in 0.01 s he: matrix structural nonzeros: 31842 (3.93e-02% fill) he: updating materials... he: solid he: ...done in 0.00 s he: ====== time 0.000000e+00 (step 1 of 11) ===== he: initial residual: 0.000000e+00 he: equation "tmp": he: dw_ul_he_neohook.3.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_ul_he_neohook.3.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_ul_he_mooney_rivlin.3.Omega( solid.kappa, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_ul_bulk_penalty.3.Omega( solid.K, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: ====== time 1.000000e-01 (step 2 of 11) ===== he: updating variables... he: ...done angle: 10.0 he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 1.541089e-02 (rel: 1.000000e+00) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 1, residual: 8.356770e-03 (rel: 5.422639e-01) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 2, residual: 5.223665e-04 (rel: 3.389593e-02) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 3, residual: 1.447311e-04 (rel: 9.391479e-03) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 4, residual: 1.698668e-05 (rel: 1.102252e-03) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 5, residual: 3.602796e-06 (rel: 2.337825e-04) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 6, residual: 8.222021e-07 (rel: 5.335201e-05) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 7, residual: 1.883167e-07 (rel: 1.221971e-05) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 8, residual: 4.324542e-08 (rel: 2.806160e-06) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 9, residual: 9.969490e-09 (rel: 6.469119e-07) he: equation "tmp": he: dw_ul_he_neohook.3.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_ul_he_neohook.3.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_ul_he_mooney_rivlin.3.Omega( solid.kappa, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_ul_bulk_penalty.3.Omega( solid.K, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: ====== time 2.000000e-01 (step 3 of 11) ===== he: updating variables... he: ...done angle: 20.0 he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 1.551383e-02 (rel: 1.000000e+00) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 1, residual: 8.253864e-03 (rel: 5.320327e-01) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 2, residual: 5.348966e-04 (rel: 3.447870e-02) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 3, residual: 1.510784e-04 (rel: 9.738304e-03) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 4, residual: 1.798043e-05 (rel: 1.158994e-03) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 5, residual: 3.844827e-06 (rel: 2.478323e-04) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 6, residual: 8.838313e-07 (rel: 5.697055e-05) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 7, residual: 2.038750e-07 (rel: 1.314150e-05) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 8, residual: 4.715045e-08 (rel: 3.039253e-06) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 9, residual: 1.094941e-08 (rel: 7.057840e-07) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 10, residual: 2.554731e-09 (rel: 1.646745e-07) he: equation "tmp": he: dw_ul_he_neohook.3.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_ul_he_neohook.3.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_ul_he_mooney_rivlin.3.Omega( solid.kappa, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_ul_bulk_penalty.3.Omega( solid.K, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: ====== time 3.000000e-01 (step 4 of 11) ===== he: updating variables... he: ...done angle: 30.0 he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 1.562693e-02 (rel: 1.000000e+00) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 1, residual: 8.142631e-03 (rel: 5.210640e-01) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 2, residual: 5.486954e-04 (rel: 3.511217e-02) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 3, residual: 1.577994e-04 (rel: 1.009792e-02) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 4, residual: 1.914623e-05 (rel: 1.225208e-03) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 5, residual: 4.145618e-06 (rel: 2.652868e-04) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 6, residual: 9.629515e-07 (rel: 6.162128e-05) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 7, residual: 2.243673e-07 (rel: 1.435773e-05) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 8, residual: 5.240844e-08 (rel: 3.353726e-06) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 9, residual: 1.229468e-08 (rel: 7.867621e-07) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 10, residual: 2.899402e-09 (rel: 1.855388e-07) he: equation "tmp": he: dw_ul_he_neohook.3.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_ul_he_neohook.3.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_ul_he_mooney_rivlin.3.Omega( solid.kappa, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_ul_bulk_penalty.3.Omega( solid.K, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: ====== time 4.000000e-01 (step 5 of 11) ===== he: updating variables... he: ...done angle: 40.0 he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 1.575110e-02 (rel: 1.000000e+00) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 1, residual: 8.024086e-03 (rel: 5.094301e-01) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 2, residual: 5.637537e-04 (rel: 3.579138e-02) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 3, residual: 1.648445e-04 (rel: 1.046558e-02) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 4, residual: 2.040012e-05 (rel: 1.295155e-03) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.06 [s] he: nls: iter: 5, residual: 4.479384e-06 (rel: 2.843855e-04) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 6, residual: 1.053415e-06 (rel: 6.687882e-05) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 7, residual: 2.485122e-07 (rel: 1.577745e-05) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 8, residual: 5.879102e-08 (rel: 3.732502e-06) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 9, residual: 1.397572e-08 (rel: 8.872851e-07) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 10, residual: 3.342260e-09 (rel: 2.121921e-07) he: equation "tmp": he: dw_ul_he_neohook.3.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_ul_he_neohook.3.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_ul_he_mooney_rivlin.3.Omega( solid.kappa, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_ul_bulk_penalty.3.Omega( solid.K, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: ====== time 5.000000e-01 (step 6 of 11) ===== he: updating variables... he: ...done angle: 50.0 he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 1.588378e-02 (rel: 1.000000e+00) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 1, residual: 7.909211e-03 (rel: 4.979425e-01) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 2, residual: 5.792393e-04 (rel: 3.646733e-02) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 3, residual: 1.717755e-04 (rel: 1.081452e-02) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 4, residual: 2.164833e-05 (rel: 1.362920e-03) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 5, residual: 4.822087e-06 (rel: 3.035855e-04) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 6, residual: 1.149961e-06 (rel: 7.239844e-05) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 7, residual: 2.753891e-07 (rel: 1.733775e-05) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 8, residual: 6.620760e-08 (rel: 4.168251e-06) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 9, residual: 1.601317e-08 (rel: 1.008146e-06) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 10, residual: 3.901070e-09 (rel: 2.456008e-07) he: equation "tmp": he: dw_ul_he_neohook.3.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_ul_he_neohook.3.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_ul_he_mooney_rivlin.3.Omega( solid.kappa, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_ul_bulk_penalty.3.Omega( solid.K, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: ====== time 6.000000e-01 (step 7 of 11) ===== he: updating variables... he: ...done angle: 60.0 he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 1.602100e-02 (rel: 1.000000e+00) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 1, residual: 7.808018e-03 (rel: 4.873615e-01) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 2, residual: 5.943190e-04 (rel: 3.709625e-02) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 3, residual: 1.782461e-04 (rel: 1.112578e-02) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 4, residual: 2.286833e-05 (rel: 1.427397e-03) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 5, residual: 5.170431e-06 (rel: 3.227284e-04) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 6, residual: 1.251845e-06 (rel: 7.813776e-05) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 7, residual: 3.048694e-07 (rel: 1.902936e-05) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 8, residual: 7.466010e-08 (rel: 4.660141e-06) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 9, residual: 1.842284e-08 (rel: 1.149918e-06) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.06 [s] he: nls: iter: 10, residual: 4.585776e-09 (rel: 2.862353e-07) he: equation "tmp": he: dw_ul_he_neohook.3.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_ul_he_neohook.3.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_ul_he_mooney_rivlin.3.Omega( solid.kappa, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_ul_bulk_penalty.3.Omega( solid.K, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: ====== time 7.000000e-01 (step 8 of 11) ===== he: updating variables... he: ...done angle: 70.0 he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 1.616034e-02 (rel: 1.000000e+00) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.06 [s] he: nls: iter: 1, residual: 7.716670e-03 (rel: 4.775065e-01) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 2, residual: 6.093173e-04 (rel: 3.770448e-02) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 3, residual: 1.844643e-04 (rel: 1.141463e-02) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 4, residual: 2.415352e-05 (rel: 1.494616e-03) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 5, residual: 5.549971e-06 (rel: 3.434315e-04) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 6, residual: 1.364782e-06 (rel: 8.445255e-05) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 7, residual: 3.380478e-07 (rel: 2.091835e-05) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 8, residual: 8.431929e-08 (rel: 5.217666e-06) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 9, residual: 2.122174e-08 (rel: 1.313199e-06) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 10, residual: 5.395472e-09 (rel: 3.338711e-07) he: equation "tmp": he: dw_ul_he_neohook.3.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_ul_he_neohook.3.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_ul_he_mooney_rivlin.3.Omega( solid.kappa, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_ul_bulk_penalty.3.Omega( solid.K, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: ====== time 8.000000e-01 (step 9 of 11) ===== he: updating variables... he: ...done angle: 80.0 he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 1.630225e-02 (rel: 1.000000e+00) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 1, residual: 7.616956e-03 (rel: 4.672335e-01) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 2, residual: 6.259202e-04 (rel: 3.839472e-02) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 3, residual: 1.911533e-04 (rel: 1.172558e-02) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 4, residual: 2.568172e-05 (rel: 1.575349e-03) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 5, residual: 6.006897e-06 (rel: 3.684705e-04) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 6, residual: 1.500309e-06 (rel: 9.203081e-05) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 7, residual: 3.776293e-07 (rel: 2.316425e-05) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 8, residual: 9.578855e-08 (rel: 5.875788e-06) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 9, residual: 2.454094e-08 (rel: 1.505371e-06) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 10, residual: 6.359093e-09 (rel: 3.900746e-07) he: equation "tmp": he: dw_ul_he_neohook.3.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_ul_he_neohook.3.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_ul_he_mooney_rivlin.3.Omega( solid.kappa, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_ul_bulk_penalty.3.Omega( solid.K, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: ====== time 9.000000e-01 (step 10 of 11) ===== he: updating variables... he: ...done angle: 90.0 he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 1.644843e-02 (rel: 1.000000e+00) he: rezidual: 0.02 [s] he: solve: 0.01 [s] he: matrix: 0.08 [s] he: nls: iter: 1, residual: 7.491609e-03 (rel: 4.554604e-01) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.07 [s] he: nls: iter: 2, residual: 6.458434e-04 (rel: 3.926474e-02) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 3, residual: 1.989611e-04 (rel: 1.209605e-02) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 4, residual: 2.761982e-05 (rel: 1.679176e-03) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 5, residual: 6.591146e-06 (rel: 4.007157e-04) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 6, residual: 1.674297e-06 (rel: 1.017907e-04) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 7, residual: 4.285124e-07 (rel: 2.605187e-05) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 8, residual: 1.105397e-07 (rel: 6.720380e-06) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 9, residual: 2.881573e-08 (rel: 1.751883e-06) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 10, residual: 7.604657e-09 (rel: 4.623332e-07) he: equation "tmp": he: dw_ul_he_neohook.3.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_ul_he_neohook.3.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_ul_he_mooney_rivlin.3.Omega( solid.kappa, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_ul_bulk_penalty.3.Omega( solid.K, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: ====== time 1.000000e+00 (step 11 of 11) ===== he: updating variables... he: ...done angle: 100.0 he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 1.659969e-02 (rel: 1.000000e+00) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 1, residual: 7.341610e-03 (rel: 4.422739e-01) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 2, residual: 6.691951e-04 (rel: 4.031370e-02) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 3, residual: 2.078760e-04 (rel: 1.252288e-02) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 4, residual: 3.002118e-05 (rel: 1.808538e-03) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 5, residual: 7.335074e-06 (rel: 4.418801e-04) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 6, residual: 1.902750e-06 (rel: 1.146256e-04) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 7, residual: 4.973201e-07 (rel: 2.995959e-05) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 8, residual: 1.310330e-07 (rel: 7.893701e-06) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 9, residual: 3.489920e-08 (rel: 2.102400e-06) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.06 [s] he: nls: iter: 10, residual: 9.414685e-09 (rel: 5.671602e-07) he: equation "tmp": he: dw_ul_he_neohook.3.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_ul_he_neohook.3.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_ul_he_mooney_rivlin.3.Omega( solid.kappa, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_ul_bulk_penalty.3.Omega( solid.K, v, u ) he: updating materials... he: solid he: ...done in 0.00 s ... c:\Python27\lib\site-packages\sfepy\tests\../examples/large_deformation/hyperelastic_ul.py solved he: left over: ['verbose', 'centre', 'rotate_yz', '__file__', '__name__', 'data_dir', 'nm', 'stress_strain', '_filename', '__package__', '__builtins__', '__doc__'] ... hyperelastic formulation: TL he: reading mesh [line2, tri3, quad4, tetra4, hexa8] (c:\Python27\lib\site-packages\sfepy\meshes\3d\cylinder.mesh)... he: ...done in 0.01 s he: creating regions... he: Right he: Omega he: Left he: ...done in 0.01 s he: equation "balance": he: dw_tl_he_neohook.i.Omega( solid.mu, v, u ) + dw_tl_he_mooney_rivlin.i.Omega( solid.kappa, v, u ) + dw_tl_bulk_penalty.i.Omega( solid.K, v, u ) = 0 he: using solvers: ts: ts nls: newton ls: ls he: updating variables... he: ...done angle: 0.0 he: setting up dof connectivities... he: ...done in 0.00 s he: matrix shape: (900L, 900L) he: assembling matrix graph... he: ...done in 0.01 s he: matrix structural nonzeros: 31842 (3.93e-02% fill) he: updating materials... he: solid he: ...done in 0.00 s he: ====== time 0.000000e+00 (step 1 of 11) ===== he: initial residual: 0.000000e+00 he: equation "tmp": he: dw_tl_he_neohook.i.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_tl_he_neohook.i.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_tl_he_mooney_rivlin.i.Omega( solid.kappa, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_tl_bulk_penalty.i.Omega( solid.K, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: ====== time 1.000000e-01 (step 2 of 11) ===== he: updating variables... he: ...done angle: 10.0 he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 1.542762e-02 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.03 [s] he: nls: iter: 1, residual: 1.292867e-02 (rel: 8.380213e-01) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.03 [s] he: nls: iter: 2, residual: 4.055858e-04 (rel: 2.628959e-02) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.03 [s] he: nls: iter: 3, residual: 2.188132e-06 (rel: 1.418321e-04) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.03 [s] he: nls: iter: 4, residual: 7.115966e-10 (rel: 4.612486e-08) he: rezidual: 0.00 [s] he: solve: 0.02 [s] he: matrix: 0.04 [s] he: nls: iter: 5, residual: 1.766202e-16 (rel: 1.144831e-14) he: equation "tmp": he: dw_tl_he_neohook.i.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_tl_he_neohook.i.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_tl_he_mooney_rivlin.i.Omega( solid.kappa, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_tl_bulk_penalty.i.Omega( solid.K, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: ====== time 2.000000e-01 (step 3 of 11) ===== he: updating variables... he: ...done angle: 20.0 he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 1.553477e-02 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.04 [s] he: nls: iter: 1, residual: 1.280028e-02 (rel: 8.239761e-01) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.03 [s] he: nls: iter: 2, residual: 4.026136e-04 (rel: 2.591694e-02) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.03 [s] he: nls: iter: 3, residual: 2.865480e-06 (rel: 1.844560e-04) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.03 [s] he: nls: iter: 4, residual: 1.323883e-09 (rel: 8.522069e-08) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.03 [s] he: nls: iter: 5, residual: 3.377884e-16 (rel: 2.174403e-14) he: equation "tmp": he: dw_tl_he_neohook.i.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_tl_he_neohook.i.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_tl_he_mooney_rivlin.i.Omega( solid.kappa, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_tl_bulk_penalty.i.Omega( solid.K, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: ====== time 3.000000e-01 (step 4 of 11) ===== he: updating variables... he: ...done angle: 30.0 he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 1.564786e-02 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.03 [s] he: nls: iter: 1, residual: 1.267511e-02 (rel: 8.100221e-01) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.03 [s] he: nls: iter: 2, residual: 4.001295e-04 (rel: 2.557088e-02) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.03 [s] he: nls: iter: 3, residual: 3.844027e-06 (rel: 2.456583e-04) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.03 [s] he: nls: iter: 4, residual: 2.074298e-09 (rel: 1.325611e-07) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.03 [s] he: nls: iter: 5, residual: 1.160738e-15 (rel: 7.417867e-14) he: equation "tmp": he: dw_tl_he_neohook.i.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_tl_he_neohook.i.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_tl_he_mooney_rivlin.i.Omega( solid.kappa, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_tl_bulk_penalty.i.Omega( solid.K, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: ====== time 4.000000e-01 (step 5 of 11) ===== he: updating variables... he: ...done angle: 40.0 he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 1.576732e-02 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.03 [s] he: nls: iter: 1, residual: 1.255322e-02 (rel: 7.961547e-01) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.03 [s] he: nls: iter: 2, residual: 3.981146e-04 (rel: 2.524935e-02) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.03 [s] he: nls: iter: 3, residual: 5.074043e-06 (rel: 3.218076e-04) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.03 [s] he: nls: iter: 4, residual: 2.931853e-09 (rel: 1.859449e-07) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.04 [s] he: nls: iter: 5, residual: 3.368796e-15 (rel: 2.136569e-13) he: equation "tmp": he: dw_tl_he_neohook.i.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_tl_he_neohook.i.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_tl_he_mooney_rivlin.i.Omega( solid.kappa, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_tl_bulk_penalty.i.Omega( solid.K, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: ====== time 5.000000e-01 (step 6 of 11) ===== he: updating variables... he: ...done angle: 50.0 he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 1.589360e-02 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.03 [s] he: nls: iter: 1, residual: 1.243465e-02 (rel: 7.823683e-01) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.03 [s] he: nls: iter: 2, residual: 3.965518e-04 (rel: 2.495041e-02) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.03 [s] he: nls: iter: 3, residual: 6.510580e-06 (rel: 4.096354e-04) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.03 [s] he: nls: iter: 4, residual: 3.869756e-09 (rel: 2.434789e-07) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.03 [s] he: nls: iter: 5, residual: 8.011953e-15 (rel: 5.040994e-13) he: equation "tmp": he: dw_tl_he_neohook.i.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_tl_he_neohook.i.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_tl_he_mooney_rivlin.i.Omega( solid.kappa, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_tl_bulk_penalty.i.Omega( solid.K, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: ====== time 6.000000e-01 (step 7 of 11) ===== he: updating variables... he: ...done angle: 60.0 he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 1.602719e-02 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.03 [s] he: nls: iter: 1, residual: 1.231941e-02 (rel: 7.686568e-01) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.03 [s] he: nls: iter: 2, residual: 3.954264e-04 (rel: 2.467222e-02) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.03 [s] he: nls: iter: 3, residual: 8.117581e-06 (rel: 5.064882e-04) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.03 [s] he: nls: iter: 4, residual: 4.865350e-09 (rel: 3.035685e-07) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.03 [s] he: nls: iter: 5, residual: 1.645089e-14 (rel: 1.026437e-12) he: equation "tmp": he: dw_tl_he_neohook.i.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_tl_he_neohook.i.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_tl_he_mooney_rivlin.i.Omega( solid.kappa, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_tl_bulk_penalty.i.Omega( solid.K, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: ====== time 7.000000e-01 (step 8 of 11) ===== he: updating variables... he: ...done angle: 70.0 he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 1.616861e-02 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.04 [s] he: nls: iter: 1, residual: 1.220754e-02 (rel: 7.550144e-01) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.03 [s] he: nls: iter: 2, residual: 3.947255e-04 (rel: 2.441307e-02) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.04 [s] he: nls: iter: 3, residual: 9.865673e-06 (rel: 6.101743e-04) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.03 [s] he: nls: iter: 4, residual: 5.899979e-09 (rel: 3.649032e-07) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.03 [s] he: nls: iter: 5, residual: 3.027046e-14 (rel: 1.872174e-12) he: equation "tmp": he: dw_tl_he_neohook.i.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_tl_he_neohook.i.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_tl_he_mooney_rivlin.i.Omega( solid.kappa, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_tl_bulk_penalty.i.Omega( solid.K, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: ====== time 8.000000e-01 (step 9 of 11) ===== he: updating variables... he: ...done angle: 80.0 he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 1.631842e-02 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.03 [s] he: nls: iter: 1, residual: 1.209908e-02 (rel: 7.414365e-01) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.03 [s] he: nls: iter: 2, residual: 3.944384e-04 (rel: 2.417136e-02) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.04 [s] he: nls: iter: 3, residual: 1.173037e-05 (rel: 7.188418e-04) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.03 [s] he: nls: iter: 4, residual: 6.958837e-09 (rel: 4.264405e-07) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.03 [s] he: nls: iter: 5, residual: 5.101730e-14 (rel: 3.126362e-12) he: equation "tmp": he: dw_tl_he_neohook.i.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_tl_he_neohook.i.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_tl_he_mooney_rivlin.i.Omega( solid.kappa, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_tl_bulk_penalty.i.Omega( solid.K, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: ====== time 9.000000e-01 (step 10 of 11) ===== he: updating variables... he: ...done angle: 90.0 he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 1.647719e-02 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.03 [s] he: nls: iter: 1, residual: 1.199407e-02 (rel: 7.279196e-01) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.03 [s] he: nls: iter: 2, residual: 3.945565e-04 (rel: 2.394561e-02) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.03 [s] he: nls: iter: 3, residual: 1.369112e-05 (rel: 8.309134e-04) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.03 [s] he: nls: iter: 4, residual: 8.030824e-09 (rel: 4.873903e-07) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.03 [s] he: nls: iter: 5, residual: 8.030043e-14 (rel: 4.873429e-12) he: equation "tmp": he: dw_tl_he_neohook.i.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_tl_he_neohook.i.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_tl_he_mooney_rivlin.i.Omega( solid.kappa, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_tl_bulk_penalty.i.Omega( solid.K, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: ====== time 1.000000e+00 (step 11 of 11) ===== he: updating variables... he: ...done angle: 100.0 he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 1.664552e-02 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.03 [s] he: nls: iter: 1, residual: 1.189259e-02 (rel: 7.144615e-01) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.03 [s] he: nls: iter: 2, residual: 3.950731e-04 (rel: 2.373449e-02) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.03 [s] he: nls: iter: 3, residual: 1.573093e-05 (rel: 9.450547e-04) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.03 [s] he: nls: iter: 4, residual: 9.108409e-09 (rel: 5.471987e-07) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.03 [s] he: nls: iter: 5, residual: 1.194705e-13 (rel: 7.177336e-12) he: equation "tmp": he: dw_tl_he_neohook.i.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_tl_he_neohook.i.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_tl_he_mooney_rivlin.i.Omega( solid.kappa, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_tl_bulk_penalty.i.Omega( solid.K, v, u ) he: updating materials... he: solid he: ...done in 0.00 s ... c:\Python27\lib\site-packages\sfepy\tests\../examples/large_deformation/hyperelastic.py solved he: left over: ['verbose', 'centre', '__builtins__', '__doc__', '__name__', 'data_dir', 'nm', 'stress_strain', '__package__', 'rotate_yz', '_filename', '__file__'] ... hyperelastic formulation: ULM he: reading mesh [line2, tri3, quad4, tetra4, hexa8] (c:\Python27\lib\site-packages\sfepy\meshes\3d\cylinder.mesh)... he: ...done in 0.01 s he: creating regions... he: Right he: Omega he: Left he: ...done in 0.01 s he: equation "volume": he: dw_ul_volume.3.Omega( q, u ) + dw_ul_compressible.3.Omega( solid.iK, q, p, u ) = 0 he: equation "balance": he: dw_ul_he_neohook.3.Omega( solid.mu, v, u ) + dw_ul_he_mooney_rivlin.3.Omega(solid.kappa, v, u) + dw_ul_bulk_pressure.3.Omega( v, u, p ) = 0 he: using solvers: ts: ts nls: newton ls: ls he: updating variables... he: ...done angle: 0.0 he: setting up dof connectivities... he: ...done in 0.00 s he: matrix shape: (2248L, 2248L) he: assembling matrix graph... he: ...done in 0.02 s he: matrix structural nonzeros: 62902 (1.24e-02% fill) he: updating materials... he: solid he: ...done in 0.00 s he: ====== time 0.000000e+00 (step 1 of 11) ===== he: initial residual: 0.000000e+00 he: equation "tmp": he: dw_ul_he_neohook.3.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_ul_he_neohook.3.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_ul_he_mooney_rivlin.3.Omega( solid.kappa, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_ul_bulk_pressure.3.Omega( v, u, p ) he: updating materials... he: ...done in 0.00 s he: ====== time 1.000000e-01 (step 2 of 11) ===== he: updating variables... he: ...done angle: 10.0 he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 2.825495e-03 (rel: 1.000000e+00) he: rezidual: 0.01 [s] he: solve: 0.14 [s] he: matrix: 0.05 [s] he: nls: iter: 1, residual: 1.280475e-03 (rel: 4.531862e-01) he: rezidual: 0.01 [s] he: solve: 0.15 [s] he: matrix: 0.06 [s] he: nls: iter: 2, residual: 2.066041e-04 (rel: 7.312140e-02) he: rezidual: 0.01 [s] he: solve: 0.14 [s] he: matrix: 0.07 [s] he: nls: iter: 3, residual: 4.065440e-05 (rel: 1.438842e-02) he: rezidual: 0.01 [s] he: solve: 0.15 [s] he: matrix: 0.07 [s] he: nls: iter: 4, residual: 8.758871e-06 (rel: 3.099943e-03) he: rezidual: 0.01 [s] he: solve: 0.14 [s] he: matrix: 0.07 [s] he: nls: iter: 5, residual: 1.921880e-06 (rel: 6.801925e-04) he: rezidual: 0.01 [s] he: solve: 0.15 [s] he: matrix: 0.06 [s] he: nls: iter: 6, residual: 4.244317e-07 (rel: 1.502150e-04) he: rezidual: 0.01 [s] he: solve: 0.14 [s] he: matrix: 0.06 [s] he: nls: iter: 7, residual: 9.430908e-08 (rel: 3.337790e-05) he: rezidual: 0.01 [s] he: solve: 0.14 [s] he: matrix: 0.06 [s] he: nls: iter: 8, residual: 2.109056e-08 (rel: 7.464377e-06) he: rezidual: 0.01 [s] he: solve: 0.13 [s] he: matrix: 0.07 [s] he: nls: iter: 9, residual: 4.746932e-09 (rel: 1.680036e-06) he: equation "tmp": he: dw_ul_he_neohook.3.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_ul_he_neohook.3.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_ul_he_mooney_rivlin.3.Omega( solid.kappa, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_ul_bulk_pressure.3.Omega( v, u, p ) he: updating materials... he: ...done in 0.00 s he: ====== time 2.000000e-01 (step 3 of 11) ===== he: updating variables... he: ...done angle: 20.0 he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 2.831297e-03 (rel: 1.000000e+00) he: rezidual: 0.01 [s] he: solve: 0.13 [s] he: matrix: 0.06 [s] he: nls: iter: 1, residual: 1.320418e-03 (rel: 4.663650e-01) he: rezidual: 0.01 [s] he: solve: 0.14 [s] he: matrix: 0.07 [s] he: nls: iter: 2, residual: 2.188046e-04 (rel: 7.728069e-02) he: rezidual: 0.01 [s] he: solve: 0.14 [s] he: matrix: 0.07 [s] he: nls: iter: 3, residual: 4.320423e-05 (rel: 1.525952e-02) he: rezidual: 0.01 [s] he: solve: 0.15 [s] he: matrix: 0.06 [s] he: nls: iter: 4, residual: 9.350833e-06 (rel: 3.302668e-03) he: rezidual: 0.01 [s] he: solve: 0.13 [s] he: matrix: 0.07 [s] he: nls: iter: 5, residual: 2.063723e-06 (rel: 7.288966e-04) he: rezidual: 0.01 [s] he: solve: 0.14 [s] he: matrix: 0.07 [s] he: nls: iter: 6, residual: 4.583352e-07 (rel: 1.618817e-04) he: rezidual: 0.01 [s] he: solve: 0.14 [s] he: matrix: 0.06 [s] he: nls: iter: 7, residual: 1.023736e-07 (rel: 3.615787e-05) he: rezidual: 0.01 [s] he: solve: 0.13 [s] he: matrix: 0.06 [s] he: nls: iter: 8, residual: 2.300609e-08 (rel: 8.125637e-06) he: rezidual: 0.01 [s] he: solve: 0.14 [s] he: matrix: 0.06 [s] he: nls: iter: 9, residual: 5.203114e-09 (rel: 1.837714e-06) he: equation "tmp": he: dw_ul_he_neohook.3.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_ul_he_neohook.3.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_ul_he_mooney_rivlin.3.Omega( solid.kappa, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_ul_bulk_pressure.3.Omega( v, u, p ) he: updating materials... he: ...done in 0.00 s he: ====== time 3.000000e-01 (step 4 of 11) ===== he: updating variables... he: ...done angle: 30.0 he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 2.846546e-03 (rel: 1.000000e+00) he: rezidual: 0.01 [s] he: solve: 0.13 [s] he: matrix: 0.07 [s] he: nls: iter: 1, residual: 1.365099e-03 (rel: 4.795633e-01) he: rezidual: 0.01 [s] he: solve: 0.14 [s] he: matrix: 0.07 [s] he: nls: iter: 2, residual: 2.323577e-04 (rel: 8.162794e-02) he: rezidual: 0.01 [s] he: solve: 0.14 [s] he: matrix: 0.06 [s] he: nls: iter: 3, residual: 4.609065e-05 (rel: 1.619178e-02) he: rezidual: 0.01 [s] he: solve: 0.13 [s] he: matrix: 0.06 [s] he: nls: iter: 4, residual: 1.003052e-05 (rel: 3.523751e-03) he: rezidual: 0.01 [s] he: solve: 0.15 [s] he: matrix: 0.07 [s] he: nls: iter: 5, residual: 2.229409e-06 (rel: 7.831978e-04) he: rezidual: 0.01 [s] he: solve: 0.13 [s] he: matrix: 0.06 [s] he: nls: iter: 6, residual: 4.987536e-07 (rel: 1.752136e-04) he: rezidual: 0.01 [s] he: solve: 0.14 [s] he: matrix: 0.07 [s] he: nls: iter: 7, residual: 1.122087e-07 (rel: 3.941926e-05) he: rezidual: 0.01 [s] he: solve: 0.14 [s] he: matrix: 0.06 [s] he: nls: iter: 8, residual: 2.539982e-08 (rel: 8.923032e-06) he: rezidual: 0.01 [s] he: solve: 0.14 [s] he: matrix: 0.07 [s] he: nls: iter: 9, residual: 5.787902e-09 (rel: 2.033307e-06) he: equation "tmp": he: dw_ul_he_neohook.3.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_ul_he_neohook.3.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_ul_he_mooney_rivlin.3.Omega( solid.kappa, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_ul_bulk_pressure.3.Omega( v, u, p ) he: updating materials... he: ...done in 0.00 s he: ====== time 4.000000e-01 (step 5 of 11) ===== he: updating variables... he: ...done angle: 40.0 he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 2.871296e-03 (rel: 1.000000e+00) he: rezidual: 0.01 [s] he: solve: 0.14 [s] he: matrix: 0.06 [s] he: nls: iter: 1, residual: 1.410227e-03 (rel: 4.911467e-01) he: rezidual: 0.01 [s] he: solve: 0.14 [s] he: matrix: 0.08 [s] he: nls: iter: 2, residual: 2.466238e-04 (rel: 8.589286e-02) he: rezidual: 0.01 [s] he: solve: 0.15 [s] he: matrix: 0.07 [s] he: nls: iter: 3, residual: 4.916863e-05 (rel: 1.712420e-02) he: rezidual: 0.01 [s] he: solve: 0.14 [s] he: matrix: 0.07 [s] he: nls: iter: 4, residual: 1.075370e-05 (rel: 3.745244e-03) he: rezidual: 0.01 [s] he: solve: 0.14 [s] he: matrix: 0.07 [s] he: nls: iter: 5, residual: 2.406194e-06 (rel: 8.380167e-04) he: rezidual: 0.01 [s] he: solve: 0.14 [s] he: matrix: 0.07 [s] he: nls: iter: 6, residual: 5.423534e-07 (rel: 1.888880e-04) he: rezidual: 0.01 [s] he: solve: 0.14 [s] he: matrix: 0.06 [s] he: nls: iter: 7, residual: 1.230088e-07 (rel: 4.284089e-05) he: rezidual: 0.01 [s] he: solve: 0.13 [s] he: matrix: 0.07 [s] he: nls: iter: 8, residual: 2.809049e-08 (rel: 9.783210e-06) he: rezidual: 0.01 [s] he: solve: 0.14 [s] he: matrix: 0.07 [s] he: nls: iter: 9, residual: 6.463575e-09 (rel: 2.251101e-06) he: equation "tmp": he: dw_ul_he_neohook.3.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_ul_he_neohook.3.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_ul_he_mooney_rivlin.3.Omega( solid.kappa, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_ul_bulk_pressure.3.Omega( v, u, p ) he: updating materials... he: ...done in 0.00 s he: ====== time 5.000000e-01 (step 6 of 11) ===== he: updating variables... he: ...done angle: 50.0 he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 2.905113e-03 (rel: 1.000000e+00) he: rezidual: 0.01 [s] he: solve: 0.14 [s] he: matrix: 0.08 [s] he: nls: iter: 1, residual: 1.452418e-03 (rel: 4.999523e-01) he: rezidual: 0.01 [s] he: solve: 0.13 [s] he: matrix: 0.06 [s] he: nls: iter: 2, residual: 2.609390e-04 (rel: 8.982058e-02) he: rezidual: 0.01 [s] he: solve: 0.15 [s] he: matrix: 0.06 [s] he: nls: iter: 3, residual: 5.229190e-05 (rel: 1.799995e-02) he: rezidual: 0.01 [s] he: solve: 0.14 [s] he: matrix: 0.06 [s] he: nls: iter: 4, residual: 1.148899e-05 (rel: 3.954748e-03) he: rezidual: 0.01 [s] he: solve: 0.13 [s] he: matrix: 0.06 [s] he: nls: iter: 5, residual: 2.587489e-06 (rel: 8.906671e-04) he: rezidual: 0.01 [s] he: solve: 0.13 [s] he: matrix: 0.07 [s] he: nls: iter: 6, residual: 5.878880e-07 (rel: 2.023632e-04) he: rezidual: 0.01 [s] he: solve: 0.15 [s] he: matrix: 0.06 [s] he: nls: iter: 7, residual: 1.345896e-07 (rel: 4.632851e-05) he: rezidual: 0.01 [s] he: solve: 0.14 [s] he: matrix: 0.06 [s] he: nls: iter: 8, residual: 3.107131e-08 (rel: 1.069539e-05) he: rezidual: 0.01 [s] he: solve: 0.14 [s] he: matrix: 0.06 [s] he: nls: iter: 9, residual: 7.240371e-09 (rel: 2.492285e-06) he: equation "tmp": he: dw_ul_he_neohook.3.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_ul_he_neohook.3.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_ul_he_mooney_rivlin.3.Omega( solid.kappa, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_ul_bulk_pressure.3.Omega( v, u, p ) he: updating materials... he: ...done in 0.00 s he: ====== time 6.000000e-01 (step 7 of 11) ===== he: updating variables... he: ...done angle: 60.0 he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 2.947090e-03 (rel: 1.000000e+00) he: rezidual: 0.01 [s] he: solve: 0.13 [s] he: matrix: 0.08 [s] he: nls: iter: 1, residual: 1.493233e-03 (rel: 5.066805e-01) he: rezidual: 0.01 [s] he: solve: 0.14 [s] he: matrix: 0.07 [s] he: nls: iter: 2, residual: 2.752354e-04 (rel: 9.339226e-02) he: rezidual: 0.01 [s] he: solve: 0.14 [s] he: matrix: 0.06 [s] he: nls: iter: 3, residual: 5.551565e-05 (rel: 1.883745e-02) he: rezidual: 0.01 [s] he: solve: 0.14 [s] he: matrix: 0.08 [s] he: nls: iter: 4, residual: 1.226910e-05 (rel: 4.163123e-03) he: rezidual: 0.01 [s] he: solve: 0.14 [s] he: matrix: 0.07 [s] he: nls: iter: 5, residual: 2.785385e-06 (rel: 9.451306e-04) he: rezidual: 0.01 [s] he: solve: 0.15 [s] he: matrix: 0.06 [s] he: nls: iter: 6, residual: 6.391432e-07 (rel: 2.168726e-04) he: rezidual: 0.01 [s] he: solve: 0.13 [s] he: matrix: 0.07 [s] he: nls: iter: 7, residual: 1.480443e-07 (rel: 5.023406e-05) he: rezidual: 0.01 [s] he: solve: 0.14 [s] he: matrix: 0.07 [s] he: nls: iter: 8, residual: 3.464537e-08 (rel: 1.175579e-05) he: rezidual: 0.01 [s] he: solve: 0.13 [s] he: matrix: 0.06 [s] he: nls: iter: 9, residual: 8.200973e-09 (rel: 2.782736e-06) he: equation "tmp": he: dw_ul_he_neohook.3.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_ul_he_neohook.3.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_ul_he_mooney_rivlin.3.Omega( solid.kappa, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_ul_bulk_pressure.3.Omega( v, u, p ) he: updating materials... he: ...done in 0.00 s he: ====== time 7.000000e-01 (step 8 of 11) ===== he: updating variables... he: ...done angle: 70.0 he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 2.996058e-03 (rel: 1.000000e+00) he: rezidual: 0.01 [s] he: solve: 0.14 [s] he: matrix: 0.06 [s] he: nls: iter: 1, residual: 1.538165e-03 (rel: 5.133961e-01) he: rezidual: 0.01 [s] he: solve: 0.15 [s] he: matrix: 0.07 [s] he: nls: iter: 2, residual: 2.904725e-04 (rel: 9.695156e-02) he: rezidual: 0.01 [s] he: solve: 0.13 [s] he: matrix: 0.08 [s] he: nls: iter: 3, residual: 5.924227e-05 (rel: 1.977340e-02) he: rezidual: 0.01 [s] he: solve: 0.14 [s] he: matrix: 0.06 [s] he: nls: iter: 4, residual: 1.322565e-05 (rel: 4.414350e-03) he: rezidual: 0.01 [s] he: solve: 0.13 [s] he: matrix: 0.07 [s] he: nls: iter: 5, residual: 3.039323e-06 (rel: 1.014441e-03) he: rezidual: 0.01 [s] he: solve: 0.14 [s] he: matrix: 0.07 [s] he: nls: iter: 6, residual: 7.073666e-07 (rel: 2.360991e-04) he: rezidual: 0.01 [s] he: solve: 0.13 [s] he: matrix: 0.07 [s] he: nls: iter: 7, residual: 1.664786e-07 (rel: 5.556587e-05) he: rezidual: 0.01 [s] he: solve: 0.14 [s] he: matrix: 0.06 [s] he: nls: iter: 8, residual: 3.965172e-08 (rel: 1.323463e-05) he: rezidual: 0.01 [s] he: solve: 0.15 [s] he: matrix: 0.06 [s] he: nls: iter: 9, residual: 9.569221e-09 (rel: 3.193937e-06) he: equation "tmp": he: dw_ul_he_neohook.3.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_ul_he_neohook.3.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_ul_he_mooney_rivlin.3.Omega( solid.kappa, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_ul_bulk_pressure.3.Omega( v, u, p ) he: updating materials... he: ...done in 0.00 s he: ====== time 8.000000e-01 (step 9 of 11) ===== he: updating variables... he: ...done angle: 80.0 he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 3.050781e-03 (rel: 1.000000e+00) he: rezidual: 0.01 [s] he: solve: 0.15 [s] he: matrix: 0.07 [s] he: nls: iter: 1, residual: 1.591299e-03 (rel: 5.216040e-01) he: rezidual: 0.01 [s] he: solve: 0.14 [s] he: matrix: 0.06 [s] he: nls: iter: 2, residual: 3.080438e-04 (rel: 1.009721e-01) he: rezidual: 0.01 [s] he: solve: 0.16 [s] he: matrix: 0.07 [s] he: nls: iter: 3, residual: 6.402806e-05 (rel: 2.098744e-02) he: rezidual: 0.01 [s] he: solve: 0.20 [s] he: matrix: 0.09 [s] he: nls: iter: 4, residual: 1.453852e-05 (rel: 4.765508e-03) he: rezidual: 0.01 [s] he: solve: 0.17 [s] he: matrix: 0.09 [s] he: nls: iter: 5, residual: 3.403540e-06 (rel: 1.115629e-03) he: rezidual: 0.01 [s] he: solve: 0.15 [s] he: matrix: 0.08 [s] he: nls: iter: 6, residual: 8.084584e-07 (rel: 2.650005e-04) he: rezidual: 0.01 [s] he: solve: 0.13 [s] he: matrix: 0.06 [s] he: nls: iter: 7, residual: 1.944950e-07 (rel: 6.375255e-05) he: rezidual: 0.01 [s] he: solve: 0.13 [s] he: matrix: 0.07 [s] he: nls: iter: 8, residual: 4.741585e-08 (rel: 1.554220e-05) he: rezidual: 0.01 [s] he: solve: 0.13 [s] he: matrix: 0.08 [s] he: nls: iter: 9, residual: 1.172689e-08 (rel: 3.843897e-06) he: rezidual: 0.01 [s] he: solve: 0.14 [s] he: matrix: 0.07 [s] he: nls: iter: 10, residual: 2.945911e-09 (rel: 9.656253e-07) he: equation "tmp": he: dw_ul_he_neohook.3.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_ul_he_neohook.3.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_ul_he_mooney_rivlin.3.Omega( solid.kappa, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_ul_bulk_pressure.3.Omega( v, u, p ) he: updating materials... he: ...done in 0.00 s he: ====== time 9.000000e-01 (step 10 of 11) ===== he: updating variables... he: ...done angle: 90.0 he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 3.109464e-03 (rel: 1.000000e+00) he: rezidual: 0.01 [s] he: solve: 0.15 [s] he: matrix: 0.06 [s] he: nls: iter: 1, residual: 1.650237e-03 (rel: 5.307142e-01) he: rezidual: 0.01 [s] he: solve: 0.14 [s] he: matrix: 0.06 [s] he: nls: iter: 2, residual: 3.283547e-04 (rel: 1.055985e-01) he: rezidual: 0.01 [s] he: solve: 0.15 [s] he: matrix: 0.06 [s] he: nls: iter: 3, residual: 7.011744e-05 (rel: 2.254968e-02) he: rezidual: 0.01 [s] he: solve: 0.14 [s] he: matrix: 0.06 [s] he: nls: iter: 4, residual: 1.630204e-05 (rel: 5.242717e-03) he: rezidual: 0.01 [s] he: solve: 0.13 [s] he: matrix: 0.07 [s] he: nls: iter: 5, residual: 3.911218e-06 (rel: 1.257843e-03) he: rezidual: 0.01 [s] he: solve: 0.14 [s] he: matrix: 0.07 [s] he: nls: iter: 6, residual: 9.537691e-07 (rel: 3.067310e-04) he: rezidual: 0.01 [s] he: solve: 0.15 [s] he: matrix: 0.07 [s] he: nls: iter: 7, residual: 2.359252e-07 (rel: 7.587327e-05) he: rezidual: 0.01 [s] he: solve: 0.13 [s] he: matrix: 0.07 [s] he: nls: iter: 8, residual: 5.921798e-08 (rel: 1.904443e-05) he: rezidual: 0.01 [s] he: solve: 0.13 [s] he: matrix: 0.07 [s] he: nls: iter: 9, residual: 1.509739e-08 (rel: 4.855303e-06) he: rezidual: 0.01 [s] he: solve: 0.14 [s] he: matrix: 0.07 [s] he: nls: iter: 10, residual: 3.913855e-09 (rel: 1.258691e-06) he: equation "tmp": he: dw_ul_he_neohook.3.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_ul_he_neohook.3.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_ul_he_mooney_rivlin.3.Omega( solid.kappa, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_ul_bulk_pressure.3.Omega( v, u, p ) he: updating materials... he: ...done in 0.00 s he: ====== time 1.000000e+00 (step 11 of 11) ===== he: updating variables... he: ...done angle: 100.0 he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 3.169908e-03 (rel: 1.000000e+00) he: rezidual: 0.01 [s] he: solve: 0.14 [s] he: matrix: 0.07 [s] he: nls: iter: 1, residual: 1.705699e-03 (rel: 5.380912e-01) he: rezidual: 0.01 [s] he: solve: 0.15 [s] he: matrix: 0.06 [s] he: nls: iter: 2, residual: 3.504169e-04 (rel: 1.105448e-01) he: rezidual: 0.01 [s] he: solve: 0.14 [s] he: matrix: 0.06 [s] he: nls: iter: 3, residual: 7.726134e-05 (rel: 2.437337e-02) he: rezidual: 0.01 [s] he: solve: 0.15 [s] he: matrix: 0.06 [s] he: nls: iter: 4, residual: 1.845776e-05 (rel: 5.822806e-03) he: rezidual: 0.01 [s] he: solve: 0.15 [s] he: matrix: 0.06 [s] he: nls: iter: 5, residual: 4.553882e-06 (rel: 1.436598e-03) he: rezidual: 0.01 [s] he: solve: 0.14 [s] he: matrix: 0.06 [s] he: nls: iter: 6, residual: 1.144430e-06 (rel: 3.610296e-04) he: rezidual: 0.01 [s] he: solve: 0.13 [s] he: matrix: 0.06 [s] he: nls: iter: 7, residual: 2.923765e-07 (rel: 9.223502e-05) he: rezidual: 0.01 [s] he: solve: 0.14 [s] he: matrix: 0.06 [s] he: nls: iter: 8, residual: 7.594018e-08 (rel: 2.395659e-05) he: rezidual: 0.01 [s] he: solve: 0.13 [s] he: matrix: 0.07 [s] he: nls: iter: 9, residual: 2.006413e-08 (rel: 6.329563e-06) he: rezidual: 0.01 [s] he: solve: 0.15 [s] he: matrix: 0.06 [s] he: nls: iter: 10, residual: 5.395421e-09 (rel: 1.702075e-06) he: equation "tmp": he: dw_ul_he_neohook.3.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_ul_he_neohook.3.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_ul_he_mooney_rivlin.3.Omega( solid.kappa, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_ul_bulk_pressure.3.Omega( v, u, p ) he: updating materials... he: ...done in 0.00 s ... c:\Python27\lib\site-packages\sfepy\tests\../examples/large_deformation/hyperelastic_ul_up.py solved ... allowed error: rel = 1.000000e-03, abs = 2.948074e-04 ... ||TLF - ULF||: 6.335181e-04 ... ||ULF - ULF_mixed||: 1.949368e-04 +++ test_solution: ok >>> all passed in 34.73 s <<< c:\Python27\lib\site-packages\sfepy\tests\test_input_acoustics.py he: left over: ['verbose', '__builtins__', '_filename', '__file__', '__package__', 'TestInput', 'input_name', '__name__', 'output_name', '__doc__'] he: left over: ['w', '__builtins__', '__doc__', 'rho', '__name__', 'v_n', 'verbose', 'c', 'data_dir', '__package__', '_filename', '__file__'] >>> test instance prepared (1 test(s)) ... solving ../examples/acoustics/acoustics.py... he: reading mesh [line2, tri3, quad4, tetra4, hexa8] (c:\Python27\lib\site-packages\sfepy\meshes\2d\special\two_rectangles.mesh)... he: ...done in 0.00 s he: creating regions... he: Gamma_out he: Omega he: Gamma_in he: ...done in 0.01 s he: equation "Acoustic pressure": he: 117649.0 * dw_laplace.i.Omega( one.one, q, p ) - 1000000.0 * dw_volume_dot.i.Omega( q, p ) - 343000j * dw_surface_dot.i.Gamma_out( q, p ) = 182355950j * dw_surface_integrate.i.Gamma_in( q ) he: using solvers: ts: no ts nls: newton ls: ls he: updating variables... he: ...done he: setting up dof connectivities... he: ...done in 0.00 s he: matrix shape: (61L, 61L) he: assembling matrix graph... he: ...done in 0.00 s he: matrix structural nonzeros: 445 (1.20e-01% fill) he: updating materials... he: one he: ...done in 0.00 s he: nls: iter: 0, residual: 1.116698e+07 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.00 [s] he: matrix: 0.00 [s] he: nls: iter: 1, residual: 1.118384e-07 (rel: 1.001510e-14) ... ../examples/acoustics/acoustics.py solved +++ test_input: ok >>> all passed in 0.05 s <<< c:\Python27\lib\site-packages\sfepy\tests\test_input_acoustics3d.py he: left over: ['verbose', '__builtins__', '_filename', '__file__', '__package__', 'TestInput', 'input_name', '__name__', 'output_name', '__doc__'] he: left over: ['verbose', 'por', 'coef_r', 'v_n', 'rhoc1', 'rhoc2', 'data_dir', 'tw', '__package__', 'k2', 'k1', '__doc__', 'dh', '__builtins__', '__file__', 'R', 'rho', '__name__', 'freq', 'Z', 'c', 'coef_k', 'w', '_filename'] >>> test instance prepared (1 test(s)) ... solving ../examples/acoustics/acoustics3d.py... he: reading mesh [line2, tri3, quad4, tetra4, hexa8] (c:\Python27\lib\site-packages\sfepy\meshes\3d\acoustics_mesh3d.mesh)... he: ...done in 0.06 s he: creating regions... he: Omega_2 he: Omega_1 he: Gamma_in he: Gamma_out he: Omega he: Gamma_12 he: Gamma_12_1 he: Gamma_12_2 he: ...done in 0.03 s he: equation "Acoustic pressure": he: 48.9591921733 * dw_volume_dot.i.Omega_1(q_1, p_1) + (62.9615348435-19.4783636886j) * dw_volume_dot.i.Omega_2(q_2, p_2) - dw_laplace.i.Omega_1(q_1, p_1) - dw_laplace.i.Omega_2(q_2, p_2) - 6.9970845481j * dw_surface_dot.i.Gamma_out(q_1, p_1) + (23.5231162815+8.79079240116j) * dw_jump.i.Gamma_12_1(q_1, p_1, tr(p_2)) + (31.4071015098+4.33837391784j) * dw_jump.i.Gamma_12_2(q_2, p_2, tr(p_1)) = 3720j * dw_surface_integrate.i.Gamma_in(q_1) he: using solvers: ts: no ts nls: newton ls: ls he: updating variables... he: ...done he: setting up dof connectivities... he: ...done in 0.01 s he: matrix shape: (5425L, 5425L) he: assembling matrix graph... he: ...done in 0.01 s he: matrix structural nonzeros: 129619 (4.40e-03% fill) he: updating materials... he: ...done in 0.00 s he: nls: iter: 0, residual: 2.926322e-01 (rel: 1.000000e+00) he: rezidual: 0.04 [s] he: solve: 0.85 [s] he: matrix: 0.06 [s] he: nls: iter: 1, residual: 2.589789e-14 (rel: 8.849981e-14) ... ../examples/acoustics/acoustics3d.py solved +++ test_input: ok >>> all passed in 1.38 s <<< c:\Python27\lib\site-packages\sfepy\tests\test_input_active_fibres.py he: left over: ['verbose', '__builtins__', '_filename', '__file__', '__package__', 'TestInputEvolutionary', 'input_name', '__doc__', '__name__', 'output_name_trunk'] he: left over: ['verbose', '__builtins__', '__file__', '__name__', 'vf_matrix', 'get_pars_fibres', 'data_dir', 'nm', '_filename', '__package__', 'stress_strain', 'vf_fibres2', 'vf_fibres1', '__doc__'] >>> test instance prepared (1 test(s)) ... solving ../examples/large_deformation/active_fibres.py... he: reading mesh [line2, tri3, quad4, tetra4, hexa8] (c:\Python27\lib\site-packages\sfepy\meshes\3d\cylinder.mesh)... he: ...done in 0.01 s he: creating regions... he: Right he: Omega he: Left he: ...done in 0.01 s he: equation "balance": he: dw_tl_he_neohook.i.Omega( solid.mu, v, u ) + dw_tl_bulk_penalty.i.Omega( solid.K, v, u ) + dw_tl_fib_a.i.Omega( f1.fmax, f1.eps_opt, f1.s, f1.fdir, f1.act, v, u ) + dw_tl_fib_a.i.Omega( f2.fmax, f2.eps_opt, f2.s, f2.fdir, f2.act, v, u ) = 0 he: using solvers: ts: ts nls: newton ls: ls he: updating variables... he: ...done he: setting up dof connectivities... he: ...done in 0.00 s he: matrix shape: (981L, 981L) he: assembling matrix graph... he: ...done in 0.01 s he: matrix structural nonzeros: 34839 (3.62e-02% fill) he: updating materials... he: f1 0.0 he: solid he: f2 1.0 he: ...done in 0.00 s he: ====== time 0.000000e+00 (step 1 of 21) ===== he: initial residual: 1.871491e-03 he: equation "tmp": he: dw_tl_he_neohook.i.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_tl_he_neohook.i.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_tl_bulk_penalty.i.Omega( solid.K, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: ====== time 5.000000e-02 (step 2 of 21) ===== he: updating variables... he: ...done he: updating materials... he: f1 0.0244717418524 he: solid he: f2 0.975528258148 he: ...done in 0.00 s he: nls: iter: 0, residual: 1.825735e-03 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: linesearch: iter 1, (2.66532e-03 < 2.00831e-03) (new ls: 1.000000e-01) he: nls: iter: 1, residual: 1.626534e-03 (rel: 8.908930e-01) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: linesearch: iter 2, (2.10643e-03 < 1.78919e-03) (new ls: 1.000000e-01) he: nls: iter: 2, residual: 1.450798e-03 (rel: 7.946380e-01) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: linesearch: iter 3, (1.67033e-03 < 1.59588e-03) (new ls: 1.000000e-01) he: nls: iter: 3, residual: 1.295386e-03 (rel: 7.095148e-01) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 4, residual: 1.328215e-03 (rel: 7.274964e-01) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 5, residual: 3.892189e-06 (rel: 2.131848e-03) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: linesearch: iter 6, (1.05986e-05 < 4.28141e-06) (new ls: 1.000000e-01) he: nls: iter: 6, residual: 3.562892e-06 (rel: 1.951484e-03) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: linesearch: iter 7, (8.59975e-06 < 3.91918e-06) (new ls: 1.000000e-01) he: nls: iter: 7, residual: 3.256841e-06 (rel: 1.783852e-03) he: equation "tmp": he: dw_tl_he_neohook.i.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_tl_he_neohook.i.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_tl_bulk_penalty.i.Omega( solid.K, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: ====== time 1.000000e-01 (step 3 of 21) ===== he: updating variables... he: ...done he: updating materials... he: f1 0.0954915028125 he: solid he: f2 0.904508497187 he: ...done in 0.00 s he: nls: iter: 0, residual: 1.290162e-04 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 1, residual: 2.119555e-05 (rel: 1.642859e-01) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 2, residual: 2.025765e-07 (rel: 1.570163e-03) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 3, residual: 5.876780e-10 (rel: 4.555071e-06) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 4, residual: 2.888251e-16 (rel: 2.238672e-12) he: equation "tmp": he: dw_tl_he_neohook.i.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_tl_he_neohook.i.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_tl_bulk_penalty.i.Omega( solid.K, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: ====== time 1.500000e-01 (step 4 of 21) ===== he: updating variables... he: ...done he: updating materials... he: f1 0.206107373854 he: solid he: f2 0.793892626146 he: ...done in 0.00 s he: nls: iter: 0, residual: 2.009435e-04 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 1, residual: 6.279918e-05 (rel: 3.125215e-01) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 2, residual: 2.040178e-08 (rel: 1.015299e-04) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 3, residual: 3.632013e-10 (rel: 1.807479e-06) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 4, residual: 8.890194e-17 (rel: 4.424225e-13) he: equation "tmp": he: dw_tl_he_neohook.i.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_tl_he_neohook.i.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_tl_bulk_penalty.i.Omega( solid.K, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: ====== time 2.000000e-01 (step 5 of 21) ===== he: updating variables... he: ...done he: updating materials... he: f1 0.345491502813 he: solid he: f2 0.654508497187 he: ...done in 0.00 s he: nls: iter: 0, residual: 2.558460e-04 (rel: 1.000000e+00) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 1, residual: 1.007349e-04 (rel: 3.937324e-01) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 2, residual: 8.769534e-08 (rel: 3.427661e-04) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 3, residual: 5.148864e-09 (rel: 2.012485e-05) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 4, residual: 4.575917e-16 (rel: 1.788544e-12) he: equation "tmp": he: dw_tl_he_neohook.i.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_tl_he_neohook.i.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_tl_bulk_penalty.i.Omega( solid.K, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: ====== time 2.500000e-01 (step 6 of 21) ===== he: updating variables... he: ...done he: updating materials... he: f1 0.5 he: solid he: f2 0.5 he: ...done in 0.00 s he: nls: iter: 0, residual: 2.873782e-04 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 1, residual: 1.257634e-04 (rel: 4.376235e-01) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 2, residual: 2.153476e-07 (rel: 7.493525e-04) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 3, residual: 2.161331e-08 (rel: 7.520862e-05) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 4, residual: 1.056494e-14 (rel: 3.676318e-11) he: equation "tmp": he: dw_tl_he_neohook.i.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_tl_he_neohook.i.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_tl_bulk_penalty.i.Omega( solid.K, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: ====== time 3.000000e-01 (step 7 of 21) ===== he: updating variables... he: ...done he: updating materials... he: f1 0.654508497187 he: solid he: f2 0.345491502813 he: ...done in 0.00 s he: nls: iter: 0, residual: 2.916571e-04 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 1, residual: 1.285152e-04 (rel: 4.406380e-01) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 2, residual: 2.618252e-07 (rel: 8.977158e-04) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 3, residual: 2.712685e-08 (rel: 9.300938e-05) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 4, residual: 1.848262e-14 (rel: 6.337105e-11) he: equation "tmp": he: dw_tl_he_neohook.i.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_tl_he_neohook.i.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_tl_bulk_penalty.i.Omega( solid.K, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: ====== time 3.500000e-01 (step 8 of 21) ===== he: updating variables... he: ...done he: updating materials... he: f1 0.793892626146 he: solid he: f2 0.206107373854 he: ...done in 0.00 s he: nls: iter: 0, residual: 2.670530e-04 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 1, residual: 1.073173e-04 (rel: 4.018576e-01) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 2, residual: 1.492691e-07 (rel: 5.589491e-04) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 3, residual: 9.930141e-09 (rel: 3.718415e-05) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 4, residual: 2.108239e-15 (rel: 7.894459e-12) he: equation "tmp": he: dw_tl_he_neohook.i.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_tl_he_neohook.i.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_tl_bulk_penalty.i.Omega( solid.K, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: ====== time 4.000000e-01 (step 9 of 21) ===== he: updating variables... he: ...done he: updating materials... he: f1 0.904508497187 he: solid he: f2 0.0954915028125 he: ...done in 0.00 s he: nls: iter: 0, residual: 2.148077e-04 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 1, residual: 6.938624e-05 (rel: 3.230155e-01) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 2, residual: 4.220623e-08 (rel: 1.964837e-04) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 3, residual: 1.006049e-09 (rel: 4.683488e-06) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 4, residual: 9.720532e-17 (rel: 4.525224e-13) he: equation "tmp": he: dw_tl_he_neohook.i.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_tl_he_neohook.i.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_tl_bulk_penalty.i.Omega( solid.K, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: ====== time 4.500000e-01 (step 10 of 21) ===== he: updating variables... he: ...done he: updating materials... he: f1 0.975528258148 he: solid he: f2 0.0244717418524 he: ...done in 0.00 s he: nls: iter: 0, residual: 1.393949e-04 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 1, residual: 2.924375e-05 (rel: 2.097906e-01) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 2, residual: 5.128380e-09 (rel: 3.679029e-05) he: rezidual: 0.00 [s] he: solve: 0.02 [s] he: matrix: 0.05 [s] he: nls: iter: 3, residual: 1.745864e-11 (rel: 1.252459e-07) he: equation "tmp": he: dw_tl_he_neohook.i.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_tl_he_neohook.i.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_tl_bulk_penalty.i.Omega( solid.K, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: ====== time 5.000000e-01 (step 11 of 21) ===== he: updating variables... he: ...done he: updating materials... he: f1 1.0 he: solid he: f2 0.0 he: ...done in 0.00 s he: nls: iter: 0, residual: 4.836127e-05 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 1, residual: 3.522194e-06 (rel: 7.283088e-02) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 2, residual: 5.688862e-11 (rel: 1.176326e-06) he: equation "tmp": he: dw_tl_he_neohook.i.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_tl_he_neohook.i.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_tl_bulk_penalty.i.Omega( solid.K, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: ====== time 5.500000e-01 (step 12 of 21) ===== he: updating variables... he: ...done he: updating materials... he: f1 0.975528258148 he: solid he: f2 0.0244717418524 he: ...done in 0.00 s he: nls: iter: 0, residual: 4.847485e-05 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 1, residual: 3.536683e-06 (rel: 7.295912e-02) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 2, residual: 4.926579e-11 (rel: 1.016316e-06) he: equation "tmp": he: dw_tl_he_neohook.i.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_tl_he_neohook.i.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_tl_bulk_penalty.i.Omega( solid.K, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: ====== time 6.000000e-01 (step 13 of 21) ===== he: updating variables... he: ...done he: updating materials... he: f1 0.904508497187 he: solid he: f2 0.0954915028125 he: ...done in 0.00 s he: nls: iter: 0, residual: 1.403499e-04 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 1, residual: 2.959040e-05 (rel: 2.108331e-01) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 2, residual: 3.244964e-09 (rel: 2.312053e-05) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 3, residual: 7.018625e-12 (rel: 5.000805e-08) he: equation "tmp": he: dw_tl_he_neohook.i.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_tl_he_neohook.i.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_tl_bulk_penalty.i.Omega( solid.K, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: ====== time 6.500000e-01 (step 14 of 21) ===== he: updating variables... he: ...done he: updating materials... he: f1 0.793892626146 he: solid he: f2 0.206107373854 he: ...done in 0.00 s he: nls: iter: 0, residual: 2.171127e-04 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 1, residual: 7.063745e-05 (rel: 3.253492e-01) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 2, residual: 1.852272e-08 (rel: 8.531386e-05) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 3, residual: 2.356467e-10 (rel: 1.085366e-06) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 4, residual: 9.471481e-17 (rel: 4.362472e-13) he: equation "tmp": he: dw_tl_he_neohook.i.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_tl_he_neohook.i.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_tl_bulk_penalty.i.Omega( solid.K, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: ====== time 7.000000e-01 (step 15 of 21) ===== he: updating variables... he: ...done he: updating materials... he: f1 0.654508497187 he: solid he: f2 0.345491502813 he: ...done in 0.00 s he: nls: iter: 0, residual: 2.706736e-04 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 1, residual: 1.096382e-04 (rel: 4.050569e-01) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 2, residual: 4.604913e-08 (rel: 1.701279e-04) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 3, residual: 1.532707e-09 (rel: 5.662568e-06) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 4, residual: 8.940401e-17 (rel: 3.303020e-13) he: equation "tmp": he: dw_tl_he_neohook.i.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_tl_he_neohook.i.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_tl_bulk_penalty.i.Omega( solid.K, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: ====== time 7.500000e-01 (step 16 of 21) ===== he: updating variables... he: ...done he: updating materials... he: f1 0.5 he: solid he: f2 0.5 he: ...done in 0.00 s he: nls: iter: 0, residual: 2.960306e-04 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 1, residual: 1.313465e-04 (rel: 4.436922e-01) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 2, residual: 6.663732e-08 (rel: 2.251028e-04) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 3, residual: 3.357193e-09 (rel: 1.134070e-05) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 4, residual: 1.666338e-16 (rel: 5.628940e-13) he: equation "tmp": he: dw_tl_he_neohook.i.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_tl_he_neohook.i.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_tl_bulk_penalty.i.Omega( solid.K, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: ====== time 8.000000e-01 (step 17 of 21) ===== he: updating variables... he: ...done he: updating materials... he: f1 0.345491502813 he: solid he: f2 0.654508497187 he: ...done in 0.00 s he: nls: iter: 0, residual: 2.916571e-04 (rel: 1.000000e+00) he: rezidual: 0.01 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 1, residual: 1.282192e-04 (rel: 4.396232e-01) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 2, residual: 5.965099e-08 (rel: 2.045244e-04) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 3, residual: 2.733531e-09 (rel: 9.372412e-06) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 4, residual: 1.326539e-16 (rel: 4.548282e-13) he: equation "tmp": he: dw_tl_he_neohook.i.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_tl_he_neohook.i.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_tl_bulk_penalty.i.Omega( solid.K, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: ====== time 8.500000e-01 (step 18 of 21) ===== he: updating variables... he: ...done he: updating materials... he: f1 0.206107373854 he: solid he: f2 0.793892626146 he: ...done in 0.00 s he: nls: iter: 0, residual: 2.592476e-04 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 1, residual: 1.022776e-04 (rel: 3.945170e-01) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 2, residual: 3.379945e-08 (rel: 1.303752e-04) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 3, residual: 8.395345e-10 (rel: 3.238350e-06) he: rezidual: 0.00 [s] he: solve: 0.02 [s] he: matrix: 0.05 [s] he: nls: iter: 4, residual: 8.441667e-17 (rel: 3.256218e-13) he: equation "tmp": he: dw_tl_he_neohook.i.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_tl_he_neohook.i.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_tl_bulk_penalty.i.Omega( solid.K, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: ====== time 9.000000e-01 (step 19 of 21) ===== he: updating variables... he: ...done he: updating materials... he: f1 0.0954915028125 he: solid he: f2 0.904508497187 he: ...done in 0.00 s he: nls: iter: 0, residual: 2.030405e-04 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 1, residual: 6.346827e-05 (rel: 3.125891e-01) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 2, residual: 1.169401e-08 (rel: 5.759444e-05) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 3, residual: 9.054352e-11 (rel: 4.459381e-07) he: equation "tmp": he: dw_tl_he_neohook.i.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_tl_he_neohook.i.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_tl_bulk_penalty.i.Omega( solid.K, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: ====== time 9.500000e-01 (step 20 of 21) ===== he: updating variables... he: ...done he: updating materials... he: f1 0.0244717418524 he: solid he: f2 0.975528258148 he: ...done in 0.00 s he: nls: iter: 0, residual: 1.290137e-04 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 1, residual: 2.590710e-05 (rel: 2.008089e-01) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 2, residual: 1.836518e-09 (rel: 1.423506e-05) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 3, residual: 2.039869e-12 (rel: 1.581126e-08) he: equation "tmp": he: dw_tl_he_neohook.i.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_tl_he_neohook.i.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_tl_bulk_penalty.i.Omega( solid.K, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: ====== time 1.000000e+00 (step 21 of 21) ===== he: updating variables... he: ...done he: updating materials... he: f1 0.0 he: solid he: f2 1.0 he: ...done in 0.00 s he: nls: iter: 0, residual: 4.416200e-05 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 1, residual: 3.058863e-06 (rel: 6.926460e-02) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.05 [s] he: nls: iter: 2, residual: 2.560147e-11 (rel: 5.797172e-07) he: equation "tmp": he: dw_tl_he_neohook.i.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_tl_he_neohook.i.Omega( solid.mu, v, u ) he: updating materials... he: solid he: ...done in 0.00 s he: equation "tmp": he: dw_tl_bulk_penalty.i.Omega( solid.K, v, u ) he: updating materials... he: solid he: ...done in 0.00 s ... ../examples/large_deformation/active_fibres.py solved +++ test_input: ok >>> all passed in 7.24 s <<< c:\Python27\lib\site-packages\sfepy\tests\test_input_biot.py he: left over: ['verbose', '__builtins__', '_filename', '__file__', '__package__', 'TestInput', 'input_name', '__name__', 'output_name', '__doc__'] he: left over: ['verbose', '__builtins__', '__doc__', '__name__', 'data_dir', 'nm', '_filename', '__package__', '__file__'] >>> test instance prepared (1 test(s)) ... solving ../examples/multi_physics/biot.py... he: reading mesh [line2, tri3, quad4, tetra4, hexa8] (c:\Python27\lib\site-packages\sfepy\meshes\3d\cube_medium_hexa.mesh)... he: ...done in 0.01 s he: creating regions... he: Top he: Left he: Omega he: Bottom he: ...done in 0.01 s he: equation "eq_2": he: dw_biot.i1.Omega( m.alpha, u, q ) + dw_diffusion.i1.Omega( m.K, q, p ) = 0 he: equation "eq_1": he: dw_lin_elastic_iso.i2.Omega( m.lam, m.mu, v, u ) - dw_biot.i1.Omega( m.alpha, v, p ) = 0 he: using solvers: ts: no ts nls: newton ls: ls_d he: updating variables... he: ...done he: setting up dof connectivities... he: ...done in 0.00 s he: matrix shape: (4719L, 4719L) he: assembling matrix graph... he: ...done in 0.07 s he: matrix structural nonzeros: 411339 (1.85e-02% fill) he: updating materials... he: m he: ...done in 0.01 s he: nls: iter: 0, residual: 1.955551e+00 (rel: 1.000000e+00) he: rezidual: 0.01 [s] he: solve: 1.12 [s] he: matrix: 0.13 [s] he: nls: iter: 1, residual: 6.416294e-15 (rel: 3.281067e-15) ... ../examples/multi_physics/biot.py solved +++ test_input: ok >>> all passed in 1.53 s <<< c:\Python27\lib\site-packages\sfepy\tests\test_input_biot_npbc.py he: left over: ['verbose', '__builtins__', '_filename', '__file__', '__package__', 'TestInput', 'input_name', '__name__', 'output_name', '__doc__'] he: left over: ['verbose', 'dim', '_filename', 'filename', 'output_dir'] >>> test instance prepared (1 test(s)) ... solving ../examples/multi_physics/biot_npbc.py... he: reading mesh [line2, tri3, quad4, tetra4, hexa8] (c:\Python27\lib\site-packages\sfepy\meshes\3d\cylinder.mesh)... he: ...done in 0.01 s he: creating regions... he: Omega he: Rigid he: Inlet he: Outlet he: Walls he: ...done in 0.02 s he: equation "eq_2": he: dw_biot.i.Omega( m.alpha, u, q ) + dw_diffusion.i.Omega( m.K, q, p ) = 0 he: equation "eq_1": he: dw_lin_elastic.i.Omega( m.D, v, u ) - dw_biot.i.Omega( m.alpha, v, p ) = 0 he: using solvers: ts: no ts nls: newton ls: ls he: updating variables... he: ...done he: setting up dof connectivities... he: ...done in 0.00 s he: lcbc: no_penetration he: lcbc: rigid he: dofs: total 1351, free 724, constrained 627, new 398 he: -> active 1122 he: matrix shape: (1351L, 1351L) he: assembling matrix graph... he: ...done in 0.03 s he: matrix structural nonzeros: 62563 (3.43e-02% fill) he: updating materials... he: m he: ...done in 0.00 s he: nls: iter: 0, residual: 6.009200e-02 (rel: 1.000000e+00) he: rezidual: 0.01 [s] he: solve: 0.02 [s] he: matrix: 0.04 [s] he: nls: iter: 1, residual: 1.227666e-16 (rel: 2.042978e-15) he: equation "tmp": he: ev_diffusion_velocity.i.Omega( m.K, p ) he: updating materials... he: m he: ...done in 0.00 s he: equation "tmp": he: ev_cauchy_stress.i.Omega( m.D, u ) he: updating materials... he: m he: ...done in 0.00 s ... ../examples/multi_physics/biot_npbc.py solved +++ test_input: ok >>> all passed in 0.30 s <<< c:\Python27\lib\site-packages\sfepy\tests\test_input_biot_npbc_lagrange.py he: left over: ['verbose', '__builtins__', '_filename', '__file__', '__package__', 'TestInput', 'input_name', '__name__', 'output_name', '__doc__'] he: left over: ['dim', 'verbose', '_filename', 'filename', 'output_dir'] >>> test instance prepared (1 test(s)) ... solving ../examples/multi_physics/biot_npbc_lagrange.py... he: reading mesh [line2, tri3, quad4, tetra4, hexa8] (c:\Python27\lib\site-packages\sfepy\meshes\3d\cylinder.mesh)... he: ...done in 0.01 s he: creating regions... he: Omega he: Rigid he: Inlet he: Outlet he: Walls he: ...done in 0.01 s he: equation "eq_3": he: dw_non_penetration.2.Walls( u, vl ) = 0 he: equation "eq_2": he: dw_biot.2.Omega( m.alpha, u, q ) + dw_diffusion.2.Omega( m.K, q, p ) = 0 he: equation "eq_1": he: dw_lin_elastic.2.Omega( m.D, v, u ) - dw_biot.2.Omega( m.alpha, v, p ) + dw_non_penetration.2.Walls( v, ul ) = 0 he: using solvers: ts: no ts nls: newton ls: ls he: updating variables... he: ...done he: setting up dof connectivities... he: ...done in 0.00 s he: lcbc: rigid he: dofs: total 1547, free 1508, constrained 39, new 6 he: -> active 1514 he: matrix shape: (1547L, 1547L) he: assembling matrix graph... he: ...done in 0.03 s he: matrix structural nonzeros: 70459 (2.94e-02% fill) he: updating materials... he: m he: ...done in 0.00 s he: nls: iter: 0, residual: 6.009200e-02 (rel: 1.000000e+00) he: rezidual: 0.01 [s] he: solve: 0.06 [s] he: matrix: 0.03 [s] he: nls: iter: 1, residual: 2.024846e-16 (rel: 3.369577e-15) he: equation "tmp": he: ev_diffusion_velocity.2.Omega( m.K, p ) he: updating materials... he: m he: ...done in 0.00 s he: equation "tmp": he: ev_cauchy_stress.2.Omega( m.D, u ) he: updating materials... he: m he: ...done in 0.00 s ... ../examples/multi_physics/biot_npbc_lagrange.py solved +++ test_input: ok >>> all passed in 0.51 s <<< c:\Python27\lib\site-packages\sfepy\tests\test_input_cube.py he: left over: ['verbose', '__builtins__', '_filename', '__file__', '__package__', 'TestInput', 'input_name', '__name__', 'output_name', '__doc__'] he: left over: ['verbose', '__file__', '__name__', 'data_dir', '__package__', '__builtins__', '_filename', '__doc__'] >>> test instance prepared (1 test(s)) ... solving ../examples/diffusion/cube.py... he: reading mesh [line2, tri3, quad4, tetra4, hexa8] (c:\Python27\lib\site-packages\sfepy\meshes\3d\cube_medium_hexa.mesh)... he: ...done in 0.01 s he: creating regions... he: Top he: Bottom he: Omega he: Surface he: Left he: ...done in 0.01 s he: equation "nice_equation": he: dw_laplace.i.Omega( coef.val, s, T ) = 0 he: using solvers: ts: no ts nls: newton ls: ls he: updating variables... he: ...done he: setting up dof connectivities... he: ...done in 0.00 s he: matrix shape: (729L, 729L) he: assembling matrix graph... he: ...done in 0.00 s he: matrix structural nonzeros: 15625 (2.94e-02% fill) he: updating materials... he: coef he: ...done in 0.00 s he: nls: iter: 0, residual: 4.783598e+00 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.00 [s] he: nls: iter: 1, residual: 4.008183e-15 (rel: 8.379013e-16) ... ../examples/diffusion/cube.py solved +++ test_input: ok >>> all passed in 0.18 s <<< c:\Python27\lib\site-packages\sfepy\tests\test_input_darcy_flow_multicomp.py he: left over: ['verbose', '__builtins__', '_filename', '__file__', '__package__', 'TestInput', 'input_name', '__name__', 'output_name', '__doc__'] he: left over: ['Struct', 'G_bar', 'mat_fun', '__file__', '__doc__', '__name__', 'verbose', 'alpha2', 'alpha1', 'data_dir', 'nm', '__package__', '__builtins__', '_filename', 'postproc'] >>> test instance prepared (1 test(s)) ... solving ../examples/diffusion/darcy_flow_multicomp.py... he: reading mesh [line2, tri3, quad4, tetra4, hexa8] (c:\Python27\lib\site-packages\sfepy\meshes\3d\cube_medium_hexa.mesh)... he: ...done in 0.01 s he: creating regions... he: Sigma_1 he: Source he: Sink he: Omega he: Omega2 he: Omega1 he: ...done in 0.01 s he: equation "komp2": he: dw_diffusion.5.Omega2(mat.K, q2, p2) + dw_volume_dot.5.Omega2(mat.G_alfa, q2, p2) - dw_volume_dot.5.Omega2(mat.G_alfa, q2, p1) = dw_volume_integrate.5.Sink(mat.f_2, q2) he: equation "komp1": he: dw_diffusion.5.Omega1(mat.K, q1, p1) + dw_volume_dot.5.Omega1(mat.G_alfa, q1, p1) - dw_volume_dot.5.Omega1(mat.G_alfa, q1, p2) = dw_volume_integrate.5.Source(mat.f_1, q1) he: using solvers: ts: no ts nls: newton ls: ls he: updating variables... he: ...done he: setting up dof connectivities... he: ...done in 0.00 s he: matrix shape: (2661L, 2661L) he: assembling matrix graph... he: ...done in 0.01 s he: matrix structural nonzeros: 119133 (1.68e-02% fill) he: updating materials... he: mat he: ...done in 0.02 s he: nls: iter: 0, residual: 1.000000e-02 (rel: 1.000000e+00) he: rezidual: 0.04 [s] he: solve: 0.13 [s] he: matrix: 0.11 [s] he: nls: iter: 1, residual: 6.188667e-17 (rel: 6.188667e-15) he: equation "tmp": he: ev_integrate_mat.5.Omega(mat.G_alfa, p1) he: updating materials... he: mat he: ...done in 0.01 s ... ../examples/diffusion/darcy_flow_multicomp.py solved +++ test_input: ok >>> all passed in 0.56 s <<< c:\Python27\lib\site-packages\sfepy\tests\test_input_elastic_contact_planes.py he: left over: ['verbose', '__builtins__', '_filename', '__file__', '__package__', 'TestInput', 'input_name', '__name__', 'output_name', '__doc__'] he: left over: ['f0', 'dn', 'verbose', '__builtins__', '__doc__', '__name__', 'az', 'ds', 'main', 'data_dir', 'k', '__package__', '_filename', '__file__'] >>> test instance prepared (1 test(s)) ... solving ../examples/linear_elasticity/elastic_contact_planes.py... he: reading mesh [line2, tri3, quad4, tetra4, hexa8] (c:\Python27\lib\site-packages\sfepy\meshes\3d\cube_medium_hexa.mesh)... he: ...done in 0.01 s he: creating regions... he: Top he: Omega he: Bottom he: ...done in 0.01 s he: equation "elasticity": he: dw_lin_elastic_iso.2.Omega(solid.lam, solid.mu, v, u) + dw_contact_plane.2.Top(cp0.f, cp0.n, cp0.a, cp0.bs, v, u) + dw_contact_plane.2.Top(cp1.f, cp1.n, cp1.a, cp1.bs, v, u) + dw_contact_plane.2.Top(cp2.f, cp2.n, cp2.a, cp2.bs, v, u) + dw_contact_plane.2.Top(cp3.f, cp3.n, cp3.a, cp3.bs, v, u) = 0 he: using solvers: ts: no ts nls: newton ls: lsd he: updating variables... he: ...done he: setting up dof connectivities... he: ...done in 0.00 s he: matrix shape: (3630L, 3630L) he: assembling matrix graph... he: ...done in 0.04 s he: matrix structural nonzeros: 242172 (1.84e-02% fill) he: updating materials... he: solid he: cp1 he: cp0 he: cp3 he: cp2 he: ...done in 0.01 s he: nls: iter: 0, residual: 3.318495e+02 (rel: 1.000000e+00) he: rezidual: 0.02 [s] he: solve: 0.57 [s] he: matrix: 0.08 [s] he: nls: iter: 1, residual: 1.417974e+00 (rel: 4.272942e-03) he: rezidual: 0.02 [s] he: solve: 0.58 [s] he: matrix: 0.09 [s] he: nls: iter: 2, residual: 3.858461e-01 (rel: 1.162714e-03) he: rezidual: 0.01 [s] he: solve: 0.65 [s] he: matrix: 0.09 [s] he: nls: iter: 3, residual: 8.571644e-02 (rel: 2.582991e-04) he: rezidual: 0.02 [s] he: solve: 0.54 [s] he: matrix: 0.08 [s] he: nls: iter: 4, residual: 1.315954e-02 (rel: 3.965515e-05) he: rezidual: 0.01 [s] he: solve: 0.54 [s] he: matrix: 0.07 [s] he: nls: iter: 5, residual: 6.076395e-04 (rel: 1.831070e-06) he: rezidual: 0.02 [s] he: solve: 0.56 [s] he: matrix: 0.07 [s] he: nls: iter: 6, residual: 1.998774e-06 (rel: 6.023133e-09) he: rezidual: 0.02 [s] he: solve: 0.58 [s] he: matrix: 0.08 [s] he: nls: iter: 7, residual: 3.923866e-11 (rel: 1.182423e-13) ... ../examples/linear_elasticity/elastic_contact_planes.py solved +++ test_input: ok >>> all passed in 4.98 s <<< c:\Python27\lib\site-packages\sfepy\tests\test_input_elastic_contact_sphere.py he: left over: ['verbose', '__builtins__', '_filename', '__file__', '__package__', 'TestInput', 'input_name', '__name__', 'output_name', '__doc__'] he: left over: ['f0', 'verbose', '__builtins__', '__file__', '__doc__', '__name__', 'data_dir', 'k', '__package__', '_filename', 'main'] >>> test instance prepared (1 test(s)) ... solving ../examples/linear_elasticity/elastic_contact_sphere.py... he: reading mesh [line2, tri3, quad4, tetra4, hexa8] (c:\Python27\lib\site-packages\sfepy\meshes\3d\cube_medium_hexa.mesh)... he: ...done in 0.02 s he: creating regions... he: Top he: Omega he: Bottom he: ...done in 0.01 s he: equation "elasticity": he: dw_lin_elastic_iso.2.Omega(solid.lam, solid.mu, v, u) + dw_contact_sphere.2.Top(cs.f, cs.c, cs.r, v, u) = 0 he: using solvers: ts: no ts nls: newton ls: ls he: updating variables... he: ...done he: setting up dof connectivities... he: ...done in 0.00 s he: matrix shape: (3630L, 3630L) he: assembling matrix graph... he: ...done in 0.05 s he: matrix structural nonzeros: 242172 (1.84e-02% fill) he: updating materials... he: cs he: solid he: ...done in 0.01 s he: nls: iter: 0, residual: 3.774543e+02 (rel: 1.000000e+00) he: rezidual: 0.01 [s] he: solve: 0.58 [s] he: matrix: 0.07 [s] he: nls: iter: 1, residual: 1.931773e+00 (rel: 5.117899e-03) he: rezidual: 0.01 [s] he: solve: 0.57 [s] he: matrix: 0.06 [s] he: linesearch: iter 2, (3.17341e+01 < 3.86355e+00) (new ls: 1.000000e-01) he: nls: iter: 2, residual: 3.100783e-01 (rel: 8.214989e-04) he: rezidual: 0.01 [s] he: solve: 0.56 [s] he: matrix: 0.06 [s] he: linesearch: iter 3, (3.03796e+01 < 6.20157e-01) (new ls: 1.000000e-01) he: linesearch: iter 3, (1.51739e+00 < 6.20157e-01) (new ls: 1.000000e-02) he: nls: iter: 3, residual: 3.069746e-01 (rel: 8.132761e-04) he: rezidual: 0.01 [s] he: solve: 0.58 [s] he: matrix: 0.06 [s] he: linesearch: iter 4, (3.03839e+01 < 6.13949e-01) (new ls: 1.000000e-01) he: linesearch: iter 4, (1.74112e+00 < 6.13949e-01) (new ls: 1.000000e-02) he: nls: iter: 4, residual: 3.039936e-01 (rel: 8.053785e-04) he: rezidual: 0.01 [s] he: solve: 0.57 [s] he: matrix: 0.07 [s] he: linesearch: iter 5, (1.94821e+01 < 6.07987e-01) (new ls: 1.000000e-01) he: linesearch: iter 5, (1.08653e+00 < 6.07987e-01) (new ls: 1.000000e-02) he: nls: iter: 5, residual: 3.008785e-01 (rel: 7.971255e-04) he: rezidual: 0.01 [s] he: solve: 0.58 [s] he: matrix: 0.06 [s] he: linesearch: iter 6, (1.78497e+01 < 6.01757e-01) (new ls: 1.000000e-01) he: linesearch: iter 6, (1.11411e+00 < 6.01757e-01) (new ls: 1.000000e-02) he: nls: iter: 6, residual: 3.013698e-01 (rel: 7.984272e-04) he: rezidual: 0.01 [s] he: solve: 0.61 [s] he: matrix: 0.07 [s] he: linesearch: iter 7, (1.09936e+01 < 6.02740e-01) (new ls: 1.000000e-01) he: linesearch: iter 7, (6.79175e-01 < 6.02740e-01) (new ls: 1.000000e-02) he: nls: iter: 7, residual: 2.968858e-01 (rel: 7.865476e-04) he: rezidual: 0.01 [s] he: solve: 0.59 [s] he: matrix: 0.06 [s] he: nls: iter: 8, residual: 5.904641e-01 (rel: 1.564333e-03) he: rezidual: 0.01 [s] he: solve: 0.59 [s] he: matrix: 0.07 [s] he: nls: iter: 9, residual: 1.319303e-01 (rel: 3.495266e-04) he: rezidual: 0.01 [s] he: solve: 0.59 [s] he: matrix: 0.07 [s] he: nls: iter: 10, residual: 1.249511e-01 (rel: 3.310365e-04) he: rezidual: 0.01 [s] he: solve: 0.57 [s] he: matrix: 0.07 [s] he: linesearch: iter 11, (2.51798e-01 < 2.49902e-01) (new ls: 1.000000e-01) he: nls: iter: 11, residual: 1.141365e-01 (rel: 3.023849e-04) he: rezidual: 0.01 [s] he: solve: 0.62 [s] he: matrix: 0.06 [s] he: nls: iter: 12, residual: 2.073594e-02 (rel: 5.493630e-05) ... ../examples/linear_elasticity/elastic_contact_sphere.py solved +++ test_input: ok >>> all passed in 8.19 s <<< c:\Python27\lib\site-packages\sfepy\tests\test_input_elastic_shifted_periodic.py he: left over: ['verbose', '__builtins__', '_filename', '__file__', '__package__', 'TestInput', 'input_name', '__name__', 'output_name', '__doc__'] he: left over: ['verbose', 'linear_tension', '__builtins__', '__file__', '__name__', 'get_von_mises_stress', 'data_dir', 'nm', 'get_shift', 'per', '__package__', 'stiffness_from_lame', '_filename', 'post_process', '__doc__'] >>> test instance prepared (1 test(s)) ... solving ../examples/linear_elasticity/elastic_shifted_periodic.py... he: reading mesh [line2, tri3, quad4, tetra4, hexa8] (c:\Python27\lib\site-packages\sfepy\meshes\3d\block.mesh)... he: ...done in 0.01 s he: creating regions... he: Right he: Bottom he: Far he: Top he: Near he: Omega he: Left he: ...done in 0.01 s he: equation "elasticity": he: dw_lin_elastic.2.Omega(solid.D, v, u) = -dw_surface_ltr.2.Bottom(load.val, v) he: using solvers: ts: no ts nls: newton ls: ls he: updating variables... he: ...done he: setting up dof connectivities... he: ...done in 0.00 s he: lcbc: shifted he: dofs: total 1160, free 928, constrained 232, new 0 he: -> active 928 he: matrix shape: (1160L, 1160L) he: assembling matrix graph... he: ...done in 0.01 s he: matrix structural nonzeros: 78312 (5.82e-02% fill) he: updating materials... he: solid he: load he: ...done in 0.01 s he: nls: iter: 0, residual: 2.576813e+01 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.03 [s] he: nls: iter: 1, residual: 1.403148e-14 (rel: 5.445286e-16) he: equation "tmp": he: ev_cauchy_stress.2.Omega(solid.D, u) he: updating materials... he: solid he: ...done in 0.00 s ... ../examples/linear_elasticity/elastic_shifted_periodic.py solved +++ test_input: ok >>> all passed in 0.18 s <<< c:\Python27\lib\site-packages\sfepy\tests\test_input_its2D_2.py he: left over: ['verbose', '__builtins__', '_filename', '__file__', '__package__', 'TestInput', 'input_name', '__name__', 'output_name', '__doc__'] he: left over: ['verbose', 'asphalt', '__builtins__', '__file__', '__name__', 'stiffness_from_youngpoisson', 'data_dir', 'refinement_level', 'young', 'lame_from_youngpoisson', '__package__', 'stress_strain', 'refine_mesh', 'output_dir', 'poisson', '_filename', '__doc__'] >>> test instance prepared (1 test(s)) ... solving ../examples/linear_elasticity/its2D_2.py... he: reading mesh [line2, tri3, quad4, tetra4, hexa8] (c:\Python27\lib\site-packages\sfepy\meshes\2d\its2D.mesh)... he: ...done in 0.00 s he: creating regions... he: Top he: Bottom he: Omega he: Left he: ...done in 0.01 s he: equation "balance_of_forces": he: dw_lin_elastic_iso.2.Omega(Asphalt.lam, Asphalt.mu, v, u ) = dw_point_load.0.Top(Load.val, v) he: using solvers: ts: no ts nls: newton ls: ls he: updating variables... he: ...done he: setting up dof connectivities... he: ...done in 0.00 s he: matrix shape: (94L, 94L) he: assembling matrix graph... he: ...done in 0.00 s he: matrix structural nonzeros: 1070 (1.21e-01% fill) he: updating materials... he: Asphalt he: Load he: ...done in 0.00 s he: nls: iter: 0, residual: 1.000000e+03 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.00 [s] he: matrix: 0.00 [s] he: nls: iter: 1, residual: 3.773833e-12 (rel: 3.773833e-15) he: equation "tmp": he: ev_cauchy_strain.2.Omega(u) he: updating materials... he: ...done in 0.00 s he: equation "tmp": he: ev_cauchy_stress.2.Omega(Asphalt.D, u) he: updating materials... he: Asphalt he: ...done in 0.00 s ... ../examples/linear_elasticity/its2D_2.py solved +++ test_input: ok >>> all passed in 0.07 s <<< c:\Python27\lib\site-packages\sfepy\tests\test_input_laplace_1d.py he: left over: ['verbose', '__builtins__', '_filename', '__file__', '__package__', 'TestInput', 'input_name', '__name__', 'output_name', '__doc__'] he: left over: ['verbose', '__builtins__', '__file__', '__doc__', 'Mesh', 'nm', 'UserMeshIO', '__name__', '__package__', '_filename', 'get_coef', 'mesh_hook'] >>> test instance prepared (1 test(s)) ... solving ../examples/diffusion/laplace_1d.py... he: reading mesh [user] (function:mesh_hook)... he: ...done in 0.00 s he: creating regions... he: Gamma_Right he: Omega he: Gamma_Left he: ...done in 0.01 s he: equation "Temperature": he: dw_laplace.i.Omega(coef.val, s, t) = 0 he: using solvers: ts: no ts nls: newton ls: ls he: updating variables... he: ...done he: setting up dof connectivities... he: ...done in 0.00 s he: matrix shape: (99L, 99L) he: assembling matrix graph... he: ...done in 0.00 s he: matrix structural nonzeros: 295 (3.01e-02% fill) he: updating materials... he: coef he: ...done in 0.00 s he: nls: iter: 0, residual: 4.298429e+00 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.00 [s] he: matrix: 0.00 [s] he: nls: iter: 1, residual: 1.734146e-14 (rel: 4.034372e-15) ... ../examples/diffusion/laplace_1d.py solved +++ test_input: ok >>> all passed in 0.04 s <<< c:\Python27\lib\site-packages\sfepy\tests\test_input_laplace_coupling_lcbcs.py he: left over: ['verbose', '__builtins__', '_filename', '__file__', '__package__', 'TestInput', 'input_name', '__name__', 'output_name', '__doc__'] he: left over: ['verbose', '__builtins__', '__file__', 'get_shift1', '__name__', 'get_shift2', 'data_dir', 'nm', 'per', '__package__', '_filename', '__doc__'] >>> test instance prepared (1 test(s)) ... solving ../examples/diffusion/laplace_coupling_lcbcs.py... he: reading mesh [line2, tri3, quad4, tetra4, hexa8] (c:\Python27\lib\site-packages\sfepy\meshes\2d\square_quad.mesh)... he: ...done in 0.00 s he: creating regions... he: Right he: Bottom he: Top he: Omega2 he: Omega1 he: Left he: Gamma he: Gamma2 he: Gamma1 he: Omega_m2 he: Omega_m1 he: Omega_m11 he: Omega_m12 he: ...done in 0.01 s he: equation "eq1": he: dw_laplace.i1.Omega1(v1, u1) = 0 he: equation "eq2": he: dw_laplace.i1.Omega2(v2, u2) = 0 he: using solvers: ts: no ts nls: newton ls: ls he: updating variables... he: ...done he: setting up dof connectivities... he: ...done in 0.00 s he: lcbc: mean2 he: lcbc: mean12 he: lcbc: shifted1 he: lcbc: shifted2 he: lcbc: mean11 he: dofs: total 110, free 28, constrained 82, new 3 he: -> active 31 he: matrix shape: (110L, 110L) he: assembling matrix graph... he: ...done in 0.00 s he: matrix structural nonzeros: 806 (6.66e-02% fill) he: updating materials... he: ...done in 0.00 s he: nls: iter: 0, residual: 9.110435e+00 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.00 [s] he: matrix: 0.00 [s] he: nls: iter: 1, residual: 2.417431e-15 (rel: 2.653474e-16) ... ../examples/diffusion/laplace_coupling_lcbcs.py solved +++ test_input: ok >>> all passed in 0.07 s <<< c:\Python27\lib\site-packages\sfepy\tests\test_input_laplace_time_ebcs.py he: left over: ['verbose', '__builtins__', '_filename', '__file__', '__package__', 'TestInputEvolutionary', 'input_name', '__doc__', '__name__', 'output_name_trunk'] he: left over: ['verbose', '__builtins__', '__doc__', '__name__', 'data_dir', '__package__', 'is_ebc', '_filename', '__file__'] >>> test instance prepared (1 test(s)) ... solving ../examples/diffusion/laplace_time_ebcs.py... he: reading mesh [line2, tri3, quad4, tetra4, hexa8] (c:\Python27\lib\site-packages\sfepy\meshes\3d\cube_medium_tetra.mesh)... he: ...done in 0.01 s he: creating regions... he: Right he: Top he: Bottom he: Omega he: Left he: ...done in 0.01 s he: equation "eq": he: dw_laplace.2.Omega( one.val, s, t ) = 0 he: using solvers: ts: ts nls: newton ls: ls he: updating variables... he: ...done he: setting up dof connectivities... he: ...done in 0.00 s he: matrix shape: (341L, 341L) he: assembling matrix graph... he: ...done in 0.00 s he: matrix structural nonzeros: 4049 (3.48e-02% fill) he: ====== time 0.000000e+00 (step 1 of 5) ===== he: updating materials... he: one he: ...done in 0.00 s he: nls: iter: 0, residual: 1.194325e+00 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.00 [s] he: matrix: 0.00 [s] he: nls: iter: 1, residual: 2.153292e-15 (rel: 1.802936e-15) he: ====== time 1.000000e+00 (step 2 of 5) ===== he: updating variables... he: ...done he: setting up dof connectivities... he: ...done in 0.00 s he: matrix shape: (341L, 341L) he: assembling matrix graph... he: ...done in 0.00 s he: matrix structural nonzeros: 4039 (3.47e-02% fill) he: updating materials... he: one he: ...done in 0.00 s he: nls: iter: 0, residual: 2.038845e+00 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.00 [s] he: matrix: 0.00 [s] he: nls: iter: 1, residual: 2.705189e-15 (rel: 1.326824e-15) he: ====== time 2.000000e+00 (step 3 of 5) ===== he: updating variables... he: ...done he: setting up dof connectivities... he: ...done in 0.00 s he: matrix shape: (334L, 334L) he: assembling matrix graph... he: ...done in 0.00 s he: matrix structural nonzeros: 3942 (3.53e-02% fill) he: updating materials... he: one he: ...done in 0.00 s he: nls: iter: 0, residual: 2.245056e+00 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.00 [s] he: matrix: 0.00 [s] he: nls: iter: 1, residual: 3.170676e-15 (rel: 1.412292e-15) he: ====== time 3.000000e+00 (step 4 of 5) ===== he: updating variables... he: ...done he: setting up dof connectivities... he: ...done in 0.00 s he: matrix shape: (341L, 341L) he: assembling matrix graph... he: ...done in 0.00 s he: matrix structural nonzeros: 4049 (3.48e-02% fill) he: updating materials... he: one he: ...done in 0.00 s he: nls: iter: 0, residual: 1.261162e+00 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.00 [s] he: matrix: 0.00 [s] he: nls: iter: 1, residual: 4.579038e-15 (rel: 3.630809e-15) he: ====== time 4.000000e+00 (step 5 of 5) ===== he: updating variables... he: ...done he: setting up dof connectivities... he: ...done in 0.00 s he: matrix shape: (334L, 334L) he: assembling matrix graph... he: ...done in 0.00 s he: matrix structural nonzeros: 3942 (3.53e-02% fill) he: updating materials... he: one he: ...done in 0.00 s he: nls: iter: 0, residual: 2.840053e+00 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.00 [s] he: matrix: 0.00 [s] he: nls: iter: 1, residual: 4.072885e-15 (rel: 1.434088e-15) ... ../examples/diffusion/laplace_time_ebcs.py solved +++ test_input: ok >>> all passed in 0.45 s <<< c:\Python27\lib\site-packages\sfepy\tests\test_input_linear_elastic.py he: left over: ['verbose', '__builtins__', '_filename', '__file__', '__package__', 'TestInput', 'input_name', '__name__', 'output_name', '__doc__'] he: left over: ['data_dir', 'verbose', '_filename', '__builtins__', '__package__', '__doc__', '__name__', '__file__'] >>> test instance prepared (1 test(s)) ... solving ../examples/linear_elasticity/linear_elastic.py... he: reading mesh [line2, tri3, quad4, tetra4, hexa8] (c:\Python27\lib\site-packages\sfepy\meshes\3d\cylinder.mesh)... he: ...done in 0.01 s he: creating regions... he: SomewhereTop he: Right he: Omega he: Left he: ...done in 0.01 s he: equation "balance_of_forces": he: dw_lin_elastic_iso.i.Omega(solid.lam, solid.mu, v, u) = 0 he: using solvers: ts: no ts nls: newton ls: ls he: updating variables... he: ...done he: setting up dof connectivities... he: ...done in 0.00 s he: matrix shape: (886L, 886L) he: assembling matrix graph... he: ...done in 0.01 s he: matrix structural nonzeros: 31052 (3.96e-02% fill) he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 4.778054e-03 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.01 [s] he: nls: iter: 1, residual: 7.193990e-18 (rel: 1.505632e-15) ... ../examples/linear_elasticity/linear_elastic.py solved +++ test_input: ok >>> all passed in 0.12 s <<< c:\Python27\lib\site-packages\sfepy\tests\test_input_linear_elastic_damping.py he: left over: ['verbose', '__builtins__', '_filename', '__file__', '__package__', 'TestInputEvolutionary', 'input_name', '__doc__', '__name__', 'output_name_trunk'] he: left over: ['verbose', '__builtins__', 'print_times', 'ebc_sin', '__file__', 'deepcopy', 'nm', '__name__', 'adapt_time_step', '__package__', 'ls', '_filename', '__doc__'] >>> test instance prepared (1 test(s)) ... solving ../examples/linear_elasticity/linear_elastic_damping.py... he: reading mesh [line2, tri3, quad4, tetra4, hexa8] (c:\Python27\lib\site-packages\sfepy\meshes\3d\cylinder.mesh)... he: ...done in 0.01 s he: creating regions... he: SomewhereTop he: Right he: Omega he: Left he: ...done in 0.01 s he: equation "balance_of_forces in time": he: dw_volume_dot.i.Omega( solid.c, v, du/dt ) + dw_lin_elastic_iso.i.Omega( solid.lam, solid.mu, v, u ) = 0 he: using solvers: ts: ts nls: newton ls: ls he: updating variables... he: ...done he: setting up dof connectivities... he: ...done in 0.00 s he: matrix shape: (872L, 872L) he: assembling matrix graph... he: ...done in 0.01 s he: matrix structural nonzeros: 30398 (4.00e-02% fill) he: updating materials... he: solid he: ...done in 0.00 s he: ====== time 0.000000e+00 (dt 1.000000e-02, wait 0, step 1 of 1) ===== he: initial residual: 4.776140e-03 he: ====== time 1.000000e-02 (dt 1.000000e-02, wait 0, step 2 of 2) ===== he: updating variables... he: ...done he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 4.778523e-03 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.02 [s] he: nls: iter: 1, residual: 3.565620e-18 (rel: 7.461762e-16) he: ====== time 2.000000e-02 (dt 1.000000e-02, wait 0, step 3 of 3) ===== he: updating variables... he: ...done he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 1.015046e-03 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.02 [s] he: nls: iter: 1, residual: 1.144800e-18 (rel: 1.127831e-15) he: ====== time 3.000000e-02 (dt 1.000000e-02, wait 0, step 4 of 4) ===== he: updating variables... he: ...done he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 4.779821e-04 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.02 [s] he: nls: iter: 1, residual: 1.048748e-18 (rel: 2.194115e-15) he: ====== time 4.000000e-02 (dt 1.000000e-02, wait 0, step 5 of 5) ===== he: updating variables... he: ...done he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 3.290553e-04 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.02 [s] he: nls: iter: 1, residual: 9.498657e-19 (rel: 2.886645e-15) he: ====== time 5.000000e-02 (dt 1.000000e-02, wait 0, step 6 of 6) ===== he: updating variables... he: ...done he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 2.711463e-04 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.02 [s] he: nls: iter: 1, residual: 8.363598e-19 (rel: 3.084533e-15) he: ====== time 6.000000e-02 (dt 1.000000e-02, wait 0, step 7 of 7) ===== he: updating variables... he: ...done he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 2.420130e-04 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.02 [s] he: nls: iter: 1, residual: 8.700386e-19 (rel: 3.595008e-15) he: ====== time 7.000000e-02 (dt 1.000000e-02, wait 0, step 8 of 8) ===== he: updating variables... he: ...done he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 2.242887e-04 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.02 [s] he: nls: iter: 1, residual: 9.136908e-19 (rel: 4.073726e-15) he: ====== time 8.000000e-02 (dt 1.000000e-02, wait 0, step 9 of 9) ===== he: updating variables... he: ...done he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 2.118625e-04 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.02 [s] he: nls: iter: 1, residual: 8.529171e-19 (rel: 4.025804e-15) he: ====== time 9.000000e-02 (dt 1.000000e-02, wait 0, step 10 of 10) ===== he: updating variables... he: ...done he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 2.021237e-04 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.02 [s] he: matrix: 0.03 [s] he: nls: iter: 1, residual: 1.023807e-18 (rel: 5.065252e-15) he: ====== time 1.000000e-01 (dt 1.000000e-02, wait 0, step 11 of 11) ===== he: updating variables... he: ...done he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 1.937955e-04 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.02 [s] he: nls: iter: 1, residual: 9.738970e-19 (rel: 5.025385e-15) he: ====== time 1.100000e-01 (dt 1.000000e-02, wait 0, step 12 of 12) ===== he: updating variables... he: ...done he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 1.861977e-04 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.02 [s] he: nls: iter: 1, residual: 9.303022e-19 (rel: 4.996315e-15) he: ====== time 1.200000e-01 (dt 1.000000e-02, wait 0, step 13 of 13) ===== he: updating variables... he: ...done he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 1.789485e-04 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.02 [s] he: nls: iter: 1, residual: 1.069501e-18 (rel: 5.976582e-15) he: ====== time 1.300000e-01 (dt 1.000000e-02, wait 0, step 14 of 14) ===== he: updating variables... he: ...done he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 1.718285e-04 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.02 [s] he: nls: iter: 1, residual: 1.120530e-18 (rel: 6.521213e-15) he: ====== time 1.400000e-01 (dt 1.000000e-02, wait 0, step 15 of 15) ===== he: updating variables... he: ...done he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 1.647120e-04 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.02 [s] he: nls: iter: 1, residual: 1.078559e-18 (rel: 6.548150e-15) he: ====== time 1.500000e-01 (dt 1.000000e-02, wait 0, step 16 of 16) ===== he: updating variables... he: ...done he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 1.575314e-04 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.02 [s] he: nls: iter: 1, residual: 1.087918e-18 (rel: 6.906038e-15) he: ====== time 1.600000e-01 (dt 1.000000e-02, wait 0, step 17 of 17) ===== he: updating variables... he: ...done he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 1.502570e-04 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.02 [s] he: nls: iter: 1, residual: 1.060205e-18 (rel: 7.055946e-15) he: ====== time 1.700000e-01 (dt 1.000000e-02, wait 0, step 18 of 18) ===== he: updating variables... he: ...done he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 1.428862e-04 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.02 [s] he: nls: iter: 1, residual: 1.225346e-18 (rel: 8.575681e-15) he: ====== time 1.800000e-01 (dt 1.000000e-02, wait 0, step 19 of 19) ===== he: updating variables... he: ...done he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 1.354368e-04 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.02 [s] he: nls: iter: 1, residual: 1.263646e-18 (rel: 9.330152e-15) he: ====== time 1.900000e-01 (dt 1.000000e-02, wait 0, step 20 of 20) ===== he: updating variables... he: ...done he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 1.279447e-04 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.02 [s] he: nls: iter: 1, residual: 1.268220e-18 (rel: 9.912252e-15) he: ====== time 2.000000e-01 (dt 1.000000e-02, wait 0, step 21 of 21) ===== he: updating variables... he: ...done he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 1.204621e-04 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.02 [s] he: matrix: 0.02 [s] he: nls: iter: 1, residual: 1.326347e-18 (rel: 1.101050e-14) he: ====== time 2.100000e-01 (dt 1.000000e-02, wait 0, step 22 of 22) ===== he: updating variables... he: ...done he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 1.130582e-04 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.02 [s] he: nls: iter: 1, residual: 1.159744e-18 (rel: 1.025794e-14) he: ====== time 2.200000e-01 (dt 1.000000e-02, wait 0, step 23 of 23) ===== he: updating variables... he: ...done he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 1.058209e-04 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.02 [s] he: nls: iter: 1, residual: 1.244014e-18 (rel: 1.175584e-14) he: ====== time 2.300000e-01 (dt 1.000000e-02, wait 0, step 24 of 24) ===== he: updating variables... he: ...done he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 9.885933e-05 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.02 [s] he: nls: iter: 1, residual: 1.212984e-18 (rel: 1.226980e-14) he: ====== time 2.400000e-01 (dt 1.000000e-02, wait 0, step 25 of 25) ===== he: updating variables... he: ...done he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 9.230605e-05 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.02 [s] he: nls: iter: 1, residual: 1.525124e-18 (rel: 1.652247e-14) he: ====== time 2.500000e-01 (dt 1.000000e-02, wait 0, step 26 of 26) ===== he: updating variables... he: ...done he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 8.631914e-05 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.03 [s] he: nls: iter: 1, residual: 1.280901e-18 (rel: 1.483913e-14) he: ====== time 2.600000e-01 (dt 1.000000e-02, wait 0, step 27 of 27) ===== he: updating variables... he: ...done he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 8.108047e-05 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.02 [s] he: nls: iter: 1, residual: 1.375438e-18 (rel: 1.696386e-14) he: ====== time 2.700000e-01 (dt 1.000000e-02, wait 0, step 28 of 28) ===== he: updating variables... he: ...done he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 7.678740e-05 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.02 [s] he: nls: iter: 1, residual: 1.449311e-18 (rel: 1.887433e-14) he: ====== time 2.800000e-01 (dt 1.000000e-02, wait 0, step 29 of 29) ===== he: updating variables... he: ...done he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 7.363404e-05 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.02 [s] he: nls: iter: 1, residual: 1.501469e-18 (rel: 2.039096e-14) he: ====== time 2.900000e-01 (dt 1.000000e-02, wait 0, step 30 of 30) ===== he: updating variables... he: ...done he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 7.178120e-05 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.02 [s] he: nls: iter: 1, residual: 1.509681e-18 (rel: 2.103171e-14) he: ====== time 3.000000e-01 (dt 1.000000e-02, wait 0, step 31 of 31) ===== he: updating variables... he: ...done he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 7.132130e-05 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.02 [s] he: nls: iter: 1, residual: 1.535769e-18 (rel: 2.153310e-14) he: ====== time 3.100000e-01 (dt 1.000000e-02, wait 0, step 32 of 32) ===== he: updating variables... he: ...done he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 7.225266e-05 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.02 [s] he: nls: iter: 1, residual: 1.365998e-18 (rel: 1.890585e-14) he: ====== time 3.200000e-01 (dt 1.000000e-02, wait 0, step 33 of 33) ===== he: updating variables... he: ...done he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 7.447709e-05 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.02 [s] he: nls: iter: 1, residual: 1.505301e-18 (rel: 2.021159e-14) he: ====== time 3.300000e-01 (dt 1.000000e-02, wait 0, step 34 of 34) ===== he: updating variables... he: ...done he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 7.782246e-05 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.02 [s] he: nls: iter: 1, residual: 1.492796e-18 (rel: 1.918208e-14) he: ====== time 3.400000e-01 (dt 1.000000e-02, wait 0, step 35 of 35) ===== he: updating variables... he: ...done he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 8.207800e-05 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.02 [s] he: nls: iter: 1, residual: 1.650069e-18 (rel: 2.010368e-14) he: ====== time 3.500000e-01 (dt 1.000000e-02, wait 0, step 36 of 36) ===== he: updating variables... he: ...done he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 8.702690e-05 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.02 [s] he: nls: iter: 1, residual: 1.448831e-18 (rel: 1.664808e-14) he: ====== time 3.600000e-01 (dt 1.000000e-02, wait 0, step 37 of 37) ===== he: updating variables... he: ...done he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 9.246750e-05 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.02 [s] he: nls: iter: 1, residual: 1.368526e-18 (rel: 1.480008e-14) he: ====== time 3.700000e-01 (dt 1.000000e-02, wait 0, step 38 of 38) ===== he: updating variables... he: ...done he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 9.822306e-05 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.02 [s] he: nls: iter: 1, residual: 1.492923e-18 (rel: 1.519931e-14) he: ====== time 3.800000e-01 (dt 1.000000e-02, wait 0, step 39 of 39) ===== he: updating variables... he: ...done he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 1.041437e-04 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.02 [s] he: nls: iter: 1, residual: 1.560039e-18 (rel: 1.497968e-14) he: ====== time 3.900000e-01 (dt 1.000000e-02, wait 0, step 40 of 40) ===== he: updating variables... he: ...done he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 1.101040e-04 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.02 [s] he: nls: iter: 1, residual: 1.514577e-18 (rel: 1.375587e-14) he: ====== time 4.000000e-01 (dt 1.000000e-02, wait 0, step 41 of 41) ===== he: updating variables... he: ...done he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 1.159999e-04 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.02 [s] he: nls: iter: 1, residual: 1.583073e-18 (rel: 1.364719e-14) he: ====== time 4.100000e-01 (dt 1.000000e-02, wait 0, step 42 of 42) ===== he: updating variables... he: ...done he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 1.217444e-04 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.02 [s] he: nls: iter: 1, residual: 1.537605e-18 (rel: 1.262978e-14) he: ====== time 4.200000e-01 (dt 1.000000e-02, wait 0, step 43 of 43) ===== he: updating variables... he: ...done he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 1.272644e-04 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.02 [s] he: nls: iter: 1, residual: 1.544117e-18 (rel: 1.213314e-14) he: ====== time 4.300000e-01 (dt 1.000000e-02, wait 0, step 44 of 44) ===== he: updating variables... he: ...done he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 1.324984e-04 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.02 [s] he: nls: iter: 1, residual: 1.524295e-18 (rel: 1.150426e-14) he: ====== time 4.400000e-01 (dt 1.000000e-02, wait 0, step 45 of 45) ===== he: updating variables... he: ...done he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 1.373940e-04 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.02 [s] he: nls: iter: 1, residual: 1.512277e-18 (rel: 1.100687e-14) he: ====== time 4.500000e-01 (dt 1.000000e-02, wait 0, step 46 of 46) ===== he: updating variables... he: ...done he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 1.419066e-04 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.02 [s] he: nls: iter: 1, residual: 1.472613e-18 (rel: 1.037734e-14) he: ====== time 4.600000e-01 (dt 1.000000e-02, wait 0, step 47 of 47) ===== he: updating variables... he: ...done he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 1.459981e-04 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.02 [s] he: matrix: 0.02 [s] he: nls: iter: 1, residual: 1.540609e-18 (rel: 1.055225e-14) he: ====== time 4.700000e-01 (dt 1.000000e-02, wait 0, step 48 of 48) ===== he: updating variables... he: ...done he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 1.496359e-04 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.02 [s] he: nls: iter: 1, residual: 1.419950e-18 (rel: 9.489363e-15) he: ====== time 4.800000e-01 (dt 1.000000e-02, wait 0, step 49 of 49) ===== he: updating variables... he: ...done he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 1.527926e-04 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.02 [s] he: nls: iter: 1, residual: 1.494727e-18 (rel: 9.782715e-15) he: ====== time 4.900000e-01 (dt 1.000000e-02, wait 0, step 50 of 50) ===== he: updating variables... he: ...done he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 1.554450e-04 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.02 [s] he: nls: iter: 1, residual: 1.566471e-18 (rel: 1.007733e-14) he: ====== time 5.000000e-01 (dt 1.000000e-02, wait 0, step 51 of 51) ===== he: updating variables... he: ...done he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 1.575740e-04 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.02 [s] he: matrix: 0.02 [s] he: nls: iter: 1, residual: 1.573126e-18 (rel: 9.983409e-15) he: ====== time 6.000000e-01 (dt 1.000000e-01, wait 0, step 52 of 52) ===== he: updating variables... he: ...done he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 1.154647e-03 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.03 [s] he: matrix: 0.02 [s] he: nls: iter: 1, residual: 2.580808e-18 (rel: 2.235148e-15) he: ====== time 7.000000e-01 (dt 1.000000e-01, wait 0, step 53 of 53) ===== he: updating variables... he: ...done he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 7.257446e-04 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.02 [s] he: nls: iter: 1, residual: 1.919508e-18 (rel: 2.644881e-15) he: ====== time 8.000000e-01 (dt 1.000000e-01, wait 0, step 54 of 54) ===== he: updating variables... he: ...done he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 7.668834e-05 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.02 [s] he: nls: iter: 1, residual: 1.623094e-18 (rel: 2.116481e-14) he: ====== time 9.000000e-01 (dt 1.000000e-01, wait 0, step 55 of 55) ===== he: updating variables... he: ...done he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 6.965644e-04 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.02 [s] he: nls: iter: 1, residual: 1.493774e-18 (rel: 2.144488e-15) he: ====== time 1.000000e+00 (dt 1.000000e-01, wait 0, step 56 of 56) ===== he: updating variables... he: ...done he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 1.140586e-03 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.01 [s] he: matrix: 0.02 [s] he: nls: iter: 1, residual: 1.850391e-18 (rel: 1.622316e-15) [ 0. 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.1 0.11 0.12 0.13 0.14 0.15 0.16 0.17 0.18 0.19 0.2 0.21 0.22 0.23 0.24 0.25 0.26 0.27 0.28 0.29 0.3 0.31 0.32 0.33 0.34 0.35 0.36 0.37 0.38 0.39 0.4 0.41 0.42 0.43 0.44 0.45 0.46 0.47 0.48 0.49 0.5 0.6 0.7 0.8 0.9 1. ] ... ../examples/linear_elasticity/linear_elastic_damping.py solved +++ test_input: ok >>> all passed in 4.42 s <<< c:\Python27\lib\site-packages\sfepy\tests\test_input_linear_elastic_iga.py he: left over: ['verbose', '__builtins__', '_filename', '__file__', '__package__', 'TestInput', 'input_name', '__name__', 'output_name', '__doc__'] he: left over: ['verbose', '__builtins__', '__doc__', '__name__', 'data_dir', '__package__', 'stiffness_from_lame', '_filename', 'get_ebcs', '__file__'] >>> test instance prepared (1 test(s)) ... solving ../examples/linear_elasticity/linear_elastic_iga.py... he: creating regions... he: Omega he: Gamma2 he: Gamma1 he: ...done in 0.01 s he: equation "balance_of_forces": he: dw_lin_elastic.i.Omega(solid.D, v, u) = 0 he: using solvers: ts: no ts nls: newton ls: ls he: updating variables... he: ...done he: setting up dof connectivities... he: ...done in 0.00 s he: matrix shape: (720L, 720L) he: assembling matrix graph... he: ...done in 0.06 s he: matrix structural nonzeros: 133056 (2.57e-01% fill) he: updating materials... he: solid he: ...done in 0.00 s he: nls: iter: 0, residual: 1.904810e-01 (rel: 1.000000e+00) he: rezidual: 0.00 [s] he: solve: 0.02 [s] he: matrix: 0.08 [s] he: nls: iter: 1, residual: 3.768813e-16 (rel: 1.978577e-15) ... ../examples/linear_elasticity/linear_elastic_iga.py solved +++ test_input: ok >>> all passed in 0.60 s <<< c:\Python27\lib\site-packages\sfepy\tests\test_input_linear_elastic_mM.py he: left over: ['verbose', '__builtins__', '_filename', '__file__', '__package__', 'TestInput', 'input_name', '__name__', 'output_name', '__doc__'] he: left over: ['get_homog_coefs_linear', 'get_elements', '__builtins__', '__file__', '__doc__', '__name__', 'verbose', 'data_dir', 'nm', 'recover_micro_hook', '__package__', 'save_recovery_region', '_filename', 'post_process', 'os', 'micro_filename'] >>> test instance prepared (1 test(s)) ... solving ../examples/homogenization/linear_elastic_mM.py... macro: reading mesh [line2, tri3, quad4, tetra4, hexa8] (c:\Python27\lib\site-packages\sfepy\meshes\3d\cylinder.mesh)... macro: ...done in 0.01 s macro: creating regions... macro: Right macro: Omega macro: Recovery macro: Left macro: ...done in 0.01 s macro: equation "balance_of_forces": macro: dw_lin_elastic.i.Omega( solid.D, v, u ) = 0 macro: using solvers: ts: no ts nls: newton ls: ls macro: updating variables... macro: ...done macro: setting up dof connectivities... macro: ...done in 0.00 s macro: matrix shape: (900L, 900L) macro: assembling matrix graph... macro: ...done in 0.01 s macro: matrix structural nonzeros: 31842 (3.93e-02% fill) macro: updating materials... macro: solid macro: ...done in 0.02 s macro: nls: iter: 0, residual: 2.204111e+07 (rel: 1.000000e+00) macro: rezidual: 0.00 [s] macro: solve: 0.01 [s] macro: matrix: 0.01 [s] macro: nls: iter: 1, residual: 3.436164e-08 (rel: 1.558979e-15) macro: equation "tmp": macro: ev_cauchy_stress.i.Omega( solid.D, u ) macro: updating materials... macro: solid macro: ...done in 0.02 s macro: equation "tmp": macro: ev_cauchy_strain.i.Omega( u ) macro: updating materials... macro: ...done in 0.00 s macro: equation "tmp": macro: ev_cauchy_strain.i.Recovery( u ) macro: updating materials... macro: ...done in 0.00 s macro: left over: ['compute_micro_u', 'recovery_le', 'Struct', 'cb', 'compute_stress_strain_u', 'verbose', 'compute_mac_stress_part', 'requirements', 'nm', '_filename', 'data_dir', 'per', '__package__', 'add_stress_p', 'get_box_volume', 'set_elastic_u', 'expr_coefs', '__doc__', 'all_periodic', 'stiffness_from_youngpoisson_mixed', 'dim', '__builtins__', '__file__', 'equation_corrs', '__name__', 'region_lbn', 'ii', 'region_rtf', 'bulk_from_youngpoisson', 'define_box_regions', 'coefs'] macro: reading mesh [line2, tri3, quad4, tetra4, hexa8] (c:\Python27\lib\site-packages\sfepy\meshes\3d\matrix_fiber.mesh)... macro: ...done in 0.02 s macro: creating regions... macro: Right macro: Ym macro: Bottom macro: Far macro: Corners macro: Yc macro: Near macro: Y macro: Top macro: Left macro: ...done in 0.03 s ii: 0, iel: 50 ii: 1, iel: 150 ii: 2, iel: 250 ii: 3, iel: 350 ii: 4, iel: 450 ii: 5, iel: 550 ii: 6, iel: 650 ii: 7, iel: 750 ii: 8, iel: 850 ii: 9, iel: 950 ii: 10, iel: 1050 ii: 11, iel: 1150 ii: 12, iel: 1250 ... ../examples/homogenization/linear_elastic_mM.py solved +++ test_input: ok >>> all passed in 3.90 s <<< c:\Python27\lib\site-packages\sfepy\tests\test_input_linear_elastic_probes.py macro: left over: ['verbose', '__builtins__', '_filename', '__file__', '__package__', 'TestInput', 'input_name', '__name__', 'output_name', '__doc__'] macro: left over: ['lam', 'plt', 'Struct', '__builtins__', '__file__', 'deepcopy', 'fm', 'verbose', 'par_fun', 'mu', 'data_dir', 'nm', 'solid', '__name__', 'Probe', '__package__', 'stiffness_from_lame', '_filename', 'post_process', 'os', '__doc__'] >>> test instance prepared (1 test(s)) ... solving ../examples/linear_elasticity/linear_elastic_probes.py... macro: reading mesh [line2, tri3, quad4, tetra4, hexa8] (c:\Python27\lib\site-packages\sfepy\meshes\3d\cylinder.mesh)... macro: ...done in 0.01 s macro: creating regions... macro: SomewhereTop macro: Right macro: Omega macro: Left macro: ...done in 0.01 s macro: equation "balance_of_forces": macro: dw_lin_elastic_iso.i.Omega(solid.lam, solid.mu, v, u) = 0 macro: using solvers: ts: no ts nls: newton ls: ls macro: updating variables... macro: ...done macro: setting up dof connectivities... macro: ...done in 0.00 s macro: matrix shape: (886L, 886L) macro: assembling matrix graph... macro: ...done in 0.01 s macro: matrix structural nonzeros: 31052 (3.96e-02% fill) macro: updating materials... macro: solid macro: ...done in 0.00 s macro: nls: iter: 0, residual: 4.778054e-03 (rel: 1.000000e+00) macro: rezidual: 0.00 [s] macro: solve: 0.01 [s] macro: matrix: 0.01 [s] macro: nls: iter: 1, residual: 7.193990e-18 (rel: 1.505632e-15) macro: equation "tmp": macro: ev_cauchy_strain.i.Omega( u ) macro: updating materials... macro: ...done in 0.00 s macro: equation "tmp": macro: ev_cauchy_stress.i.Omega( solid.D, u ) macro: updating materials... macro: solid macro: ...done in 0.00 s >>> Traceback (most recent call last): File "run_tests.py", line 206, in main() File "run_tests.py", line 200, in main op.walk(options.test_dir, run_tests, stats) File "c:\Python27\lib\ntpath.py", line 268, in walk func(arg, top, names) File "run_tests.py", line 131, in run_tests n_fail, n_total, test_time = run_test(conf_name, options) File "run_tests.py", line 99, in run_test ok, n_fail, n_total = test.run(options.debug) File "c:\Python27\lib\site-packages\sfepy\base\testing.py", line 58, in run ret = test_method() File "c:\Python27\lib\site-packages\sfepy\tests\tests_basic.py", line 73, in t est_input post_process_hook_final=self.post_process_hook_final) File "c:\Python27\lib\site-packages\sfepy\applications\pde_solver_app.py", lin e 55, in solve_pde return app(nls_status=nls_status) File "c:\Python27\lib\site-packages\sfepy\applications\application.py", line 2 9, in call_basic return self.call(**kwargs) File "c:\Python27\lib\site-packages\sfepy\applications\pde_solver_app.py", lin e 219, in call post_process_hook=self.post_process_hook): File "c:\Python27\lib\site-packages\sfepy\solvers\ts_solvers.py", line 37, in __call__ file_per_var=None) File "c:\Python27\lib\site-packages\sfepy\discrete\problem.py", line 736, in s ave_state out = post_process_hook(out, self, state, extend=extend) File "c:\Python27\lib\site-packages\sfepy\examples\linear_elasticity\linear_el astic_probes.py", line 99, in post_process probe = Probe(out, mesh) File "c:\Python27\lib\site-packages\sfepy\postprocess\probes_vtk.py", line 36, in __init__ self.vtkprobe.SetSource(self.vtkdata) AttributeError: SetSource