<div>Thanks for your guide. </div><div>By monitoring the solutions during SNES, until it was diverged, I noticed that only at, same, one point, located on boundary, and at, same, one field,  repeatedly the solution went out of my specified constraint.</div><br><div>Amir</div><br><div class="gmail_quote_attribution">On Aug 19 2018, at 10:26 pm, Matthew Knepley <knepley@gmail.com> wrote:</div><blockquote><br><div><div><div><div><div>On Sun, Aug 19, 2018 at 10:19 AM Amir <<a href="mailto:mhbaghaei@mail.sjtu.edu.cn" title="mailto:mhbaghaei@mail.sjtu.edu.cn">mhbaghaei@mail.sjtu.edu.cn</a>> wrote:</div></div><blockquote><div>This is what I get after I modify the constraints, and using lu:</div><div>0 SNES Function norm 6.839975279451e+01</div><div>0 SNES VI Function norm 68.3997 Active lower constraints 0/0 upper constraints 0/0 Percent of total 0. Percent of bounded 0.</div><div>Linear solve converged due to CONVERGED_RTOL iterations 1</div><div>Line search: gnorm after quadratic fit 1.167382139621e+02</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 9.661035212791e+01 lambda=1.4620823686291085e-02</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 7.872462383177e+01 lambda=1.4620823686291085e-03</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 7.010735352812e+01 lambda=1.4620823686291087e-04</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.845164828370e+01 lambda=1.4620823686291088e-05</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.840024692820e+01 lambda=1.4620823686291089e-06</div><div>Line search: Cubically determined step, current gnorm 6.839974867692e+01 lambda=1.4620823686291089e-07</div></blockquote><div><br></div><div>So clearly the direction coming out of your Newton equation is not a descent direction when combined with the constraints. Since</div><div>it solves fine without the constraints, I suspect that your Jacobian is alright and your constraints are inconsistent with the solution.</div><div>Step 1, look at the solution you get without constraints. How does it violate your constraints?</div><div><br></div><div>   Matt</div><div> </div><blockquote><div>1 SNES Function norm 6.839974867692e+01</div><div>1 SNES VI Function norm 68.3997 Active lower constraints 0/0 upper constraints 0/0 Percent of total 0. Percent of bounded 0.</div><div>Linear solve converged due to CONVERGED_RTOL iterations 1</div><div>Line search: gnorm after quadratic fit 1.313376791907e+02</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 1.102614855398e+02 lambda=1.4115774027030809e-02</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 9.096517666290e+01 lambda=1.4115774027030810e-03</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 7.587134084601e+01 lambda=1.4115774027030810e-04</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.877273920434e+01 lambda=1.4115774027030810e-05</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.840408582565e+01 lambda=1.4115774027030811e-06</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.839977811504e+01 lambda=1.4115774027030813e-07</div><div>Line search: Cubically determined step, current gnorm 6.839974761445e+01 lambda=1.7645101306134412e-08</div><div>2 SNES Function norm 6.839974761445e+01</div><div>2 SNES VI Function norm 68.3997 Active lower constraints 0/0 upper constraints 0/0 Percent of total 0. Percent of bounded 0.</div><div>Linear solve converged due to CONVERGED_RTOL iterations 2</div><div>Line search: gnorm after quadratic fit 1.053816835919e+02</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 1.064875244204e+02 lambda=2.9620424834498244e-02</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 8.928160190384e+01 lambda=2.9620424834498247e-03</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 7.462621000022e+01 lambda=2.9620424834498249e-04</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.862607236334e+01 lambda=2.9620424834498250e-05</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.840219410268e+01 lambda=2.9620424834498251e-06</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.839975697009e+01 lambda=2.9620424834498251e-07</div><div>Line search: Cubically determined step, current gnorm 6.839974486019e+01 lambda=1.0073837635308560e-07</div><div>3 SNES Function norm 6.839974486019e+01</div><div>3 SNES VI Function norm 68.3997 Active lower constraints 0/0 upper constraints 0/0 Percent of total 0. Percent of bounded 0.</div><div>Linear solve converged due to CONVERGED_RTOL iterations 2</div><div>Line search: gnorm after quadratic fit 1.235990792724e+02</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 1.013971718389e+02 lambda=2.3488863243422971e-02</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 7.781731557824e+01 lambda=2.3488863243422975e-03</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.872972483427e+01 lambda=2.3488863243422976e-04</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.840207967049e+01 lambda=2.3488863243422978e-05</div><div>Line search: Cubically determined step, current gnorm 6.839957781288e+01 lambda=4.7431867318372446e-06</div><div>4 SNES Function norm 6.839957781288e+01</div><div>4 SNES VI Function norm 68.4009 Active lower constraints 0/0 upper constraints 0/0 Percent of total 0. Percent of bounded 0.</div><div>Linear solve converged due to CONVERGED_RTOL iterations 2</div><div>Line search: gnorm after quadratic fit 1.076412782767e+02</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 8.067259042361e+01 lambda=2.9654955573060381e-02</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.991099677071e+01 lambda=2.9654955573060384e-03</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.841319422751e+01 lambda=2.9654955573060385e-04</div><div>Line search: Cubically determined step, current gnorm 6.839661391957e+01 lambda=8.6072251116571315e-05</div><div>5 SNES Function norm 6.839661391957e+01</div><div>5 SNES VI Function norm 68.4308 Active lower constraints 0/0 upper constraints 0/0 Percent of total 0. Percent of bounded 0.</div><div>Linear solve converged due to CONVERGED_RTOL iterations 2</div><div>Line search: gnorm after quadratic fit 9.167521781305e+01</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 8.101349837507e+01 lambda=6.4300240519536755e-02</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 7.041475585659e+01 lambda=6.6547102741568651e-03</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.844705599568e+01 lambda=6.6547102741568653e-04</div><div>Line search: Cubically determined step, current gnorm 6.839121015135e+01 lambda=1.4965127147397144e-04</div><div>6 SNES Function norm 6.839121015135e+01</div><div>6 SNES VI Function norm 68.5745 Active lower constraints 0/0 upper constraints 0/0 Percent of total 0. Percent of bounded 0.</div><div>Linear solve converged due to CONVERGED_RTOL iterations 3</div><div>Line search: gnorm after quadratic fit 9.587731188484e+01</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 8.197287845440e+01 lambda=5.2775823101139679e-02</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.940305546680e+01 lambda=5.2775823101139686e-03</div><div>Line search: Cubically determined step, current gnorm 6.837177518619e+01 lambda=6.3969948101878912e-04</div><div>7 SNES Function norm 6.837177518619e+01</div><div>7 SNES VI Function norm 71.0652 Active lower constraints 0/0 upper constraints 0/0 Percent of total 0. Percent of bounded 0.</div><div>Linear solve converged due to CONVERGED_RTOL iterations 3</div><div>Line search: gnorm after quadratic fit 9.548865181502e+01</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 7.876974019161e+01 lambda=5.0026957605696178e-02</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.891985531663e+01 lambda=5.1260581886742491e-03</div><div>Line search: Cubically determined step, current gnorm 6.832012833110e+01 lambda=9.3300874565775075e-04</div><div>8 SNES Function norm 6.832012833110e+01</div><div>8 SNES VI Function norm 79.8569 Active lower constraints 0/0 upper constraints 0/0 Percent of total 0. Percent of bounded 0.</div><div>Linear solve converged due to CONVERGED_RTOL iterations 4</div><div>Line search: gnorm after quadratic fit 9.541909192352e+01</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 7.493349067288e+01 lambda=5.1388541796399699e-02</div><div>Line search: Cubically determined step, current gnorm 6.781311229195e+01 lambda=7.6704381452756282e-03</div><div>9 SNES Function norm 6.781311229195e+01</div><div>9 SNES VI Function norm 206.467 Active lower constraints 0/0 upper constraints 0/0 Percent of total 0. Percent of bounded 0.</div><div>Linear solve converged due to CONVERGED_RTOL iterations 6</div><div>Line search: gnorm after quadratic fit 9.394254710662e+01</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 7.306526091645e+01 lambda=5.3539445292607950e-02</div><div>Line search: Cubically determined step, current gnorm 6.721835250392e+01 lambda=9.5561464381222188e-03</div><div>10 SNES Function norm 6.721835250392e+01</div><div>10 SNES VI Function norm 368.277 Active lower constraints 0/0 upper constraints 0/0 Percent of total 0. Percent of bounded 0.</div><div>Linear solve converged due to CONVERGED_RTOL iterations 7</div><div>Line search: gnorm after quadratic fit 9.387716261564e+01</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 7.070868333837e+01 lambda=4.9443994240769853e-02</div><div>Line search: Cubically determined step, current gnorm 6.660967991379e+01 lambda=1.0810468882260445e-02</div><div>11 SNES Function norm 6.660967991379e+01</div><div>11 SNES VI Function norm 567.089 Active lower constraints 0/0 upper constraints 0/0 Percent of total 0. Percent of bounded 0.</div><div>Linear solve converged due to CONVERGED_RTOL iterations 9</div><div>Line search: gnorm after quadratic fit 9.346306901087e+01</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 7.048649174104e+01 lambda=4.9192073733008407e-02</div><div>Line search: Cubically determined step, current gnorm 6.608645472517e+01 lambda=1.0044112192046832e-02</div><div>12 SNES Function norm 6.608645472517e+01</div><div>12 SNES VI Function norm 734.752 Active lower constraints 0/0 upper constraints 0/0 Percent of total 0. Percent of bounded 0.</div><div>Linear solve converged due to CONVERGED_RTOL iterations 8</div><div>Line search: gnorm after quadratic fit 9.452595067062e+01</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 7.208837086000e+01 lambda=4.5128845776024190e-02</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.800731174084e+01 lambda=6.5132735322282484e-03</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.634685717588e+01 lambda=6.5132735322282484e-04</div><div>Line search: Cubically determined step, current gnorm 6.608637269301e+01 lambda=6.5132735322282487e-05</div><div>13 SNES Function norm 6.608637269301e+01</div><div>13 SNES VI Function norm 735.814 Active lower constraints 0/0 upper constraints 0/0 Percent of total 0. Percent of bounded 0.</div><div>Linear solve converged due to CONVERGED_RTOL iterations 8</div><div>Line search: gnorm after quadratic fit 9.318693996551e+01</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 7.139943393009e+01 lambda=4.8118075850982413e-02</div><div>Line search: Cubically determined step, current gnorm 6.568600823304e+01 lambda=7.8756632431255542e-03</div><div>14 SNES Function norm 6.568600823304e+01</div><div>14 SNES VI Function norm 881.479 Active lower constraints 0/0 upper constraints 0/0 Percent of total 0. Percent of bounded 0.</div><div>Linear solve converged due to CONVERGED_RTOL iterations 9</div><div>Line search: gnorm after quadratic fit 9.242424276014e+01</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 7.226307852468e+01 lambda=4.8021995541678457e-02</div><div>Line search: Cubically determined step, current gnorm 6.540144178550e+01 lambda=6.6931701304969058e-03</div><div>15 SNES Function norm 6.540144178550e+01</div><div>15 SNES VI Function norm 1006.28 Active lower constraints 0/0 upper constraints 0/0 Percent of total 0. Percent of bounded 0.</div><div>Linear solve converged due to CONVERGED_RTOL iterations 10</div><div>Line search: gnorm after quadratic fit 9.154992244384e+01</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 7.122067564574e+01 lambda=5.0277415210907568e-02</div><div>Line search: Cubically determined step, current gnorm 6.506969399500e+01 lambda=7.8660529387648906e-03</div><div>16 SNES Function norm 6.506969399500e+01</div><div>16 SNES VI Function norm 1131.4 Active lower constraints 0/0 upper constraints 0/0 Percent of total 0. Percent of bounded 0.</div><div>Linear solve converged due to CONVERGED_RTOL iterations 10</div><div>Line search: gnorm after quadratic fit 9.263744099217e+01</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 7.045825559590e+01 lambda=4.3976587783298791e-02</div><div>Line search: Cubically determined step, current gnorm 6.477740872258e+01 lambda=6.6653314620139965e-03</div><div>17 SNES Function norm 6.477740872258e+01</div><div>17 SNES VI Function norm 1252.29 Active lower constraints 0/0 upper constraints 0/0 Percent of total 0. Percent of bounded 0.</div><div>Linear solve converged due to CONVERGED_RTOL iterations 10</div><div>Line search: gnorm after quadratic fit 9.034042193429e+01</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 7.097176730586e+01 lambda=4.9082849797732621e-02</div><div>Line search: Cubically determined step, current gnorm 6.448770943479e+01 lambda=7.1730414269226593e-03</div><div>18 SNES Function norm 6.448770943479e+01</div><div>18 SNES VI Function norm 1371.72 Active lower constraints 0/0 upper constraints 0/0 Percent of total 0. Percent of bounded 0.</div><div>Linear solve converged due to CONVERGED_RTOL iterations 10</div><div>Line search: gnorm after quadratic fit 9.028256513212e+01</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.990381886066e+01 lambda=4.4193294185915478e-02</div><div>Line search: Cubically determined step, current gnorm 6.422123220407e+01 lambda=6.6250183714544201e-03</div><div>19 SNES Function norm 6.422123220407e+01</div><div>19 SNES VI Function norm 1490.87 Active lower constraints 0/0 upper constraints 0/0 Percent of total 0. Percent of bounded 0.</div><div>Linear solve converged due to CONVERGED_RTOL iterations 10</div><div>Line search: gnorm after quadratic fit 8.908334154738e+01</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 7.060161630260e+01 lambda=4.6452357120490895e-02</div><div>Line search: Cubically determined step, current gnorm 6.411445689017e+01 lambda=6.3417192629449783e-03</div><div>20 SNES Function norm 6.411445689017e+01</div><div>20 SNES VI Function norm 1597.29 Active lower constraints 0/0 upper constraints 0/0 Percent of total 0. Percent of bounded 0.</div><div>Linear solve converged due to CONVERGED_RTOL iterations 11</div><div>Line search: gnorm after quadratic fit 8.460267552084e+01</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 7.051195181220e+01 lambda=5.8933530173091400e-02</div><div>Line search: Cubically determined step, current gnorm 6.379276529842e+01 lambda=9.2451193719573169e-03</div><div>21 SNES Function norm 6.379276529842e+01</div><div>21 SNES VI Function norm 1729.35 Active lower constraints 0/0 upper constraints 0/0 Percent of total 0. Percent of bounded 0.</div><div>Linear solve converged due to CONVERGED_RTOL iterations 10</div><div>Line search: gnorm after quadratic fit 8.847710713735e+01</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 7.007434740347e+01 lambda=4.8620669909360585e-02</div><div>Line search: Cubically determined step, current gnorm 6.358026444281e+01 lambda=6.8922394760267177e-03</div><div>22 SNES Function norm 6.358026444281e+01</div><div>22 SNES VI Function norm 1849.77 Active lower constraints 0/0 upper constraints 0/0 Percent of total 0. Percent of bounded 0.</div><div>Linear solve converged due to CONVERGED_RTOL iterations 10</div><div>Line search: gnorm after quadratic fit 9.435581241699e+01</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 7.335475468652e+01 lambda=3.2814238215953598e-02</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.652843583618e+01 lambda=3.2814238215953598e-03</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.529985391051e+01 lambda=3.2814238215953602e-04</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.366726270581e+01 lambda=3.2814238215953603e-05</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.358115196460e+01 lambda=3.2814238215953603e-06</div><div>Line search: Cubically determined step, current gnorm 6.358025861312e+01 lambda=3.2814238215953605e-07</div><div>23 SNES Function norm 6.358025861312e+01</div><div>23 SNES VI Function norm 1849.78 Active lower constraints 0/0 upper constraints 0/0 Percent of total 0. Percent of bounded 0.</div><div>Linear solve converged due to CONVERGED_RTOL iterations 10</div><div>Line search: gnorm after quadratic fit 8.812900895448e+01</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.822426303212e+01 lambda=4.8421685266110105e-02</div><div>Line search: Cubically determined step, current gnorm 6.335842666537e+01 lambda=8.4588034499471741e-03</div><div>24 SNES Function norm 6.335842666537e+01</div><div>24 SNES VI Function norm 1993.6 Active lower constraints 0/0 upper constraints 0/0 Percent of total 0. Percent of bounded 0.</div><div>Linear solve converged due to CONVERGED_RTOL iterations 11</div><div>Line search: gnorm after quadratic fit 9.094159193064e+01</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.742173830047e+01 lambda=4.3243453252109078e-02</div><div>Line search: Cubically determined step, current gnorm 6.316236273224e+01 lambda=7.7494584282214512e-03</div><div>25 SNES Function norm 6.316236273224e+01</div><div>25 SNES VI Function norm 2126.62 Active lower constraints 0/0 upper constraints 0/0 Percent of total 0. Percent of bounded 0.</div><div>Linear solve converged due to CONVERGED_RTOL iterations 11</div><div>Line search: gnorm after quadratic fit 9.160540289460e+01</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.722886652974e+01 lambda=4.1792497999647044e-02</div><div>Line search: Cubically determined step, current gnorm 6.300020419956e+01 lambda=7.3335939016311632e-03</div><div>26 SNES Function norm 6.300020419956e+01</div><div>26 SNES VI Function norm 2252.12 Active lower constraints 0/0 upper constraints 0/0 Percent of total 0. Percent of bounded 0.</div><div>Linear solve converged due to CONVERGED_RTOL iterations 10</div><div>Line search: gnorm after quadratic fit 8.890109786397e+01</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.686097246249e+01 lambda=4.7041389395739201e-02</div><div>Line search: Cubically determined step, current gnorm 6.283084254493e+01 lambda=9.0792221408297433e-03</div><div>27 SNES Function norm 6.283084254493e+01</div><div>27 SNES VI Function norm 2404.84 Active lower constraints 0/0 upper constraints 0/0 Percent of total 0. Percent of bounded 0.</div><div>Linear solve converged due to CONVERGED_RTOL iterations 11</div><div>Line search: gnorm after quadratic fit 8.998798828268e+01</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.576469464321e+01 lambda=4.2887406731046039e-02</div><div>Line search: Cubically determined step, current gnorm 6.266968273817e+01 lambda=9.2948647727201898e-03</div><div>28 SNES Function norm 6.266968273817e+01</div><div>28 SNES VI Function norm 2558.11 Active lower constraints 0/0 upper constraints 0/0 Percent of total 0. Percent of bounded 0.</div><div>Linear solve converged due to CONVERGED_RTOL iterations 11</div><div>Line search: gnorm after quadratic fit 9.199146621649e+01</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.564298288247e+01 lambda=4.0193768728643060e-02</div><div>Line search: Cubically determined step, current gnorm 6.252861090423e+01 lambda=8.3932894622200105e-03</div><div>29 SNES Function norm 6.252861090423e+01</div><div>29 SNES VI Function norm 2698.56 Active lower constraints 0/0 upper constraints 0/0 Percent of total 0. Percent of bounded 0.</div><div>Linear solve converged due to CONVERGED_RTOL iterations 11</div><div>Line search: gnorm after quadratic fit 9.382191898563e+01</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.772162800128e+01 lambda=3.5914901407691102e-02</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.485302824556e+01 lambda=4.7901702628153379e-03</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.285215446174e+01 lambda=4.7901702628153379e-04</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.253029194340e+01 lambda=4.7901702628153380e-05</div><div>Line search: Cubically determined step, current gnorm 6.252812967035e+01 lambda=1.5127710264417741e-05</div><div>30 SNES Function norm 6.252812967035e+01</div><div>30 SNES VI Function norm 2698.83 Active lower constraints 0/0 upper constraints 0/0 Percent of total 0. Percent of bounded 0.</div><div>Linear solve converged due to CONVERGED_RTOL iterations 11</div><div>Line search: gnorm after quadratic fit 9.118652411357e+01</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.495202040381e+01 lambda=4.0596953253201326e-02</div><div>Line search: Cubically determined step, current gnorm 6.239041273754e+01 lambda=9.5128837682561859e-03</div><div>31 SNES Function norm 6.239041273754e+01</div><div>31 SNES VI Function norm 2854.38 Active lower constraints 0/0 upper constraints 0/0 Percent of total 0. Percent of bounded 0.</div><div>Linear solve converged due to CONVERGED_RTOL iterations 11</div><div>Line search: gnorm after quadratic fit 9.077825361440e+01</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.467284344645e+01 lambda=4.0583557291813177e-02</div><div>Line search: Cubically determined step, current gnorm 6.225884047363e+01 lambda=9.7846838097028124e-03</div><div>32 SNES Function norm 6.225884047363e+01</div><div>32 SNES VI Function norm 3012.55 Active lower constraints 0/0 upper constraints 0/0 Percent of total 0. Percent of bounded 0.</div><div>Linear solve converged due to CONVERGED_RTOL iterations 11</div><div>Line search: gnorm after quadratic fit 8.920999190500e+01</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.466103609666e+01 lambda=4.2498592612658738e-02</div><div>Line search: Cubically determined step, current gnorm 6.213582712328e+01 lambda=1.0171220001350099e-02</div><div>33 SNES Function norm 6.213582712328e+01</div><div>33 SNES VI Function norm 3175.23 Active lower constraints 0/0 upper constraints 0/0 Percent of total 0. Percent of bounded 0.</div><div>Linear solve converged due to CONVERGED_RTOL iterations 11</div><div>Line search: gnorm after quadratic fit 9.011924620812e+01</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.495110751676e+01 lambda=4.2453491229013124e-02</div><div>Line search: Cubically determined step, current gnorm 6.201636649414e+01 lambda=9.3685443832871688e-03</div><div>34 SNES Function norm 6.201636649414e+01</div><div>34 SNES VI Function norm 3322.22 Active lower constraints 0/0 upper constraints 0/0 Percent of total 0. Percent of bounded 0.</div><div>Linear solve converged due to CONVERGED_RTOL iterations 11</div><div>Line search: gnorm after quadratic fit 9.006024690051e+01</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.414041989575e+01 lambda=4.0940226638739433e-02</div><div>Line search: Cubically determined step, current gnorm 6.190218472271e+01 lambda=1.0273994032172384e-02</div><div>35 SNES Function norm 6.190218472271e+01</div><div>35 SNES VI Function norm 3484.25 Active lower constraints 0/0 upper constraints 0/0 Percent of total 0. Percent of bounded 0.</div><div>Linear solve converged due to CONVERGED_RTOL iterations 11</div><div>Line search: gnorm after quadratic fit 8.898991234613e+01</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.426189506521e+01 lambda=4.3426020819675164e-02</div><div>Line search: Cubically determined step, current gnorm 6.178678684779e+01 lambda=1.0559984512509757e-02</div><div>36 SNES Function norm 6.178678684779e+01</div><div>36 SNES VI Function norm 3648.36 Active lower constraints 0/0 upper constraints 0/0 Percent of total 0. Percent of bounded 0.</div><div>Linear solve converged due to CONVERGED_RTOL iterations 11</div><div>Line search: gnorm after quadratic fit 8.648025040907e+01</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.404650680722e+01 lambda=4.6749921483990771e-02</div><div>Line search: Cubically determined step, current gnorm 6.165099494681e+01 lambda=1.1866483746353381e-02</div><div>37 SNES Function norm 6.165099494681e+01</div><div>37 SNES VI Function norm 3822.88 Active lower constraints 0/0 upper constraints 0/0 Percent of total 0. Percent of bounded 0.</div><div>Linear solve converged due to CONVERGED_RTOL iterations 11</div><div>Line search: gnorm after quadratic fit 8.858076523292e+01</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.520416756669e+01 lambda=4.2063577748779259e-02</div><div>Line search: Cubically determined step, current gnorm 6.152598711257e+01 lambda=7.9475995139781722e-03</div><div>38 SNES Function norm 6.152598711257e+01</div><div>38 SNES VI Function norm 3944.22 Active lower constraints 0/0 upper constraints 0/0 Percent of total 0. Percent of bounded 0.</div><div>Linear solve converged due to CONVERGED_RTOL iterations 11</div><div>Line search: gnorm after quadratic fit 8.742192022113e+01</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.343055909293e+01 lambda=3.8810591937350140e-02</div><div>Line search: Cubically determined step, current gnorm 6.135131913134e+01 lambda=9.8750504088198892e-03</div><div>39 SNES Function norm 6.135131913134e+01</div><div>39 SNES VI Function norm 4066.15 Active lower constraints 0/1 upper constraints 0/0 Percent of total 0. Percent of bounded 0.</div><div>Linear solve converged due to CONVERGED_RTOL iterations 11</div><div>Line search: gnorm after quadratic fit 9.286112632409e+01</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 7.061265007302e+01 lambda=2.9694426041134040e-02</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.436749315348e+01 lambda=2.9694426041134044e-03</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.299789903123e+01 lambda=2.9694426041134044e-04</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.141298285103e+01 lambda=2.9694426041134044e-05</div><div>Line search: Cubically determined step, current gnorm 6.135130530880e+01 lambda=2.9694426041134044e-06</div><div>40 SNES Function norm 6.135130530880e+01</div><div>40 SNES VI Function norm 4066.23 Active lower constraints 0/1 upper constraints 0/0 Percent of total 0. Percent of bounded 0.</div><div>Linear solve converged due to CONVERGED_RTOL iterations 11</div><div>Line search: gnorm after quadratic fit 8.660062754407e+01</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.674944119178e+01 lambda=4.3170190544802138e-02</div><div>Line search: Cubically determined step, current gnorm 6.120973541057e+01 lambda=6.1770273873522300e-03</div><div>41 SNES Function norm 6.120973541057e+01</div><div>41 SNES VI Function norm 4152.83 Active lower constraints 0/1 upper constraints 0/0 Percent of total 0. Percent of bounded 0.</div><div>Linear solve converged due to CONVERGED_RTOL iterations 11</div><div>Line search: gnorm after quadratic fit 8.564645091492e+01</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.676153862982e+01 lambda=4.4220462778047379e-02</div><div>Line search: Cubically determined step, current gnorm 6.104393345352e+01 lambda=6.3589031499342100e-03</div><div>42 SNES Function norm 6.104393345352e+01</div><div>42 SNES VI Function norm 4233.82 Active lower constraints 0/1 upper constraints 0/0 Percent of total 0. Percent of bounded 0.</div><div>Linear solve converged due to CONVERGED_RTOL iterations 11</div><div>Line search: gnorm after quadratic fit 8.368658054045e+01</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.715427252993e+01 lambda=4.7003756543045970e-02</div><div>Line search: Cubically determined step, current gnorm 6.101753259623e+01 lambda=6.3787492486959132e-03</div><div>43 SNES Function norm 6.101753259623e+01</div><div>43 SNES VI Function norm 4314.49 Active lower constraints 0/1 upper constraints 0/0 Percent of total 0. Percent of bounded 0.</div><div>Linear solve converged due to CONVERGED_RTOL iterations 11</div><div>Line search: gnorm after quadratic fit 8.442307716417e+01</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.496496980051e+01 lambda=4.5726539795380154e-02</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.142635208964e+01 lambda=8.4157732414286580e-03</div><div>Line search: Cubically determined step, current gnorm 6.095430498328e+01 lambda=2.1294466610363495e-03</div><div>44 SNES Function norm 6.095430498328e+01</div><div>44 SNES VI Function norm 4341.47 Active lower constraints 0/1 upper constraints 0/0 Percent of total 0. Percent of bounded 0.</div><div>Linear solve converged due to CONVERGED_RTOL iterations 12</div><div>Line search: gnorm after quadratic fit 8.439891204437e+01</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.240254752159e+01 lambda=4.3289599290802970e-02</div><div>Line search: Cubically determined step, current gnorm 6.081365387079e+01 lambda=1.2916446928212236e-02</div><div>45 SNES Function norm 6.081365387079e+01</div><div>45 SNES VI Function norm 4500.53 Active lower constraints 0/1 upper constraints 0/0 Percent of total 0. Percent of bounded 0.</div><div>Linear solve converged due to CONVERGED_RTOL iterations 11</div><div>Line search: gnorm after quadratic fit 9.664056729635e+01</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.944851102419e+01 lambda=1.8635858087631831e-02</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.394837998027e+01 lambda=1.8635858087631832e-03</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.278442438314e+01 lambda=1.8635858087631834e-04</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.093400706607e+01 lambda=1.8635858087631834e-05</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.081489139942e+01 lambda=1.8635858087631835e-06</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.081365630503e+01 lambda=1.8635858087631836e-07</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.081365482998e+01 lambda=7.3438624052552081e-08</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.081365418254e+01 lambda=2.3870801624167493e-08</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.081365397577e+01 lambda=8.0330117961387493e-09</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.081365390593e+01 lambda=2.6874115312641255e-09</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.081365388261e+01 lambda=8.9982494953058537e-10</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.081365387467e+01 lambda=3.0083351315759278e-10</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.081365387219e+01 lambda=1.0110298276678126e-10</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.081365387129e+01 lambda=3.3292703628714462e-11</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.081365387093e+01 lambda=1.0739379834615402e-11</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.081365387083e+01 lambda=3.6237724732408161e-12</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.081365387082e+01 lambda=1.2456833192766385e-12</div><div>Line search: Cubically determined step, current gnorm 6.081365387079e+01 lambda=3.5247146389212346e-13</div><div>46 SNES Function norm 6.081365387079e+01</div><div>46 SNES VI Function norm 4500.53 Active lower constraints 0/1 upper constraints 0/0 Percent of total 0. Percent of bounded 0.</div><div>Linear solve converged due to CONVERGED_RTOL iterations 11</div><div>Line search: gnorm after quadratic fit 1.433547268242e+03</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 7.106740636849e+01 lambda=2.4859604812085913e-02</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.525497375887e+01 lambda=7.2498301889276232e-03</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.367730849091e+01 lambda=7.2498301889276236e-04</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.174333248068e+01 lambda=7.2498301889276242e-05</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.083400168850e+01 lambda=7.2498301889276242e-06</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.081378891362e+01 lambda=7.2498301889276242e-07</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.081365527548e+01 lambda=1.0753239548189211e-07</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.081365445012e+01 lambda=4.4350670800939511e-08</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.081365407923e+01 lambda=1.5953160232199837e-08</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.081365394762e+01 lambda=5.8822684037182187e-09</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.081365389903e+01 lambda=2.1599843490971948e-09</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.081365388111e+01 lambda=7.9357831777119289e-10</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.081365387457e+01 lambda=2.9181178721013368e-10</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.081365387223e+01 lambda=1.0733281860745436e-10</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.081365387132e+01 lambda=3.9236995040615546e-11</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.081365387106e+01 lambda=1.4351792294805992e-11</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.081365387081e+01 lambda=4.8683006404131606e-12</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.081365387087e+01 lambda=2.0863028416245319e-12</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.081365387080e+01 lambda=5.3058618006527772e-13</div><div>Line search: unable to find good step length! After 19 tries</div><div>Line search: fnorm=6.0813653870787959e+01, gnorm=6.0813653870799314e+01, ynorm=7.0504021074247237e+02, minlambda=9.9999999999999998e-13, lambda=5.3058618006527772e-13, initial slope=-4.7232673282579080e+03</div><div>47 SNES Function norm 6.081365387079e+01</div><div>47 SNES VI Function norm 4500.53 Active lower constraints 0/1 upper constraints 0/0 Percent of total 0. Percent of bounded 0.</div><div>Linear solve converged due to CONVERGED_RTOL iterations 11</div><div>Line search: gnorm after quadratic fit 1.280890413731e+02</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 7.053050135324e+01 lambda=2.2733898438250671e-02</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.399979736493e+01 lambda=2.2733898438250673e-03</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.297299206174e+01 lambda=2.2733898438250674e-04</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.098427084652e+01 lambda=2.2733898438250677e-05</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.081555475063e+01 lambda=2.2733898438250679e-06</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.081365745105e+01 lambda=2.2733898438250681e-07</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.081365507363e+01 lambda=9.2077672155667702e-08</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.081365429323e+01 lambda=3.2329991301020441e-08</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.081365402243e+01 lambda=1.1595563414967544e-08</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.081365392495e+01 lambda=4.1432803927719340e-09</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.081365389024e+01 lambda=1.4815603373453140e-09</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.081365387781e+01 lambda=5.2922825349649990e-10</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.081365387324e+01 lambda=1.8866396505817014e-10</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.081365387173e+01 lambda=6.7568116386381729e-11</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.081365387107e+01 lambda=2.3800109593117919e-11</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.081365387094e+01 lambda=8.6880387507986266e-12</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.081365387083e+01 lambda=2.8858967256610185e-12</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.081365387082e+01 lambda=1.0348935437400166e-12</div><div>Line search: Cubically determined step, current gnorm 6.081365387077e+01 lambda=2.7070538024707468e-13</div><div>48 SNES Function norm 6.081365387077e+01</div><div>48 SNES VI Function norm 4500.53 Active lower constraints 0/1 upper constraints 0/0 Percent of total 0. Percent of bounded 0.</div><div>Linear solve converged due to CONVERGED_RTOL iterations 11</div><div>Line search: gnorm after quadratic fit 1.170943716019e+02</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 7.038597782185e+01 lambda=2.2207476459146142e-02</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.399318571001e+01 lambda=2.2207476459146142e-03</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.295164029113e+01 lambda=2.2207476459146143e-04</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.097743654776e+01 lambda=2.2207476459146143e-05</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.081546079136e+01 lambda=2.2207476459146144e-06</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.081365719872e+01 lambda=2.2207476459146145e-07</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.081365503919e+01 lambda=9.0117449908413481e-08</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.081365427749e+01 lambda=3.1364523245356123e-08</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.081365401596e+01 lambda=1.1183735694685806e-08</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.081365392232e+01 lambda=3.9706025565052754e-09</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.081365388909e+01 lambda=1.4107330444508480e-09</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.081365387724e+01 lambda=5.0110371464381371e-10</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.081365387308e+01 lambda=1.7820909495605996e-10</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.081365387164e+01 lambda=6.3314737242538209e-11</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.081365387114e+01 lambda=2.2207910898013729e-11</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.081365387085e+01 lambda=7.4515329000007742e-12</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.081365387082e+01 lambda=2.7872933392432575e-12</div><div>Line search: Cubic step no good, shrinking lambda, current gnorm 6.081365387083e+01 lambda=8.8269934321345058e-13</div><div>Line search: unable to find good step length! After 18 tries</div><div>Line search: fnorm=6.0813653870766274e+01, gnorm=6.0813653870831025e+01, ynorm=7.0578565675503091e+02, minlambda=9.9999999999999998e-13, lambda=8.8269934321345058e-13, initial slope=-3.8123379550451641e+03</div><div>Nonlinear solve did not converge due to DIVERGED_LINE_SEARCH iterations 48</div><br><div>Amir</div><br><div>On Aug 18 2018, at 11:42 pm, Smith, Barry F. <<a href="mailto:bsmith@mcs.anl.gov" title="mailto:bsmith@mcs.anl.gov">bsmith@mcs.anl.gov</a>> wrote:</div><blockquote><br><div><br><div>-snes_monitor -snes_linesearch_monitor -pc_type lu</div><br><div>produces what?</div><br><div>Barry</div><br><div>The line search failing means that the direction computed is not a descent direction which could be due to a bad Jacobian or inconsistent constraints.</div><br><br><blockquote><div>On Aug 17, 2018, at 11:56 PM, Amir <<a href="mailto:mhbaghaei@mail.sjtu.edu.cn" title="mailto:mhbaghaei@mail.sjtu.edu.cn">mhbaghaei@mail.sjtu.edu.cn</a>> wrote:</div><br><div>Hi</div><div>I need to apply constraint to primitive variables. I am troubling to set up SNESVI. After running the code, I feel SNES is not actually solving. After one iteration of SNES, I got DIVERGED_LINESEARCH, with NO changing in function norm.</div><div>I think the way I set VI maybe wrong. Sorry for my simple explanation.</div><div>Very thankful for your time.</div><div>Amir</div><div>SNESCreate(PETSC_COMM_WORLD,&snes);</div><div>SNESSetApplicationContext(snes,&user);</div><div>SNESSetDM(snes,<a href="http://user.dm" title="http://user.dm">user.dm</a>);</div><div>SNESSetFunction(snes,NULL,FormFunction,&user);</div><div>FormInitialSolution(user.X,&user);</div><div>SNESSetSolution(snes,user.X);</div><div>SNESMonitorSet(snes,MySNESMonitor,&user,NULL);</div><div>SNESSetType(snes,SNESVINEWTONSSLS);</div><div>SNESSetNormSchedule(snes, SNES_NORM_ALWAYS);</div><div>SNESVISetVariableBounds(snes,user.XLOWER,user.XUPPER);</div><div>SNESGetLineSearch(snes,&linesearch);</div><div>SNESLineSearchSetType(linesearch,SNESLINESEARCHBT);</div><div>SNESLineSearchSetVIFunctions(linesearch, NULL, SNESNormFunction);</div><div>SNESSetFromOptions(snes);</div><div>SNESSolve(snes,NULL,user.X);</div></blockquote></div></blockquote><div><img alt="Open Tracking" width="0" height="0" style="border:0;width:0;height:0" src="https://link.getmailspring.com/open/1534687923.local-d7e5039a-f013-v1.4.2-f587b7b7@getmailspring.com?recipient=cGV0c2MtdXNlcnNAbWNzLmFubC5nb3Y%3D"></div></blockquote></div><br><div><br></div><div>--</div><div><div><div><div><div><div>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.</div><div>-- Norbert Wiener</div></div><div><br></div><div><div><a href="http://www.caam.rice.edu/~mk51/" title="http://www.caam.rice.edu/~mk51/">https://www.cse.buffalo.edu/~knepley/</a></div></div></div></div></div></div></div></div></blockquote>