On Mon, Jun 20, 2011 at 10:19 PM, Tian(ICT) <rongtian at ncic.ac.cn> wrote:
> ** > Thanks for quick answering. > I looked into SNESSetFromOptions() in the debugger and MF was not chosen. > Still the wired thing is the code never breaks at the point I set inside > FormJacobian when running in the debug mode. > I will leave this issue alone again and look at it later. > > Actually I am struggling with a SNES solve. The followings are some > outputs. > -snes_ratiomonitor shows that residual norm almost equals to one but the > SNES converged, > is this a good convergence? > No. Are you sure this system is consistent, meaning it has a solution? Have you tried putting in an exact solution for a simpler problem? Matt > > atol=1e-050, rtol=1e-008, stol=1e-008, maxit=50, maxf=10000 > function................. > residual force................. > 0 > 0 > 0 > 0 > 0 > 0 > 0 > 0 > 0 > 0 > 0 > 0 > 0 > 40000 <== initial value > 0 > 0 > 40000 > 0 > 0 > 40000 > 0 > 0 > 40000 > 0 > 0 SNES Function norm 8.000000000000e+004 > 0 SNES Function norm 8.000000000000e+004 > ------------------->Jacobian > function................. > residual force................. > -8.2245e+007 > -5.12049e+006 > 1.03104e+008 > 2.44985e+007 > -1.11118e+007 > -3.144e+007 > -2.37583e+007 > -535905 > 2.79206e+007 > 5.63514e+007 > -2.26332e+007 > -7.54711e+007 > -5.76723e+007 > 1.544e+007 > 7.56104e+007 > 2.34587e+007 > 9.57426e+006 > -2.56718e+007 > 2.7089e+007 > 3.32996e+007 > -3.09772e+007 > 3.2278e+007 > -1.87524e+007 > -4.30752e+007 > function................. > residual force................. > -8.2245e+007 > -5.12049e+006 > 1.03104e+008 > 2.44985e+007 > -1.11118e+007 > -3.144e+007 > -2.37583e+007 > -535905 > 2.79206e+007 > 5.63514e+007 > -2.26332e+007 > -7.54711e+007 > -5.76723e+007 > 1.544e+007 > 7.56104e+007 > 2.34587e+007 > 9.57426e+006 > -2.56718e+007 > 2.7089e+007 > 3.32996e+007 > -3.09772e+007 > 3.2278e+007 > -1.87524e+007 > -4.30752e+007 > function................. > residual force................. > -2.05272e+007 > 3.44926e+006 > 2.5202e+007 > -1.55108e+006 > -9.56839e+006 > 2.17527e+006 > 3.56516e+006 > 6.12324e+006 > -5.9437e+006 > 1.41997e+007 > -1.21796e+007 > -1.92224e+007 > -1.63165e+007 > 1.00523e+007 > 2.14974e+007 > -1.87655e+006 > -3.44175e+006 > 4.3047e+006 > 1.50918e+007 > 1.66532e+007 > -1.8379e+007 > 7.41472e+006 > -1.09283e+007 > -9.63423e+006 > function................. > residual force................. > -4.83477e+006 > 1.9056e+006 > 5.06631e+006 > -1.74833e+006 > -3.81725e+006 > 3.0505e+006 > 1.93196e+006 > 2.614e+006 > -3.9812e+006 > 3.71686e+006 > -4.00712e+006 > -4.46187e+006 > -4.24801e+006 > 3.70415e+006 > 5.07447e+006 > -1.59935e+006 > -2.18778e+006 > 3.63332e+006 > 3.74435e+006 > 5.0913e+006 > -5.46539e+006 > 3.03729e+006 > -3.14292e+006 > -2.91613e+006 > function................. > residual force................. > -1.62518e+006 > 635442 > 1.14687e+006 > -395026 > -1.35624e+006 > 1.30091e+006 > 258576 > 922003 > -1.54824e+006 > 1.34306e+006 > -1.18937e+006 > -1.19905e+006 > -1.44957e+006 > 1.18631e+006 > 1.32948e+006 > -233476 > -848651 > 1.52147e+006 > 630822 > 1.59179e+006 > -1.62993e+006 > 1.4708e+006 > -781281 > -921506 > function................. > residual force................. > -680584 > 150805 > 298173 > -4235.39 > -556100 > 495315 > -83367.9 > 401881 > -609708 > 583119 > -325632 > -382676 > -608984 > 356434 > 416726 > 79947.2 > -362413 > 611645 > 24823.1 > 640620 > -552863 > 689283 > -145594 > -276612 > function................. > residual force................. > -252778 > -43702.9 > 107996 > 39434.6 > -252471 > 108972 > -56021.1 > 202917 > -167131 > 237714 > -13235.2 > -160064 > -245478 > 46386.1 > 170214 > 53958.6 > -164661 > 167576 > -41053 > 309826 > -119435 > 264223 > 74941.4 > -108128 > function................. > residual force................. > -37722.5 > -30354.3 > 12055.9 > 19055.1 > -45688.3 > -1869.42 > -17961.6 > 41506 > -11175 > 39667.1 > 24371.9 > -24644.8 > -40629.2 > 12856 > 25503.9 > 17749.7 > -1976.73 > 10689.7 > -20792.3 > 89625.4 > 2237.53 > 40633.6 > 69659.9 > -12797.9 > function................. > residual force................. > -3191.46 > -3795.77 > 508.943 > 2612.15 > -4047.01 > -1039.77 > -2394.34 > 3957.38 > -341.114 > 3526.13 > 3531.34 > -1843.35 > -3574.59 > 36099.4 > 1775.2 > 2391.73 > 35967.4 > 179.561 > -2883.93 > 44445.2 > 1264.25 > 3514.31 > 43842.1 > -503.71 > function................. > residual force................. > -268.146 > -347.447 > 46.5142 > 257.768 > -335.67 > -97.6787 > -234.991 > 341.388 > -35.3256 > 299.924 > 338.987 > -161.929 > -306.308 > 39615 > 132.745 > 236.496 > 39642.5 > -7.34212 > -285.009 > 40378.6 > 141.125 > 300.268 > 40366.6 > -18.108 > function................. > residual force................. > -22.758 > -27.2585 > 11.605 > 23.791 > -26.2078 > -3.38911 > -21.3068 > 28.26 > -10.762 > 25.2783 > 28.3454 > -19.8467 > -27.0516 > 39963.9 > 11.1365 > 21.8233 > 39967.5 > -1.27698 > -26.2771 > 40032.7 > 13.3758 > 26.5009 > 40032.8 > -0.842484 > 1 SNES Function norm 7.999847670241e+004 > 1 SNES Function norm 7.999847670241e+004 0.999981 > ------------------->Jacobian > function................. > residual force................. > -223.361 > -326.156 > 36.1377 > 229.267 > -321.922 > -82.1057 > -210.332 > 329.183 > -52.6793 > 249.815 > 322.469 > -153.907 > -255.672 > 39637.8 > 127.173 > 211.072 > 39655.1 > 12.6163 > -250.93 > 40359 > 119.951 > 250.141 > 40344.5 > -7.18602 > function................. > residual force................. > -69.9794 > -97.0147 > 18.7512 > 72.7704 > -94.1566 > -22.1285 > -66.1397 > 98.02 > -20.4693 > 78.3769 > 97.6066 > -52.1529 > -81.3498 > 39886.4 > 37.8587 > 66.7162 > 39894.4 > 0.222982 > -79.8688 > 40108.4 > 40.0624 > 79.4743 > 40106.3 > -2.14464 > 2 SNES Function norm 7.999823436070e+004 > 2 SNES Function norm 7.999823436070e+004 0.999997 > ------------------->Jacobian > function................. > residual force................. > -85.9273 > -122.698 > 20.3604 > 89.7785 > -119.971 > -28.2927 > -81.9324 > 124.468 > -25.2852 > 96.1213 > 123.187 > -63.7311 > -99.4852 > 39858.1 > 47.741 > 82.4929 > 39866.9 > 2.51439 > -98.3175 > 40136.8 > 48.6046 > 97.2698 > 40133.2 <== almost no change > -1.91145 > 3 SNES Function norm 7.999822446658e+004 > 3 SNES Function norm 7.999822446658e+004 1 > ------------------->Jacobian > function................. > residual force................. > -90.9631 > -130.918 > 20.7647 > 95.1669 > -128.284 > -30.1982 > -86.9654 > 132.973 > -26.9162 > 101.696 > 131.372 > -67.3645 > -105.178 > 39849.1 > 50.8645 > 87.519 > 39858 > 3.38427 > -104.154 > 40145.9 > 51.2208 > 102.878 > 40141.8 > -1.7555 > function................. > residual force................. > -87.4304 > -125.151 > 20.4826 > 91.387 > -122.451 > -28.8623 > -83.4344 > 127.006 > -25.7711 > 97.7856 > 125.63 > -64.8168 > -101.185 > 39855.4 > 48.6736 > 83.9928 > 39864.2 > 2.77216 > -100.06 > 40139.5 > 49.3875 > 98.944 > 40135.8 > -1.86551 > function................. > residual force................. > -86.3776 > -123.432 > 20.3972 > 90.2603 > -120.714 > -28.4634 > -82.3823 > 125.228 > -25.4307 > 96.6198 > 123.919 > -64.0564 > -99.9944 > 39857.3 > 48.0204 > 82.9421 > 39866.1 > 2.59144 > -98.8394 > 40137.6 > 48.8393 > 97.7713 > 40134 > -1.89774 > function................. > residual force................. > -86.0623 > -122.918 > 20.3714 > 89.923 > -120.194 > -28.3439 > -82.0673 > 124.696 > -25.3288 > 96.2708 > 123.407 > -63.8287 > -99.6379 > 39857.9 > 47.8248 > 82.6276 > 39866.6 > 2.53748 > -98.474 > 40137.1 > 48.675 > 97.4202 > 40133.5 > -1.90734 > function................. > residual force................. > -85.9678 > -122.764 > 20.3637 > 89.8218 > -120.038 > -28.3081 > -81.9729 > 124.536 > -25.2983 > 96.1661 > 123.253 > -63.7604 > -99.531 > 39858 > 47.7661 > 82.5333 > 39866.8 > 2.52132 > -98.3644 > 40136.9 > 48.6258 > 97.3149 > 40133.3 > -1.91022 > function................. > residual force................. > -85.9395 > -122.717 > 20.3614 > 89.7915 > -119.991 > -28.2973 > -81.9445 > 124.489 > -25.2891 > 96.1347 > 123.207 > -63.7399 > -99.4989 > 39858.1 > 47.7485 > 82.505 > 39866.9 > 2.51647 > -98.3316 > 40136.8 > 48.611 > 97.2833 > 40133.2 > -1.91108 > function................. > residual force................. > -85.931 > -122.704 > 20.3607 > 89.7824 > -119.977 > -28.2941 > -81.936 > 124.474 > -25.2863 > 96.1253 > 123.193 > -63.7337 > -99.4893 > 39858.1 > 47.7432 > 82.4965 > 39866.9 > 2.51502 > -98.3217 > 40136.8 > 48.6065 > 97.2738 > 40133.2 > -1.91134 > function................. > residual force................. > -85.9284 > -122.699 > 20.3605 > 89.7796 > -119.973 > -28.2931 > -81.9335 > 124.47 > -25.2855 > 96.1225 > 123.189 > -63.7319 > -99.4864 > 39858.1 > 47.7416 > 82.494 > 39866.9 > 2.51458 > -98.3187 > 40136.8 > 48.6052 > 97.271 > 40133.2 > -1.91142 > function................. > residual force................. > -85.9277 > -122.698 > 20.3604 > 89.7788 > -119.972 > -28.2928 > -81.9327 > 124.469 > -25.2853 > 96.1216 > 123.188 > -63.7313 > -99.4856 > 39858.1 > 47.7412 > 82.4932 > 39866.9 > 2.51445 > -98.3179 > 40136.8 > 48.6048 > 97.2702 > 40133.2 > -1.91144 > function................. > residual force................. > -85.9274 > -122.698 > 20.3604 > 89.7786 > -119.972 > -28.2927 > -81.9325 > 124.468 > -25.2852 > 96.1214 > 123.187 > -63.7312 > -99.4853 > 39858.1 > 47.741 > 82.493 > 39866.9 > 2.51441 > -98.3176 > 40136.8 > 48.6047 > 97.2699 > 40133.2 > -1.91145 > function................. > residual force................. > -85.9274 > -122.698 > 20.3604 > 89.7785 > -119.971 > -28.2927 > -81.9324 > 124.468 > -25.2852 > 96.1213 > 123.187 > -63.7311 > -99.4852 > 39858.1 > 47.741 > 82.4929 > 39866.9 > 2.5144 > -98.3175 > 40136.8 > 48.6047 > 97.2698 > 40133.2 > -1.91145 > function................. > residual force................. > -85.9273 > -122.698 > 20.3604 > 89.7785 > -119.971 > -28.2927 > -81.9324 > 124.468 > -25.2852 > 96.1213 > 123.187 > -63.7311 > -99.4852 > 39858.1 > 47.741 > 82.4929 > 39866.9 > 2.5144 > -98.3175 > 40136.8 > 48.6046 > 97.2698 > 40133.2 > -1.91145 > function................. > residual force................. > -85.9273 > -122.698 > 20.3604 > 89.7785 > -119.971 > -28.2927 > -81.9324 > 124.468 > -25.2852 > 96.1213 > 123.187 > -63.7311 > -99.4852 > 39858.1 > 47.741 > 82.4929 > 39866.9 > 2.5144 > -98.3175 > 40136.8 > 48.6046 > 97.2698 > 40133.2 > -1.91145 > function................. > residual force................. > -85.9273 > -122.698 > 20.3604 > 89.7785 > -119.971 > -28.2927 > -81.9324 > 124.468 > -25.2852 > 96.1213 > 123.187 > -63.7311 > -99.4852 > 39858.1 > 47.741 > 82.4929 > 39866.9 > 2.51439 > -98.3175 > 40136.8 > 48.6046 > 97.2698 > 40133.2 > -1.91145 > function................. > residual force................. > -85.9273 > -122.698 > 20.3604 > 89.7785 > -119.971 > -28.2927 > -81.9324 > 124.468 > -25.2852 > 96.1213 > 123.187 > -63.7311 > -99.4852 > 39858.1 > 47.741 > 82.4929 > 39866.9 > 2.51439 > -98.3175 > 40136.8 > 48.6046 > 97.2698 > 40133.2 > -1.91145 > function................. > residual force................. > -85.9273 > -122.698 > 20.3604 > 89.7785 > -119.971 > -28.2927 > -81.9324 > 124.468 > -25.2852 > 96.1213 > 123.187 > -63.7311 > -99.4852 > 39858.1 > 47.741 > 82.4929 > 39866.9 > 2.51439 > -98.3175 > 40136.8 > 48.6046 > 97.2698 > 40133.2 > -1.91145 > function................. > residual force................. > -85.9273 > -122.698 > 20.3604 > 89.7785 > -119.971 > -28.2927 > -81.9324 > 124.468 > -25.2852 > 96.1213 > 123.187 > -63.7311 > -99.4852 > 39858.1 > 47.741 > 82.4929 > 39866.9 > 2.51439 > -98.3175 > 40136.8 > 48.6046 > 97.2698 > 40133.2 > -1.91145 > function................. > residual force................. > -85.9273 > -122.698 > 20.3604 > 89.7785 > -119.971 > -28.2927 > -81.9324 > 124.468 > -25.2852 > 96.1213 > 123.187 > -63.7311 > -99.4852 > 39858.1 > 47.741 > 82.4929 > 39866.9 > 2.51439 > -98.3175 > 40136.8 > 48.6046 > 97.2698 > 40133.2 > -1.91145 > function................. > residual force................. > -85.9273 > -122.698 > 20.3604 > 89.7785 > -119.971 > -28.2927 > -81.9324 > 124.468 > -25.2852 > 96.1213 > 123.187 > -63.7311 > -99.4852 > 39858.1 > 47.741 > 82.4929 > 39866.9 > 2.51439 > -98.3175 > 40136.8 > 48.6046 > 97.2698 > 40133.2 > -1.91145 > function................. > residual force................. > -85.9273 > -122.698 > 20.3604 > 89.7785 > -119.971 > -28.2927 > -81.9324 > 124.468 > -25.2852 > 96.1213 > 123.187 > -63.7311 > -99.4852 > 39858.1 > 47.741 > 82.4929 > 39866.9 > 2.51439 > -98.3175 > 40136.8 > 48.6046 > 97.2698 > 40133.2 > -1.91145 > function................. > residual force................. > -85.9273 > -122.698 > 20.3604 > 89.7785 > -119.971 > -28.2927 > -81.9324 > 124.468 > -25.2852 > 96.1213 > 123.187 > -63.7311 > -99.4852 > 39858.1 > 47.741 > 82.4929 > 39866.9 > 2.51439 > -98.3175 > 40136.8 > 48.6046 > 97.2698 > 40133.2 > -1.91145 > 4 SNES Function norm 7.999822446658e+004 > 4 SNES Function norm 7.999822446658e+004 1 > SNES Object: > type: tr > mu=0.25, eta=0.75, sigma=0.0001 > delta0=0.2, delta1=0.3, delta2=0.75, delta3=2 > maximum iterations=50, maximum function evaluations=10000 > tolerances: relative=1e-008, absolute=1e-050, solution=1e-008 > total number of linear solver iterations=4 > total number of function evaluations=36 > KSP Object: > type: gmres > GMRES: restart=30, using Classical (unmodified) Gram-Schmidt > Orthogonalization with no iterative refinement > GMRES: happy breakdown tolerance 1e-030 > maximum iterations=10000, initial guess is zero > tolerances: relative=1e-005, absolute=1e-050, divergence=10000 > left preconditioning > PC Object: > type: ilu > ILU: 0 levels of fill > ILU: factor fill ratio allocated 1 > ILU: tolerance for zero pivot 1e-012 > out-of-place factorization > matrix ordering: natural > ILU: factor fill ratio needed 1 > Factored matrix follows > Matrix Object: > type=seqaij, rows=24, cols=24 > total: nonzeros=300, allocated nonzeros=840 > using I-node routines: found 15 nodes, limit used is 5 > linear system matrix = precond matrix: > Matrix Object: > type=seqaij, rows=24, cols=24 > total: nonzeros=300, allocated nonzeros=1920 > using I-node routines: found 15 nodes, limit used is 5 > Nonlinear solve converged due to CONVERGED_TR_DELTA > STEP 0 (Newton iterations: 3) > > WARNING! There are options you set that were not used! > WARNING! could be spelling mistake, etc! > Option left: name:-options_dump no value > > > ----- Original Message ----- > *From:* Matthew Knepley <knepley at gmail.com> > *To:* PETSc users list <petsc-users at mcs.anl.gov> > *Sent:* Tuesday, June 21, 2011 5:46 AM > *Subject:* Re: [petsc-users] MSVC debug runs in a different branch. > > On Mon, Jun 20, 2011 at 9:40 PM, Tian(ICT) <rongtian at ncic.ac.cn> wrote: > >> Hi, >> Good morning/evening! >> I have this issue for quite a while. Using MSVC with Petsc, when debuging >> the code, the FormJacobian was never called, instead calling >> "MatMFFDComputeJacobian(SNES snes,Vec x,Mat *jac,Mat *B,MatStructure >> *flag,void *dummy)", though I donot use matrix free option. >> >> However, when runing the same thing in command line (still a debug >> compilation), I can see FormJacobian was called correctly. >> >> Any hint? >> > > My guess is that you are getting the matrix-free option from somewhere > (maybe a .petscrc file?). You can use -options_dump to see > everything passed in. If not, look at SNESSetFromOptions() in the debugger, > and see why it is choosing MF. > > Matt > > >> Rong >> >> > > > -- > What most experimenters take for granted before they begin their > experiments is infinitely more interesting than any results to which their > experiments lead. > -- Norbert Wiener > > -- What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead. -- Norbert Wiener -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110620/a24e48e4/attachment-0001.htm>
