using TopOpt
using WGLMakie
WGLMakie.activate!(; resize_to=:parent)
using Bonito
if haskey(ENV, "QUARTO_PROJECT_DIR")
Bonito.Page(exportable=true, offline=true)
else
Bonito.browser_display()
end
display_app(app) = display(app)Heat Conduction: Conductivity Tree Benchmark
Description
This tutorial demonstrates heat conduction topology optimization using the classic “conductivity tree” benchmark from Bendsøe & Sigmund’s Topology Optimization textbook (§1.3). The problem minimizes thermal compliance (temperature work) subject to a volume constraint.
The setup: distributed heat flux enters the top edge, the bottom edge is held at T = 0 (heat sink), and the sides are insulated. The optimal design forms a branching tree structure that efficiently conducts heat from the top to the bottom boundary — analogous to structural force flow paths.
Setup
WGLMakie.activate!(; resize_to=:parent) selects the browser renderer and fills the Quarto output column. Bonito.Page(exportable=true, offline=true) embeds the assets needed by visualize(...; static=true) in the Quarto output, so the visualization does not require a running Julia process.
Define the Problem
The HeatTree problem sets up the standard benchmark with heat flux on top and fixed temperature at the bottom:
nels = (80, 40) # 80×40 rectangular mesh
k = 1.0 # thermal conductivity
q = 1.0 # heat flux on top boundary
problem = HeatTree(nels, (1.0, 1.0), k; q=q)The boundary conditions:
- Top edge: Distributed heat flux q (Neumann BC)
- Bottom edge: Fixed temperature T = 0 (Dirichlet BC)
- Left/Right edges: Insulated (zero heat flux)
Parameter Settings
V = 0.3 # volume fraction (30% high-conductivity material)
xmin = 1e-3 # minimum density
rmin = 2.0 # density filter radius
x0 = fill(V, getncells(problem)) # start from uniform 30% density
solver = FEASolver(DirectSolver, problem; xmin=xmin, penalty=PowerPenaltyFun(1.0))
comp = ThermalComplianceFun(solver)
filter = DensityFilterFun(solver; rmin=rmin)The 30% volume constraint forces the optimizer to strategically place high- conductivity material for maximum heat dissipation.
Objective and Constraint
obj = x -> comp(filter(PseudoDensities(x))) # thermal compliance
volfrac = VolumeFun(solver)
constr = x -> volfrac(filter(PseudoDensities(x))) - V # volume constraintThermal compliance measures the work done by the heat flux against the temperature field: J = QᵀT. Minimizing it drives heat efficiently toward the cold boundary.
Optimization Setup
model = Model(obj)
addvar!(model, zeros(length(x0)), ones(length(x0)))
add_ineq_constraint!(model, constr)
alg = MMA87()
convcriteria = Nonconvex.KKTCriteria()Continuation Strategy
A fixed penalty often leaves intermediate (gray) densities. We use a continuation strategy that ramps the penalty exponent through 1.0, 2.0, and 3.0, with exponentially decaying KKT tolerances so each subproblem is solved accurately before the penalty increases:
ps = [1.0, 2.0, 3.0]
tols = [1e-2, 1e-3, 1e-4]
x = x0Run Optimization
Each continuation step:
- Set the penalty exponent
p - Set the convergence tolerance
tol - Solve the subproblem with MMA87
- Warm-start the next step from the current solution
for j in eachindex(ps)
global x, convcriteria, res
p = ps[j]
tol = tols[j]
setpenalty!(solver, p)
options = MMAOptions(; tol=Tolerance(; kkt=tol), maxiter=1000, convcriteria)
res = optimize(model, alg, x; options)
global x = res.minimizer
println("Step $j: penalty=$p, tol=$tol, obj=$(obj(x))")
end[ Info: iter obj Δobj violation kkt_residual [ Info: 0 1.1e+04 Inf 0.0e+00 1.1e+01 [ Info: 1 1.1e+04 4.5e-06 0.0e+00 4.5e-06 Step 1: penalty=1.0, tol=0.01, obj=10641.835721160478 [ Info: iter obj Δobj violation kkt_residual [ Info: 0 3.5e+04 Inf 0.0e+00 7.5e+01 [ Info: 1 3.5e+04 2.5e-05 0.0e+00 4.5e-06 Step 2: penalty=2.0, tol=0.001, obj=35199.6480080139 [ Info: iter obj Δobj violation kkt_residual [ Info: 0 1.1e+05 Inf 0.0e+00 3.6e+02 [ Info: 1 1.1e+05 5.5e-03 0.0e+00 5.5e-03 [ Info: 2 1.1e+05 3.2e-03 0.0e+00 2.4e-03 [ Info: 3 1.1e+05 1.3e-03 0.0e+00 3.6e-03 [ Info: 4 1.1e+05 9.4e-05 0.0e+00 3.5e-03 [ Info: 5 1.1e+05 1.5e-05 0.0e+00 3.5e-03 [ Info: 6 1.1e+05 8.2e-06 0.0e+00 3.5e-03 [ Info: 7 1.1e+05 3.5e-06 0.0e+00 3.5e-03 [ Info: 8 1.1e+05 9.2e-07 0.0e+00 3.5e-03 [ Info: 9 1.1e+05 7.3e-08 0.0e+00 3.5e-03 [ Info: 10 1.1e+05 1.0e-08 0.0e+00 3.5e-03 [ Info: 11 1.1e+05 1.9e-08 0.0e+00 3.5e-03 [ Info: 12 1.1e+05 8.3e-09 0.0e+00 3.5e-03 [ Info: 13 1.1e+05 2.2e-08 0.0e+00 3.5e-03 [ Info: 14 1.1e+05 4.1e-08 0.0e+00 3.5e-03 [ Info: 15 1.1e+05 1.2e-07 0.0e+00 3.5e-03 [ Info: 16 1.1e+05 2.7e-07 0.0e+00 3.5e-03 [ Info: 17 1.1e+05 8.6e-07 0.0e+00 3.5e-03 [ Info: 18 1.1e+05 2.6e-06 0.0e+00 3.5e-03 [ Info: 19 1.1e+05 1.0e-05 0.0e+00 3.5e-03 [ Info: 20 1.1e+05 4.1e-05 0.0e+00 3.5e-03 [ Info: 21 1.1e+05 2.1e-04 0.0e+00 3.3e-03 [ Info: 22 1.1e+05 1.1e-03 0.0e+00 2.2e-03 [ Info: 23 1.1e+05 7.6e-03 4.9e-09 5.4e-03 [ Info: 24 1.1e+05 5.5e-03 0.0e+00 1.1e-02 [ Info: 25 1.1e+05 9.5e-04 0.0e+00 9.7e-02 [ Info: 26 1.1e+05 2.2e-02 0.0e+00 1.0e+00 [ Info: 27 1.1e+05 3.0e+00 0.0e+00 1.2e+01 [ Info: 28 1.1e+05 5.2e+02 0.0e+00 1.6e+02 [ Info: 29 8.1e+04 3.3e+04 0.0e+00 2.6e+03 [ Info: 30 4.6e+04 3.6e+04 0.0e+00 1.2e+04 [ Info: 31 8.0e+04 3.5e+04 0.0e+00 1.2e+04 [ Info: 32 4.9e+04 3.2e+04 0.0e+00 7.8e+03 [ Info: 33 3.5e+04 1.4e+04 0.0e+00 1.0e+04 [ Info: 34 2.1e+04 1.4e+04 0.0e+00 1.8e+03 [ Info: 35 2.1e+04 2.6e+02 0.0e+00 2.5e+03 [ Info: 36 1.8e+04 2.2e+03 0.0e+00 9.8e+02 [ Info: 37 1.7e+04 1.0e+03 0.0e+00 1.9e+02 [ Info: 38 1.7e+04 3.6e+02 0.0e+00 5.8e+01 [ Info: 39 1.7e+04 1.6e+02 0.0e+00 2.4e+01 [ Info: 40 1.7e+04 7.9e+01 0.0e+00 1.9e+01 [ Info: 41 1.7e+04 6.1e+01 0.0e+00 1.7e+01 [ Info: 42 1.7e+04 6.4e+01 0.0e+00 1.9e+01 [ Info: 43 1.7e+04 4.4e+01 0.0e+00 2.0e+01 [ Info: 44 1.7e+04 3.1e+01 0.0e+00 2.1e+01 [ Info: 45 1.7e+04 2.6e+01 0.0e+00 2.0e+01 [ Info: 46 1.7e+04 2.3e+01 0.0e+00 1.4e+01 [ Info: 47 1.7e+04 1.7e+01 0.0e+00 9.5e+00 [ Info: 48 1.7e+04 1.5e+01 0.0e+00 9.9e+00 [ Info: 49 1.7e+04 1.8e+01 0.0e+00 1.1e+01 [ Info: 50 1.7e+04 1.3e+01 0.0e+00 1.0e+01 [ Info: 51 1.7e+04 5.9e+00 0.0e+00 9.6e+00 [ Info: 52 1.7e+04 3.8e+00 0.0e+00 5.1e+00 [ Info: 53 1.7e+04 3.2e+00 0.0e+00 4.8e+00 [ Info: 54 1.7e+04 2.7e+00 0.0e+00 4.9e+00 [ Info: 55 1.7e+04 2.0e+00 0.0e+00 4.9e+00 [ Info: 56 1.7e+04 2.0e+00 0.0e+00 5.1e+00 [ Info: 57 1.7e+04 2.8e+00 0.0e+00 5.1e+00 [ Info: 58 1.7e+04 3.0e+00 0.0e+00 5.0e+00 [ Info: 59 1.7e+04 2.4e+00 0.0e+00 5.2e+00 [ Info: 60 1.7e+04 2.2e+00 0.0e+00 5.1e+00 [ Info: 61 1.7e+04 2.8e+00 0.0e+00 5.1e+00 [ Info: 62 1.7e+04 2.7e+00 0.0e+00 5.0e+00 [ Info: 63 1.7e+04 1.7e+00 0.0e+00 4.7e+00 [ Info: 64 1.7e+04 1.5e+00 0.0e+00 3.6e+00 [ Info: 65 1.7e+04 1.2e+00 0.0e+00 3.4e+00 [ Info: 66 1.7e+04 1.1e+00 0.0e+00 2.2e+00 [ Info: 67 1.7e+04 8.3e-01 0.0e+00 1.8e+00 [ Info: 68 1.7e+04 7.7e-01 0.0e+00 1.7e+00 [ Info: 69 1.7e+04 8.1e-01 0.0e+00 1.5e+00 [ Info: 70 1.7e+04 8.8e-01 0.0e+00 1.5e+00 [ Info: 71 1.7e+04 8.1e-01 0.0e+00 2.2e+00 [ Info: 72 1.7e+04 5.0e-01 0.0e+00 3.6e+00 [ Info: 73 1.7e+04 5.6e-02 0.0e+00 5.8e+00 [ Info: 74 1.7e+04 2.2e-01 0.0e+00 3.7e+00 [ Info: 75 1.7e+04 1.6e-01 0.0e+00 4.8e+00 [ Info: 76 1.7e+04 2.9e-01 0.0e+00 5.5e+00 [ Info: 77 1.7e+04 7.8e-01 0.0e+00 3.3e+00 [ Info: 78 1.7e+04 1.2e+00 0.0e+00 1.1e+00 [ Info: 79 1.7e+04 4.9e-01 0.0e+00 8.5e-01 [ Info: 80 1.7e+04 3.0e-01 0.0e+00 8.9e-01 [ Info: 81 1.7e+04 2.7e-01 0.0e+00 8.9e-01 [ Info: 82 1.7e+04 2.3e-01 0.0e+00 8.2e-01 [ Info: 83 1.7e+04 1.7e-01 0.0e+00 7.7e-01 [ Info: 84 1.7e+04 1.8e-01 0.0e+00 7.2e-01 [ Info: 85 1.7e+04 2.0e-01 0.0e+00 7.1e-01 [ Info: 86 1.7e+04 2.2e-01 0.0e+00 7.5e-01 [ Info: 87 1.7e+04 2.5e-01 0.0e+00 7.9e-01 [ Info: 88 1.7e+04 2.9e-01 0.0e+00 8.2e-01 [ Info: 89 1.7e+04 3.2e-01 0.0e+00 8.1e-01 [ Info: 90 1.7e+04 3.5e-01 0.0e+00 8.8e-01 [ Info: 91 1.7e+04 3.7e-01 0.0e+00 8.2e-01 [ Info: 92 1.6e+04 2.9e-01 0.0e+00 1.3e+00 [ Info: 93 1.6e+04 2.6e-02 0.0e+00 4.3e+00 [ Info: 94 1.7e+04 9.5e-01 0.0e+00 6.0e+00 [ Info: 95 1.7e+04 8.2e-01 0.0e+00 8.3e+00 [ Info: 96 1.7e+04 1.1e+00 0.0e+00 3.4e+00 [ Info: 97 1.6e+04 1.7e+00 0.0e+00 2.4e+00 [ Info: 98 1.6e+04 7.4e-01 0.0e+00 1.9e+00 [ Info: 99 1.6e+04 4.5e-01 0.0e+00 1.0e+00 [ Info: 100 1.6e+04 3.7e-01 0.0e+00 9.0e-01 [ Info: 101 1.6e+04 2.8e-01 0.0e+00 8.2e-01 [ Info: 102 1.6e+04 2.5e-01 0.0e+00 7.9e-01 [ Info: 103 1.6e+04 2.7e-01 0.0e+00 7.4e-01 [ Info: 104 1.6e+04 3.1e-01 0.0e+00 8.4e-01 [ Info: 105 1.6e+04 3.5e-01 0.0e+00 9.8e-01 [ Info: 106 1.6e+04 3.8e-01 0.0e+00 9.8e-01 [ Info: 107 1.6e+04 3.9e-01 0.0e+00 9.9e-01 [ Info: 108 1.6e+04 4.5e-01 0.0e+00 1.0e+00 [ Info: 109 1.6e+04 5.3e-01 0.0e+00 1.1e+00 [ Info: 110 1.6e+04 5.6e-01 0.0e+00 1.1e+00 [ Info: 111 1.6e+04 4.5e-01 0.0e+00 9.7e-01 [ Info: 112 1.6e+04 3.9e-01 0.0e+00 9.4e-01 [ Info: 113 1.6e+04 3.5e-01 0.0e+00 8.7e-01 [ Info: 114 1.6e+04 3.0e-01 0.0e+00 1.0e+00 [ Info: 115 1.6e+04 2.0e-01 0.0e+00 2.1e+00 [ Info: 116 1.6e+04 4.5e-01 0.0e+00 5.1e+00 [ Info: 117 1.6e+04 8.7e-01 0.0e+00 4.8e+00 [ Info: 118 1.6e+04 2.4e+00 0.0e+00 1.2e+01 [ Info: 119 1.6e+04 9.9e-01 0.0e+00 7.8e+00 [ Info: 120 1.6e+04 5.3e-01 0.0e+00 1.4e+01 [ Info: 121 1.6e+04 2.0e+00 0.0e+00 6.5e+00 [ Info: 122 1.6e+04 3.3e+00 0.0e+00 2.6e+00 [ Info: 123 1.6e+04 9.1e-01 0.0e+00 1.1e+00 [ Info: 124 1.6e+04 1.3e-01 0.0e+00 1.1e+00 [ Info: 125 1.6e+04 9.3e-02 0.0e+00 1.1e+00 [ Info: 126 1.6e+04 1.0e-01 0.0e+00 1.1e+00 [ Info: 127 1.6e+04 9.8e-02 0.0e+00 1.0e+00 [ Info: 128 1.6e+04 1.1e-01 0.0e+00 1.0e+00 [ Info: 129 1.6e+04 1.2e-01 0.0e+00 9.9e-01 [ Info: 130 1.6e+04 1.3e-01 0.0e+00 7.8e-01 [ Info: 131 1.6e+04 9.0e-02 0.0e+00 5.7e-01 [ Info: 132 1.6e+04 7.0e-02 0.0e+00 5.3e-01 [ Info: 133 1.6e+04 6.6e-02 0.0e+00 5.0e-01 [ Info: 134 1.6e+04 6.7e-02 0.0e+00 4.8e-01 [ Info: 135 1.6e+04 6.9e-02 0.0e+00 4.6e-01 [ Info: 136 1.6e+04 6.9e-02 0.0e+00 4.4e-01 [ Info: 137 1.6e+04 5.7e-02 0.0e+00 4.1e-01 [ Info: 138 1.6e+04 5.1e-02 0.0e+00 4.1e-01 [ Info: 139 1.6e+04 5.1e-02 0.0e+00 4.2e-01 [ Info: 140 1.6e+04 5.1e-02 0.0e+00 4.4e-01 [ Info: 141 1.6e+04 5.1e-02 0.0e+00 4.8e-01 [ Info: 142 1.6e+04 3.1e-02 0.0e+00 1.0e+00 [ Info: 143 1.6e+04 2.8e-02 0.0e+00 1.5e+00 [ Info: 144 1.6e+04 2.1e-02 0.0e+00 1.7e+00 [ Info: 145 1.6e+04 1.2e-01 0.0e+00 9.5e-01 [ Info: 146 1.6e+04 1.2e-01 0.0e+00 5.1e-01 [ Info: 147 1.6e+04 5.4e-02 0.0e+00 4.1e-01 [ Info: 148 1.6e+04 2.8e-02 0.0e+00 3.7e-01 [ Info: 149 1.6e+04 2.8e-02 0.0e+00 3.7e-01 [ Info: 150 1.6e+04 2.9e-02 0.0e+00 3.7e-01 [ Info: 151 1.6e+04 2.9e-02 0.0e+00 3.7e-01 [ Info: 152 1.6e+04 3.0e-02 0.0e+00 3.7e-01 [ Info: 153 1.6e+04 3.2e-02 0.0e+00 3.7e-01 [ Info: 154 1.6e+04 3.4e-02 0.0e+00 3.6e-01 [ Info: 155 1.6e+04 3.5e-02 0.0e+00 3.4e-01 [ Info: 156 1.6e+04 3.7e-02 0.0e+00 3.2e-01 [ Info: 157 1.6e+04 3.9e-02 0.0e+00 3.0e-01 [ Info: 158 1.6e+04 4.1e-02 0.0e+00 2.9e-01 [ Info: 159 1.6e+04 4.4e-02 0.0e+00 3.0e-01 [ Info: 160 1.6e+04 3.9e-02 0.0e+00 2.9e-01 [ Info: 161 1.6e+04 3.1e-02 0.0e+00 2.9e-01 [ Info: 162 1.6e+04 3.0e-02 0.0e+00 3.0e-01 [ Info: 163 1.6e+04 2.8e-02 0.0e+00 3.0e-01 [ Info: 164 1.6e+04 1.1e-02 0.0e+00 9.3e-01 [ Info: 165 1.6e+04 1.1e-01 0.0e+00 1.9e+00 [ Info: 166 1.6e+04 7.9e-01 0.0e+00 6.1e+00 [ Info: 167 1.6e+04 7.4e-01 0.0e+00 4.9e+00 [ Info: 168 1.6e+04 3.2e-01 0.0e+00 5.7e+00 [ Info: 169 1.6e+04 8.0e-01 0.0e+00 2.5e+00 [ Info: 170 1.6e+04 5.6e-01 0.0e+00 5.7e-01 [ Info: 171 1.6e+04 8.0e-02 0.0e+00 4.9e-01 [ Info: 172 1.6e+04 1.6e-02 0.0e+00 3.7e-01 [ Info: 173 1.6e+04 1.6e-02 0.0e+00 3.6e-01 [ Info: 174 1.6e+04 1.9e-02 0.0e+00 3.6e-01 [ Info: 175 1.6e+04 2.1e-02 0.0e+00 3.6e-01 [ Info: 176 1.6e+04 2.4e-02 0.0e+00 3.5e-01 [ Info: 177 1.6e+04 2.7e-02 0.0e+00 3.5e-01 [ Info: 178 1.6e+04 3.1e-02 0.0e+00 3.8e-01 [ Info: 179 1.6e+04 3.6e-02 0.0e+00 4.1e-01 [ Info: 180 1.6e+04 3.9e-02 0.0e+00 4.3e-01 [ Info: 181 1.6e+04 3.6e-02 0.0e+00 3.1e-01 [ Info: 182 1.6e+04 2.3e-02 0.0e+00 3.0e-01 [ Info: 183 1.6e+04 2.1e-02 0.0e+00 3.2e-01 [ Info: 184 1.6e+04 2.0e-02 0.0e+00 3.2e-01 [ Info: 185 1.6e+04 1.7e-02 0.0e+00 3.3e-01 [ Info: 186 1.6e+04 1.7e-02 0.0e+00 3.3e-01 [ Info: 187 1.6e+04 1.6e-02 0.0e+00 3.3e-01 [ Info: 188 1.6e+04 1.5e-02 0.0e+00 3.2e-01 [ Info: 189 1.6e+04 1.7e-02 0.0e+00 3.2e-01 [ Info: 190 1.6e+04 1.7e-02 0.0e+00 3.0e-01 [ Info: 191 1.6e+04 1.8e-02 0.0e+00 2.9e-01 [ Info: 192 1.6e+04 1.8e-02 0.0e+00 2.7e-01 [ Info: 193 1.6e+04 1.8e-02 0.0e+00 2.6e-01 [ Info: 194 1.6e+04 1.7e-02 0.0e+00 2.5e-01 [ Info: 195 1.6e+04 1.7e-02 0.0e+00 2.5e-01 [ Info: 196 1.6e+04 1.7e-02 0.0e+00 2.4e-01 [ Info: 197 1.6e+04 1.7e-02 0.0e+00 2.4e-01 [ Info: 198 1.6e+04 1.8e-02 0.0e+00 2.4e-01 [ Info: 199 1.6e+04 1.9e-02 0.0e+00 2.4e-01 [ Info: 200 1.6e+04 1.9e-02 0.0e+00 2.4e-01 [ Info: 201 1.6e+04 1.7e-02 0.0e+00 2.3e-01 [ Info: 202 1.6e+04 1.6e-02 0.0e+00 2.6e-01 [ Info: 203 1.6e+04 1.2e-02 0.0e+00 5.6e-01 [ Info: 204 1.6e+04 7.9e-04 0.0e+00 8.6e-01 [ Info: 205 1.6e+04 2.3e-02 0.0e+00 7.0e-01 [ Info: 206 1.6e+04 3.7e-02 0.0e+00 2.7e-01 [ Info: 207 1.6e+04 2.9e-02 0.0e+00 2.5e-01 [ Info: 208 1.6e+04 2.2e-02 0.0e+00 2.4e-01 [ Info: 209 1.6e+04 2.0e-02 0.0e+00 2.4e-01 [ Info: 210 1.6e+04 2.1e-02 0.0e+00 2.5e-01 [ Info: 211 1.6e+04 2.1e-02 0.0e+00 2.6e-01 [ Info: 212 1.6e+04 2.2e-02 0.0e+00 2.6e-01 [ Info: 213 1.6e+04 2.3e-02 0.0e+00 2.7e-01 [ Info: 214 1.6e+04 2.4e-02 0.0e+00 2.8e-01 [ Info: 215 1.6e+04 2.6e-02 0.0e+00 2.9e-01 [ Info: 216 1.6e+04 2.7e-02 0.0e+00 3.0e-01 [ Info: 217 1.6e+04 3.0e-02 0.0e+00 3.2e-01 [ Info: 218 1.6e+04 3.2e-02 0.0e+00 3.3e-01 [ Info: 219 1.6e+04 3.5e-02 0.0e+00 3.5e-01 [ Info: 220 1.6e+04 3.8e-02 0.0e+00 3.7e-01 [ Info: 221 1.6e+04 4.1e-02 0.0e+00 3.9e-01 [ Info: 222 1.6e+04 3.8e-02 0.0e+00 2.9e-01 [ Info: 223 1.6e+04 2.5e-02 0.0e+00 5.2e-01 [ Info: 224 1.6e+04 2.9e-02 0.0e+00 4.6e-01 [ Info: 225 1.6e+04 3.7e-02 0.0e+00 3.5e-01 [ Info: 226 1.6e+04 3.6e-02 0.0e+00 3.7e-01 [ Info: 227 1.6e+04 3.4e-02 0.0e+00 3.9e-01 [ Info: 228 1.6e+04 3.5e-02 0.0e+00 4.2e-01 [ Info: 229 1.6e+04 3.7e-02 0.0e+00 4.4e-01 [ Info: 230 1.6e+04 4.0e-02 0.0e+00 4.7e-01 [ Info: 231 1.6e+04 3.4e-02 0.0e+00 2.0e-01 [ Info: 232 1.6e+04 1.2e-02 0.0e+00 1.9e-01 [ Info: 233 1.6e+04 9.9e-03 0.0e+00 1.8e-01 [ Info: 234 1.6e+04 9.3e-03 0.0e+00 1.8e-01 [ Info: 235 1.6e+04 8.8e-03 0.0e+00 1.8e-01 [ Info: 236 1.6e+04 8.3e-03 0.0e+00 1.7e-01 [ Info: 237 1.6e+04 7.9e-03 0.0e+00 1.7e-01 [ Info: 238 1.6e+04 7.6e-03 0.0e+00 1.7e-01 [ Info: 239 1.6e+04 7.3e-03 0.0e+00 1.6e-01 [ Info: 240 1.6e+04 7.1e-03 0.0e+00 1.6e-01 [ Info: 241 1.6e+04 6.6e-03 0.0e+00 1.6e-01 [ Info: 242 1.6e+04 5.4e-03 0.0e+00 2.7e-01 [ Info: 243 1.6e+04 6.0e-03 0.0e+00 7.5e-01 [ Info: 244 1.6e+04 3.5e-02 0.0e+00 1.4e+00 [ Info: 245 1.6e+04 1.6e-02 0.0e+00 1.4e+00 [ Info: 246 1.6e+04 5.3e-02 0.0e+00 5.7e-01 [ Info: 247 1.6e+04 3.6e-02 0.0e+00 1.6e-01 [ Info: 248 1.6e+04 8.3e-03 0.0e+00 1.5e-01 [ Info: 249 1.6e+04 5.9e-03 0.0e+00 1.4e-01 [ Info: 250 1.6e+04 5.7e-03 0.0e+00 1.4e-01 [ Info: 251 1.6e+04 5.7e-03 0.0e+00 1.4e-01 [ Info: 252 1.6e+04 5.8e-03 0.0e+00 1.4e-01 [ Info: 253 1.6e+04 5.9e-03 0.0e+00 1.4e-01 [ Info: 254 1.6e+04 5.9e-03 0.0e+00 1.4e-01 [ Info: 255 1.6e+04 6.1e-03 0.0e+00 1.4e-01 [ Info: 256 1.6e+04 6.2e-03 0.0e+00 1.5e-01 [ Info: 257 1.6e+04 6.3e-03 0.0e+00 1.5e-01 [ Info: 258 1.6e+04 6.5e-03 0.0e+00 1.6e-01 [ Info: 259 1.6e+04 6.7e-03 0.0e+00 1.6e-01 [ Info: 260 1.6e+04 6.9e-03 0.0e+00 1.7e-01 [ Info: 261 1.6e+04 7.2e-03 0.0e+00 1.7e-01 [ Info: 262 1.6e+04 7.0e-03 0.0e+00 1.7e-01 [ Info: 263 1.6e+04 5.2e-03 0.0e+00 3.5e-01 [ Info: 264 1.6e+04 1.9e-03 0.0e+00 6.2e-01 [ Info: 265 1.6e+04 1.1e-03 0.0e+00 6.9e-01 [ Info: 266 1.6e+04 1.7e-02 0.0e+00 3.2e-01 [ Info: 267 1.6e+04 1.9e-02 0.0e+00 2.1e-01 [ Info: 268 1.6e+04 1.2e-02 0.0e+00 2.1e-01 [ Info: 269 1.6e+04 1.2e-02 0.0e+00 2.2e-01 [ Info: 270 1.6e+04 1.2e-02 0.0e+00 2.4e-01 [ Info: 271 1.6e+04 1.3e-02 0.0e+00 2.5e-01 [ Info: 272 1.6e+04 1.5e-02 0.0e+00 2.7e-01 [ Info: 273 1.6e+04 1.6e-02 0.0e+00 2.9e-01 [ Info: 274 1.6e+04 1.8e-02 0.0e+00 2.2e-01 [ Info: 275 1.6e+04 1.2e-02 0.0e+00 2.1e-01 [ Info: 276 1.6e+04 9.8e-03 0.0e+00 2.0e-01 [ Info: 277 1.6e+04 1.1e-02 0.0e+00 2.0e-01 [ Info: 278 1.6e+04 9.8e-03 0.0e+00 1.9e-01 [ Info: 279 1.6e+04 1.0e-02 0.0e+00 1.9e-01 [ Info: 280 1.6e+04 9.8e-03 0.0e+00 1.8e-01 [ Info: 281 1.6e+04 7.1e-03 0.0e+00 3.2e-01 [ Info: 282 1.6e+04 7.7e-03 0.0e+00 9.2e-01 [ Info: 283 1.6e+04 5.1e-02 0.0e+00 1.6e+00 [ Info: 284 1.6e+04 3.1e-02 0.0e+00 1.8e+00 [ Info: 285 1.6e+04 6.8e-02 0.0e+00 7.7e-01 [ Info: 286 1.6e+04 5.9e-02 0.0e+00 1.6e-01 [ Info: 287 1.6e+04 1.4e-02 0.0e+00 1.6e-01 [ Info: 288 1.6e+04 7.3e-03 0.0e+00 1.5e-01 [ Info: 289 1.6e+04 7.0e-03 0.0e+00 1.5e-01 [ Info: 290 1.6e+04 6.8e-03 0.0e+00 1.5e-01 [ Info: 291 1.6e+04 6.7e-03 0.0e+00 1.5e-01 [ Info: 292 1.6e+04 6.6e-03 0.0e+00 1.4e-01 [ Info: 293 1.6e+04 6.6e-03 0.0e+00 1.4e-01 [ Info: 294 1.6e+04 6.5e-03 0.0e+00 1.4e-01 [ Info: 295 1.6e+04 6.5e-03 0.0e+00 1.4e-01 [ Info: 296 1.6e+04 6.6e-03 0.0e+00 1.4e-01 [ Info: 297 1.6e+04 6.5e-03 0.0e+00 1.3e-01 [ Info: 298 1.6e+04 6.6e-03 0.0e+00 1.3e-01 [ Info: 299 1.6e+04 6.5e-03 0.0e+00 1.3e-01 [ Info: 300 1.6e+04 5.7e-03 0.0e+00 2.0e-01 [ Info: 301 1.6e+04 1.0e-03 0.0e+00 5.5e-01 [ Info: 302 1.6e+04 1.5e-02 0.0e+00 9.7e-01 [ Info: 303 1.6e+04 9.3e-03 0.0e+00 1.1e+00 [ Info: 304 1.6e+04 2.5e-02 0.0e+00 5.1e-01 [ Info: 305 1.6e+04 2.6e-02 0.0e+00 1.3e-01 [ Info: 306 1.6e+04 1.1e-02 0.0e+00 1.3e-01 [ Info: 307 1.6e+04 7.4e-03 0.0e+00 1.4e-01 [ Info: 308 1.6e+04 7.4e-03 0.0e+00 1.5e-01 [ Info: 309 1.6e+04 8.0e-03 0.0e+00 1.6e-01 [ Info: 310 1.6e+04 9.1e-03 0.0e+00 1.7e-01 [ Info: 311 1.6e+04 1.0e-02 0.0e+00 1.9e-01 [ Info: 312 1.6e+04 1.2e-02 0.0e+00 2.2e-01 [ Info: 313 1.6e+04 1.5e-02 0.0e+00 2.5e-01 [ Info: 314 1.6e+04 1.7e-02 0.0e+00 2.8e-01 [ Info: 315 1.6e+04 2.0e-02 0.0e+00 3.1e-01 [ Info: 316 1.6e+04 2.5e-02 0.0e+00 3.6e-01 [ Info: 317 1.6e+04 3.0e-02 0.0e+00 3.9e-01 [ Info: 318 1.6e+04 3.6e-02 0.0e+00 4.8e-01 [ Info: 319 1.6e+04 4.0e-02 0.0e+00 5.5e-01 [ Info: 320 1.6e+04 7.7e-03 0.0e+00 1.8e+00 [ Info: 321 1.6e+04 4.9e-01 0.0e+00 4.5e+00 [ Info: 322 1.6e+04 1.8e+00 0.0e+00 1.2e+01 [ Info: 323 1.6e+04 1.6e+00 0.0e+00 8.2e+00 [ Info: 324 1.6e+04 2.6e+00 0.0e+00 3.6e+00 [ Info: 325 1.6e+04 1.6e+00 0.0e+00 4.7e-01 [ Info: 326 1.6e+04 2.1e-01 0.0e+00 3.5e-01 [ Info: 327 1.6e+04 1.4e-02 0.0e+00 3.5e-01 [ Info: 328 1.6e+04 1.4e-02 0.0e+00 3.3e-01 [ Info: 329 1.6e+04 1.4e-02 0.0e+00 3.2e-01 [ Info: 330 1.6e+04 1.4e-02 0.0e+00 3.0e-01 [ Info: 331 1.6e+04 1.5e-02 0.0e+00 2.9e-01 [ Info: 332 1.6e+04 1.6e-02 0.0e+00 2.9e-01 [ Info: 333 1.6e+04 1.7e-02 0.0e+00 2.9e-01 [ Info: 334 1.6e+04 1.8e-02 0.0e+00 2.4e-01 [ Info: 335 1.6e+04 1.4e-02 0.0e+00 2.3e-01 [ Info: 336 1.6e+04 1.2e-02 0.0e+00 2.0e-01 [ Info: 337 1.6e+04 1.1e-02 0.0e+00 2.1e-01 [ Info: 338 1.6e+04 1.1e-02 0.0e+00 2.0e-01 [ Info: 339 1.6e+04 1.1e-02 0.0e+00 2.7e-01 [ Info: 340 1.6e+04 8.5e-03 0.0e+00 4.1e-01 [ Info: 341 1.6e+04 2.8e-04 0.0e+00 7.9e-01 [ Info: 342 1.6e+04 4.8e-03 0.0e+00 9.3e-01 [ Info: 343 1.6e+04 1.6e-02 0.0e+00 6.1e-01 [ Info: 344 1.6e+04 2.5e-02 0.0e+00 2.4e-01 [ Info: 345 1.6e+04 1.5e-02 0.0e+00 1.9e-01 [ Info: 346 1.6e+04 9.1e-03 0.0e+00 2.0e-01 [ Info: 347 1.6e+04 9.9e-03 0.0e+00 2.0e-01 [ Info: 348 1.6e+04 9.2e-03 0.0e+00 2.0e-01 [ Info: 349 1.6e+04 9.7e-03 0.0e+00 2.0e-01 [ Info: 350 1.6e+04 1.0e-02 0.0e+00 2.0e-01 [ Info: 351 1.6e+04 1.1e-02 0.0e+00 2.1e-01 [ Info: 352 1.6e+04 1.1e-02 0.0e+00 2.2e-01 [ Info: 353 1.6e+04 1.2e-02 0.0e+00 2.3e-01 [ Info: 354 1.6e+04 9.5e-03 0.0e+00 1.9e-01 [ Info: 355 1.6e+04 6.6e-03 0.0e+00 1.8e-01 [ Info: 356 1.6e+04 7.1e-03 0.0e+00 1.7e-01 [ Info: 357 1.6e+04 5.5e-03 0.0e+00 1.6e-01 [ Info: 358 1.6e+04 4.6e-03 0.0e+00 1.7e-01 [ Info: 359 1.6e+04 2.1e-03 0.0e+00 2.9e-01 [ Info: 360 1.6e+04 5.1e-03 0.0e+00 5.8e-01 [ Info: 361 1.6e+04 2.0e-02 0.0e+00 1.0e+00 [ Info: 362 1.6e+04 6.4e-03 0.0e+00 1.1e+00 [ Info: 363 1.6e+04 3.1e-02 0.0e+00 4.6e-01 [ Info: 364 1.6e+04 2.2e-02 0.0e+00 1.2e-01 [ Info: 365 1.6e+04 4.5e-03 0.0e+00 1.1e-01 [ Info: 366 1.6e+04 2.9e-03 0.0e+00 1.1e-01 [ Info: 367 1.6e+04 2.7e-03 0.0e+00 1.0e-01 [ Info: 368 1.6e+04 2.7e-03 0.0e+00 9.7e-02 [ Info: 369 1.6e+04 2.7e-03 0.0e+00 9.3e-02 [ Info: 370 1.6e+04 2.9e-03 0.0e+00 9.1e-02 [ Info: 371 1.6e+04 3.0e-03 0.0e+00 8.9e-02 [ Info: 372 1.6e+04 3.1e-03 0.0e+00 8.7e-02 [ Info: 373 1.6e+04 3.3e-03 0.0e+00 9.1e-02 [ Info: 374 1.6e+04 3.6e-03 0.0e+00 1.0e-01 [ Info: 375 1.6e+04 3.6e-03 0.0e+00 1.1e-01 [ Info: 376 1.6e+04 3.1e-03 0.0e+00 1.2e-01 [ Info: 377 1.6e+04 3.4e-03 0.0e+00 1.2e-01 [ Info: 378 1.6e+04 3.4e-03 0.0e+00 1.3e-01 [ Info: 379 1.6e+04 3.4e-03 0.0e+00 1.3e-01 [ Info: 380 1.6e+04 3.2e-03 0.0e+00 1.4e-01 [ Info: 381 1.6e+04 9.8e-05 0.0e+00 3.6e-01 [ Info: 382 1.6e+04 5.9e-03 0.0e+00 6.7e-01 [ Info: 383 1.6e+04 4.2e-03 0.0e+00 6.9e-01 [ Info: 384 1.6e+04 1.1e-02 0.0e+00 3.5e-01 [ Info: 385 1.6e+04 1.5e-02 0.0e+00 2.1e-01 [ Info: 386 1.6e+04 8.2e-03 0.0e+00 1.4e-01 [ Info: 387 1.6e+04 4.5e-03 0.0e+00 1.5e-01 [ Info: 388 1.6e+04 4.6e-03 0.0e+00 1.5e-01 [ Info: 389 1.6e+04 4.8e-03 0.0e+00 1.5e-01 [ Info: 390 1.6e+04 5.0e-03 0.0e+00 1.5e-01 [ Info: 391 1.6e+04 5.4e-03 0.0e+00 1.5e-01 [ Info: 392 1.6e+04 5.9e-03 0.0e+00 1.5e-01 [ Info: 393 1.6e+04 6.4e-03 0.0e+00 1.5e-01 [ Info: 394 1.6e+04 7.1e-03 0.0e+00 1.5e-01 [ Info: 395 1.6e+04 7.9e-03 0.0e+00 1.5e-01 [ Info: 396 1.6e+04 8.8e-03 0.0e+00 1.6e-01 [ Info: 397 1.6e+04 9.6e-03 0.0e+00 1.8e-01 [ Info: 398 1.6e+04 1.1e-02 0.0e+00 1.9e-01 [ Info: 399 1.6e+04 1.2e-02 0.0e+00 2.2e-01 [ Info: 400 1.6e+04 1.4e-02 0.0e+00 2.2e-01 [ Info: 401 1.6e+04 9.5e-03 0.0e+00 2.9e-01 [ Info: 402 1.6e+04 5.2e-03 0.0e+00 5.3e-01 [ Info: 403 1.6e+04 8.7e-02 0.0e+00 1.7e+00 [ Info: 404 1.6e+04 1.6e-01 0.0e+00 2.5e+00 [ Info: 405 1.6e+04 6.5e-02 0.0e+00 2.8e+00 [ Info: 406 1.6e+04 2.0e-01 0.0e+00 1.2e+00 [ Info: 407 1.6e+04 1.6e-01 0.0e+00 2.2e-01 [ Info: 408 1.6e+04 3.1e-02 0.0e+00 2.4e-01 [ Info: 409 1.6e+04 1.5e-02 0.0e+00 2.7e-01 [ Info: 410 1.6e+04 1.7e-02 0.0e+00 2.9e-01 [ Info: 411 1.6e+04 1.3e-02 0.0e+00 2.5e-01 [ Info: 412 1.6e+04 9.5e-03 0.0e+00 2.3e-01 [ Info: 413 1.6e+04 8.8e-03 0.0e+00 2.2e-01 [ Info: 414 1.6e+04 8.0e-03 0.0e+00 2.1e-01 [ Info: 415 1.6e+04 7.5e-03 0.0e+00 2.0e-01 [ Info: 416 1.6e+04 7.1e-03 0.0e+00 1.8e-01 [ Info: 417 1.6e+04 6.7e-03 0.0e+00 1.7e-01 [ Info: 418 1.6e+04 6.4e-03 0.0e+00 1.6e-01 [ Info: 419 1.6e+04 6.2e-03 0.0e+00 1.5e-01 [ Info: 420 1.6e+04 6.1e-03 0.0e+00 1.5e-01 [ Info: 421 1.6e+04 5.9e-03 0.0e+00 1.4e-01 [ Info: 422 1.6e+04 5.8e-03 0.0e+00 1.3e-01 [ Info: 423 1.6e+04 5.5e-03 0.0e+00 1.3e-01 [ Info: 424 1.6e+04 5.5e-03 0.0e+00 1.4e-01 [ Info: 425 1.6e+04 4.9e-03 0.0e+00 1.4e-01 [ Info: 426 1.6e+04 2.6e-03 0.0e+00 2.7e-01 [ Info: 427 1.6e+04 1.1e-02 0.0e+00 6.4e-01 [ Info: 428 1.6e+04 4.4e-02 0.0e+00 1.1e+00 [ Info: 429 1.6e+04 3.0e-02 0.0e+00 1.4e+00 [ Info: 430 1.6e+04 4.7e-02 0.0e+00 8.1e-01 [ Info: 431 1.6e+04 5.5e-02 0.0e+00 1.5e-01 [ Info: 432 1.6e+04 1.3e-02 0.0e+00 1.9e-01 [ Info: 433 1.6e+04 5.0e-03 0.0e+00 1.8e-01 [ Info: 434 1.6e+04 4.6e-03 0.0e+00 1.7e-01 [ Info: 435 1.6e+04 4.7e-03 0.0e+00 1.7e-01 [ Info: 436 1.6e+04 5.2e-03 0.0e+00 1.8e-01 [ Info: 437 1.6e+04 5.6e-03 0.0e+00 1.9e-01 [ Info: 438 1.6e+04 6.1e-03 0.0e+00 2.0e-01 [ Info: 439 1.6e+04 6.7e-03 0.0e+00 2.1e-01 [ Info: 440 1.6e+04 7.5e-03 0.0e+00 2.2e-01 [ Info: 441 1.6e+04 7.0e-03 0.0e+00 1.0e-01 [ Info: 442 1.6e+04 3.4e-03 0.0e+00 9.7e-02 [ Info: 443 1.6e+04 3.4e-03 0.0e+00 9.8e-02 [ Info: 444 1.6e+04 3.4e-03 0.0e+00 9.5e-02 [ Info: 445 1.6e+04 2.9e-03 0.0e+00 9.5e-02 [ Info: 446 1.6e+04 2.7e-03 0.0e+00 1.1e-01 [ Info: 447 1.6e+04 1.6e-03 0.0e+00 2.1e-01 [ Info: 448 1.6e+04 3.5e-03 0.0e+00 4.2e-01 [ Info: 449 1.6e+04 3.4e-02 0.0e+00 7.8e-01 [ Info: 450 1.6e+04 5.9e-02 0.0e+00 1.3e+00 [ Info: 451 1.6e+04 4.7e-03 0.0e+00 1.4e+00 [ Info: 452 1.6e+04 7.1e-02 0.0e+00 6.1e-01 [ Info: 453 1.6e+04 3.8e-02 0.0e+00 8.5e-02 [ Info: 454 1.6e+04 3.8e-03 0.0e+00 6.5e-02 [ Info: 455 1.6e+04 1.2e-03 0.0e+00 6.4e-02 [ Info: 456 1.6e+04 8.7e-04 0.0e+00 6.5e-02 [ Info: 457 1.6e+04 8.3e-04 0.0e+00 6.7e-02 [ Info: 458 1.6e+04 8.3e-04 0.0e+00 6.8e-02 [ Info: 459 1.6e+04 1.2e-03 1.7e-09 6.7e-02 [ Info: 460 1.6e+04 4.9e-04 0.0e+00 6.7e-02 [ Info: 461 1.6e+04 8.4e-04 0.0e+00 6.7e-02 [ Info: 462 1.6e+04 8.7e-04 0.0e+00 6.6e-02 [ Info: 463 1.6e+04 9.0e-04 0.0e+00 6.6e-02 [ Info: 464 1.6e+04 9.4e-04 0.0e+00 6.5e-02 [ Info: 465 1.6e+04 1.0e-03 0.0e+00 6.4e-02 [ Info: 466 1.6e+04 1.0e-03 0.0e+00 6.1e-02 [ Info: 467 1.6e+04 1.1e-03 0.0e+00 5.4e-02 [ Info: 468 1.6e+04 9.6e-04 0.0e+00 6.9e-02 [ Info: 469 1.6e+04 8.0e-04 0.0e+00 6.6e-02 [ Info: 470 1.6e+04 5.5e-04 0.0e+00 1.1e-01 [ Info: 471 1.6e+04 5.0e-04 0.0e+00 1.5e-01 [ Info: 472 1.6e+04 6.6e-04 0.0e+00 2.3e-01 [ Info: 473 1.6e+04 1.2e-04 0.0e+00 2.9e-01 [ Info: 474 1.6e+04 1.1e-03 0.0e+00 2.7e-01 [ Info: 475 1.6e+04 1.9e-03 0.0e+00 1.4e-01 [ Info: 476 1.6e+04 6.0e-04 0.0e+00 8.0e-02 [ Info: 477 1.6e+04 7.7e-04 0.0e+00 5.4e-02 [ Info: 478 1.6e+04 4.8e-04 0.0e+00 5.4e-02 [ Info: 479 1.6e+04 5.8e-05 0.0e+00 5.4e-02 [ Info: 480 1.6e+04 8.0e-04 0.0e+00 5.2e-02 [ Info: 481 1.6e+04 4.4e-04 0.0e+00 5.2e-02 [ Info: 482 1.6e+04 4.9e-04 0.0e+00 5.2e-02 [ Info: 483 1.6e+04 5.5e-04 0.0e+00 5.2e-02 [ Info: 484 1.6e+04 6.2e-04 0.0e+00 5.3e-02 [ Info: 485 1.6e+04 6.6e-04 0.0e+00 5.9e-02 [ Info: 486 1.6e+04 3.0e-04 0.0e+00 5.8e-02 [ Info: 487 1.6e+04 2.0e-04 0.0e+00 5.1e-02 [ Info: 488 1.6e+04 2.6e-04 0.0e+00 4.8e-02 [ Info: 489 1.6e+04 4.2e-04 0.0e+00 4.9e-02 [ Info: 490 1.6e+04 5.0e-04 0.0e+00 4.3e-02 [ Info: 491 1.6e+04 1.5e-03 6.6e-09 3.8e-02 [ Info: 492 1.6e+04 6.0e-04 0.0e+00 2.9e-02 [ Info: 493 1.6e+04 5.8e-04 3.0e-09 2.4e-02 [ Info: 494 1.6e+04 2.5e-04 0.0e+00 1.9e-02 [ Info: 495 1.6e+04 1.3e-04 0.0e+00 2.1e-02 [ Info: 496 1.6e+04 2.5e-04 0.0e+00 1.8e-02 [ Info: 497 1.6e+04 4.4e-04 0.0e+00 1.9e-02 [ Info: 498 1.6e+04 7.0e-05 0.0e+00 1.7e-02 [ Info: 499 1.6e+04 7.5e-05 0.0e+00 1.7e-02 [ Info: 500 1.6e+04 5.9e-05 0.0e+00 1.7e-02 [ Info: 501 1.6e+04 5.8e-05 0.0e+00 1.7e-02 [ Info: 502 1.6e+04 6.0e-05 0.0e+00 1.6e-02 [ Info: 503 1.6e+04 4.9e-05 0.0e+00 1.6e-02 [ Info: 504 1.6e+04 5.3e-05 0.0e+00 1.6e-02 [ Info: 505 1.6e+04 1.9e-05 0.0e+00 2.9e-02 [ Info: 506 1.6e+04 1.0e-04 0.0e+00 6.6e-02 [ Info: 507 1.6e+04 2.2e-04 0.0e+00 9.9e-02 [ Info: 508 1.6e+04 1.5e-06 0.0e+00 8.3e-02 [ Info: 509 1.6e+04 3.2e-04 0.0e+00 2.7e-02 [ Info: 510 1.6e+04 2.8e-04 0.0e+00 1.7e-02 [ Info: 511 1.6e+04 5.9e-05 0.0e+00 1.7e-02 [ Info: 512 1.6e+04 5.9e-05 0.0e+00 1.7e-02 [ Info: 513 1.6e+04 5.9e-05 0.0e+00 1.7e-02 [ Info: 514 1.6e+04 6.1e-05 0.0e+00 1.7e-02 [ Info: 515 1.6e+04 6.7e-05 0.0e+00 1.7e-02 [ Info: 516 1.6e+04 7.3e-05 0.0e+00 1.8e-02 [ Info: 517 1.6e+04 7.9e-05 0.0e+00 2.0e-02 [ Info: 518 1.6e+04 8.6e-05 0.0e+00 2.1e-02 [ Info: 519 1.6e+04 9.3e-05 0.0e+00 2.2e-02 [ Info: 520 1.6e+04 1.0e-04 0.0e+00 2.3e-02 [ Info: 521 1.6e+04 1.1e-04 0.0e+00 2.5e-02 [ Info: 522 1.6e+04 1.2e-04 0.0e+00 2.6e-02 [ Info: 523 1.6e+04 1.1e-04 0.0e+00 2.8e-02 [ Info: 524 1.6e+04 1.2e-04 0.0e+00 2.9e-02 [ Info: 525 1.6e+04 1.1e-04 0.0e+00 3.0e-02 [ Info: 526 1.6e+04 1.4e-04 0.0e+00 3.2e-02 [ Info: 527 1.6e+04 8.9e-05 0.0e+00 3.2e-02 [ Info: 528 1.6e+04 1.7e-04 0.0e+00 3.4e-02 [ Info: 529 1.6e+04 1.6e-04 0.0e+00 3.5e-02 [ Info: 530 1.6e+04 9.0e-05 0.0e+00 3.4e-02 [ Info: 531 1.6e+04 2.1e-04 0.0e+00 3.8e-02 [ Info: 532 1.6e+04 2.2e-04 0.0e+00 3.8e-02 [ Info: 533 1.6e+04 2.4e-04 0.0e+00 3.9e-02 [ Info: 534 1.6e+04 2.0e-04 0.0e+00 4.0e-02 [ Info: 535 1.6e+04 2.2e-04 0.0e+00 4.1e-02 [ Info: 536 1.6e+04 2.3e-04 0.0e+00 4.2e-02 [ Info: 537 1.6e+04 2.4e-04 0.0e+00 4.3e-02 [ Info: 538 1.6e+04 2.5e-04 0.0e+00 4.5e-02 [ Info: 539 1.6e+04 2.6e-04 0.0e+00 4.6e-02 [ Info: 540 1.6e+04 2.8e-04 0.0e+00 4.7e-02 [ Info: 541 1.6e+04 2.9e-04 0.0e+00 4.8e-02 [ Info: 542 1.6e+04 3.1e-04 0.0e+00 4.9e-02 [ Info: 543 1.6e+04 3.2e-04 0.0e+00 5.1e-02 [ Info: 544 1.6e+04 3.4e-04 0.0e+00 5.2e-02 [ Info: 545 1.6e+04 3.5e-04 0.0e+00 5.3e-02 [ Info: 546 1.6e+04 3.7e-04 0.0e+00 5.4e-02 [ Info: 547 1.6e+04 3.9e-04 0.0e+00 5.6e-02 [ Info: 548 1.6e+04 4.1e-04 0.0e+00 5.7e-02 [ Info: 549 1.6e+04 4.3e-04 0.0e+00 5.8e-02 [ Info: 550 1.6e+04 4.4e-04 0.0e+00 5.9e-02 [ Info: 551 1.6e+04 3.0e-04 0.0e+00 8.1e-02 [ Info: 552 1.6e+04 6.9e-06 0.0e+00 1.4e-01 [ Info: 553 1.6e+04 1.3e-04 0.0e+00 1.5e-01 [ Info: 554 1.6e+04 1.1e-03 0.0e+00 6.8e-02 [ Info: 555 1.6e+04 8.6e-04 0.0e+00 6.4e-02 [ Info: 556 1.6e+04 6.5e-04 0.0e+00 6.4e-02 [ Info: 557 1.6e+04 5.8e-04 0.0e+00 6.5e-02 [ Info: 558 1.6e+04 5.9e-04 0.0e+00 6.6e-02 [ Info: 559 1.6e+04 6.0e-04 0.0e+00 6.7e-02 [ Info: 560 1.6e+04 6.2e-04 0.0e+00 6.8e-02 [ Info: 561 1.6e+04 6.4e-04 0.0e+00 6.9e-02 [ Info: 562 1.6e+04 6.6e-04 0.0e+00 7.0e-02 [ Info: 563 1.6e+04 6.8e-04 0.0e+00 7.1e-02 [ Info: 564 1.6e+04 7.0e-04 0.0e+00 7.1e-02 [ Info: 565 1.6e+04 7.2e-04 0.0e+00 7.2e-02 [ Info: 566 1.6e+04 7.4e-04 0.0e+00 7.3e-02 [ Info: 567 1.6e+04 7.5e-04 0.0e+00 7.4e-02 [ Info: 568 1.6e+04 7.7e-04 0.0e+00 7.4e-02 [ Info: 569 1.6e+04 7.8e-04 0.0e+00 7.5e-02 [ Info: 570 1.6e+04 8.3e-04 0.0e+00 7.6e-02 [ Info: 571 1.6e+04 6.9e-04 0.0e+00 7.5e-02 [ Info: 572 1.6e+04 3.6e-04 0.0e+00 1.4e-01 [ Info: 573 1.6e+04 8.3e-04 0.0e+00 2.6e-01 [ Info: 574 1.6e+04 1.6e-04 0.0e+00 2.7e-01 [ Info: 575 1.6e+04 2.4e-03 0.0e+00 1.1e-01 [ Info: 576 1.6e+04 2.3e-03 0.0e+00 7.8e-02 [ Info: 577 1.6e+04 1.2e-03 0.0e+00 7.5e-02 [ Info: 578 1.6e+04 9.8e-04 0.0e+00 7.5e-02 [ Info: 579 1.6e+04 9.2e-04 0.0e+00 7.6e-02 [ Info: 580 1.6e+04 9.3e-04 0.0e+00 7.6e-02 [ Info: 581 1.6e+04 9.5e-04 0.0e+00 7.6e-02 [ Info: 582 1.6e+04 9.6e-04 0.0e+00 7.7e-02 [ Info: 583 1.6e+04 9.7e-04 0.0e+00 7.7e-02 [ Info: 584 1.6e+04 9.8e-04 0.0e+00 7.7e-02 [ Info: 585 1.6e+04 9.9e-04 0.0e+00 7.7e-02 [ Info: 586 1.6e+04 1.0e-03 0.0e+00 7.7e-02 [ Info: 587 1.6e+04 1.0e-03 0.0e+00 7.7e-02 [ Info: 588 1.6e+04 1.0e-03 0.0e+00 7.7e-02 [ Info: 589 1.6e+04 1.0e-03 0.0e+00 7.7e-02 [ Info: 590 1.6e+04 1.0e-03 0.0e+00 7.6e-02 [ Info: 591 1.6e+04 1.1e-03 0.0e+00 7.6e-02 [ Info: 592 1.6e+04 1.1e-03 0.0e+00 7.6e-02 [ Info: 593 1.6e+04 1.1e-03 0.0e+00 7.6e-02 [ Info: 594 1.6e+04 8.0e-04 0.0e+00 9.9e-02 [ Info: 595 1.6e+04 2.2e-04 0.0e+00 1.8e-01 [ Info: 596 1.6e+04 3.3e-04 0.0e+00 2.0e-01 [ Info: 597 1.6e+04 1.7e-03 0.0e+00 1.3e-01 [ Info: 598 1.6e+04 1.9e-03 0.0e+00 7.8e-02 [ Info: 599 1.6e+04 1.6e-03 0.0e+00 7.9e-02 [ Info: 600 1.6e+04 1.3e-03 0.0e+00 7.8e-02 [ Info: 601 1.6e+04 1.2e-03 0.0e+00 7.9e-02 [ Info: 602 1.6e+04 1.2e-03 0.0e+00 8.1e-02 [ Info: 603 1.6e+04 1.2e-03 0.0e+00 8.2e-02 [ Info: 604 1.6e+04 1.2e-03 0.0e+00 8.4e-02 [ Info: 605 1.6e+04 1.2e-03 0.0e+00 8.5e-02 [ Info: 606 1.6e+04 1.3e-03 0.0e+00 8.7e-02 [ Info: 607 1.6e+04 1.3e-03 0.0e+00 8.9e-02 [ Info: 608 1.6e+04 1.4e-03 0.0e+00 9.0e-02 [ Info: 609 1.6e+04 1.4e-03 0.0e+00 9.3e-02 [ Info: 610 1.6e+04 1.5e-03 0.0e+00 9.5e-02 [ Info: 611 1.6e+04 1.5e-03 0.0e+00 9.8e-02 [ Info: 612 1.6e+04 1.6e-03 0.0e+00 1.0e-01 [ Info: 613 1.6e+04 1.7e-03 0.0e+00 1.0e-01 [ Info: 614 1.6e+04 1.8e-03 0.0e+00 1.1e-01 [ Info: 615 1.6e+04 2.0e-03 0.0e+00 1.1e-01 [ Info: 616 1.6e+04 2.1e-03 0.0e+00 1.2e-01 [ Info: 617 1.6e+04 2.3e-03 0.0e+00 1.2e-01 [ Info: 618 1.6e+04 2.3e-03 0.0e+00 1.3e-01 [ Info: 619 1.6e+04 2.7e-03 0.0e+00 1.3e-01 [ Info: 620 1.6e+04 2.9e-03 0.0e+00 1.4e-01 [ Info: 621 1.6e+04 3.1e-03 0.0e+00 1.4e-01 [ Info: 622 1.6e+04 3.4e-03 0.0e+00 1.5e-01 [ Info: 623 1.6e+04 3.7e-03 0.0e+00 1.6e-01 [ Info: 624 1.6e+04 4.1e-03 0.0e+00 1.7e-01 [ Info: 625 1.6e+04 2.6e-03 0.0e+00 9.5e-02 [ Info: 626 1.6e+04 1.4e-03 0.0e+00 8.7e-02 [ Info: 627 1.6e+04 1.6e-03 0.0e+00 8.3e-02 [ Info: 628 1.6e+04 1.3e-03 0.0e+00 7.8e-02 [ Info: 629 1.6e+04 1.2e-03 0.0e+00 7.6e-02 [ Info: 630 1.6e+04 5.7e-04 0.0e+00 1.3e-01 [ Info: 631 1.6e+04 2.5e-04 0.0e+00 1.9e-01 [ Info: 632 1.6e+04 6.5e-04 0.0e+00 1.4e-01 [ Info: 633 1.6e+04 1.8e-03 0.0e+00 8.6e-02 [ Info: 634 1.6e+04 1.3e-03 0.0e+00 6.3e-02 [ Info: 635 1.6e+04 8.2e-04 0.0e+00 5.8e-02 [ Info: 636 1.6e+04 7.6e-04 0.0e+00 5.7e-02 [ Info: 637 1.6e+04 6.5e-04 0.0e+00 5.5e-02 [ Info: 638 1.6e+04 5.9e-04 0.0e+00 5.2e-02 [ Info: 639 1.6e+04 9.2e-04 3.0e-09 5.0e-02 [ Info: 640 1.6e+04 1.3e-04 0.0e+00 4.8e-02 [ Info: 641 1.6e+04 4.7e-04 0.0e+00 4.7e-02 [ Info: 642 1.6e+04 4.4e-04 0.0e+00 4.5e-02 [ Info: 643 1.6e+04 4.0e-04 0.0e+00 4.3e-02 [ Info: 644 1.6e+04 3.8e-04 0.0e+00 4.1e-02 [ Info: 645 1.6e+04 3.5e-04 0.0e+00 4.0e-02 [ Info: 646 1.6e+04 2.1e-04 0.0e+00 3.9e-02 [ Info: 647 1.6e+04 1.3e-03 4.9e-09 3.8e-02 [ Info: 648 1.6e+04 1.9e-04 0.0e+00 3.6e-02 [ Info: 649 1.6e+04 2.5e-04 0.0e+00 3.5e-02 [ Info: 650 1.6e+04 1.8e-04 0.0e+00 3.6e-02 [ Info: 651 1.6e+04 9.2e-05 0.0e+00 6.3e-02 [ Info: 652 1.6e+04 1.3e-04 0.0e+00 6.3e-02 [ Info: 653 1.6e+04 3.2e-04 0.0e+00 3.2e-02 [ Info: 654 1.6e+04 2.9e-04 0.0e+00 2.9e-02 [ Info: 655 1.6e+04 2.1e-04 0.0e+00 2.8e-02 [ Info: 656 1.6e+04 1.6e-04 0.0e+00 2.7e-02 [ Info: 657 1.6e+04 6.3e-05 0.0e+00 2.6e-02 [ Info: 658 1.6e+04 1.3e-04 0.0e+00 2.6e-02 [ Info: 659 1.6e+04 3.0e-04 7.2e-10 2.5e-02 [ Info: 660 1.6e+04 5.7e-04 6.5e-09 2.4e-02 [ Info: 661 1.6e+04 4.4e-04 0.0e+00 2.3e-02 [ Info: 662 1.6e+04 6.8e-04 7.0e-09 2.3e-02 [ Info: 663 1.6e+04 4.9e-04 0.0e+00 2.2e-02 [ Info: 664 1.6e+04 1.9e-04 0.0e+00 2.1e-02 [ Info: 665 1.6e+04 3.5e-04 0.0e+00 2.0e-02 [ Info: 666 1.6e+04 4.6e-06 0.0e+00 2.0e-02 [ Info: 667 1.6e+04 1.5e-04 0.0e+00 1.9e-02 [ Info: 668 1.6e+04 6.5e-05 0.0e+00 1.9e-02 [ Info: 669 1.6e+04 9.0e-05 1.6e-10 1.8e-02 [ Info: 670 1.6e+04 4.9e-04 5.8e-09 2.2e-02 [ Info: 671 1.6e+04 5.4e-04 0.0e+00 3.9e-02 [ Info: 672 1.6e+04 3.2e-05 0.0e+00 4.0e-02 [ Info: 673 1.6e+04 1.2e-04 0.0e+00 2.2e-02 [ Info: 674 1.6e+04 4.6e-05 0.0e+00 1.4e-02 [ Info: 675 1.6e+04 5.6e-04 5.7e-09 1.5e-02 [ Info: 676 1.6e+04 4.3e-04 0.0e+00 1.4e-02 [ Info: 677 1.6e+04 4.3e-05 0.0e+00 1.3e-02 [ Info: 678 1.6e+04 4.3e-04 0.0e+00 1.3e-02 [ Info: 679 1.6e+04 4.9e-04 0.0e+00 1.3e-02 [ Info: 680 1.6e+04 2.4e-04 0.0e+00 1.2e-02 [ Info: 681 1.6e+04 2.9e-04 0.0e+00 1.2e-02 [ Info: 682 1.6e+04 4.8e-04 0.0e+00 1.1e-02 [ Info: 683 1.6e+04 5.3e-04 0.0e+00 1.1e-02 [ Info: 684 1.6e+04 2.1e-05 0.0e+00 1.1e-02 [ Info: 685 1.6e+04 6.8e-05 0.0e+00 1.0e-02 [ Info: 686 1.6e+04 1.0e-04 0.0e+00 1.0e-02 [ Info: 687 1.6e+04 1.2e-04 1.3e-09 9.7e-03 [ Info: 688 1.6e+04 9.0e-05 0.0e+00 9.4e-03 [ Info: 689 1.6e+04 1.6e-05 0.0e+00 9.1e-03 [ Info: 690 1.6e+04 7.1e-06 0.0e+00 8.7e-03 [ Info: 691 1.6e+04 6.2e-04 7.6e-09 1.4e-02 [ Info: 692 1.6e+04 6.1e-04 0.0e+00 2.1e-02 [ Info: 693 1.6e+04 6.9e-07 0.0e+00 1.6e-02 [ Info: 694 1.6e+04 2.5e-05 0.0e+00 7.8e-03 [ Info: 695 1.6e+04 1.9e-05 0.0e+00 7.3e-03 [ Info: 696 1.6e+04 3.1e-04 0.0e+00 7.1e-03 [ Info: 697 1.6e+04 3.3e-04 0.0e+00 6.8e-03 [ Info: 698 1.6e+04 9.6e-05 0.0e+00 6.6e-03 [ Info: 699 1.6e+04 4.8e-04 0.0e+00 6.4e-03 [ Info: 700 1.6e+04 6.0e-04 0.0e+00 6.2e-03 [ Info: 701 1.6e+04 5.2e-04 0.0e+00 6.0e-03 [ Info: 702 1.6e+04 5.3e-04 0.0e+00 5.8e-03 [ Info: 703 1.6e+04 6.0e-06 0.0e+00 5.6e-03 [ Info: 704 1.6e+04 4.2e-04 0.0e+00 5.5e-03 [ Info: 705 1.6e+04 4.3e-04 0.0e+00 5.3e-03 [ Info: 706 1.6e+04 5.6e-04 0.0e+00 5.1e-03 [ Info: 707 1.6e+04 5.7e-04 0.0e+00 5.0e-03 [ Info: 708 1.6e+04 4.9e-04 0.0e+00 4.8e-03 [ Info: 709 1.6e+04 1.2e-04 0.0e+00 5.3e-03 [ Info: 710 1.6e+04 3.7e-04 0.0e+00 1.4e-02 [ Info: 711 1.6e+04 2.9e-05 0.0e+00 2.4e-02 [ Info: 712 1.6e+04 1.7e-05 0.0e+00 2.3e-02 [ Info: 713 1.6e+04 3.3e-05 0.0e+00 8.7e-03 [ Info: 714 1.6e+04 3.1e-05 0.0e+00 3.8e-03 [ Info: 715 1.6e+04 4.8e-04 0.0e+00 3.7e-03 [ Info: 716 1.6e+04 2.9e-04 0.0e+00 3.5e-03 [ Info: 717 1.6e+04 7.8e-04 0.0e+00 3.4e-03 [ Info: 718 1.6e+04 6.6e-04 8.3e-09 3.3e-03 [ Info: 719 1.6e+04 6.5e-04 0.0e+00 3.3e-03 [ Info: 720 1.6e+04 8.5e-06 0.0e+00 3.2e-03 [ Info: 721 1.6e+04 1.7e-04 2.0e-09 3.1e-03 [ Info: 722 1.6e+04 3.9e-04 7.0e-09 3.0e-03 [ Info: 723 1.6e+04 5.5e-04 0.0e+00 2.9e-03 [ Info: 724 1.6e+04 2.9e-04 3.7e-09 2.8e-03 [ Info: 725 1.6e+04 9.9e-04 0.0e+00 2.8e-03 [ Info: 726 1.6e+04 7.8e-04 9.8e-10 2.7e-03 [ Info: 727 1.6e+04 8.0e-05 0.0e+00 2.5e-03 [ Info: 728 1.6e+04 1.7e-06 0.0e+00 3.9e-03 [ Info: 729 1.6e+04 1.2e-05 0.0e+00 1.1e-02 [ Info: 730 1.6e+04 2.0e-05 0.0e+00 1.9e-02 [ Info: 731 1.6e+04 1.0e-05 0.0e+00 1.9e-02 [ Info: 732 1.6e+04 7.6e-05 3.2e-10 8.3e-03 [ Info: 733 1.6e+04 3.4e-05 0.0e+00 1.9e-03 [ Info: 734 1.6e+04 6.4e-04 8.0e-09 1.9e-03 [ Info: 735 1.6e+04 6.3e-04 0.0e+00 1.9e-03 [ Info: 736 1.6e+04 3.6e-04 0.0e+00 1.8e-03 [ Info: 737 1.6e+04 2.0e-04 0.0e+00 1.8e-03 [ Info: 738 1.6e+04 6.4e-05 0.0e+00 1.7e-03 [ Info: 739 1.6e+04 5.6e-05 0.0e+00 1.7e-03 [ Info: 740 1.6e+04 1.7e-04 0.0e+00 1.6e-03 [ Info: 741 1.6e+04 2.9e-04 0.0e+00 1.6e-03 [ Info: 742 1.6e+04 3.3e-04 2.9e-09 1.5e-03 [ Info: 743 1.6e+04 1.2e-04 4.4e-09 1.5e-03 [ Info: 744 1.6e+04 6.9e-05 5.3e-09 1.5e-03 [ Info: 745 1.6e+04 2.3e-04 2.4e-09 1.4e-03 [ Info: 746 1.6e+04 4.7e-07 2.4e-09 1.4e-03 [ Info: 747 1.6e+04 1.9e-04 0.0e+00 1.3e-03 [ Info: 748 1.6e+04 7.3e-04 9.3e-09 1.3e-03 [ Info: 749 1.6e+04 7.3e-04 0.0e+00 1.2e-03 [ Info: 750 1.6e+04 1.6e-06 0.0e+00 1.7e-03 [ Info: 751 1.6e+04 4.9e-07 0.0e+00 2.4e-03 [ Info: 752 1.6e+04 3.3e-06 0.0e+00 6.9e-03 [ Info: 753 1.6e+04 1.9e-05 0.0e+00 1.2e-02 [ Info: 754 1.6e+04 1.2e-05 0.0e+00 1.2e-02 [ Info: 755 1.6e+04 2.1e-04 2.1e-09 4.6e-03 [ Info: 756 1.6e+04 1.8e-04 0.0e+00 9.0e-04 [ Info: 757 1.6e+04 1.3e-05 0.0e+00 8.9e-04 [ Info: 758 1.6e+04 1.3e-04 1.6e-09 8.7e-04 [ Info: 759 1.6e+04 1.6e-04 3.7e-09 8.4e-04 [ Info: 760 1.6e+04 2.2e-04 6.4e-09 8.2e-04 [ Info: 761 1.6e+04 2.7e-04 9.8e-09 8.0e-04 [ Info: 762 1.6e+04 7.7e-04 0.0e+00 8.1e-04 [ Info: 763 1.6e+04 3.9e-04 0.0e+00 7.9e-04 [ Info: 764 1.6e+04 3.1e-04 0.0e+00 7.7e-04 [ Info: 765 1.6e+04 7.0e-04 0.0e+00 7.2e-04 [ Info: 766 1.6e+04 1.6e-04 2.0e-09 7.0e-04 [ Info: 767 1.6e+04 5.6e-04 9.1e-09 7.2e-04 [ Info: 768 1.6e+04 8.6e-04 0.0e+00 7.6e-04 [ Info: 769 1.6e+04 1.4e-04 0.0e+00 1.7e-03 [ Info: 770 1.6e+04 3.0e-06 0.0e+00 2.5e-03 [ Info: 771 1.6e+04 3.5e-07 0.0e+00 2.0e-03 [ Info: 772 1.6e+04 1.5e-06 0.0e+00 6.5e-04 [ Info: 773 1.6e+04 3.0e-06 0.0e+00 5.7e-04 [ Info: 774 1.6e+04 2.9e-04 0.0e+00 5.5e-04 [ Info: 775 1.6e+04 2.5e-04 0.0e+00 5.4e-04 [ Info: 776 1.6e+04 2.2e-04 0.0e+00 7.7e-04 [ Info: 777 1.6e+04 7.6e-04 0.0e+00 4.7e-04 [ Info: 778 1.6e+04 4.6e-04 5.9e-09 4.6e-04 [ Info: 779 1.6e+04 4.6e-04 0.0e+00 4.7e-04 [ Info: 780 1.6e+04 2.1e-04 0.0e+00 4.5e-04 [ Info: 781 1.6e+04 2.2e-04 0.0e+00 4.4e-04 [ Info: 782 1.6e+04 2.0e-04 0.0e+00 6.3e-04 [ Info: 783 1.6e+04 6.3e-04 0.0e+00 3.9e-04 [ Info: 784 1.6e+04 1.4e-04 0.0e+00 3.8e-04 [ Info: 785 1.6e+04 7.0e-04 7.2e-09 5.6e-04 [ Info: 786 1.6e+04 1.3e-03 0.0e+00 7.2e-04 [ Info: 787 1.6e+04 7.2e-04 0.0e+00 1.0e-03 [ Info: 788 1.6e+04 2.5e-06 0.0e+00 1.5e-03 [ Info: 789 1.6e+04 9.6e-07 0.0e+00 1.3e-03 [ Info: 790 1.6e+04 1.7e-06 0.0e+00 4.3e-04 [ Info: 791 1.6e+04 3.4e-04 4.3e-09 3.4e-04 [ Info: 792 1.6e+04 3.7e-04 9.1e-09 7.1e-04 [ Info: 793 1.6e+04 8.1e-05 8.0e-09 6.3e-04 [ Info: 794 1.6e+04 6.3e-04 0.0e+00 3.1e-04 [ Info: 795 1.6e+04 4.0e-07 0.0e+00 2.7e-04 [ Info: 796 1.6e+04 3.2e-04 4.0e-09 3.2e-04 [ Info: 797 1.6e+04 5.3e-04 0.0e+00 2.7e-04 [ Info: 798 1.6e+04 9.8e-04 9.7e-09 7.6e-04 [ Info: 799 1.6e+04 1.2e-03 0.0e+00 3.9e-04 [ Info: 800 1.6e+04 7.2e-04 4.1e-09 3.2e-04 [ Info: 801 1.6e+04 2.0e-04 1.5e-09 2.2e-04 [ Info: 802 1.6e+04 4.8e-05 9.2e-10 2.2e-04 [ Info: 803 1.6e+04 7.5e-05 0.0e+00 2.1e-04 [ Info: 804 1.6e+04 6.7e-05 0.0e+00 2.1e-04 [ Info: 805 1.6e+04 7.0e-05 0.0e+00 2.0e-04 [ Info: 806 1.6e+04 8.6e-05 0.0e+00 2.3e-04 [ Info: 807 1.6e+04 6.5e-05 0.0e+00 2.9e-04 [ Info: 808 1.6e+04 1.7e-04 0.0e+00 4.6e-04 [ Info: 809 1.6e+04 1.9e-04 0.0e+00 2.7e-04 [ Info: 810 1.6e+04 2.6e-04 0.0e+00 6.1e-04 [ Info: 811 1.6e+04 3.5e-07 0.0e+00 1.1e-03 [ Info: 812 1.6e+04 2.2e-07 0.0e+00 1.1e-03 [ Info: 813 1.6e+04 7.5e-07 0.0e+00 5.1e-04 [ Info: 814 1.6e+04 1.4e-04 1.7e-09 2.4e-04 [ Info: 815 1.6e+04 5.5e-04 0.0e+00 4.1e-04 [ Info: 816 1.6e+04 1.4e-04 0.0e+00 5.6e-04 [ Info: 817 1.6e+04 1.8e-04 0.0e+00 7.3e-04 [ Info: 818 1.6e+04 7.3e-04 0.0e+00 1.1e-04 [ Info: 819 1.6e+04 2.7e-04 3.5e-09 2.7e-04 [ Info: 820 1.6e+04 2.5e-04 6.7e-09 5.2e-04 [ Info: 821 1.6e+04 2.0e-04 9.1e-09 7.2e-04 [ Info: 822 1.6e+04 7.2e-04 0.0e+00 1.2e-04 [ Info: 823 1.6e+04 6.7e-05 0.0e+00 1.2e-04 [ Info: 824 1.6e+04 3.2e-04 0.0e+00 3.9e-04 [ Info: 825 1.6e+04 7.4e-05 0.0e+00 3.1e-04 [ Info: 826 1.6e+04 3.1e-04 0.0e+00 1.1e-04 [ Info: 827 1.6e+04 4.1e-04 5.2e-09 4.1e-04 [ Info: 828 1.6e+04 3.5e-05 4.7e-09 3.7e-04 [ Info: 829 1.6e+04 1.3e-04 6.4e-09 5.1e-04 [ Info: 830 1.6e+04 5.8e-05 7.2e-09 5.6e-04 [ Info: 831 1.6e+04 5.2e-06 7.1e-09 5.6e-04 [ Info: 832 1.6e+04 1.4e-04 8.9e-09 7.0e-04 [ Info: 833 1.6e+04 3.2e-04 4.8e-09 3.7e-04 [ Info: 834 1.6e+04 3.7e-04 0.0e+00 3.4e-04 [ Info: 835 1.6e+04 2.0e-07 0.0e+00 2.8e-04 [ Info: 836 1.6e+04 1.9e-04 2.4e-09 1.9e-04 [ Info: 837 1.6e+04 4.1e-06 2.4e-09 1.9e-04 [ Info: 838 1.6e+04 5.4e-05 1.7e-09 1.3e-04 [ Info: 839 1.6e+04 3.2e-06 1.7e-09 1.3e-04 [ Info: 840 1.6e+04 4.4e-06 1.7e-09 1.3e-04 [ Info: 841 1.6e+04 9.6e-06 1.5e-09 1.2e-04 [ Info: 842 1.6e+04 3.0e-06 1.5e-09 1.2e-04 [ Info: 843 1.6e+04 9.8e-06 1.4e-09 1.1e-04 [ Info: 844 1.6e+04 5.9e-06 1.3e-09 1.0e-04 [ Info: 845 1.6e+04 1.3e-05 1.1e-09 8.9e-05 Step 3: penalty=3.0, tol=0.0001, obj=16484.54898674185
The MMA87 optimizer iteratively updates the density field to minimize thermal compliance while satisfying the volume constraint. Ramping the penalty drives the final design toward a crisp 0/1 distribution.
Results
@show obj(res.minimizer)
@show constr(res.minimizer)obj(res.minimizer) = 16484.54898674185
constr(res.minimizer) = 1.1325817506069313e-91.1325817506069313e-9The final design should show the characteristic branching “tree” pattern that efficiently conducts heat from the top to the bottom boundary.
Visualization
fig = visualize(problem; static=true, topology=res.minimizer)
display_app(fig)undeformed mesh
load arrows
support arrows
The visualization reveals the tree-like structure — high-conductivity material forms branching paths similar to how structural members carry loads in compliance minimization. This analogy between heat conduction and structural mechanics is a fundamental insight in topology optimization.