[petsc-users] Setting SNESVI

Sun Aug 19 13:00:01 CDT 2018

   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 <mhbaghaei at mail.sjtu.edu.cn> 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. <bsmith at mcs.anl.gov> 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 <mhbaghaei at mail.sjtu.edu.cn> 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);
> 