Can you share the code with [email protected]
This is very odd and unexpected and may indicate an error in PETSc. Barry > On Aug 19, 2018, at 10:12 PM, Amir <[email protected]> wrote: > > I just put the constraint at the PETSC_INFINITY, this is what I get, which is > the same as when I use NEWTONLS: > > 0 SNES Function norm 6.839972031896e+01 > 0 SNES VI Function norm 68.3997 Active lower constraints 0/0 upper > constraints 0/0 Percent of total 0. Percent of bounded 0. > Linear solve converged due to CONVERGED_RTOL iterations 2 > Line search: gnorm after quadratic fit 6.231294332518e+01 > Line search: Quadratically determined step, lambda=1.0000000000000001e-01 > 1 SNES Function norm 6.231294332518e+01 > 1 SNES VI Function norm 62.3129 Active lower constraints 0/0 upper > constraints 0/0 Percent of total 0. Percent of bounded 0. > Linear solve converged due to CONVERGED_RTOL iterations 2 > Line search: gnorm after quadratic fit 5.677271460653e+01 > Line search: Quadratically determined step, lambda=1.0000000000000001e-01 > 2 SNES Function norm 5.677271460653e+01 > 2 SNES VI Function norm 56.7727 Active lower constraints 0/0 upper > constraints 0/0 Percent of total 0. Percent of bounded 0. > Linear solve converged due to CONVERGED_RTOL iterations 2 > Line search: gnorm after quadratic fit 4.945867427788e+01 > Line search: Quadratically determined step, lambda=1.5314449290174556e-01 > 3 SNES Function norm 4.945867427788e+01 > 3 SNES VI Function norm 49.4587 Active lower constraints 0/0 upper > constraints 0/0 Percent of total 0. Percent of bounded 0. > Linear solve converged due to CONVERGED_RTOL iterations 2 > Line search: gnorm after quadratic fit 3.912546059379e+01 > Line search: Quadratically determined step, lambda=2.8355714983096392e-01 > 4 SNES Function norm 3.912546059379e+01 > 4 SNES VI Function norm 39.1255 Active lower constraints 0/0 upper > constraints 0/0 Percent of total 0. Percent of bounded 0. > Linear solve converged due to CONVERGED_RTOL iterations 2 > Line search: Using full step: fnorm 3.912546059379e+01 gnorm > 3.171372000187e+01 > 5 SNES Function norm 3.171372000187e+01 > 5 SNES VI Function norm 31.7137 Active lower constraints 0/0 upper > constraints 0/0 Percent of total 0. Percent of bounded 0. > Linear solve converged due to CONVERGED_RTOL iterations 2 > Line search: Using full step: fnorm 3.171372000187e+01 gnorm > 1.036076366647e-01 > 6 SNES Function norm 1.036076366647e-01 > 6 SNES VI Function norm 0.103608 Active lower constraints 0/0 upper > constraints 0/0 Percent of total 0. Percent of bounded 0. > Linear solve converged due to CONVERGED_RTOL iterations 2 > Line search: Using full step: fnorm 1.036076366647e-01 gnorm > 1.913739253850e-06 > 7 SNES Function norm 1.913739253850e-06 > 7 SNES VI Function norm 1.91374e-06 Active lower constraints 0/0 upper > constraints 0/0 Percent of total 0. Percent of bounded 0. > Nonlinear solve converged due to CONVERGED_FNORM_RELATIVE iterations 7 > > > On Aug 20 2018, at 2:00 am, Smith, Barry F. <[email protected]> wrote: > > > Note that no constraints are reached at the first iterations despite the > terrible convergence. I suspect something is wrong with either the > formulation or the Jacobian (that is nothing to do with the constraints). > What do you get with the same options if you do not set any constraints? The > same terrible convergence or something different (send it). In PETSc 3.9 you > can run with the option -snes_test_jacobian to see if the Jacobians are > accurate or out of whack. > > Barry > > > On Aug 19, 2018, at 9:19 AM, Amir <[email protected]> wrote: > > This is what I get after I modify the constraints, and using lu: > 0 SNES Function norm 6.839975279451e+01 > 0 SNES VI Function norm 68.3997 Active lower constraints 0/0 upper > constraints 0/0 Percent of total 0. Percent of bounded 0. > Linear solve converged due to CONVERGED_RTOL iterations 1 > Line search: gnorm after quadratic fit 1.167382139621e+02 > Line search: Cubic step no good, shrinking lambda, current gnorm > 9.661035212791e+01 lambda=1.4620823686291085e-02 > Line search: Cubic step no good, shrinking lambda, current gnorm > 7.872462383177e+01 lambda=1.4620823686291085e-03 > Line search: Cubic step no good, shrinking lambda, current gnorm > 7.010735352812e+01 lambda=1.4620823686291087e-04 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.845164828370e+01 lambda=1.4620823686291088e-05 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.840024692820e+01 lambda=1.4620823686291089e-06 > Line search: Cubically determined step, current gnorm 6.839974867692e+01 > lambda=1.4620823686291089e-07 > 1 SNES Function norm 6.839974867692e+01 > 1 SNES VI Function norm 68.3997 Active lower constraints 0/0 upper > constraints 0/0 Percent of total 0. Percent of bounded 0. > Linear solve converged due to CONVERGED_RTOL iterations 1 > Line search: gnorm after quadratic fit 1.313376791907e+02 > Line search: Cubic step no good, shrinking lambda, current gnorm > 1.102614855398e+02 lambda=1.4115774027030809e-02 > Line search: Cubic step no good, shrinking lambda, current gnorm > 9.096517666290e+01 lambda=1.4115774027030810e-03 > Line search: Cubic step no good, shrinking lambda, current gnorm > 7.587134084601e+01 lambda=1.4115774027030810e-04 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.877273920434e+01 lambda=1.4115774027030810e-05 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.840408582565e+01 lambda=1.4115774027030811e-06 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.839977811504e+01 lambda=1.4115774027030813e-07 > Line search: Cubically determined step, current gnorm 6.839974761445e+01 > lambda=1.7645101306134412e-08 > 2 SNES Function norm 6.839974761445e+01 > 2 SNES VI Function norm 68.3997 Active lower constraints 0/0 upper > constraints 0/0 Percent of total 0. Percent of bounded 0. > Linear solve converged due to CONVERGED_RTOL iterations 2 > Line search: gnorm after quadratic fit 1.053816835919e+02 > Line search: Cubic step no good, shrinking lambda, current gnorm > 1.064875244204e+02 lambda=2.9620424834498244e-02 > Line search: Cubic step no good, shrinking lambda, current gnorm > 8.928160190384e+01 lambda=2.9620424834498247e-03 > Line search: Cubic step no good, shrinking lambda, current gnorm > 7.462621000022e+01 lambda=2.9620424834498249e-04 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.862607236334e+01 lambda=2.9620424834498250e-05 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.840219410268e+01 lambda=2.9620424834498251e-06 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.839975697009e+01 lambda=2.9620424834498251e-07 > Line search: Cubically determined step, current gnorm 6.839974486019e+01 > lambda=1.0073837635308560e-07 > 3 SNES Function norm 6.839974486019e+01 > 3 SNES VI Function norm 68.3997 Active lower constraints 0/0 upper > constraints 0/0 Percent of total 0. Percent of bounded 0. > Linear solve converged due to CONVERGED_RTOL iterations 2 > Line search: gnorm after quadratic fit 1.235990792724e+02 > Line search: Cubic step no good, shrinking lambda, current gnorm > 1.013971718389e+02 lambda=2.3488863243422971e-02 > Line search: Cubic step no good, shrinking lambda, current gnorm > 7.781731557824e+01 lambda=2.3488863243422975e-03 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.872972483427e+01 lambda=2.3488863243422976e-04 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.840207967049e+01 lambda=2.3488863243422978e-05 > Line search: Cubically determined step, current gnorm 6.839957781288e+01 > lambda=4.7431867318372446e-06 > 4 SNES Function norm 6.839957781288e+01 > 4 SNES VI Function norm 68.4009 Active lower constraints 0/0 upper > constraints 0/0 Percent of total 0. Percent of bounded 0. > Linear solve converged due to CONVERGED_RTOL iterations 2 > Line search: gnorm after quadratic fit 1.076412782767e+02 > Line search: Cubic step no good, shrinking lambda, current gnorm > 8.067259042361e+01 lambda=2.9654955573060381e-02 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.991099677071e+01 lambda=2.9654955573060384e-03 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.841319422751e+01 lambda=2.9654955573060385e-04 > Line search: Cubically determined step, current gnorm 6.839661391957e+01 > lambda=8.6072251116571315e-05 > 5 SNES Function norm 6.839661391957e+01 > 5 SNES VI Function norm 68.4308 Active lower constraints 0/0 upper > constraints 0/0 Percent of total 0. Percent of bounded 0. > Linear solve converged due to CONVERGED_RTOL iterations 2 > Line search: gnorm after quadratic fit 9.167521781305e+01 > Line search: Cubic step no good, shrinking lambda, current gnorm > 8.101349837507e+01 lambda=6.4300240519536755e-02 > Line search: Cubic step no good, shrinking lambda, current gnorm > 7.041475585659e+01 lambda=6.6547102741568651e-03 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.844705599568e+01 lambda=6.6547102741568653e-04 > Line search: Cubically determined step, current gnorm 6.839121015135e+01 > lambda=1.4965127147397144e-04 > 6 SNES Function norm 6.839121015135e+01 > 6 SNES VI Function norm 68.5745 Active lower constraints 0/0 upper > constraints 0/0 Percent of total 0. Percent of bounded 0. > Linear solve converged due to CONVERGED_RTOL iterations 3 > Line search: gnorm after quadratic fit 9.587731188484e+01 > Line search: Cubic step no good, shrinking lambda, current gnorm > 8.197287845440e+01 lambda=5.2775823101139679e-02 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.940305546680e+01 lambda=5.2775823101139686e-03 > Line search: Cubically determined step, current gnorm 6.837177518619e+01 > lambda=6.3969948101878912e-04 > 7 SNES Function norm 6.837177518619e+01 > 7 SNES VI Function norm 71.0652 Active lower constraints 0/0 upper > constraints 0/0 Percent of total 0. Percent of bounded 0. > Linear solve converged due to CONVERGED_RTOL iterations 3 > Line search: gnorm after quadratic fit 9.548865181502e+01 > Line search: Cubic step no good, shrinking lambda, current gnorm > 7.876974019161e+01 lambda=5.0026957605696178e-02 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.891985531663e+01 lambda=5.1260581886742491e-03 > Line search: Cubically determined step, current gnorm 6.832012833110e+01 > lambda=9.3300874565775075e-04 > 8 SNES Function norm 6.832012833110e+01 > 8 SNES VI Function norm 79.8569 Active lower constraints 0/0 upper > constraints 0/0 Percent of total 0. Percent of bounded 0. > Linear solve converged due to CONVERGED_RTOL iterations 4 > Line search: gnorm after quadratic fit 9.541909192352e+01 > Line search: Cubic step no good, shrinking lambda, current gnorm > 7.493349067288e+01 lambda=5.1388541796399699e-02 > Line search: Cubically determined step, current gnorm 6.781311229195e+01 > lambda=7.6704381452756282e-03 > 9 SNES Function norm 6.781311229195e+01 > 9 SNES VI Function norm 206.467 Active lower constraints 0/0 upper > constraints 0/0 Percent of total 0. Percent of bounded 0. > Linear solve converged due to CONVERGED_RTOL iterations 6 > Line search: gnorm after quadratic fit 9.394254710662e+01 > Line search: Cubic step no good, shrinking lambda, current gnorm > 7.306526091645e+01 lambda=5.3539445292607950e-02 > Line search: Cubically determined step, current gnorm 6.721835250392e+01 > lambda=9.5561464381222188e-03 > 10 SNES Function norm 6.721835250392e+01 > 10 SNES VI Function norm 368.277 Active lower constraints 0/0 upper > constraints 0/0 Percent of total 0. Percent of bounded 0. > Linear solve converged due to CONVERGED_RTOL iterations 7 > Line search: gnorm after quadratic fit 9.387716261564e+01 > Line search: Cubic step no good, shrinking lambda, current gnorm > 7.070868333837e+01 lambda=4.9443994240769853e-02 > Line search: Cubically determined step, current gnorm 6.660967991379e+01 > lambda=1.0810468882260445e-02 > 11 SNES Function norm 6.660967991379e+01 > 11 SNES VI Function norm 567.089 Active lower constraints 0/0 upper > constraints 0/0 Percent of total 0. Percent of bounded 0. > Linear solve converged due to CONVERGED_RTOL iterations 9 > Line search: gnorm after quadratic fit 9.346306901087e+01 > Line search: Cubic step no good, shrinking lambda, current gnorm > 7.048649174104e+01 lambda=4.9192073733008407e-02 > Line search: Cubically determined step, current gnorm 6.608645472517e+01 > lambda=1.0044112192046832e-02 > 12 SNES Function norm 6.608645472517e+01 > 12 SNES VI Function norm 734.752 Active lower constraints 0/0 upper > constraints 0/0 Percent of total 0. Percent of bounded 0. > Linear solve converged due to CONVERGED_RTOL iterations 8 > Line search: gnorm after quadratic fit 9.452595067062e+01 > Line search: Cubic step no good, shrinking lambda, current gnorm > 7.208837086000e+01 lambda=4.5128845776024190e-02 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.800731174084e+01 lambda=6.5132735322282484e-03 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.634685717588e+01 lambda=6.5132735322282484e-04 > Line search: Cubically determined step, current gnorm 6.608637269301e+01 > lambda=6.5132735322282487e-05 > 13 SNES Function norm 6.608637269301e+01 > 13 SNES VI Function norm 735.814 Active lower constraints 0/0 upper > constraints 0/0 Percent of total 0. Percent of bounded 0. > Linear solve converged due to CONVERGED_RTOL iterations 8 > Line search: gnorm after quadratic fit 9.318693996551e+01 > Line search: Cubic step no good, shrinking lambda, current gnorm > 7.139943393009e+01 lambda=4.8118075850982413e-02 > Line search: Cubically determined step, current gnorm 6.568600823304e+01 > lambda=7.8756632431255542e-03 > 14 SNES Function norm 6.568600823304e+01 > 14 SNES VI Function norm 881.479 Active lower constraints 0/0 upper > constraints 0/0 Percent of total 0. Percent of bounded 0. > Linear solve converged due to CONVERGED_RTOL iterations 9 > Line search: gnorm after quadratic fit 9.242424276014e+01 > Line search: Cubic step no good, shrinking lambda, current gnorm > 7.226307852468e+01 lambda=4.8021995541678457e-02 > Line search: Cubically determined step, current gnorm 6.540144178550e+01 > lambda=6.6931701304969058e-03 > 15 SNES Function norm 6.540144178550e+01 > 15 SNES VI Function norm 1006.28 Active lower constraints 0/0 upper > constraints 0/0 Percent of total 0. Percent of bounded 0. > Linear solve converged due to CONVERGED_RTOL iterations 10 > Line search: gnorm after quadratic fit 9.154992244384e+01 > Line search: Cubic step no good, shrinking lambda, current gnorm > 7.122067564574e+01 lambda=5.0277415210907568e-02 > Line search: Cubically determined step, current gnorm 6.506969399500e+01 > lambda=7.8660529387648906e-03 > 16 SNES Function norm 6.506969399500e+01 > 16 SNES VI Function norm 1131.4 Active lower constraints 0/0 upper > constraints 0/0 Percent of total 0. Percent of bounded 0. > Linear solve converged due to CONVERGED_RTOL iterations 10 > Line search: gnorm after quadratic fit 9.263744099217e+01 > Line search: Cubic step no good, shrinking lambda, current gnorm > 7.045825559590e+01 lambda=4.3976587783298791e-02 > Line search: Cubically determined step, current gnorm 6.477740872258e+01 > lambda=6.6653314620139965e-03 > 17 SNES Function norm 6.477740872258e+01 > 17 SNES VI Function norm 1252.29 Active lower constraints 0/0 upper > constraints 0/0 Percent of total 0. Percent of bounded 0. > Linear solve converged due to CONVERGED_RTOL iterations 10 > Line search: gnorm after quadratic fit 9.034042193429e+01 > Line search: Cubic step no good, shrinking lambda, current gnorm > 7.097176730586e+01 lambda=4.9082849797732621e-02 > Line search: Cubically determined step, current gnorm 6.448770943479e+01 > lambda=7.1730414269226593e-03 > 18 SNES Function norm 6.448770943479e+01 > 18 SNES VI Function norm 1371.72 Active lower constraints 0/0 upper > constraints 0/0 Percent of total 0. Percent of bounded 0. > Linear solve converged due to CONVERGED_RTOL iterations 10 > Line search: gnorm after quadratic fit 9.028256513212e+01 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.990381886066e+01 lambda=4.4193294185915478e-02 > Line search: Cubically determined step, current gnorm 6.422123220407e+01 > lambda=6.6250183714544201e-03 > 19 SNES Function norm 6.422123220407e+01 > 19 SNES VI Function norm 1490.87 Active lower constraints 0/0 upper > constraints 0/0 Percent of total 0. Percent of bounded 0. > Linear solve converged due to CONVERGED_RTOL iterations 10 > Line search: gnorm after quadratic fit 8.908334154738e+01 > Line search: Cubic step no good, shrinking lambda, current gnorm > 7.060161630260e+01 lambda=4.6452357120490895e-02 > Line search: Cubically determined step, current gnorm 6.411445689017e+01 > lambda=6.3417192629449783e-03 > 20 SNES Function norm 6.411445689017e+01 > 20 SNES VI Function norm 1597.29 Active lower constraints 0/0 upper > constraints 0/0 Percent of total 0. Percent of bounded 0. > Linear solve converged due to CONVERGED_RTOL iterations 11 > Line search: gnorm after quadratic fit 8.460267552084e+01 > Line search: Cubic step no good, shrinking lambda, current gnorm > 7.051195181220e+01 lambda=5.8933530173091400e-02 > Line search: Cubically determined step, current gnorm 6.379276529842e+01 > lambda=9.2451193719573169e-03 > 21 SNES Function norm 6.379276529842e+01 > 21 SNES VI Function norm 1729.35 Active lower constraints 0/0 upper > constraints 0/0 Percent of total 0. Percent of bounded 0. > Linear solve converged due to CONVERGED_RTOL iterations 10 > Line search: gnorm after quadratic fit 8.847710713735e+01 > Line search: Cubic step no good, shrinking lambda, current gnorm > 7.007434740347e+01 lambda=4.8620669909360585e-02 > Line search: Cubically determined step, current gnorm 6.358026444281e+01 > lambda=6.8922394760267177e-03 > 22 SNES Function norm 6.358026444281e+01 > 22 SNES VI Function norm 1849.77 Active lower constraints 0/0 upper > constraints 0/0 Percent of total 0. Percent of bounded 0. > Linear solve converged due to CONVERGED_RTOL iterations 10 > Line search: gnorm after quadratic fit 9.435581241699e+01 > Line search: Cubic step no good, shrinking lambda, current gnorm > 7.335475468652e+01 lambda=3.2814238215953598e-02 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.652843583618e+01 lambda=3.2814238215953598e-03 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.529985391051e+01 lambda=3.2814238215953602e-04 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.366726270581e+01 lambda=3.2814238215953603e-05 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.358115196460e+01 lambda=3.2814238215953603e-06 > Line search: Cubically determined step, current gnorm 6.358025861312e+01 > lambda=3.2814238215953605e-07 > 23 SNES Function norm 6.358025861312e+01 > 23 SNES VI Function norm 1849.78 Active lower constraints 0/0 upper > constraints 0/0 Percent of total 0. Percent of bounded 0. > Linear solve converged due to CONVERGED_RTOL iterations 10 > Line search: gnorm after quadratic fit 8.812900895448e+01 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.822426303212e+01 lambda=4.8421685266110105e-02 > Line search: Cubically determined step, current gnorm 6.335842666537e+01 > lambda=8.4588034499471741e-03 > 24 SNES Function norm 6.335842666537e+01 > 24 SNES VI Function norm 1993.6 Active lower constraints 0/0 upper > constraints 0/0 Percent of total 0. Percent of bounded 0. > Linear solve converged due to CONVERGED_RTOL iterations 11 > Line search: gnorm after quadratic fit 9.094159193064e+01 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.742173830047e+01 lambda=4.3243453252109078e-02 > Line search: Cubically determined step, current gnorm 6.316236273224e+01 > lambda=7.7494584282214512e-03 > 25 SNES Function norm 6.316236273224e+01 > 25 SNES VI Function norm 2126.62 Active lower constraints 0/0 upper > constraints 0/0 Percent of total 0. Percent of bounded 0. > Linear solve converged due to CONVERGED_RTOL iterations 11 > Line search: gnorm after quadratic fit 9.160540289460e+01 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.722886652974e+01 lambda=4.1792497999647044e-02 > Line search: Cubically determined step, current gnorm 6.300020419956e+01 > lambda=7.3335939016311632e-03 > 26 SNES Function norm 6.300020419956e+01 > 26 SNES VI Function norm 2252.12 Active lower constraints 0/0 upper > constraints 0/0 Percent of total 0. Percent of bounded 0. > Linear solve converged due to CONVERGED_RTOL iterations 10 > Line search: gnorm after quadratic fit 8.890109786397e+01 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.686097246249e+01 lambda=4.7041389395739201e-02 > Line search: Cubically determined step, current gnorm 6.283084254493e+01 > lambda=9.0792221408297433e-03 > 27 SNES Function norm 6.283084254493e+01 > 27 SNES VI Function norm 2404.84 Active lower constraints 0/0 upper > constraints 0/0 Percent of total 0. Percent of bounded 0. > Linear solve converged due to CONVERGED_RTOL iterations 11 > Line search: gnorm after quadratic fit 8.998798828268e+01 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.576469464321e+01 lambda=4.2887406731046039e-02 > Line search: Cubically determined step, current gnorm 6.266968273817e+01 > lambda=9.2948647727201898e-03 > 28 SNES Function norm 6.266968273817e+01 > 28 SNES VI Function norm 2558.11 Active lower constraints 0/0 upper > constraints 0/0 Percent of total 0. Percent of bounded 0. > Linear solve converged due to CONVERGED_RTOL iterations 11 > Line search: gnorm after quadratic fit 9.199146621649e+01 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.564298288247e+01 lambda=4.0193768728643060e-02 > Line search: Cubically determined step, current gnorm 6.252861090423e+01 > lambda=8.3932894622200105e-03 > 29 SNES Function norm 6.252861090423e+01 > 29 SNES VI Function norm 2698.56 Active lower constraints 0/0 upper > constraints 0/0 Percent of total 0. Percent of bounded 0. > Linear solve converged due to CONVERGED_RTOL iterations 11 > Line search: gnorm after quadratic fit 9.382191898563e+01 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.772162800128e+01 lambda=3.5914901407691102e-02 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.485302824556e+01 lambda=4.7901702628153379e-03 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.285215446174e+01 lambda=4.7901702628153379e-04 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.253029194340e+01 lambda=4.7901702628153380e-05 > Line search: Cubically determined step, current gnorm 6.252812967035e+01 > lambda=1.5127710264417741e-05 > 30 SNES Function norm 6.252812967035e+01 > 30 SNES VI Function norm 2698.83 Active lower constraints 0/0 upper > constraints 0/0 Percent of total 0. Percent of bounded 0. > Linear solve converged due to CONVERGED_RTOL iterations 11 > Line search: gnorm after quadratic fit 9.118652411357e+01 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.495202040381e+01 lambda=4.0596953253201326e-02 > Line search: Cubically determined step, current gnorm 6.239041273754e+01 > lambda=9.5128837682561859e-03 > 31 SNES Function norm 6.239041273754e+01 > 31 SNES VI Function norm 2854.38 Active lower constraints 0/0 upper > constraints 0/0 Percent of total 0. Percent of bounded 0. > Linear solve converged due to CONVERGED_RTOL iterations 11 > Line search: gnorm after quadratic fit 9.077825361440e+01 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.467284344645e+01 lambda=4.0583557291813177e-02 > Line search: Cubically determined step, current gnorm 6.225884047363e+01 > lambda=9.7846838097028124e-03 > 32 SNES Function norm 6.225884047363e+01 > 32 SNES VI Function norm 3012.55 Active lower constraints 0/0 upper > constraints 0/0 Percent of total 0. Percent of bounded 0. > Linear solve converged due to CONVERGED_RTOL iterations 11 > Line search: gnorm after quadratic fit 8.920999190500e+01 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.466103609666e+01 lambda=4.2498592612658738e-02 > Line search: Cubically determined step, current gnorm 6.213582712328e+01 > lambda=1.0171220001350099e-02 > 33 SNES Function norm 6.213582712328e+01 > 33 SNES VI Function norm 3175.23 Active lower constraints 0/0 upper > constraints 0/0 Percent of total 0. Percent of bounded 0. > Linear solve converged due to CONVERGED_RTOL iterations 11 > Line search: gnorm after quadratic fit 9.011924620812e+01 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.495110751676e+01 lambda=4.2453491229013124e-02 > Line search: Cubically determined step, current gnorm 6.201636649414e+01 > lambda=9.3685443832871688e-03 > 34 SNES Function norm 6.201636649414e+01 > 34 SNES VI Function norm 3322.22 Active lower constraints 0/0 upper > constraints 0/0 Percent of total 0. Percent of bounded 0. > Linear solve converged due to CONVERGED_RTOL iterations 11 > Line search: gnorm after quadratic fit 9.006024690051e+01 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.414041989575e+01 lambda=4.0940226638739433e-02 > Line search: Cubically determined step, current gnorm 6.190218472271e+01 > lambda=1.0273994032172384e-02 > 35 SNES Function norm 6.190218472271e+01 > 35 SNES VI Function norm 3484.25 Active lower constraints 0/0 upper > constraints 0/0 Percent of total 0. Percent of bounded 0. > Linear solve converged due to CONVERGED_RTOL iterations 11 > Line search: gnorm after quadratic fit 8.898991234613e+01 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.426189506521e+01 lambda=4.3426020819675164e-02 > Line search: Cubically determined step, current gnorm 6.178678684779e+01 > lambda=1.0559984512509757e-02 > 36 SNES Function norm 6.178678684779e+01 > 36 SNES VI Function norm 3648.36 Active lower constraints 0/0 upper > constraints 0/0 Percent of total 0. Percent of bounded 0. > Linear solve converged due to CONVERGED_RTOL iterations 11 > Line search: gnorm after quadratic fit 8.648025040907e+01 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.404650680722e+01 lambda=4.6749921483990771e-02 > Line search: Cubically determined step, current gnorm 6.165099494681e+01 > lambda=1.1866483746353381e-02 > 37 SNES Function norm 6.165099494681e+01 > 37 SNES VI Function norm 3822.88 Active lower constraints 0/0 upper > constraints 0/0 Percent of total 0. Percent of bounded 0. > Linear solve converged due to CONVERGED_RTOL iterations 11 > Line search: gnorm after quadratic fit 8.858076523292e+01 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.520416756669e+01 lambda=4.2063577748779259e-02 > Line search: Cubically determined step, current gnorm 6.152598711257e+01 > lambda=7.9475995139781722e-03 > 38 SNES Function norm 6.152598711257e+01 > 38 SNES VI Function norm 3944.22 Active lower constraints 0/0 upper > constraints 0/0 Percent of total 0. Percent of bounded 0. > Linear solve converged due to CONVERGED_RTOL iterations 11 > Line search: gnorm after quadratic fit 8.742192022113e+01 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.343055909293e+01 lambda=3.8810591937350140e-02 > Line search: Cubically determined step, current gnorm 6.135131913134e+01 > lambda=9.8750504088198892e-03 > 39 SNES Function norm 6.135131913134e+01 > 39 SNES VI Function norm 4066.15 Active lower constraints 0/1 upper > constraints 0/0 Percent of total 0. Percent of bounded 0. > Linear solve converged due to CONVERGED_RTOL iterations 11 > Line search: gnorm after quadratic fit 9.286112632409e+01 > Line search: Cubic step no good, shrinking lambda, current gnorm > 7.061265007302e+01 lambda=2.9694426041134040e-02 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.436749315348e+01 lambda=2.9694426041134044e-03 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.299789903123e+01 lambda=2.9694426041134044e-04 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.141298285103e+01 lambda=2.9694426041134044e-05 > Line search: Cubically determined step, current gnorm 6.135130530880e+01 > lambda=2.9694426041134044e-06 > 40 SNES Function norm 6.135130530880e+01 > 40 SNES VI Function norm 4066.23 Active lower constraints 0/1 upper > constraints 0/0 Percent of total 0. Percent of bounded 0. > Linear solve converged due to CONVERGED_RTOL iterations 11 > Line search: gnorm after quadratic fit 8.660062754407e+01 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.674944119178e+01 lambda=4.3170190544802138e-02 > Line search: Cubically determined step, current gnorm 6.120973541057e+01 > lambda=6.1770273873522300e-03 > 41 SNES Function norm 6.120973541057e+01 > 41 SNES VI Function norm 4152.83 Active lower constraints 0/1 upper > constraints 0/0 Percent of total 0. Percent of bounded 0. > Linear solve converged due to CONVERGED_RTOL iterations 11 > Line search: gnorm after quadratic fit 8.564645091492e+01 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.676153862982e+01 lambda=4.4220462778047379e-02 > Line search: Cubically determined step, current gnorm 6.104393345352e+01 > lambda=6.3589031499342100e-03 > 42 SNES Function norm 6.104393345352e+01 > 42 SNES VI Function norm 4233.82 Active lower constraints 0/1 upper > constraints 0/0 Percent of total 0. Percent of bounded 0. > Linear solve converged due to CONVERGED_RTOL iterations 11 > Line search: gnorm after quadratic fit 8.368658054045e+01 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.715427252993e+01 lambda=4.7003756543045970e-02 > Line search: Cubically determined step, current gnorm 6.101753259623e+01 > lambda=6.3787492486959132e-03 > 43 SNES Function norm 6.101753259623e+01 > 43 SNES VI Function norm 4314.49 Active lower constraints 0/1 upper > constraints 0/0 Percent of total 0. Percent of bounded 0. > Linear solve converged due to CONVERGED_RTOL iterations 11 > Line search: gnorm after quadratic fit 8.442307716417e+01 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.496496980051e+01 lambda=4.5726539795380154e-02 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.142635208964e+01 lambda=8.4157732414286580e-03 > Line search: Cubically determined step, current gnorm 6.095430498328e+01 > lambda=2.1294466610363495e-03 > 44 SNES Function norm 6.095430498328e+01 > 44 SNES VI Function norm 4341.47 Active lower constraints 0/1 upper > constraints 0/0 Percent of total 0. Percent of bounded 0. > Linear solve converged due to CONVERGED_RTOL iterations 12 > Line search: gnorm after quadratic fit 8.439891204437e+01 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.240254752159e+01 lambda=4.3289599290802970e-02 > Line search: Cubically determined step, current gnorm 6.081365387079e+01 > lambda=1.2916446928212236e-02 > 45 SNES Function norm 6.081365387079e+01 > 45 SNES VI Function norm 4500.53 Active lower constraints 0/1 upper > constraints 0/0 Percent of total 0. Percent of bounded 0. > Linear solve converged due to CONVERGED_RTOL iterations 11 > Line search: gnorm after quadratic fit 9.664056729635e+01 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.944851102419e+01 lambda=1.8635858087631831e-02 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.394837998027e+01 lambda=1.8635858087631832e-03 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.278442438314e+01 lambda=1.8635858087631834e-04 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.093400706607e+01 lambda=1.8635858087631834e-05 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.081489139942e+01 lambda=1.8635858087631835e-06 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.081365630503e+01 lambda=1.8635858087631836e-07 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.081365482998e+01 lambda=7.3438624052552081e-08 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.081365418254e+01 lambda=2.3870801624167493e-08 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.081365397577e+01 lambda=8.0330117961387493e-09 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.081365390593e+01 lambda=2.6874115312641255e-09 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.081365388261e+01 lambda=8.9982494953058537e-10 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.081365387467e+01 lambda=3.0083351315759278e-10 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.081365387219e+01 lambda=1.0110298276678126e-10 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.081365387129e+01 lambda=3.3292703628714462e-11 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.081365387093e+01 lambda=1.0739379834615402e-11 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.081365387083e+01 lambda=3.6237724732408161e-12 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.081365387082e+01 lambda=1.2456833192766385e-12 > Line search: Cubically determined step, current gnorm 6.081365387079e+01 > lambda=3.5247146389212346e-13 > 46 SNES Function norm 6.081365387079e+01 > 46 SNES VI Function norm 4500.53 Active lower constraints 0/1 upper > constraints 0/0 Percent of total 0. Percent of bounded 0. > Linear solve converged due to CONVERGED_RTOL iterations 11 > Line search: gnorm after quadratic fit 1.433547268242e+03 > Line search: Cubic step no good, shrinking lambda, current gnorm > 7.106740636849e+01 lambda=2.4859604812085913e-02 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.525497375887e+01 lambda=7.2498301889276232e-03 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.367730849091e+01 lambda=7.2498301889276236e-04 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.174333248068e+01 lambda=7.2498301889276242e-05 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.083400168850e+01 lambda=7.2498301889276242e-06 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.081378891362e+01 lambda=7.2498301889276242e-07 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.081365527548e+01 lambda=1.0753239548189211e-07 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.081365445012e+01 lambda=4.4350670800939511e-08 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.081365407923e+01 lambda=1.5953160232199837e-08 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.081365394762e+01 lambda=5.8822684037182187e-09 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.081365389903e+01 lambda=2.1599843490971948e-09 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.081365388111e+01 lambda=7.9357831777119289e-10 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.081365387457e+01 lambda=2.9181178721013368e-10 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.081365387223e+01 lambda=1.0733281860745436e-10 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.081365387132e+01 lambda=3.9236995040615546e-11 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.081365387106e+01 lambda=1.4351792294805992e-11 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.081365387081e+01 lambda=4.8683006404131606e-12 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.081365387087e+01 lambda=2.0863028416245319e-12 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.081365387080e+01 lambda=5.3058618006527772e-13 > Line search: unable to find good step length! After 19 tries > 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 > 47 SNES Function norm 6.081365387079e+01 > 47 SNES VI Function norm 4500.53 Active lower constraints 0/1 upper > constraints 0/0 Percent of total 0. Percent of bounded 0. > Linear solve converged due to CONVERGED_RTOL iterations 11 > Line search: gnorm after quadratic fit 1.280890413731e+02 > Line search: Cubic step no good, shrinking lambda, current gnorm > 7.053050135324e+01 lambda=2.2733898438250671e-02 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.399979736493e+01 lambda=2.2733898438250673e-03 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.297299206174e+01 lambda=2.2733898438250674e-04 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.098427084652e+01 lambda=2.2733898438250677e-05 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.081555475063e+01 lambda=2.2733898438250679e-06 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.081365745105e+01 lambda=2.2733898438250681e-07 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.081365507363e+01 lambda=9.2077672155667702e-08 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.081365429323e+01 lambda=3.2329991301020441e-08 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.081365402243e+01 lambda=1.1595563414967544e-08 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.081365392495e+01 lambda=4.1432803927719340e-09 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.081365389024e+01 lambda=1.4815603373453140e-09 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.081365387781e+01 lambda=5.2922825349649990e-10 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.081365387324e+01 lambda=1.8866396505817014e-10 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.081365387173e+01 lambda=6.7568116386381729e-11 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.081365387107e+01 lambda=2.3800109593117919e-11 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.081365387094e+01 lambda=8.6880387507986266e-12 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.081365387083e+01 lambda=2.8858967256610185e-12 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.081365387082e+01 lambda=1.0348935437400166e-12 > Line search: Cubically determined step, current gnorm 6.081365387077e+01 > lambda=2.7070538024707468e-13 > 48 SNES Function norm 6.081365387077e+01 > 48 SNES VI Function norm 4500.53 Active lower constraints 0/1 upper > constraints 0/0 Percent of total 0. Percent of bounded 0. > Linear solve converged due to CONVERGED_RTOL iterations 11 > Line search: gnorm after quadratic fit 1.170943716019e+02 > Line search: Cubic step no good, shrinking lambda, current gnorm > 7.038597782185e+01 lambda=2.2207476459146142e-02 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.399318571001e+01 lambda=2.2207476459146142e-03 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.295164029113e+01 lambda=2.2207476459146143e-04 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.097743654776e+01 lambda=2.2207476459146143e-05 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.081546079136e+01 lambda=2.2207476459146144e-06 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.081365719872e+01 lambda=2.2207476459146145e-07 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.081365503919e+01 lambda=9.0117449908413481e-08 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.081365427749e+01 lambda=3.1364523245356123e-08 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.081365401596e+01 lambda=1.1183735694685806e-08 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.081365392232e+01 lambda=3.9706025565052754e-09 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.081365388909e+01 lambda=1.4107330444508480e-09 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.081365387724e+01 lambda=5.0110371464381371e-10 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.081365387308e+01 lambda=1.7820909495605996e-10 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.081365387164e+01 lambda=6.3314737242538209e-11 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.081365387114e+01 lambda=2.2207910898013729e-11 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.081365387085e+01 lambda=7.4515329000007742e-12 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.081365387082e+01 lambda=2.7872933392432575e-12 > Line search: Cubic step no good, shrinking lambda, current gnorm > 6.081365387083e+01 lambda=8.8269934321345058e-13 > Line search: unable to find good step length! After 18 tries > 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 > Nonlinear solve did not converge due to DIVERGED_LINE_SEARCH iterations 48 > > Amir > > On Aug 18 2018, at 11:42 pm, Smith, Barry F. <[email protected]> wrote: > > > -snes_monitor -snes_linesearch_monitor -pc_type lu > > produces what? > > Barry > > 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. > > > On Aug 17, 2018, at 11:56 PM, Amir <[email protected]> wrote: > > Hi > 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. > I think the way I set VI maybe wrong. Sorry for my simple explanation. > Very thankful for your time. > Amir > SNESCreate(PETSC_COMM_WORLD,&snes); > SNESSetApplicationContext(snes,&user); > SNESSetDM(snes,user.dm); > SNESSetFunction(snes,NULL,FormFunction,&user); > FormInitialSolution(user.X,&user); > SNESSetSolution(snes,user.X); > SNESMonitorSet(snes,MySNESMonitor,&user,NULL); > SNESSetType(snes,SNESVINEWTONSSLS); > SNESSetNormSchedule(snes, SNES_NORM_ALWAYS); > SNESVISetVariableBounds(snes,user.XLOWER,user.XUPPER); > SNESGetLineSearch(snes,&linesearch); > SNESLineSearchSetType(linesearch,SNESLINESEARCHBT); > SNESLineSearchSetVIFunctions(linesearch, NULL, SNESNormFunction); > SNESSetFromOptions(snes); > SNESSolve(snes,NULL,user.X);
