[petsc-users] MSVC debug runs in a different branch.
Tian(ICT)
rongtian at ncic.ac.cn
Mon Jun 20 17:19:27 CDT 2011
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?
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
To: PETSc users list
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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110621/95f10813/attachment.htm>
More information about the petsc-users
mailing list