Skip to content

Commit

Permalink
Minor change for testing
Browse files Browse the repository at this point in the history
  • Loading branch information
zjwegert committed Nov 15, 2024
1 parent dd9f060 commit 467ec50
Showing 1 changed file with 3 additions and 3 deletions.
6 changes: 3 additions & 3 deletions scripts/Embedded/Examples/2d_thermal.jl
Original file line number Diff line number Diff line change
Expand Up @@ -5,11 +5,11 @@ using GridapEmbedded, GridapEmbedded.LevelSetCutters
using GridapTopOpt: StateParamIntegrandWithMeasure

path="./results/UnfittedFEM_thermal_compliance_ALM/"
rm(path,force=true,recursive=true)
mkpath(path)
n = 101
n = 51
order = 1
γ = 0.1
γ_reinit = 0.5
max_steps = floor(Int,order*minimum(n)/10)
vf = 0.4
α_coeff = 4max_steps*γ
Expand Down Expand Up @@ -99,7 +99,7 @@ vel_ext = VelocityExtension(a_hilb,U_reg,V_reg)

## Optimiser
optimiser = AugmentedLagrangian(pcfs,ls_evo,vel_ext,φh;debug=true,
γ,γ_reinit,verbose=true,constraint_names=[:Vol])
γ,verbose=true,constraint_names=[:Vol])
for (it,uh,φh,state) in optimiser
if iszero(it % iter_mod)
writevtk(Ω,path*"Omega$it",cellfields=["φ"=>φh,"|∇(φ)|"=>(norm (φh)),"uh"=>uh,"velh"=>FEFunction(V_φ,state.vel)])
Expand Down

0 comments on commit 467ec50

Please sign in to comment.