<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>
</head>
<body dir="ltr">
<div id="divtagdefaultwrapper" style="font-size:12pt;color:#000000;font-family:Calibri,Helvetica,sans-serif;" dir="ltr">
<p></p>
<pre>Hello,</pre>
<pre>I do topological optimization on a 3D problem</pre>
<pre>using the SIMP method whit p=4.</pre>
<pre><br></pre>
<pre><pre style="font-size: 16px;">I want a merely 0.001 (black) and 1 (white) design,</pre><pre style="font-size: 16px;">but it doesn't happen.
</pre><pre style="font-size: 16px;"><br></pre><pre style="font-size: 16px;">I report <u>output:</u></pre><pre style="font-size: 16px;">----------</pre></pre>
<div>IPOPT Note : No hessian given ==> LBFGS hessian approximation enabled</div>
<div>This is Ipopt version 3.10.2, running with linear solver mumps.</div>
<div><br>
</div>
<div>Number of nonzeros in equality constraint Jacobian...: 0</div>
<div>Number of nonzeros in inequality constraint Jacobian.: 13202</div>
<div>Number of nonzeros in Lagrangian Hessian.............: 0</div>
<div><br>
</div>
<div>Total number of variables............................: 13202</div>
<div> variables with only lower bounds: 0</div>
<div> variables with lower and upper bounds: 13202</div>
<div> variables with only upper bounds: 0</div>
<div>Total number of equality constraints.................: 0</div>
<div>Total number of inequality constraints...............: 1</div>
<div> inequality constraints with only lower bounds: 0</div>
<div> inequality constraints with lower and upper bounds: 1</div>
<div> inequality constraints with only upper bounds: 0</div>
<div><br>
</div>
<div> -- FESpace: Nb of Nodes 13202 Nb of DoF 39606</div>
<div> -- Solve : </div>
<div> min -0.064877 max 0.093047</div>
<div>iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls</div>
<div> 0 1.1217395e+04 0.00e+00 8.89e+01 0.0 0.00e+00 - 0.00e+00 0.00e+00 0</div>
<div> -- Solve : </div>
<div> min -0.0660702 max 0.0942643</div>
<div> 1 1.1708250e+04 0.00e+00 1.40e+02 -0.3 5.25e+01 - 1.00e+00 1.19e-01f 1</div>
<div> -- Solve : </div>
<div> min -0.0649783 max 0.0931706</div>
<br>
<p></p>
<p><br>
</p>
<p>... and so on</p>
<p><br>
</p>
<p></p>
<div> 299 8.5618728e+03 0.00e+00 1.63e+01 -0.9 3.08e-02 - 8.76e-01 4.78e-01f 2</div>
<div> -- Solve : </div>
<div> min -0.0491874 max 0.066941</div>
<div> -- Solve : </div>
<div> min -0.0491929 max 0.0669537</div>
<div>iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls</div>
<div> 300 8.5570186e+03 0.00e+00 2.06e+01 -0.9 3.96e-02 - 1.00e+00 4.13e-01f 2</div>
<div><br>
</div>
<div>Number of Iterations....: 300</div>
<div><br>
</div>
<div> (scaled) (unscaled)</div>
<div>Objective...............: 2.7318990891333351e+03 8.5570186106938800e+03</div>
<div>Dual infeasibility......: 2.0563058857838509e+01 6.4408849521282519e+01</div>
<div>Constraint violation....: 0.0000000000000000e+00 0.0000000000000000e+00</div>
<div>Complementarity.........: 1.4337684852579810e-01 4.4909358696960705e-01</div>
<div>Overall NLP error.......: 2.0563058857838509e+01 6.4408849521282519e+01</div>
<div><br>
</div>
<div><br>
</div>
<div>Number of objective function evaluations = 557</div>
<div>Number of objective gradient evaluations = 301</div>
<div>Number of equality constraint evaluations = 0</div>
<div>Number of inequality constraint evaluations = 558</div>
<div>Number of equality constraint Jacobian evaluations = 0</div>
<div>Number of inequality constraint Jacobian evaluations = 301</div>
<div>Number of Lagrangian Hessian evaluations = 0</div>
<div>Total CPU secs in IPOPT (w/o function evaluations) = 62.984</div>
<div>Total CPU secs in NLP function evaluations = 31406.348</div>
<div><br>
</div>
<div>EXIT: Maximum Number of Iterations Exceeded.</div>
<br>
<p></p>
<p>---------------------</p>
<p><br>
</p>
<p>I don't understand why (at the end) there are parts of for example 0.2 density; I want black or white, not grey. ( I have used the reference toll (10^-8) and I think that the problem is not the "<span style="font-family: Calibri, Helvetica, sans-serif, EmojiFont, "Apple Color Emoji", "Segoe UI Emoji", NotoColorEmoji, "Segoe UI Symbol", "Android Emoji", EmojiSymbols; font-size: 16px;">Maximum
Number of Iterations Exceeded" )</span></p>
<p><span style="font-family: Calibri, Helvetica, sans-serif, EmojiFont, "Apple Color Emoji", "Segoe UI Emoji", NotoColorEmoji, "Segoe UI Symbol", "Android Emoji", EmojiSymbols; font-size: 16px;"><br>
</span></p>
<p><span style="font-family: Calibri, Helvetica, sans-serif, EmojiFont, "Apple Color Emoji", "Segoe UI Emoji", NotoColorEmoji, "Segoe UI Symbol", "Android Emoji", EmojiSymbols; font-size: 16px;">Thank you</span></p>
<p><span style="font-family: Calibri, Helvetica, sans-serif, EmojiFont, "Apple Color Emoji", "Segoe UI Emoji", NotoColorEmoji, "Segoe UI Symbol", "Android Emoji", EmojiSymbols; font-size: 16px;"><br>
</span></p>
<p><span style="font-family: Calibri, Helvetica, sans-serif, EmojiFont, "Apple Color Emoji", "Segoe UI Emoji", NotoColorEmoji, "Segoe UI Symbol", "Android Emoji", EmojiSymbols; font-size: 16px;">Monica</span></p>
</div>
</body>
</html>