Load up necessary packages
using Cthoniosusing Meshes
import Meshes: SimpleMesh, vizFile management
mesh_file = "../../../examples/hole_array/hole_array.exo""../../../examples/hole_array/hole_array.exo"Functions for BCs
func_1(x, t) = -5. * t#(0., -5. * t)
func_2(x, t) = 0.0
func_3(x, t) = @SVector [0., -0.025 * t]func_3 (generic function with 1 method)Boundary conditions
disp_bcs = [
DirichletBC("yminus_nodeset", [1, 2], func_2)
DirichletBC("yplus_nodeset", [1], func_2)
DirichletBC("yplus_nodeset", [2], func_1)
]
traction_bcs = [
]Any[]Sections
sections = Section[
Section(
"Block1", 2,
SolidMechanics(NeoHookean(), PlaneStrain()),
MaterialProperties(
"bulk modulus" => 0.833,
"shear modulus" => 1.0
)
)
]1-element Vector{Section}:
Section{SolidMechanics{2, 2, 0, NeoHookean, PlaneStrain, Cthonios.StressDivergence{2, 2, 0, NeoHookean, PlaneStrain}}, MaterialProperties{Dict{Symbol, Any}}}("Block1", 2, SolidMechanics{2, 2, 0, NeoHookean, PlaneStrain, Cthonios.StressDivergence{2, 2, 0, NeoHookean, PlaneStrain}}(NeoHookean(), PlaneStrain(), Cthonios.StressDivergence{2, 2, 0, NeoHookean, PlaneStrain}(NeoHookean(), PlaneStrain())), MaterialProperties{Dict{Symbol, Any}}(Dict{Symbol, Any}(Symbol("bulk modulus") => 0.833, Symbol("shear modulus") => 1.0)))Domain setup
domain = Domain(mesh_file, sections, disp_bcs, traction_bcs)Cthonios.Domain{FiniteElementContainers.VectorizedNodalField{Float64, 2, 2, 8416, Vector{Float64}}, DofManager{Int64, 2, 8416, Vector{Float64}, Vector{Int64}}, @NamedTuple{Block1::Cthonios.SectionInternal{NonAllocatedFunctionSpace{2, Vector{Int64}, FiniteElementContainers.VectorizedElementField{Int64, 2, 3, 14768, Vector{Int64}}, FiniteElementContainers.VectorizedElementField{Int64, 2, 6, 14768, Vector{Int64}}, ReferenceFiniteElements.ReferenceFE{Int64, Float64, ReferenceFiniteElements.Tri3{SpecialPolynomials.Lagrange, 2}, ReferenceFiniteElements.Edge{2, SpecialPolynomials.Lagrange, 1, 2}, ReferenceFiniteElements.ArrayBackend{SArray}, Vector{SVector{2, Int64}}, Vector{SVector{0, Union{}}}, Vector{Int64}, Vector{SVector{2, Float64}}, ReferenceFiniteElements.CellInterpolants{StructArrays.StructVector{ReferenceFiniteElements.Interpolants, @NamedTuple{ξ::Vector{SVector{2, Float64}}, w::Vector{Float64}, N::Vector{SVector{3, Float64}}, ∇N_ξ::Vector{SMatrix{3, 2, Float64, 6}}, ∇∇N_ξ::Vector{SArray{Tuple{3, 2, 2}, Float64, 3, 12}}}, Int64}}, Vector{Vector{SVector{2, Float64}}}, ReferenceFiniteElements.SurfaceInterpolants{StructArrays.StructArray{ReferenceFiniteElements.Interpolants, 2, @NamedTuple{ξ::Matrix{SVector{2, Float64}}, w::Matrix{Float64}, N::Matrix{SVector{3, Float64}}, ∇N_ξ::Matrix{SMatrix{3, 2, Float64, 6}}, ∇∇N_ξ::Matrix{SArray{Tuple{3, 2, 2}, Float64, 3, 12}}}, Int64}}}}, SolidMechanics{2, 2, 0, NeoHookean, PlaneStrain, Cthonios.StressDivergence{2, 2, 0, NeoHookean, PlaneStrain}}, SVector{2, Float64}}}, @NamedTuple{yplus_nodeset_2::Cthonios.DirichletBCInternal{Vector{Int32}, Vector{Int64}, typeof(Main.func_1)}, yminus_nodeset_1_2::Cthonios.DirichletBCInternal{Vector{Int32}, Vector{Int64}, typeof(Main.func_2)}, yplus_nodeset_1::Cthonios.DirichletBCInternal{Vector{Int32}, Vector{Int64}, typeof(Main.func_2)}}, Vector{Int64}, @NamedTuple{}, @NamedTuple{}}([44.29728315483824 44.09375 … 40.010176197600856 20.023250115402412; 50.65587225363096 51.0 … 30.205590373069363 30.2388316572343], DofManager
Number of nodes = 8416
Number of dofs per node = 2
Storage type = Vector{Int64}, (Block1 = Cthonios.SectionInternal{NonAllocatedFunctionSpace{2, Vector{Int64}, FiniteElementContainers.VectorizedElementField{Int64, 2, 3, 14768, Vector{Int64}}, FiniteElementContainers.VectorizedElementField{Int64, 2, 6, 14768, Vector{Int64}}, ReferenceFiniteElements.ReferenceFE{Int64, Float64, ReferenceFiniteElements.Tri3{SpecialPolynomials.Lagrange, 2}, ReferenceFiniteElements.Edge{2, SpecialPolynomials.Lagrange, 1, 2}, ReferenceFiniteElements.ArrayBackend{SArray}, Vector{SVector{2, Int64}}, Vector{SVector{0, Union{}}}, Vector{Int64}, Vector{SVector{2, Float64}}, ReferenceFiniteElements.CellInterpolants{StructArrays.StructVector{ReferenceFiniteElements.Interpolants, @NamedTuple{ξ::Vector{SVector{2, Float64}}, w::Vector{Float64}, N::Vector{SVector{3, Float64}}, ∇N_ξ::Vector{SMatrix{3, 2, Float64, 6}}, ∇∇N_ξ::Vector{SArray{Tuple{3, 2, 2}, Float64, 3, 12}}}, Int64}}, Vector{Vector{SVector{2, Float64}}}, ReferenceFiniteElements.SurfaceInterpolants{StructArrays.StructArray{ReferenceFiniteElements.Interpolants, 2, @NamedTuple{ξ::Matrix{SVector{2, Float64}}, w::Matrix{Float64}, N::Matrix{SVector{3, Float64}}, ∇N_ξ::Matrix{SMatrix{3, 2, Float64, 6}}, ∇∇N_ξ::Matrix{SArray{Tuple{3, 2, 2}, Float64, 3, 12}}}, Int64}}}}, SolidMechanics{2, 2, 0, NeoHookean, PlaneStrain, Cthonios.StressDivergence{2, 2, 0, NeoHookean, PlaneStrain}}, SVector{2, Float64}}(1, 2, NonAllocatedFunctionSpace
Reference finite element = ReferenceFE:
Element Type = ReferenceFiniteElements.Tri3{SpecialPolynomials.Lagrange, 2}
Surface Element Type = ReferenceFiniteElements.Edge{2, SpecialPolynomials.Lagrange, 1, 2}
Interpolation Type = SpecialPolynomials.Lagrange
Polynomial Degree = 1
Quadrature Degree = 2
Array Backend Type = ReferenceFiniteElements.ArrayBackend{SArray}()
Edge Nodes Storage Type = Vector{SVector{2, Int64}}
Face Nodes Storage Type = Vector{SVector{0, Union{}}}
Interior Nodes Storage Type = Vector{Int64}
Nodal Coordinates Storage Type = Vector{SVector{2, Float64}}
CellInterpolants:
Quadrature Point Storage Type = Vector{SVector{2, Float64}}
Quadrature Weights Storage Type = Vector{Float64}
Shape Function Values Storage Type = Vector{SVector{3, Float64}}
Shape Function Gradients Storage Type = Vector{SMatrix{3, 2, Float64, 6}}
Shape Function Hessians Storage Type = Vector{SArray{Tuple{3, 2, 2}, Float64, 3, 12}}
Surface Nodal Coordinates Storage Type = Vector{Vector{SVector{2, Float64}}}
CellInterpolants:
Quadrature Point Storage Type = Matrix{SVector{2, Float64}}
Quadrature Weights Storage Type = Matrix{Float64}
Shape Function Values Storage Type = Matrix{SVector{3, Float64}}
Shape Function Gradients Storage Type = Matrix{SMatrix{3, 2, Float64, 6}}
Shape Function Hessians Storage Type = Matrix{SArray{Tuple{3, 2, 2}, Float64, 3, 12}}
, SolidMechanics{2, 2, 0, NeoHookean, PlaneStrain, Cthonios.StressDivergence{2, 2, 0, NeoHookean, PlaneStrain}}(NeoHookean(), PlaneStrain(), Cthonios.StressDivergence{2, 2, 0, NeoHookean, PlaneStrain}(NeoHookean(), PlaneStrain())), [0.833, 1.0]),), (yplus_nodeset_2 = DirichletBCInternal:
Numbe of active dofs = 129
Function = func_1
, yminus_nodeset_1_2 = DirichletBCInternal:
Numbe of active dofs = 258
Function = func_2
, yplus_nodeset_1 = DirichletBCInternal:
Numbe of active dofs = 129
Function = func_2
), Int64[], NamedTuple(), NamedTuple())Objective setup
objective = Objective(domain, Cthonios.energy)Objective{Cthonios.Domain{FiniteElementContainers.VectorizedNodalField{Float64, 2, 2, 8416, Vector{Float64}}, DofManager{Int64, 2, 8416, Vector{Float64}, Vector{Int64}}, @NamedTuple{Block1::Cthonios.SectionInternal{NonAllocatedFunctionSpace{2, Vector{Int64}, FiniteElementContainers.VectorizedElementField{Int64, 2, 3, 14768, Vector{Int64}}, FiniteElementContainers.VectorizedElementField{Int64, 2, 6, 14768, Vector{Int64}}, ReferenceFiniteElements.ReferenceFE{Int64, Float64, ReferenceFiniteElements.Tri3{SpecialPolynomials.Lagrange, 2}, ReferenceFiniteElements.Edge{2, SpecialPolynomials.Lagrange, 1, 2}, ReferenceFiniteElements.ArrayBackend{SArray}, Vector{SVector{2, Int64}}, Vector{SVector{0, Union{}}}, Vector{Int64}, Vector{SVector{2, Float64}}, ReferenceFiniteElements.CellInterpolants{StructArrays.StructVector{ReferenceFiniteElements.Interpolants, @NamedTuple{ξ::Vector{SVector{2, Float64}}, w::Vector{Float64}, N::Vector{SVector{3, Float64}}, ∇N_ξ::Vector{SMatrix{3, 2, Float64, 6}}, ∇∇N_ξ::Vector{SArray{Tuple{3, 2, 2}, Float64, 3, 12}}}, Int64}}, Vector{Vector{SVector{2, Float64}}}, ReferenceFiniteElements.SurfaceInterpolants{StructArrays.StructArray{ReferenceFiniteElements.Interpolants, 2, @NamedTuple{ξ::Matrix{SVector{2, Float64}}, w::Matrix{Float64}, N::Matrix{SVector{3, Float64}}, ∇N_ξ::Matrix{SMatrix{3, 2, Float64, 6}}, ∇∇N_ξ::Matrix{SArray{Tuple{3, 2, 2}, Float64, 3, 12}}}, Int64}}}}, SolidMechanics{2, 2, 0, NeoHookean, PlaneStrain, Cthonios.StressDivergence{2, 2, 0, NeoHookean, PlaneStrain}}, SVector{2, Float64}}}, @NamedTuple{yplus_nodeset_2::Cthonios.DirichletBCInternal{Vector{Int32}, Vector{Int64}, typeof(Main.func_1)}, yminus_nodeset_1_2::Cthonios.DirichletBCInternal{Vector{Int32}, Vector{Int64}, typeof(Main.func_2)}, yplus_nodeset_1::Cthonios.DirichletBCInternal{Vector{Int32}, Vector{Int64}, typeof(Main.func_2)}}, Vector{Int64}, @NamedTuple{}, @NamedTuple{}}, typeof(Cthonios.energy), typeof(Cthonios.gradient), typeof(Cthonios.hessian), typeof(Cthonios.neumann_energy), typeof(Cthonios.neumann_gradient), typeof(Cthonios.neumann_hessian), TimerOutput}(Cthonios.Domain{FiniteElementContainers.VectorizedNodalField{Float64, 2, 2, 8416, Vector{Float64}}, DofManager{Int64, 2, 8416, Vector{Float64}, Vector{Int64}}, @NamedTuple{Block1::Cthonios.SectionInternal{NonAllocatedFunctionSpace{2, Vector{Int64}, FiniteElementContainers.VectorizedElementField{Int64, 2, 3, 14768, Vector{Int64}}, FiniteElementContainers.VectorizedElementField{Int64, 2, 6, 14768, Vector{Int64}}, ReferenceFiniteElements.ReferenceFE{Int64, Float64, ReferenceFiniteElements.Tri3{SpecialPolynomials.Lagrange, 2}, ReferenceFiniteElements.Edge{2, SpecialPolynomials.Lagrange, 1, 2}, ReferenceFiniteElements.ArrayBackend{SArray}, Vector{SVector{2, Int64}}, Vector{SVector{0, Union{}}}, Vector{Int64}, Vector{SVector{2, Float64}}, ReferenceFiniteElements.CellInterpolants{StructArrays.StructVector{ReferenceFiniteElements.Interpolants, @NamedTuple{ξ::Vector{SVector{2, Float64}}, w::Vector{Float64}, N::Vector{SVector{3, Float64}}, ∇N_ξ::Vector{SMatrix{3, 2, Float64, 6}}, ∇∇N_ξ::Vector{SArray{Tuple{3, 2, 2}, Float64, 3, 12}}}, Int64}}, Vector{Vector{SVector{2, Float64}}}, ReferenceFiniteElements.SurfaceInterpolants{StructArrays.StructArray{ReferenceFiniteElements.Interpolants, 2, @NamedTuple{ξ::Matrix{SVector{2, Float64}}, w::Matrix{Float64}, N::Matrix{SVector{3, Float64}}, ∇N_ξ::Matrix{SMatrix{3, 2, Float64, 6}}, ∇∇N_ξ::Matrix{SArray{Tuple{3, 2, 2}, Float64, 3, 12}}}, Int64}}}}, SolidMechanics{2, 2, 0, NeoHookean, PlaneStrain, Cthonios.StressDivergence{2, 2, 0, NeoHookean, PlaneStrain}}, SVector{2, Float64}}}, @NamedTuple{yplus_nodeset_2::Cthonios.DirichletBCInternal{Vector{Int32}, Vector{Int64}, typeof(Main.func_1)}, yminus_nodeset_1_2::Cthonios.DirichletBCInternal{Vector{Int32}, Vector{Int64}, typeof(Main.func_2)}, yplus_nodeset_1::Cthonios.DirichletBCInternal{Vector{Int32}, Vector{Int64}, typeof(Main.func_2)}}, Vector{Int64}, @NamedTuple{}, @NamedTuple{}}([44.29728315483824 44.09375 … 40.010176197600856 20.023250115402412; 50.65587225363096 51.0 … 30.205590373069363 30.2388316572343], DofManager
Number of nodes = 8416
Number of dofs per node = 2
Storage type = Vector{Int64}, (Block1 = Cthonios.SectionInternal{NonAllocatedFunctionSpace{2, Vector{Int64}, FiniteElementContainers.VectorizedElementField{Int64, 2, 3, 14768, Vector{Int64}}, FiniteElementContainers.VectorizedElementField{Int64, 2, 6, 14768, Vector{Int64}}, ReferenceFiniteElements.ReferenceFE{Int64, Float64, ReferenceFiniteElements.Tri3{SpecialPolynomials.Lagrange, 2}, ReferenceFiniteElements.Edge{2, SpecialPolynomials.Lagrange, 1, 2}, ReferenceFiniteElements.ArrayBackend{SArray}, Vector{SVector{2, Int64}}, Vector{SVector{0, Union{}}}, Vector{Int64}, Vector{SVector{2, Float64}}, ReferenceFiniteElements.CellInterpolants{StructArrays.StructVector{ReferenceFiniteElements.Interpolants, @NamedTuple{ξ::Vector{SVector{2, Float64}}, w::Vector{Float64}, N::Vector{SVector{3, Float64}}, ∇N_ξ::Vector{SMatrix{3, 2, Float64, 6}}, ∇∇N_ξ::Vector{SArray{Tuple{3, 2, 2}, Float64, 3, 12}}}, Int64}}, Vector{Vector{SVector{2, Float64}}}, ReferenceFiniteElements.SurfaceInterpolants{StructArrays.StructArray{ReferenceFiniteElements.Interpolants, 2, @NamedTuple{ξ::Matrix{SVector{2, Float64}}, w::Matrix{Float64}, N::Matrix{SVector{3, Float64}}, ∇N_ξ::Matrix{SMatrix{3, 2, Float64, 6}}, ∇∇N_ξ::Matrix{SArray{Tuple{3, 2, 2}, Float64, 3, 12}}}, Int64}}}}, SolidMechanics{2, 2, 0, NeoHookean, PlaneStrain, Cthonios.StressDivergence{2, 2, 0, NeoHookean, PlaneStrain}}, SVector{2, Float64}}(1, 2, NonAllocatedFunctionSpace
Reference finite element = ReferenceFE:
Element Type = ReferenceFiniteElements.Tri3{SpecialPolynomials.Lagrange, 2}
Surface Element Type = ReferenceFiniteElements.Edge{2, SpecialPolynomials.Lagrange, 1, 2}
Interpolation Type = SpecialPolynomials.Lagrange
Polynomial Degree = 1
Quadrature Degree = 2
Array Backend Type = ReferenceFiniteElements.ArrayBackend{SArray}()
Edge Nodes Storage Type = Vector{SVector{2, Int64}}
Face Nodes Storage Type = Vector{SVector{0, Union{}}}
Interior Nodes Storage Type = Vector{Int64}
Nodal Coordinates Storage Type = Vector{SVector{2, Float64}}
CellInterpolants:
Quadrature Point Storage Type = Vector{SVector{2, Float64}}
Quadrature Weights Storage Type = Vector{Float64}
Shape Function Values Storage Type = Vector{SVector{3, Float64}}
Shape Function Gradients Storage Type = Vector{SMatrix{3, 2, Float64, 6}}
Shape Function Hessians Storage Type = Vector{SArray{Tuple{3, 2, 2}, Float64, 3, 12}}
Surface Nodal Coordinates Storage Type = Vector{Vector{SVector{2, Float64}}}
CellInterpolants:
Quadrature Point Storage Type = Matrix{SVector{2, Float64}}
Quadrature Weights Storage Type = Matrix{Float64}
Shape Function Values Storage Type = Matrix{SVector{3, Float64}}
Shape Function Gradients Storage Type = Matrix{SMatrix{3, 2, Float64, 6}}
Shape Function Hessians Storage Type = Matrix{SArray{Tuple{3, 2, 2}, Float64, 3, 12}}
, SolidMechanics{2, 2, 0, NeoHookean, PlaneStrain, Cthonios.StressDivergence{2, 2, 0, NeoHookean, PlaneStrain}}(NeoHookean(), PlaneStrain(), Cthonios.StressDivergence{2, 2, 0, NeoHookean, PlaneStrain}(NeoHookean(), PlaneStrain())), [0.833, 1.0]),), (yplus_nodeset_2 = DirichletBCInternal:
Numbe of active dofs = 129
Function = func_1
, yminus_nodeset_1_2 = DirichletBCInternal:
Numbe of active dofs = 258
Function = func_2
, yplus_nodeset_1 = DirichletBCInternal:
Numbe of active dofs = 129
Function = func_2
), Int64[], NamedTuple(), NamedTuple()), Cthonios.energy, Cthonios.gradient, Cthonios.hessian, Cthonios.neumann_energy, Cthonios.neumann_gradient, Cthonios.neumann_hessian, ────────────────────────────────────────────────────────────────────
Time Allocations
─────────────────────── ────────────────────────
Tot / % measured: 174ms / 0.0% 10.6MiB / 0.0%
Section ncalls time %tot avg alloc %tot avg
────────────────────────────────────────────────────────────────────
────────────────────────────────────────────────────────────────────)Integrator setup
integrator = QuasiStatic(0.0, 1.0, 1. / 20)QuasiStatic:
Start time = 0.0
End time = 1.0
Time step = [0.05]
Post-processor
pp = ExodusPostProcessor(mesh_file, "output.e", ["displ_x", "displ_y"])
problem = Problem(objective, integrator, TrustRegionSolver, pp)Problem{QuasiStatic{Vector{Int64}, Float64, Vector{Float64}}, TrustRegionSolver{CholeskyPreconditioner{FiniteElementContainers.StaticAssembler{Float64, Int64, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}, FiniteElementContainers.VectorizedNodalField{Float64, 2, 2, 8416, Vector{Float64}}, Vector{Float64}, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Float64}, Vector{Int64}, Vector{Int64}, Vector{Float64}}, SparseArrays.CHOLMOD.Factor{Float64, Int64}, TimerOutput}, Objective{Cthonios.Domain{FiniteElementContainers.VectorizedNodalField{Float64, 2, 2, 8416, Vector{Float64}}, DofManager{Int64, 2, 8416, Vector{Float64}, Vector{Int64}}, @NamedTuple{Block1::Cthonios.SectionInternal{NonAllocatedFunctionSpace{2, Vector{Int64}, FiniteElementContainers.VectorizedElementField{Int64, 2, 3, 14768, Vector{Int64}}, FiniteElementContainers.VectorizedElementField{Int64, 2, 6, 14768, Vector{Int64}}, ReferenceFiniteElements.ReferenceFE{Int64, Float64, ReferenceFiniteElements.Tri3{SpecialPolynomials.Lagrange, 2}, ReferenceFiniteElements.Edge{2, SpecialPolynomials.Lagrange, 1, 2}, ReferenceFiniteElements.ArrayBackend{SArray}, Vector{SVector{2, Int64}}, Vector{SVector{0, Union{}}}, Vector{Int64}, Vector{SVector{2, Float64}}, ReferenceFiniteElements.CellInterpolants{StructArrays.StructVector{ReferenceFiniteElements.Interpolants, @NamedTuple{ξ::Vector{SVector{2, Float64}}, w::Vector{Float64}, N::Vector{SVector{3, Float64}}, ∇N_ξ::Vector{SMatrix{3, 2, Float64, 6}}, ∇∇N_ξ::Vector{SArray{Tuple{3, 2, 2}, Float64, 3, 12}}}, Int64}}, Vector{Vector{SVector{2, Float64}}}, ReferenceFiniteElements.SurfaceInterpolants{StructArrays.StructArray{ReferenceFiniteElements.Interpolants, 2, @NamedTuple{ξ::Matrix{SVector{2, Float64}}, w::Matrix{Float64}, N::Matrix{SVector{3, Float64}}, ∇N_ξ::Matrix{SMatrix{3, 2, Float64, 6}}, ∇∇N_ξ::Matrix{SArray{Tuple{3, 2, 2}, Float64, 3, 12}}}, Int64}}}}, SolidMechanics{2, 2, 0, NeoHookean, PlaneStrain, Cthonios.StressDivergence{2, 2, 0, NeoHookean, PlaneStrain}}, SVector{2, Float64}}}, @NamedTuple{yplus_nodeset_2::Cthonios.DirichletBCInternal{Vector{Int32}, Vector{Int64}, typeof(Main.func_1)}, yminus_nodeset_1_2::Cthonios.DirichletBCInternal{Vector{Int32}, Vector{Int64}, typeof(Main.func_2)}, yplus_nodeset_1::Cthonios.DirichletBCInternal{Vector{Int32}, Vector{Int64}, typeof(Main.func_2)}}, Vector{Int64}, @NamedTuple{}, @NamedTuple{}}, typeof(Cthonios.energy), typeof(Cthonios.gradient), typeof(Cthonios.hessian), typeof(Cthonios.neumann_energy), typeof(Cthonios.neumann_gradient), typeof(Cthonios.neumann_hessian), TimerOutput}, Vector{Float64}, Cthonios.WarmStart{FiniteElementContainers.VectorizedNodalField{Float64, 2, 2, 8416, Vector{Float64}}, Vector{Float64}, ObjectiveParameters{FiniteElementContainers.VectorizedNodalField{Float64, 2, 2, 8416, Vector{Float64}}, QuasiStatic{Vector{Int64}, Float64, Vector{Float64}}, Vector{Float64}, Vector{SVector{2, Float64}}, ComponentArrays.ComponentVector{Float64, Vector{Float64}, Tuple{ComponentArrays.Axis{(Block1 = ViewAxis(1:0, ShapedAxis((0, 3, 14768))),)}}}, ComponentArrays.ComponentVector{Float64, Vector{Float64}, Tuple{ComponentArrays.Axis{(Block1 = 1:2,)}}}, FiniteElementContainers.VectorizedNodalField{Float64, 2, 2, 8416, Vector{Float64}}, FiniteElementContainers.VectorizedNodalField{Float64, 2, 2, 8416, Vector{Float64}}, Vector{Float64}, FiniteElementContainers.VectorizedNodalField{Float64, 2, 2, 8416, Vector{Float64}}, ComponentArrays.ComponentVector{Float64, Vector{Float64}, Tuple{ComponentArrays.Axis{(Block1 = ViewAxis(1:44304, ShapedAxis((3, 14768))),)}}}}}, TimerOutput, FiniteElementContainers.VectorizedNodalField{Float64, 2, 2, 8416, Vector{Float64}}, Cthonios.TrustRegionSolverSettings}, ExodusPostProcessor{ExodusDatabase{Int32, Int32, Int32, Float64, Initialization{2, 8416, 14768, 1, 2, 3}}}, TimerOutput}(QuasiStatic:
Start time = 0.0
End time = 1.0
Time step = [0.05]
, TrustRegionSolver{CholeskyPreconditioner{FiniteElementContainers.StaticAssembler{Float64, Int64, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}, FiniteElementContainers.VectorizedNodalField{Float64, 2, 2, 8416, Vector{Float64}}, Vector{Float64}, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Float64}, Vector{Int64}, Vector{Int64}, Vector{Float64}}, SparseArrays.CHOLMOD.Factor{Float64, Int64}, TimerOutput}, Objective{Cthonios.Domain{FiniteElementContainers.VectorizedNodalField{Float64, 2, 2, 8416, Vector{Float64}}, DofManager{Int64, 2, 8416, Vector{Float64}, Vector{Int64}}, @NamedTuple{Block1::Cthonios.SectionInternal{NonAllocatedFunctionSpace{2, Vector{Int64}, FiniteElementContainers.VectorizedElementField{Int64, 2, 3, 14768, Vector{Int64}}, FiniteElementContainers.VectorizedElementField{Int64, 2, 6, 14768, Vector{Int64}}, ReferenceFiniteElements.ReferenceFE{Int64, Float64, ReferenceFiniteElements.Tri3{SpecialPolynomials.Lagrange, 2}, ReferenceFiniteElements.Edge{2, SpecialPolynomials.Lagrange, 1, 2}, ReferenceFiniteElements.ArrayBackend{SArray}, Vector{SVector{2, Int64}}, Vector{SVector{0, Union{}}}, Vector{Int64}, Vector{SVector{2, Float64}}, ReferenceFiniteElements.CellInterpolants{StructArrays.StructVector{ReferenceFiniteElements.Interpolants, @NamedTuple{ξ::Vector{SVector{2, Float64}}, w::Vector{Float64}, N::Vector{SVector{3, Float64}}, ∇N_ξ::Vector{SMatrix{3, 2, Float64, 6}}, ∇∇N_ξ::Vector{SArray{Tuple{3, 2, 2}, Float64, 3, 12}}}, Int64}}, Vector{Vector{SVector{2, Float64}}}, ReferenceFiniteElements.SurfaceInterpolants{StructArrays.StructArray{ReferenceFiniteElements.Interpolants, 2, @NamedTuple{ξ::Matrix{SVector{2, Float64}}, w::Matrix{Float64}, N::Matrix{SVector{3, Float64}}, ∇N_ξ::Matrix{SMatrix{3, 2, Float64, 6}}, ∇∇N_ξ::Matrix{SArray{Tuple{3, 2, 2}, Float64, 3, 12}}}, Int64}}}}, SolidMechanics{2, 2, 0, NeoHookean, PlaneStrain, Cthonios.StressDivergence{2, 2, 0, NeoHookean, PlaneStrain}}, SVector{2, Float64}}}, @NamedTuple{yplus_nodeset_2::Cthonios.DirichletBCInternal{Vector{Int32}, Vector{Int64}, typeof(Main.func_1)}, yminus_nodeset_1_2::Cthonios.DirichletBCInternal{Vector{Int32}, Vector{Int64}, typeof(Main.func_2)}, yplus_nodeset_1::Cthonios.DirichletBCInternal{Vector{Int32}, Vector{Int64}, typeof(Main.func_2)}}, Vector{Int64}, @NamedTuple{}, @NamedTuple{}}, typeof(Cthonios.energy), typeof(Cthonios.gradient), typeof(Cthonios.hessian), typeof(Cthonios.neumann_energy), typeof(Cthonios.neumann_gradient), typeof(Cthonios.neumann_hessian), TimerOutput}, Vector{Float64}, Cthonios.WarmStart{FiniteElementContainers.VectorizedNodalField{Float64, 2, 2, 8416, Vector{Float64}}, Vector{Float64}, ObjectiveParameters{FiniteElementContainers.VectorizedNodalField{Float64, 2, 2, 8416, Vector{Float64}}, QuasiStatic{Vector{Int64}, Float64, Vector{Float64}}, Vector{Float64}, Vector{SVector{2, Float64}}, ComponentArrays.ComponentVector{Float64, Vector{Float64}, Tuple{ComponentArrays.Axis{(Block1 = ViewAxis(1:0, ShapedAxis((0, 3, 14768))),)}}}, ComponentArrays.ComponentVector{Float64, Vector{Float64}, Tuple{ComponentArrays.Axis{(Block1 = 1:2,)}}}, FiniteElementContainers.VectorizedNodalField{Float64, 2, 2, 8416, Vector{Float64}}, FiniteElementContainers.VectorizedNodalField{Float64, 2, 2, 8416, Vector{Float64}}, Vector{Float64}, FiniteElementContainers.VectorizedNodalField{Float64, 2, 2, 8416, Vector{Float64}}, ComponentArrays.ComponentVector{Float64, Vector{Float64}, Tuple{ComponentArrays.Axis{(Block1 = ViewAxis(1:44304, ShapedAxis((3, 14768))),)}}}}}, TimerOutput, FiniteElementContainers.VectorizedNodalField{Float64, 2, 2, 8416, Vector{Float64}}, Cthonios.TrustRegionSolverSettings}(CholeskyPreconditioner{FiniteElementContainers.StaticAssembler{Float64, Int64, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}, FiniteElementContainers.VectorizedNodalField{Float64, 2, 2, 8416, Vector{Float64}}, Vector{Float64}, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Float64}, Vector{Int64}, Vector{Int64}, Vector{Float64}}, SparseArrays.CHOLMOD.Factor{Float64, Int64}, TimerOutput}(FiniteElementContainers.StaticAssembler{Float64, Int64, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}, FiniteElementContainers.VectorizedNodalField{Float64, 2, 2, 8416, Vector{Float64}}, Vector{Float64}, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Float64}, Vector{Int64}, Vector{Int64}, Vector{Float64}}([1, 2, 3, 4, 1, 2, 3, 4, 1, 2 … 15213, 15214, 14713, 14714, 15715, 15716, 15213, 15214, 14713, 14714], [1, 1, 1, 1, 2, 2, 2, 2, 3, 3 … 14713, 14713, 14713, 14713, 14714, 14714, 14714, 14714, 14714, 14714], [1, 2, 5, 6, 7, 8, 11, 12, 25, 26 … 531639, 531640, 531641, 531642, 531643, 531644, 531645, 531646, 531647, 531648], [14768], [36], [0.0 0.0 … 0.0 0.0; 0.0 0.0 … 0.0 0.0], [1.0936717307567239, -0.292890957074626, -0.30874880569639573, 0.08916014301262334, -0.7849229250603279, 0.20373081406200264, -0.292890957074626, 0.7500337522693039, 0.5059934763459568, -0.6357183841377304 … -0.7268771191494259, 0.08076579504917981, 0.9979500765131166, -0.18077748700999746, 0.516845025294151, -0.4374812602032614, -0.3360675382841535, -0.11729029295148521, -0.18077748700999746, 0.5547715531547467], [144643, 144644, 118497, 118498, 144635, 144636, 144619, 144620, 94955, 94956 … 212007, 212008, 212035, 212036, 212063, 212064, 212091, 212092, 212117, 212118], [1, 11, 21, 31, 41, 51, 61, 71, 81, 91 … 212007, 212021, 212035, 212049, 212063, 212077, 212091, 212105, 212117, 212129], [1, 2, 3, 4, 5, 6, 9623, 9624, 11471, 11472 … 16251, 16252, 16315, 16316, 15729, 15730, 16279, 16280, 16315, 16316], [5.4007533141565585, -0.12274042142879321, -1.5186858384289632, -0.06110082032272626, -1.2926637279292954, 0.012201741061584614, -0.46358225309950574, -0.49453493243157615, -0.9498052563957434, 0.6572070814006334 … 0.06184611651062613, -0.8785394483537146, 0.03552764141376441, 1.5006309589952131, 0.36180915966271, -0.22120364003548657, -0.5826270613111169, -0.5956993044027834, 0.22081790164840692, 0.81690294443827], [1, 11, 21, 31, 41, 51, 61, 71, 81, 91 … 212007, 212021, 212035, 212049, 212063, 212077, 212091, 212105, 212117, 212129], [1, 2, 3, 4, 5, 6, 9623, 9624, 11471, 11472 … 15729, 15730, 16089, 16090, 16251, 16252, 16279, 16280, 16315, 16316], [5.4007533141565585, -0.12274042142879321, -1.518685838428963, -0.06110082032272626, -1.2926637279292952, 0.012201741061584559, -0.46358225309950574, -0.49453493243157604, -0.9498052563957434, 0.6572070814006336 … 0.49343249162560143, -0.6682757065971714, -0.20353841920379576, -2.3339441745135883, 0.6479072940312096, -1.7785439159854481, -0.6800008192355074, -1.2177908150442818, -0.06791756643505412, 6.213529049365272]), SparseArrays.CHOLMOD.Factor{Float64, Int64}
type: LLt
method: supernodal
maxnnz: 0
nnz: 566054
success: true
, ────────────────────────────────────────────────────────────────────────────────
Time Allocations
─────────────────────── ────────────────────────
Tot / % measured: 4.91s / 11.2% 625MiB / 19.4%
Section ncalls time %tot avg alloc %tot avg
────────────────────────────────────────────────────────────────────────────────
TrustRegionSolver ... 1 550ms 100.0% 550ms 121MiB 100.0% 121MiB
CholeskyPrecondi... 1 171ms 31.2% 171ms 56.6MiB 46.7% 56.6MiB
Objective - he... 1 34.8ms 6.3% 34.8ms 8.54MiB 7.0% 8.54MiB
────────────────────────────────────────────────────────────────────────────────), Objective{Cthonios.Domain{FiniteElementContainers.VectorizedNodalField{Float64, 2, 2, 8416, Vector{Float64}}, DofManager{Int64, 2, 8416, Vector{Float64}, Vector{Int64}}, @NamedTuple{Block1::Cthonios.SectionInternal{NonAllocatedFunctionSpace{2, Vector{Int64}, FiniteElementContainers.VectorizedElementField{Int64, 2, 3, 14768, Vector{Int64}}, FiniteElementContainers.VectorizedElementField{Int64, 2, 6, 14768, Vector{Int64}}, ReferenceFiniteElements.ReferenceFE{Int64, Float64, ReferenceFiniteElements.Tri3{SpecialPolynomials.Lagrange, 2}, ReferenceFiniteElements.Edge{2, SpecialPolynomials.Lagrange, 1, 2}, ReferenceFiniteElements.ArrayBackend{SArray}, Vector{SVector{2, Int64}}, Vector{SVector{0, Union{}}}, Vector{Int64}, Vector{SVector{2, Float64}}, ReferenceFiniteElements.CellInterpolants{StructArrays.StructVector{ReferenceFiniteElements.Interpolants, @NamedTuple{ξ::Vector{SVector{2, Float64}}, w::Vector{Float64}, N::Vector{SVector{3, Float64}}, ∇N_ξ::Vector{SMatrix{3, 2, Float64, 6}}, ∇∇N_ξ::Vector{SArray{Tuple{3, 2, 2}, Float64, 3, 12}}}, Int64}}, Vector{Vector{SVector{2, Float64}}}, ReferenceFiniteElements.SurfaceInterpolants{StructArrays.StructArray{ReferenceFiniteElements.Interpolants, 2, @NamedTuple{ξ::Matrix{SVector{2, Float64}}, w::Matrix{Float64}, N::Matrix{SVector{3, Float64}}, ∇N_ξ::Matrix{SMatrix{3, 2, Float64, 6}}, ∇∇N_ξ::Matrix{SArray{Tuple{3, 2, 2}, Float64, 3, 12}}}, Int64}}}}, SolidMechanics{2, 2, 0, NeoHookean, PlaneStrain, Cthonios.StressDivergence{2, 2, 0, NeoHookean, PlaneStrain}}, SVector{2, Float64}}}, @NamedTuple{yplus_nodeset_2::Cthonios.DirichletBCInternal{Vector{Int32}, Vector{Int64}, typeof(Main.func_1)}, yminus_nodeset_1_2::Cthonios.DirichletBCInternal{Vector{Int32}, Vector{Int64}, typeof(Main.func_2)}, yplus_nodeset_1::Cthonios.DirichletBCInternal{Vector{Int32}, Vector{Int64}, typeof(Main.func_2)}}, Vector{Int64}, @NamedTuple{}, @NamedTuple{}}, typeof(Cthonios.energy), typeof(Cthonios.gradient), typeof(Cthonios.hessian), typeof(Cthonios.neumann_energy), typeof(Cthonios.neumann_gradient), typeof(Cthonios.neumann_hessian), TimerOutput}(Cthonios.Domain{FiniteElementContainers.VectorizedNodalField{Float64, 2, 2, 8416, Vector{Float64}}, DofManager{Int64, 2, 8416, Vector{Float64}, Vector{Int64}}, @NamedTuple{Block1::Cthonios.SectionInternal{NonAllocatedFunctionSpace{2, Vector{Int64}, FiniteElementContainers.VectorizedElementField{Int64, 2, 3, 14768, Vector{Int64}}, FiniteElementContainers.VectorizedElementField{Int64, 2, 6, 14768, Vector{Int64}}, ReferenceFiniteElements.ReferenceFE{Int64, Float64, ReferenceFiniteElements.Tri3{SpecialPolynomials.Lagrange, 2}, ReferenceFiniteElements.Edge{2, SpecialPolynomials.Lagrange, 1, 2}, ReferenceFiniteElements.ArrayBackend{SArray}, Vector{SVector{2, Int64}}, Vector{SVector{0, Union{}}}, Vector{Int64}, Vector{SVector{2, Float64}}, ReferenceFiniteElements.CellInterpolants{StructArrays.StructVector{ReferenceFiniteElements.Interpolants, @NamedTuple{ξ::Vector{SVector{2, Float64}}, w::Vector{Float64}, N::Vector{SVector{3, Float64}}, ∇N_ξ::Vector{SMatrix{3, 2, Float64, 6}}, ∇∇N_ξ::Vector{SArray{Tuple{3, 2, 2}, Float64, 3, 12}}}, Int64}}, Vector{Vector{SVector{2, Float64}}}, ReferenceFiniteElements.SurfaceInterpolants{StructArrays.StructArray{ReferenceFiniteElements.Interpolants, 2, @NamedTuple{ξ::Matrix{SVector{2, Float64}}, w::Matrix{Float64}, N::Matrix{SVector{3, Float64}}, ∇N_ξ::Matrix{SMatrix{3, 2, Float64, 6}}, ∇∇N_ξ::Matrix{SArray{Tuple{3, 2, 2}, Float64, 3, 12}}}, Int64}}}}, SolidMechanics{2, 2, 0, NeoHookean, PlaneStrain, Cthonios.StressDivergence{2, 2, 0, NeoHookean, PlaneStrain}}, SVector{2, Float64}}}, @NamedTuple{yplus_nodeset_2::Cthonios.DirichletBCInternal{Vector{Int32}, Vector{Int64}, typeof(Main.func_1)}, yminus_nodeset_1_2::Cthonios.DirichletBCInternal{Vector{Int32}, Vector{Int64}, typeof(Main.func_2)}, yplus_nodeset_1::Cthonios.DirichletBCInternal{Vector{Int32}, Vector{Int64}, typeof(Main.func_2)}}, Vector{Int64}, @NamedTuple{}, @NamedTuple{}}([44.29728315483824 44.09375 … 40.010176197600856 20.023250115402412; 50.65587225363096 51.0 … 30.205590373069363 30.2388316572343], DofManager
Number of nodes = 8416
Number of dofs per node = 2
Storage type = Vector{Int64}, (Block1 = Cthonios.SectionInternal{NonAllocatedFunctionSpace{2, Vector{Int64}, FiniteElementContainers.VectorizedElementField{Int64, 2, 3, 14768, Vector{Int64}}, FiniteElementContainers.VectorizedElementField{Int64, 2, 6, 14768, Vector{Int64}}, ReferenceFiniteElements.ReferenceFE{Int64, Float64, ReferenceFiniteElements.Tri3{SpecialPolynomials.Lagrange, 2}, ReferenceFiniteElements.Edge{2, SpecialPolynomials.Lagrange, 1, 2}, ReferenceFiniteElements.ArrayBackend{SArray}, Vector{SVector{2, Int64}}, Vector{SVector{0, Union{}}}, Vector{Int64}, Vector{SVector{2, Float64}}, ReferenceFiniteElements.CellInterpolants{StructArrays.StructVector{ReferenceFiniteElements.Interpolants, @NamedTuple{ξ::Vector{SVector{2, Float64}}, w::Vector{Float64}, N::Vector{SVector{3, Float64}}, ∇N_ξ::Vector{SMatrix{3, 2, Float64, 6}}, ∇∇N_ξ::Vector{SArray{Tuple{3, 2, 2}, Float64, 3, 12}}}, Int64}}, Vector{Vector{SVector{2, Float64}}}, ReferenceFiniteElements.SurfaceInterpolants{StructArrays.StructArray{ReferenceFiniteElements.Interpolants, 2, @NamedTuple{ξ::Matrix{SVector{2, Float64}}, w::Matrix{Float64}, N::Matrix{SVector{3, Float64}}, ∇N_ξ::Matrix{SMatrix{3, 2, Float64, 6}}, ∇∇N_ξ::Matrix{SArray{Tuple{3, 2, 2}, Float64, 3, 12}}}, Int64}}}}, SolidMechanics{2, 2, 0, NeoHookean, PlaneStrain, Cthonios.StressDivergence{2, 2, 0, NeoHookean, PlaneStrain}}, SVector{2, Float64}}(1, 2, NonAllocatedFunctionSpace
Reference finite element = ReferenceFE:
Element Type = ReferenceFiniteElements.Tri3{SpecialPolynomials.Lagrange, 2}
Surface Element Type = ReferenceFiniteElements.Edge{2, SpecialPolynomials.Lagrange, 1, 2}
Interpolation Type = SpecialPolynomials.Lagrange
Polynomial Degree = 1
Quadrature Degree = 2
Array Backend Type = ReferenceFiniteElements.ArrayBackend{SArray}()
Edge Nodes Storage Type = Vector{SVector{2, Int64}}
Face Nodes Storage Type = Vector{SVector{0, Union{}}}
Interior Nodes Storage Type = Vector{Int64}
Nodal Coordinates Storage Type = Vector{SVector{2, Float64}}
CellInterpolants:
Quadrature Point Storage Type = Vector{SVector{2, Float64}}
Quadrature Weights Storage Type = Vector{Float64}
Shape Function Values Storage Type = Vector{SVector{3, Float64}}
Shape Function Gradients Storage Type = Vector{SMatrix{3, 2, Float64, 6}}
Shape Function Hessians Storage Type = Vector{SArray{Tuple{3, 2, 2}, Float64, 3, 12}}
Surface Nodal Coordinates Storage Type = Vector{Vector{SVector{2, Float64}}}
CellInterpolants:
Quadrature Point Storage Type = Matrix{SVector{2, Float64}}
Quadrature Weights Storage Type = Matrix{Float64}
Shape Function Values Storage Type = Matrix{SVector{3, Float64}}
Shape Function Gradients Storage Type = Matrix{SMatrix{3, 2, Float64, 6}}
Shape Function Hessians Storage Type = Matrix{SArray{Tuple{3, 2, 2}, Float64, 3, 12}}
, SolidMechanics{2, 2, 0, NeoHookean, PlaneStrain, Cthonios.StressDivergence{2, 2, 0, NeoHookean, PlaneStrain}}(NeoHookean(), PlaneStrain(), Cthonios.StressDivergence{2, 2, 0, NeoHookean, PlaneStrain}(NeoHookean(), PlaneStrain())), [0.833, 1.0]),), (yplus_nodeset_2 = DirichletBCInternal:
Numbe of active dofs = 129
Function = func_1
, yminus_nodeset_1_2 = DirichletBCInternal:
Numbe of active dofs = 258
Function = func_2
, yplus_nodeset_1 = DirichletBCInternal:
Numbe of active dofs = 129
Function = func_2
), [4, 10, 14, 18, 22, 46, 48, 52, 60, 80 … 8247, 8267, 8269, 8277, 8281, 8287, 8295, 8305, 8399, 8411], NamedTuple(), NamedTuple()), Cthonios.energy, Cthonios.gradient, Cthonios.hessian, Cthonios.neumann_energy, Cthonios.neumann_gradient, Cthonios.neumann_hessian, ────────────────────────────────────────────────────────────────────────────────
Time Allocations
─────────────────────── ────────────────────────
Tot / % measured: 4.94s / 11.1% 629MiB / 19.3%
Section ncalls time %tot avg alloc %tot avg
────────────────────────────────────────────────────────────────────────────────
TrustRegionSolver ... 1 550ms 100.0% 550ms 121MiB 100.0% 121MiB
CholeskyPrecondi... 1 171ms 31.2% 171ms 56.6MiB 46.7% 56.6MiB
Objective - he... 1 34.8ms 6.3% 34.8ms 8.54MiB 7.0% 8.54MiB
────────────────────────────────────────────────────────────────────────────────), [0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 … 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0], Cthonios.WarmStart{FiniteElementContainers.VectorizedNodalField{Float64, 2, 2, 8416, Vector{Float64}}, Vector{Float64}, ObjectiveParameters{FiniteElementContainers.VectorizedNodalField{Float64, 2, 2, 8416, Vector{Float64}}, QuasiStatic{Vector{Int64}, Float64, Vector{Float64}}, Vector{Float64}, Vector{SVector{2, Float64}}, ComponentArrays.ComponentVector{Float64, Vector{Float64}, Tuple{ComponentArrays.Axis{(Block1 = ViewAxis(1:0, ShapedAxis((0, 3, 14768))),)}}}, ComponentArrays.ComponentVector{Float64, Vector{Float64}, Tuple{ComponentArrays.Axis{(Block1 = 1:2,)}}}, FiniteElementContainers.VectorizedNodalField{Float64, 2, 2, 8416, Vector{Float64}}, FiniteElementContainers.VectorizedNodalField{Float64, 2, 2, 8416, Vector{Float64}}, Vector{Float64}, FiniteElementContainers.VectorizedNodalField{Float64, 2, 2, 8416, Vector{Float64}}, ComponentArrays.ComponentVector{Float64, Vector{Float64}, Tuple{ComponentArrays.Axis{(Block1 = ViewAxis(1:44304, ShapedAxis((3, 14768))),)}}}}}([0.0 0.0 … 0.0 0.0; 0.0 0.0 … 0.0 0.0], [0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 … 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0], ObjectiveParameters{FiniteElementContainers.VectorizedNodalField{Float64, 2, 2, 8416, Vector{Float64}}, QuasiStatic{Vector{Int64}, Float64, Vector{Float64}}, Vector{Float64}, Vector{SVector{2, Float64}}, ComponentArrays.ComponentVector{Float64, Vector{Float64}, Tuple{ComponentArrays.Axis{(Block1 = ViewAxis(1:0, ShapedAxis((0, 3, 14768))),)}}}, ComponentArrays.ComponentVector{Float64, Vector{Float64}, Tuple{ComponentArrays.Axis{(Block1 = 1:2,)}}}, FiniteElementContainers.VectorizedNodalField{Float64, 2, 2, 8416, Vector{Float64}}, FiniteElementContainers.VectorizedNodalField{Float64, 2, 2, 8416, Vector{Float64}}, Vector{Float64}, FiniteElementContainers.VectorizedNodalField{Float64, 2, 2, 8416, Vector{Float64}}, ComponentArrays.ComponentVector{Float64, Vector{Float64}, Tuple{ComponentArrays.Axis{(Block1 = ViewAxis(1:44304, ShapedAxis((3, 14768))),)}}}}([0.0 0.0 … 0.0 0.0; 0.0 0.0 … 0.0 0.0], QuasiStatic:
Start time = 0.0
End time = 0.0
Time step = [0.0]
, [0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 … 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0], SVector{2, Float64}[], (Block1 = Array{Float64, 3}(undef, 0, 3, 14768)), (Block1 = Array{Float64, 3}(undef, 0, 3, 14768)), (Block1 = [0.0, 0.0]), [0.0 0.0 … 0.0 0.0; 0.0 0.0 … 0.0 0.0], [0.0 0.0 … 0.0 0.0; 0.0 0.0 … 0.0 0.0], [0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 … 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0], [0.0 0.0 … 0.0 0.0; 0.0 0.0 … 0.0 0.0], (Block1 = [0.0 0.0 … 0.0 0.0; 0.0 0.0 … 0.0 0.0; 0.0 0.0 … 0.0 0.0])), [0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 … 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0]), ────────────────────────────────────────────────────────────────────────────────
Time Allocations
─────────────────────── ────────────────────────
Tot / % measured: 5.32s / 10.3% 685MiB / 17.7%
Section ncalls time %tot avg alloc %tot avg
────────────────────────────────────────────────────────────────────────────────
TrustRegionSolver ... 1 550ms 100.0% 550ms 121MiB 100.0% 121MiB
CholeskyPrecondi... 1 171ms 31.2% 171ms 56.6MiB 46.7% 56.6MiB
Objective - he... 1 34.8ms 6.3% 34.8ms 8.54MiB 7.0% 8.54MiB
────────────────────────────────────────────────────────────────────────────────, [0.0], [0.0 0.0 … 0.0 0.0; 0.0 0.0 … 0.0 0.0], [0.0 0.0 … 0.0 0.0; 0.0 0.0 … 0.0 0.0], [0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 … 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0], [0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 … 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0], [0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 … 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0], [0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 … 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0], [0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 … 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0], [0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 … 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0], [0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 … 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0], true, TrustRegionSolverSettings
t_1 = 0.25
t_2 = 1.75
η_1 = 1.0e-10
η_2 = 0.1
η_3 = 0.5
max_trust_iters = 100
tol = 1.0e-8
max_cg_iters = 50
max_cumulative_cg_iters = 1000
cg_tol = 2.0e-10
cg_inexact_solve_ratio = 1.0e-5
tr_size = 2.0
min_tr_size = 1.0e-8
check_stability = false
use_preconditioned_inner_product_for_cg = true
use_incremental_objective = false
debug_info = true
over_iter = 0
), ExodusPostProcessor{ExodusDatabase{Int32, Int32, Int32, Float64, Initialization{2, 8416, 14768, 1, 2, 3}}}(ExodusDatabase:
File name = output.e
Mode = rw
Initialization:
Number of dim = 2
Number of nodes = 8416
Number of elem = 14768
Number of blocks = 1
Number of node sets = 2
Number of side sets = 3
Block:
Block1
NodeSet:
yminus_nodeset yplus_nodeset
SideSet:
remaining_surface_sideset yminus_sideset yplus_sideset
NodalVariable:
displ_x displ_y
, ["displ_x", "displ_y"]), ────────────────────────────────────────────────────────────────────────────────
Time Allocations
─────────────────────── ────────────────────────
Tot / % measured: 5.54s / 9.9% 698MiB / 17.4%
Section ncalls time %tot avg alloc %tot avg
────────────────────────────────────────────────────────────────────────────────
TrustRegionSolver ... 1 550ms 100.0% 550ms 121MiB 100.0% 121MiB
CholeskyPrecondi... 1 171ms 31.2% 171ms 56.6MiB 46.7% 56.6MiB
Objective - he... 1 34.8ms 6.3% 34.8ms 8.54MiB 7.0% 8.54MiB
────────────────────────────────────────────────────────────────────────────────)Finally, solve the problem
Uu, p = Cthonios.create_unknowns_and_parameters(problem)
Cthonios.solve!(problem, Uu, p)[ Info: ================================================================================================
[ Info: = Load step 1
[ Info: = Old Time 0.0
[ Info: = New Time 0.05
[ Info: = End Time 1.0
[ Info: = % Completete 5.0
[ Info: ================================================================================================
[ Info: Warm start
[ Info: Initial objective = 1.715676e-02
[ Info: Initial residual = 1.522898e-02
[ Info: Updating preconditioner, attempt = 1
[ Info: ==============================================================================================================================
[ Info: Objective Model Objective Residual Model Residual CG TR size
[ Info: ==============================================================================================================================
[ Info: -5.063125e-05 -5.071297e-05 3.952473e-05 1.526497e-02 2 2.000000e+00 interior true
[ Info: -3.493087e-10 -3.493053e-10 5.006832e-10 3.952430e-05 3 2.000000e+00 interior true
[ Info: Converged
[ Info: ================================================================================================
[ Info: = Load step 2
[ Info: = Old Time 0.05
[ Info: = New Time 0.1
[ Info: = End Time 1.0
[ Info: = % Completete 10.0
[ Info: ================================================================================================
[ Info: Warm start
[ Info: Initial objective = 6.864668e-02
[ Info: Initial residual = 1.540703e-02
[ Info: Updating preconditioner, attempt = 1
[ Info: ==============================================================================================================================
[ Info: Objective Model Objective Residual Model Residual CG TR size
[ Info: ==============================================================================================================================
[ Info: -5.162061e-05 -5.170523e-05 4.065076e-05 1.544400e-02 2 2.000000e+00 interior true
[ Info: -3.681313e-10 -3.681308e-10 5.349619e-10 4.065030e-05 3 2.000000e+00 interior true
[ Info: Converged
[ Info: ================================================================================================
[ Info: = Load step 3
[ Info: = Old Time 0.1
[ Info: = New Time 0.15000000000000002
[ Info: = End Time 1.0
[ Info: = % Completete 15.000000000000002
[ Info: ================================================================================================
[ Info: Warm start
[ Info: Initial objective = 1.547735e-01
[ Info: Initial residual = 1.558793e-02
[ Info: Updating preconditioner, attempt = 1
[ Info: ==============================================================================================================================
[ Info: Objective Model Objective Residual Model Residual CG TR size
[ Info: ==============================================================================================================================
[ Info: -5.263328e-05 -5.272091e-05 4.181887e-05 1.562593e-02 2 2.000000e+00 interior true
[ Info: -3.881529e-10 -3.881516e-10 5.719933e-10 4.181837e-05 3 2.000000e+00 interior true
[ Info: Converged
[ Info: ================================================================================================
[ Info: = Load step 4
[ Info: = Old Time 0.15000000000000002
[ Info: = New Time 0.2
[ Info: = End Time 1.0
[ Info: = % Completete 20.0
[ Info: ================================================================================================
[ Info: Warm start
[ Info: Initial objective = 2.757886e-01
[ Info: Initial residual = 1.577201e-02
[ Info: Updating preconditioner, attempt = 1
[ Info: ==============================================================================================================================
[ Info: Objective Model Objective Residual Model Residual CG TR size
[ Info: ==============================================================================================================================
[ Info: -5.367406e-05 -5.376484e-05 4.303537e-05 1.581106e-02 2 2.000000e+00 interior true
[ Info: -4.095595e-10 -4.095596e-10 6.122325e-10 4.303484e-05 3 2.000000e+00 interior true
[ Info: Converged
[ Info: ================================================================================================
[ Info: = Load step 5
[ Info: = Old Time 0.2
[ Info: = New Time 0.25
[ Info: = End Time 1.0
[ Info: = % Completete 25.0
[ Info: ================================================================================================
[ Info: Warm start
[ Info: Initial objective = 4.319400e-01
[ Info: Initial residual = 1.595983e-02
[ Info: Updating preconditioner, attempt = 1
[ Info: ==============================================================================================================================
[ Info: Objective Model Objective Residual Model Residual CG TR size
[ Info: ==============================================================================================================================
[ Info: -5.475386e-05 -5.484797e-05 4.431244e-05 1.599998e-02 2 2.000000e+00 interior true
[ Info: -4.327220e-10 -4.327128e-10 6.564758e-10 4.431187e-05 3 2.000000e+00 interior true
[ Info: Converged
[ Info: ================================================================================================
[ Info: = Load step 6
[ Info: = Old Time 0.25
[ Info: = New Time 0.3
[ Info: = End Time 1.0
[ Info: = % Completete 30.0
[ Info: ================================================================================================
[ Info: Warm start
[ Info: Initial objective = 6.234715e-01
[ Info: Initial residual = 1.615261e-02
[ Info: Updating preconditioner, attempt = 1
[ Info: ==============================================================================================================================
[ Info: Objective Model Objective Residual Model Residual CG TR size
[ Info: ==============================================================================================================================
[ Info: -5.590182e-05 -5.599952e-05 4.567699e-05 1.619393e-02 2 2.000000e+00 interior true
[ Info: -4.584322e-10 -4.584345e-10 7.064555e-10 4.567638e-05 3 2.000000e+00 interior true
[ Info: Converged
[ Info: ================================================================================================
[ Info: = Load step 7
[ Info: = Old Time 0.3
[ Info: = New Time 0.35
[ Info: = End Time 1.0
[ Info: = % Completete 35.0
[ Info: ================================================================================================
[ Info: Warm start
[ Info: Initial objective = 8.506217e-01
[ Info: Initial residual = 1.635319e-02
[ Info: Updating preconditioner, attempt = 1
[ Info: ==============================================================================================================================
[ Info: Objective Model Objective Residual Model Residual CG TR size
[ Info: ==============================================================================================================================
[ Info: -5.721805e-05 -5.731986e-05 4.720246e-05 1.639578e-02 2 2.000000e+00 interior true
[ Info: -4.894383e-10 -4.894354e-10 7.686091e-10 4.720180e-05 3 2.000000e+00 interior true
[ Info: Converged
[ Info: ================================================================================================
[ Info: = Load step 8
[ Info: = Old Time 0.35
[ Info: = New Time 0.39999999999999997
[ Info: = End Time 1.0
[ Info: = % Completete 40.0
[ Info: ================================================================================================
[ Info: Warm start
[ Info: Initial objective = 1.113620e+00
[ Info: Initial residual = 1.656977e-02
[ Info: Updating preconditioner, attempt = 1
[ Info: ==============================================================================================================================
[ Info: Objective Model Objective Residual Model Residual CG TR size
[ Info: ==============================================================================================================================
[ Info: -5.921809e-05 -5.932575e-05 4.925456e-05 1.661388e-02 2 2.000000e+00 interior true
[ Info: -5.511760e-10 -5.511644e-10 8.219089e-10 4.925383e-05 4 2.000000e+00 interior true
[ Info: Converged
[ Info: ================================================================================================
[ Info: = Load step 9
[ Info: = Old Time 0.39999999999999997
[ Info: = New Time 0.44999999999999996
[ Info: = End Time 1.0
[ Info: = % Completete 44.99999999999999
[ Info: ================================================================================================
[ Info: Warm start
[ Info: Initial objective = 1.412672e+00
[ Info: Initial residual = 1.683543e-02
[ Info: Updating preconditioner, attempt = 1
[ Info: ==============================================================================================================================
[ Info: Objective Model Objective Residual Model Residual CG TR size
[ Info: ==============================================================================================================================
[ Info: -6.746924e-05 -6.760428e-05 7.631099e-05 1.688134e-02 2 2.000000e+00 interior true
[ Info: -4.448139e-09 -4.448067e-09 1.330728e-09 7.631001e-05 4 2.000000e+00 interior true
[ Info: Converged
[ Info: ================================================================================================
[ Info: = Load step 10
[ Info: = Old Time 0.44999999999999996
[ Info: = New Time 0.49999999999999994
[ Info: = End Time 1.0
[ Info: = % Completete 49.99999999999999
[ Info: ================================================================================================
[ Info: Warm start
[ Info: Initial objective = 1.747763e+00
[ Info: Initial residual = 1.744863e-02
[ Info: Updating preconditioner, attempt = 1
[ Info: ==============================================================================================================================
[ Info: Objective Model Objective Residual Model Residual CG TR size
[ Info: ==============================================================================================================================
[ Info: -2.777264e-04 -2.830970e-04 7.609987e-03 2.291558e-02 2 2.000000e+00 interior true
[ Info: -1.720672e-04 -1.765163e-04 3.618123e-03 9.746703e-03 6 2.000000e+00 interior true
[ Info: -6.082553e-05 -6.011008e-05 1.111181e-03 3.129326e-03 7 2.000000e+00 interior true
[ Info: -1.508474e-06 -1.508459e-06 4.952984e-08 1.111190e-03 7 2.000000e+00 interior true
[ Info: -8.824053e-13 -8.745459e-13 1.666394e-10 4.933149e-08 6 2.000000e+00 interior true
[ Info: Converged
[ Info: ================================================================================================
[ Info: = Load step 11
[ Info: = Old Time 0.49999999999999994
[ Info: = New Time 0.5499999999999999
[ Info: = End Time 1.0
[ Info: = % Completete 54.99999999999999
[ Info: ================================================================================================
[ Info: Warm start
[ Info: Initial objective = 2.136345e+00
[ Info: Initial residual = 1.703130e-01
[ Info: Updating preconditioner, attempt = 1
[ Info: ==============================================================================================================================
[ Info: Objective Model Objective Residual Model Residual CG TR size
[ Info: ==============================================================================================================================
[ Info: -2.596735e-02 -2.595844e-02 8.356207e-02 2.386198e-01 2 2.000000e+00 interior true
[ Info: -8.868150e-03 -8.793829e-03 6.472778e-03 8.461356e-02 5 2.000000e+00 interior true
[ Info: -1.891442e-04 -1.801108e-04 3.334573e-04 6.238576e-03 6 2.000000e+00 interior true
[ Info: -1.495072e-06 -1.485033e-06 5.543318e-06 3.279539e-04 7 2.000000e+00 interior true
[ Info: -1.621570e-10 -1.621636e-10 5.379263e-10 5.542785e-06 6 2.000000e+00 interior true
[ Info: Converged
[ Info: ================================================================================================
[ Info: = Load step 12
[ Info: = Old Time 0.5499999999999999
[ Info: = New Time 0.6
[ Info: = End Time 1.0
[ Info: = % Completete 60.0
[ Info: ================================================================================================
[ Info: Warm start
[ Info: Initial objective = 2.464661e+00
[ Info: Initial residual = 7.433010e-02
[ Info: Updating preconditioner, attempt = 1
[ Info: ==============================================================================================================================
[ Info: Objective Model Objective Residual Model Residual CG TR size
[ Info: ==============================================================================================================================
[ Info: -7.088957e-03 -6.923074e-03 1.023488e-03 7.353444e-02 2 2.000000e+00 interior true
[ Info: -7.306785e-05 -7.178960e-05 1.463838e-04 9.105397e-04 6 2.000000e+00 interior true
[ Info: -7.491317e-08 -7.485676e-08 1.047690e-07 1.462800e-04 6 2.000000e+00 interior true
[ Info: -1.358913e-13 -1.040181e-13 3.106095e-11 1.047605e-07 4 2.000000e+00 interior true
[ Info: Converged
[ Info: ================================================================================================
[ Info: = Load step 13
[ Info: = Old Time 0.6
[ Info: = New Time 0.65
[ Info: = End Time 1.0
[ Info: = % Completete 65.0
[ Info: ================================================================================================
[ Info: Warm start
[ Info: Initial objective = 2.815981e+00
[ Info: Initial residual = 4.065052e-02
[ Info: Updating preconditioner, attempt = 1
[ Info: ==============================================================================================================================
[ Info: Objective Model Objective Residual Model Residual CG TR size
[ Info: ==============================================================================================================================
[ Info: -2.017834e-03 -1.990687e-03 3.076176e-04 4.041749e-02 2 2.000000e+00 interior true
[ Info: -4.571913e-06 -4.556908e-06 6.635017e-06 3.022532e-04 5 2.000000e+00 interior true
[ Info: -1.490608e-10 -1.490617e-10 1.414513e-10 6.634878e-06 5 2.000000e+00 interior true
[ Info: Converged
[ Info: ================================================================================================
[ Info: = Load step 14
[ Info: = Old Time 0.65
[ Info: = New Time 0.7000000000000001
[ Info: = End Time 1.0
[ Info: = % Completete 70.0
[ Info: ================================================================================================
[ Info: Warm start
[ Info: Initial objective = 3.171121e+00
[ Info: Initial residual = 2.839646e-02
[ Info: Updating preconditioner, attempt = 1
[ Info: ==============================================================================================================================
[ Info: Objective Model Objective Residual Model Residual CG TR size
[ Info: ==============================================================================================================================
[ Info: -9.361653e-04 -9.286702e-04 1.271555e-04 2.830527e-02 2 2.000000e+00 interior true
[ Info: -6.622107e-07 -6.615778e-07 7.560913e-07 1.265616e-04 5 2.000000e+00 interior true
[ Info: -1.823874e-12 -1.823735e-12 1.564895e-10 7.560467e-07 3 2.000000e+00 interior true
[ Info: Converged
[ Info: ================================================================================================
[ Info: = Load step 15
[ Info: = Old Time 0.7000000000000001
[ Info: = New Time 0.7500000000000001
[ Info: = End Time 1.0
[ Info: = % Completete 75.00000000000001
[ Info: ================================================================================================
[ Info: Warm start
[ Info: Initial objective = 3.526472e+00
[ Info: Initial residual = 2.203658e-02
[ Info: Updating preconditioner, attempt = 1
[ Info: ==============================================================================================================================
[ Info: Objective Model Objective Residual Model Residual CG TR size
[ Info: ==============================================================================================================================
[ Info: -5.390718e-04 -5.362673e-04 6.471699e-05 2.199292e-02 2 2.000000e+00 interior true
[ Info: -1.532046e-07 -1.531478e-07 1.476949e-07 6.460553e-05 4 2.000000e+00 interior true
[ Info: -8.038015e-14 -6.462940e-14 1.783371e-11 1.476852e-07 3 2.000000e+00 interior true
[ Info: Converged
[ Info: ================================================================================================
[ Info: = Load step 16
[ Info: = Old Time 0.7500000000000001
[ Info: = New Time 0.8000000000000002
[ Info: = End Time 1.0
[ Info: = % Completete 80.00000000000001
[ Info: ================================================================================================
[ Info: Warm start
[ Info: Initial objective = 3.881439e+00
[ Info: Initial residual = 1.812743e-02
[ Info: Updating preconditioner, attempt = 1
[ Info: ==============================================================================================================================
[ Info: Objective Model Objective Residual Model Residual CG TR size
[ Info: ==============================================================================================================================
[ Info: -3.501610e-04 -3.488915e-04 3.782159e-05 1.810363e-02 2 2.000000e+00 interior true
[ Info: -4.735849e-08 -4.735027e-08 4.067117e-08 3.779242e-05 4 2.000000e+00 interior true
[ Info: -1.332268e-15 -4.510369e-15 6.096810e-12 4.066447e-08 3 2.000000e+00 interior true
[ Info: Converged
[ Info: ================================================================================================
[ Info: = Load step 17
[ Info: = Old Time 0.8000000000000002
[ Info: = New Time 0.8500000000000002
[ Info: = End Time 1.0
[ Info: = % Completete 85.00000000000001
[ Info: ================================================================================================
[ Info: Warm start
[ Info: Initial objective = 4.235851e+00
[ Info: Initial residual = 1.548259e-02
[ Info: Updating preconditioner, attempt = 1
[ Info: ==============================================================================================================================
[ Info: Objective Model Objective Residual Model Residual CG TR size
[ Info: ==============================================================================================================================
[ Info: -2.456573e-04 -2.450039e-04 2.438562e-05 1.546849e-02 2 2.000000e+00 interior true
[ Info: -1.780080e-08 -1.779918e-08 1.431892e-08 2.437602e-05 4 2.000000e+00 interior true
[ Info: -1.154632e-14 -5.109267e-16 4.180022e-12 1.431224e-08 3 2.000000e+00 interior true
[ Info: Converged
[ Info: ================================================================================================
[ Info: = Load step 18
[ Info: = Old Time 0.8500000000000002
[ Info: = New Time 0.9000000000000002
[ Info: = End Time 1.0
[ Info: = % Completete 90.00000000000003
[ Info: ================================================================================================
[ Info: Warm start
[ Info: Initial objective = 4.589662e+00
[ Info: Initial residual = 1.358120e-02
[ Info: Updating preconditioner, attempt = 1
[ Info: ==============================================================================================================================
[ Info: Objective Model Objective Residual Model Residual CG TR size
[ Info: ==============================================================================================================================
[ Info: -1.818643e-04 -1.814957e-04 1.693851e-05 1.357236e-02 2 2.000000e+00 interior true
[ Info: -7.709119e-09 -7.708664e-09 6.180563e-09 1.693473e-05 4 2.000000e+00 interior true
[ Info: Converged
[ Info: ================================================================================================
[ Info: = Load step 19
[ Info: = Old Time 0.9000000000000002
[ Info: = New Time 0.9500000000000003
[ Info: = End Time 1.0
[ Info: = % Completete 95.00000000000003
[ Info: ================================================================================================
[ Info: Warm start
[ Info: Initial objective = 4.942876e+00
[ Info: Initial residual = 1.215645e-02
[ Info: Updating preconditioner, attempt = 1
[ Info: ==============================================================================================================================
[ Info: Objective Model Objective Residual Model Residual CG TR size
[ Info: ==============================================================================================================================
[ Info: -1.401277e-04 -1.399048e-04 1.248875e-05 1.215070e-02 2 2.000000e+00 interior true
[ Info: -3.726756e-09 -3.726616e-09 3.246852e-09 1.248698e-05 4 2.000000e+00 interior true
[ Info: Converged
[ Info: ================================================================================================
[ Info: = Load step 20
[ Info: = Old Time 0.9500000000000003
[ Info: = New Time 1.0000000000000002
[ Info: = End Time 1.0
[ Info: = % Completete 100.00000000000003
[ Info: ================================================================================================
[ Info: Warm start
[ Info: Initial objective = 5.295526e+00
[ Info: Initial residual = 1.105628e-02
[ Info: Updating preconditioner, attempt = 1
[ Info: ==============================================================================================================================
[ Info: Objective Model Objective Residual Model Residual CG TR size
[ Info: ==============================================================================================================================
[ Info: -1.113782e-04 -1.112359e-04 9.681376e-06 1.105242e-02 2 2.000000e+00 interior true
[ Info: -1.977335e-09 -1.977250e-09 2.105721e-09 9.680381e-06 4 2.000000e+00 interior true
[ Info: Converged
────────────────────────────────────────────────────────────────────────────────
Time Allocations
─────────────────────── ────────────────────────
Tot / % measured: 56.8s / 35.0% 4.38GiB / 34.6%
Section ncalls time %tot avg alloc %tot avg
────────────────────────────────────────────────────────────────────────────────
Problem - solve! 1 19.3s 97.2% 19.3s 1.40GiB 92.2% 1.40GiB
Integrator step 20 19.3s 97.2% 966ms 1.39GiB 92.0% 71.4MiB
TrustRegionSol... 20 19.3s 97.2% 966ms 1.38GiB 91.0% 70.6MiB
TrustRegionS... 53 9.78s 49.2% 185ms 238MiB 15.3% 4.49MiB
Objective ... 266 9.50s 47.8% 35.7ms 34.2MiB 2.2% 132KiB
Objectiv... 266 9.50s 47.8% 35.7ms 0.00B 0.0% 0.00B
TrustRegionS... 53 2.06s 10.4% 38.9ms 19.0MiB 1.2% 366KiB
Objective ... 53 1.91s 9.6% 35.9ms 6.81MiB 0.4% 132KiB
Objectiv... 53 1.90s 9.6% 35.9ms 0.00B 0.0% 0.00B
Objective - hvp 53 1.87s 9.4% 35.3ms 6.81MiB 0.4% 132KiB
Objective ... 53 1.87s 9.4% 35.2ms 0.00B 0.0% 0.00B
TrustRegionS... 20 1.74s 8.8% 87.1ms 539MiB 34.8% 27.0MiB
Warm start... 20 1.74s 8.8% 87.1ms 539MiB 34.8% 27.0MiB
Objectiv... 20 827ms 4.2% 41.3ms 104MiB 6.7% 5.18MiB
WarmStar... 20 500ms 2.5% 25.0ms 431MiB 27.8% 21.5MiB
WarmStar... 20 331ms 1.7% 16.6ms 0.00B 0.0% 0.00B
CholeskyPrec... 20 1.21s 6.1% 60.7ms 527MiB 34.0% 26.3MiB
Objective ... 20 808ms 4.1% 40.4ms 104MiB 6.7% 5.18MiB
Objective - ... 106 930ms 4.7% 8.78ms 7.36KiB 0.0% 71.1B
Objective ... 106 930ms 4.7% 8.78ms 0.00B 0.0% 0.00B
Objective - ... 73 700ms 3.5% 9.58ms 752B 0.0% 10.3B
Objectives... 73 698ms 3.5% 9.56ms 0.00B 0.0% 0.00B
ExodusPostProces... 21 8.36ms 0.0% 398μs 2.70MiB 0.2% 132KiB
TrustRegionSolver ... 1 550ms 2.8% 550ms 121MiB 7.8% 121MiB
CholeskyPrecondi... 1 171ms 0.9% 171ms 56.6MiB 3.7% 56.6MiB
Objective - he... 1 34.8ms 0.2% 34.8ms 8.54MiB 0.6% 8.54MiB
────────────────────────────────────────────────────────────────────────────────Interactive plotting
exo = ExodusDatabase("output.e", "r")
mesh = SimpleMesh(exo, 20)
viz(mesh)
This page was generated using Literate.jl.