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