Skip to content

Generalizes injectivity line search β parameter#128

Merged
jmartfrut merged 1 commit into
mainfrom
optinjectivity
Mar 5, 2026
Merged

Generalizes injectivity line search β parameter#128
jmartfrut merged 1 commit into
mainfrom
optinjectivity

Conversation

@jmartfrut

Copy link
Copy Markdown
Collaborator

Transforms the β parameter in the injectivity-preserving line search from a constant Float64 to a CellField. This change allows β to potentially vary across the computational domain, enhancing the flexibility of the injectivity preservation criterion. The constructor now ensures β is stored as a CellField, converting scalar inputs as needed, and update_cellstate! adapts its usage accordingly.

Transforms the `β` parameter in the injectivity-preserving line search from a constant `Float64` to a `CellField`. This change allows `β` to potentially vary across the computational domain, enhancing the flexibility of the injectivity preservation criterion. The constructor now ensures `β` is stored as a `CellField`, converting scalar inputs as needed, and `update_cellstate!` adapts its usage accordingly.
@jmartfrut jmartfrut merged commit c1ba357 into main Mar 5, 2026
3 of 4 checks passed
@github-actions

github-actions Bot commented Mar 5, 2026

Copy link
Copy Markdown

Benchmark Results (Julia v1)

Time benchmarks
main 6ce05ef... main / 6ce05ef...
Constitutive models/Visco-elastic Ψ 0.047 ± 0.0036 ms 0.0475 ± 0.0054 ms 0.989 ± 0.14
Constitutive models/Visco-elastic ∂Ψu 0.0614 ± 0.0017 ms 0.0628 ± 0.0023 ms 0.978 ± 0.044
Constitutive models/Visco-elastic ∂Ψuu 0.112 ± 0.016 ms 0.112 ± 0.015 ms 0.999 ± 0.2
Simulations/StaticMechanicalDirichlet 0.131 ± 0.0066 s 0.137 ± 0.0088 s 0.961 ± 0.079
Simulations/StaticMechanicalNeumann 0.129 ± 0.011 s 0.132 ± 0.013 s 0.983 ± 0.13
Simulations/ViscoElastic 20.4 s 21.2 s 0.963
Tensor algebra/Cofactor 0.07 ± 0.009 μs 0.07 ± 0.01 μs 1 ± 0.19
Tensor algebra/Det(A)Inv(A') 0.13 ± 0.01 μs 0.14 ± 0.01 μs 0.929 ± 0.097
Tensor algebra/δδ_λ_2d 30 ± 0 ns 30 ± 0 ns 1 ± 0
Tensor algebra/δδ_μ_2d 30 ± 0 ns 30 ± 0 ns 1 ± 0
time_to_load 0.0212 ± 2e-05 h 0.022 ± 0.00019 h 0.961 ± 0.0085
Memory benchmarks
main 6ce05ef... main / 6ce05ef...
Constitutive models/Visco-elastic Ψ 1.24 k allocs: 0.0777 MB 1.24 k allocs: 0.0777 MB 1
Constitutive models/Visco-elastic ∂Ψu 1.36 k allocs: 0.0898 MB 1.36 k allocs: 0.0898 MB 1
Constitutive models/Visco-elastic ∂Ψuu 2.32 k allocs: 0.155 MB 2.32 k allocs: 0.155 MB 1
Simulations/StaticMechanicalDirichlet 1.51 M allocs: 0.109 GB 1.51 M allocs: 0.109 GB 1
Simulations/StaticMechanicalNeumann 1.48 M allocs: 0.0933 GB 1.48 M allocs: 0.0933 GB 1
Simulations/ViscoElastic 0.343 G allocs: 22 GB 0.343 G allocs: 22 GB 1
Tensor algebra/Cofactor 1 allocs: 0.0781 kB 1 allocs: 0.0781 kB 1
Tensor algebra/Det(A)Inv(A') 4 allocs: 0.25 kB 4 allocs: 0.25 kB 1
Tensor algebra/δδ_λ_2d 0 allocs: 0 B 0 allocs: 0 B
Tensor algebra/δδ_μ_2d 0 allocs: 0 B 0 allocs: 0 B
time_to_load 0.141 k allocs: 10.8 kB 0.141 k allocs: 10.8 kB 1

@jmartfrut jmartfrut deleted the optinjectivity branch March 9, 2026 11:12
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

1 participant