NUTILS 8.1 "Idiyappam"
arguments
nrefine=1
poisson=0.3
diffusivity=0.01
thermal_expansion=0.01
timestep=0.1
endtime=0.2
start 2023-08-15 02:38:11
project
projection type: lsqr
solve
solving 128 dof system to machine precision using arnoldi solver
constructing direct preconditioner
peak memory use 688k
residual decreased by 15.5 orders using 0 krylov vectors
residual decreased by 0.4 orders using 1 krylov vectors
solver returned with residual 5e-17
constrained 128/128 dofs, error 2.18e-02/area
optimize
solve
solving 96 dof system to machine precision using arnoldi solver
constructing symmetric direct preconditioner
peak memory use 296k
residual decreased by 15.7 orders using 0 krylov vectors
residual decreased by 0.2 orders using 1 krylov vectors
solver returned with residual 6e-15
constrained 96/128 dofs
optimum value 3.08e+00
project
projection type: lsqr
solve
solving 96 dof system to machine precision using arnoldi solver
constructing direct preconditioner
peak memory use 332k
residual decreased by 15.8 orders using 0 krylov vectors
residual decreased by 0.2 orders using 1 krylov vectors
solver returned with residual 3e-15
constrained 96/128 dofs, error 4.47e-02/area
project
projection type: lsqr
constrained 128/128 dofs, error 0.00e+00/area
timestep 1
newton
iter 1 (0%)
solve
solving 32 dof system to tolerance 9e-04 using arnoldi solver
constructing direct preconditioner
peak memory use 187k
residual decreased by 15.9 orders using 0 krylov vectors
solver returned with residual 1e-16
estimated residual minimum at 100% of update vector
update accepted at relaxation 1.0
converged in 1 iterations to residual 1.2e-15
timestep 2
newton
iter 1 (0%)
solve
solving 32 dof system to tolerance 8e-04 using arnoldi solver
constructing direct preconditioner
peak memory use 187k
residual decreased by 15.9 orders using 0 krylov vectors
solver returned with residual 1e-16
estimated residual minimum at 100% of update vector
update accepted at relaxation 1.0
converged in 1 iterations to residual 1.8e-15
optimize
constrained 24/384 dofs
optimum value 0.00e+00
solve
solving 360 dof system to machine precision using arnoldi solver
constructing direct preconditioner
peak memory use 4,088k
residual decreased by 15.0 orders using 0 krylov vectors
residual decreased by 0.1 orders using 1 krylov vectors
solver returned with residual 4e-16
NutilsDeprecationWarning: function.norm2 is deprecated; use numpy.linalg.norm instead In /usr/local/lib/python3.11/dist-packages/nutils/function.py:1226
deformed_cylinder.vtk
MatplotlibDeprecationWarning: Unable to determine Axes to steal space for Colorbar. Using gca(), but will raise in the future. Either provide the *cax* argument to use as the Axes for the Colorbar, provide the *ax* argument to steal space from it, or add *mappable* to an Axes. In /app/thermoelastic-example.py:217
displacement.png
finish 2023-08-15 02:38:18, elapsed 0:06