1Timestep: 0000 time = 0 L_2 Error: 1.0821e-16 20 TS dt 0.1 time 0. 3 0 SNES Function norm 0.38487 4 0 KSP unpreconditioned resid norm 3.848701079585e-01 true resid norm 3.848701079585e-01 ||r(i)||/||b|| 1.000000000000e+00 5 1 KSP unpreconditioned resid norm 6.306958768405e-02 true resid norm 6.306958768405e-02 ||r(i)||/||b|| 1.638723984531e-01 6 2 KSP unpreconditioned resid norm 5.760676204313e-03 true resid norm 5.760676204313e-03 ||r(i)||/||b|| 1.496784521633e-02 7 3 KSP unpreconditioned resid norm 4.612023997603e-04 true resid norm 4.612023997603e-04 ||r(i)||/||b|| 1.198332606829e-03 8 4 KSP unpreconditioned resid norm 7.800377626583e-05 true resid norm 7.800377626586e-05 ||r(i)||/||b|| 2.026755901611e-04 9 5 KSP unpreconditioned resid norm 5.796001277060e-06 true resid norm 5.796001277043e-06 ||r(i)||/||b|| 1.505962961839e-05 10 6 KSP unpreconditioned resid norm 6.294785394574e-07 true resid norm 6.294785394633e-07 ||r(i)||/||b|| 1.635561002132e-06 11 7 KSP unpreconditioned resid norm 4.711800149142e-09 true resid norm 4.711800136142e-09 ||r(i)||/||b|| 1.224257233469e-08 12 8 KSP unpreconditioned resid norm 1.845375931471e-11 true resid norm 1.845373359363e-11 ||r(i)||/||b|| 4.794795234038e-11 13 1 SNES Function norm 1.845e-11 14 Nonlinear solve converged due to CONVERGED_FNORM_RELATIVE iterations 1 15SNES Object: 1 MPI processes 16 type: newtonls 17 maximum iterations=50, maximum function evaluations=10000 18 tolerances: relative=1e-08, absolute=1e-50, solution=1e-08 19 total number of linear solver iterations=8 20 total number of function evaluations=2 21 norm schedule ALWAYS 22 SNESLineSearch Object: 1 MPI processes 23 type: bt 24 interpolation: cubic 25 alpha=1.000000e-04 26 maxstep=1.000000e+08, minlambda=1.000000e-12 27 tolerances: relative=1.000000e-08, absolute=1.000000e-15, lambda=1.000000e-08 28 maximum iterations=40 29 KSP Object: 1 MPI processes 30 type: fgmres 31 restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement 32 happy breakdown tolerance 1e-30 33 maximum iterations=10000, initial guess is zero 34 tolerances: relative=1e-09, absolute=1e-50, divergence=10000. 35 right preconditioning 36 using UNPRECONDITIONED norm type for convergence test 37 PC Object: 1 MPI processes 38 type: mg 39 type is MULTIPLICATIVE, levels=2 cycles=v 40 Cycles per PCApply=1 41 Not using Galerkin computed coarse grid matrices 42 Coarse grid solver -- level ------------------------------- 43 KSP Object: (mg_coarse_) 1 MPI processes 44 type: preonly 45 maximum iterations=10000, initial guess is zero 46 tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 47 left preconditioning 48 using NONE norm type for convergence test 49 PC Object: (mg_coarse_) 1 MPI processes 50 type: lu 51 out-of-place factorization 52 tolerance for zero pivot 2.22045e-14 53 using diagonal shift on blocks to prevent zero pivot [INBLOCKS] 54 matrix ordering: nd 55 factor fill ratio given 5., needed 1. 56 Factored matrix follows: 57 Mat Object: 1 MPI processes 58 type: seqaij 59 rows=1, cols=1 60 package used to perform factorization: petsc 61 total: nonzeros=1, allocated nonzeros=1 62 not using I-node routines 63 linear system matrix = precond matrix: 64 Mat Object: 1 MPI processes 65 type: seqaij 66 rows=1, cols=1 67 total: nonzeros=1, allocated nonzeros=1 68 total number of mallocs used during MatSetValues calls=0 69 not using I-node routines 70 Down solver (pre-smoother) on level 1 ------------------------------- 71 KSP Object: (mg_levels_1_) 1 MPI processes 72 type: chebyshev 73 eigenvalue estimates used: min = 0.151841, max = 1.67025 74 eigenvalues estimate via gmres min 0.994613, max 1.51841 75 eigenvalues estimated using gmres with translations [0. 0.1; 0. 1.1] 76 KSP Object: (mg_levels_1_esteig_) 1 MPI processes 77 type: gmres 78 restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement 79 happy breakdown tolerance 1e-30 80 maximum iterations=10, initial guess is zero 81 tolerances: relative=1e-12, absolute=1e-50, divergence=10000. 82 left preconditioning 83 using PRECONDITIONED norm type for convergence test 84 estimating eigenvalues using noisy right hand side 85 maximum iterations=2, nonzero initial guess 86 tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 87 left preconditioning 88 using NONE norm type for convergence test 89 PC Object: (mg_levels_1_) 1 MPI processes 90 type: sor 91 type = local_symmetric, iterations = 1, local iterations = 1, omega = 1. 92 linear system matrix = precond matrix: 93 Mat Object: 1 MPI processes 94 type: seqaij 95 rows=9, cols=9 96 total: nonzeros=41, allocated nonzeros=41 97 total number of mallocs used during MatSetValues calls=0 98 not using I-node routines 99 Up solver (post-smoother) same as down solver (pre-smoother) 100 linear system matrix = precond matrix: 101 Mat Object: 1 MPI processes 102 type: seqaij 103 rows=9, cols=9 104 total: nonzeros=41, allocated nonzeros=41 105 total number of mallocs used during MatSetValues calls=0 106 not using I-node routines 107Timestep: 0001 time = 0.1 L_2 Error: 0.042609 1081 TS dt 0.1 time 0.1 109 0 SNES Function norm 0.383649 110 0 KSP unpreconditioned resid norm 3.836491353067e-01 true resid norm 3.836491353067e-01 ||r(i)||/||b|| 1.000000000000e+00 111 1 KSP unpreconditioned resid norm 6.087119748495e-02 true resid norm 6.087119748495e-02 ||r(i)||/||b|| 1.586637160965e-01 112 2 KSP unpreconditioned resid norm 5.740698672644e-03 true resid norm 5.740698672644e-03 ||r(i)||/||b|| 1.496340834459e-02 113 3 KSP unpreconditioned resid norm 4.605229592736e-04 true resid norm 4.605229592736e-04 ||r(i)||/||b|| 1.200375334889e-03 114 4 KSP unpreconditioned resid norm 7.848908845461e-05 true resid norm 7.848908845462e-05 ||r(i)||/||b|| 2.045855997874e-04 115 5 KSP unpreconditioned resid norm 5.648624058193e-06 true resid norm 5.648624058220e-06 ||r(i)||/||b|| 1.472341141523e-05 116 6 KSP unpreconditioned resid norm 6.188836190140e-07 true resid norm 6.188836190145e-07 ||r(i)||/||b|| 1.613150042733e-06 117 7 KSP unpreconditioned resid norm 4.800059529470e-09 true resid norm 4.800059540549e-09 ||r(i)||/||b|| 1.251158701742e-08 118 8 KSP unpreconditioned resid norm 1.774392929439e-11 true resid norm 1.774395612425e-11 ||r(i)||/||b|| 4.625047860478e-11 119 1 SNES Function norm 1.774e-11 120 Nonlinear solve converged due to CONVERGED_FNORM_RELATIVE iterations 1 121SNES Object: 1 MPI processes 122 type: newtonls 123 maximum iterations=50, maximum function evaluations=10000 124 tolerances: relative=1e-08, absolute=1e-50, solution=1e-08 125 total number of linear solver iterations=8 126 total number of function evaluations=2 127 norm schedule ALWAYS 128 SNESLineSearch Object: 1 MPI processes 129 type: bt 130 interpolation: cubic 131 alpha=1.000000e-04 132 maxstep=1.000000e+08, minlambda=1.000000e-12 133 tolerances: relative=1.000000e-08, absolute=1.000000e-15, lambda=1.000000e-08 134 maximum iterations=40 135 KSP Object: 1 MPI processes 136 type: fgmres 137 restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement 138 happy breakdown tolerance 1e-30 139 maximum iterations=10000, initial guess is zero 140 tolerances: relative=1e-09, absolute=1e-50, divergence=10000. 141 right preconditioning 142 using UNPRECONDITIONED norm type for convergence test 143 PC Object: 1 MPI processes 144 type: mg 145 type is MULTIPLICATIVE, levels=2 cycles=v 146 Cycles per PCApply=1 147 Not using Galerkin computed coarse grid matrices 148 Coarse grid solver -- level ------------------------------- 149 KSP Object: (mg_coarse_) 1 MPI processes 150 type: preonly 151 maximum iterations=10000, initial guess is zero 152 tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 153 left preconditioning 154 using NONE norm type for convergence test 155 PC Object: (mg_coarse_) 1 MPI processes 156 type: lu 157 out-of-place factorization 158 tolerance for zero pivot 2.22045e-14 159 using diagonal shift on blocks to prevent zero pivot [INBLOCKS] 160 matrix ordering: nd 161 factor fill ratio given 5., needed 1. 162 Factored matrix follows: 163 Mat Object: 1 MPI processes 164 type: seqaij 165 rows=1, cols=1 166 package used to perform factorization: petsc 167 total: nonzeros=1, allocated nonzeros=1 168 not using I-node routines 169 linear system matrix = precond matrix: 170 Mat Object: 1 MPI processes 171 type: seqaij 172 rows=1, cols=1 173 total: nonzeros=1, allocated nonzeros=1 174 total number of mallocs used during MatSetValues calls=0 175 not using I-node routines 176 Down solver (pre-smoother) on level 1 ------------------------------- 177 KSP Object: (mg_levels_1_) 1 MPI processes 178 type: chebyshev 179 eigenvalue estimates used: min = 0.151841, max = 1.67025 180 eigenvalues estimate via gmres min 0.994613, max 1.51841 181 eigenvalues estimated using gmres with translations [0. 0.1; 0. 1.1] 182 KSP Object: (mg_levels_1_esteig_) 1 MPI processes 183 type: gmres 184 restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement 185 happy breakdown tolerance 1e-30 186 maximum iterations=10, initial guess is zero 187 tolerances: relative=1e-12, absolute=1e-50, divergence=10000. 188 left preconditioning 189 using PRECONDITIONED norm type for convergence test 190 estimating eigenvalues using noisy right hand side 191 maximum iterations=2, nonzero initial guess 192 tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 193 left preconditioning 194 using NONE norm type for convergence test 195 PC Object: (mg_levels_1_) 1 MPI processes 196 type: sor 197 type = local_symmetric, iterations = 1, local iterations = 1, omega = 1. 198 linear system matrix = precond matrix: 199 Mat Object: 1 MPI processes 200 type: seqaij 201 rows=9, cols=9 202 total: nonzeros=41, allocated nonzeros=41 203 total number of mallocs used during MatSetValues calls=0 204 not using I-node routines 205 Up solver (post-smoother) same as down solver (pre-smoother) 206 linear system matrix = precond matrix: 207 Mat Object: 1 MPI processes 208 type: seqaij 209 rows=9, cols=9 210 total: nonzeros=41, allocated nonzeros=41 211 total number of mallocs used during MatSetValues calls=0 212 not using I-node routines 213Timestep: 0002 time = 0.2 L_2 Error: 0.084372 2142 TS dt 0.1 time 0.2 215 0 SNES Function norm 0.380922 216 0 KSP unpreconditioned resid norm 3.809222270566e-01 true resid norm 3.809222270566e-01 ||r(i)||/||b|| 1.000000000000e+00 217 1 KSP unpreconditioned resid norm 5.968841254678e-02 true resid norm 5.968841254678e-02 ||r(i)||/||b|| 1.566944859270e-01 218 2 KSP unpreconditioned resid norm 5.800384745162e-03 true resid norm 5.800384745162e-03 ||r(i)||/||b|| 1.522721525069e-02 219 3 KSP unpreconditioned resid norm 4.622462756028e-04 true resid norm 4.622462756028e-04 ||r(i)||/||b|| 1.213492526216e-03 220 4 KSP unpreconditioned resid norm 7.972477625283e-05 true resid norm 7.972477625284e-05 ||r(i)||/||b|| 2.092941041243e-04 221 5 KSP unpreconditioned resid norm 5.651496202770e-06 true resid norm 5.651496202787e-06 ||r(i)||/||b|| 1.483635188856e-05 222 6 KSP unpreconditioned resid norm 6.078212012816e-07 true resid norm 6.078212013076e-07 ||r(i)||/||b|| 1.595656956026e-06 223 7 KSP unpreconditioned resid norm 4.860650279105e-09 true resid norm 4.860650279805e-09 ||r(i)||/||b|| 1.276021700641e-08 224 8 KSP unpreconditioned resid norm 1.731140125257e-11 true resid norm 1.731140182079e-11 ||r(i)||/||b|| 4.544602701334e-11 225 1 SNES Function norm 1.731e-11 226 Nonlinear solve converged due to CONVERGED_FNORM_RELATIVE iterations 1 227SNES Object: 1 MPI processes 228 type: newtonls 229 maximum iterations=50, maximum function evaluations=10000 230 tolerances: relative=1e-08, absolute=1e-50, solution=1e-08 231 total number of linear solver iterations=8 232 total number of function evaluations=2 233 norm schedule ALWAYS 234 SNESLineSearch Object: 1 MPI processes 235 type: bt 236 interpolation: cubic 237 alpha=1.000000e-04 238 maxstep=1.000000e+08, minlambda=1.000000e-12 239 tolerances: relative=1.000000e-08, absolute=1.000000e-15, lambda=1.000000e-08 240 maximum iterations=40 241 KSP Object: 1 MPI processes 242 type: fgmres 243 restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement 244 happy breakdown tolerance 1e-30 245 maximum iterations=10000, initial guess is zero 246 tolerances: relative=1e-09, absolute=1e-50, divergence=10000. 247 right preconditioning 248 using UNPRECONDITIONED norm type for convergence test 249 PC Object: 1 MPI processes 250 type: mg 251 type is MULTIPLICATIVE, levels=2 cycles=v 252 Cycles per PCApply=1 253 Not using Galerkin computed coarse grid matrices 254 Coarse grid solver -- level ------------------------------- 255 KSP Object: (mg_coarse_) 1 MPI processes 256 type: preonly 257 maximum iterations=10000, initial guess is zero 258 tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 259 left preconditioning 260 using NONE norm type for convergence test 261 PC Object: (mg_coarse_) 1 MPI processes 262 type: lu 263 out-of-place factorization 264 tolerance for zero pivot 2.22045e-14 265 using diagonal shift on blocks to prevent zero pivot [INBLOCKS] 266 matrix ordering: nd 267 factor fill ratio given 5., needed 1. 268 Factored matrix follows: 269 Mat Object: 1 MPI processes 270 type: seqaij 271 rows=1, cols=1 272 package used to perform factorization: petsc 273 total: nonzeros=1, allocated nonzeros=1 274 not using I-node routines 275 linear system matrix = precond matrix: 276 Mat Object: 1 MPI processes 277 type: seqaij 278 rows=1, cols=1 279 total: nonzeros=1, allocated nonzeros=1 280 total number of mallocs used during MatSetValues calls=0 281 not using I-node routines 282 Down solver (pre-smoother) on level 1 ------------------------------- 283 KSP Object: (mg_levels_1_) 1 MPI processes 284 type: chebyshev 285 eigenvalue estimates used: min = 0.151841, max = 1.67025 286 eigenvalues estimate via gmres min 0.994613, max 1.51841 287 eigenvalues estimated using gmres with translations [0. 0.1; 0. 1.1] 288 KSP Object: (mg_levels_1_esteig_) 1 MPI processes 289 type: gmres 290 restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement 291 happy breakdown tolerance 1e-30 292 maximum iterations=10, initial guess is zero 293 tolerances: relative=1e-12, absolute=1e-50, divergence=10000. 294 left preconditioning 295 using PRECONDITIONED norm type for convergence test 296 estimating eigenvalues using noisy right hand side 297 maximum iterations=2, nonzero initial guess 298 tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 299 left preconditioning 300 using NONE norm type for convergence test 301 PC Object: (mg_levels_1_) 1 MPI processes 302 type: sor 303 type = local_symmetric, iterations = 1, local iterations = 1, omega = 1. 304 linear system matrix = precond matrix: 305 Mat Object: 1 MPI processes 306 type: seqaij 307 rows=9, cols=9 308 total: nonzeros=41, allocated nonzeros=41 309 total number of mallocs used during MatSetValues calls=0 310 not using I-node routines 311 Up solver (post-smoother) same as down solver (pre-smoother) 312 linear system matrix = precond matrix: 313 Mat Object: 1 MPI processes 314 type: seqaij 315 rows=9, cols=9 316 total: nonzeros=41, allocated nonzeros=41 317 total number of mallocs used during MatSetValues calls=0 318 not using I-node routines 319Timestep: 0003 time = 0.3 L_2 Error: 0.12562 3203 TS dt 0.1 time 0.3 321 0 SNES Function norm 0.377457 322 0 KSP unpreconditioned resid norm 3.774573494154e-01 true resid norm 3.774573494154e-01 ||r(i)||/||b|| 1.000000000000e+00 323 1 KSP unpreconditioned resid norm 5.954853714304e-02 true resid norm 5.954853714304e-02 ||r(i)||/||b|| 1.577622935022e-01 324 2 KSP unpreconditioned resid norm 5.848939409168e-03 true resid norm 5.848939409168e-03 ||r(i)||/||b|| 1.549562995190e-02 325 3 KSP unpreconditioned resid norm 4.630907205003e-04 true resid norm 4.630907205002e-04 ||r(i)||/||b|| 1.226868999153e-03 326 4 KSP unpreconditioned resid norm 8.053803471345e-05 true resid norm 8.053803471342e-05 ||r(i)||/||b|| 2.133698942096e-04 327 5 KSP unpreconditioned resid norm 5.701869061872e-06 true resid norm 5.701869061870e-06 ||r(i)||/||b|| 1.510599560639e-05 328 6 KSP unpreconditioned resid norm 6.025388075754e-07 true resid norm 6.025388075944e-07 ||r(i)||/||b|| 1.596309645388e-06 329 7 KSP unpreconditioned resid norm 4.883261060044e-09 true resid norm 4.883261063505e-09 ||r(i)||/||b|| 1.293725256924e-08 330 8 KSP unpreconditioned resid norm 1.729375636379e-11 true resid norm 1.729375104729e-11 ||r(i)||/||b|| 4.581643747055e-11 331 1 SNES Function norm 1.729e-11 332 Nonlinear solve converged due to CONVERGED_FNORM_RELATIVE iterations 1 333SNES Object: 1 MPI processes 334 type: newtonls 335 maximum iterations=50, maximum function evaluations=10000 336 tolerances: relative=1e-08, absolute=1e-50, solution=1e-08 337 total number of linear solver iterations=8 338 total number of function evaluations=2 339 norm schedule ALWAYS 340 SNESLineSearch Object: 1 MPI processes 341 type: bt 342 interpolation: cubic 343 alpha=1.000000e-04 344 maxstep=1.000000e+08, minlambda=1.000000e-12 345 tolerances: relative=1.000000e-08, absolute=1.000000e-15, lambda=1.000000e-08 346 maximum iterations=40 347 KSP Object: 1 MPI processes 348 type: fgmres 349 restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement 350 happy breakdown tolerance 1e-30 351 maximum iterations=10000, initial guess is zero 352 tolerances: relative=1e-09, absolute=1e-50, divergence=10000. 353 right preconditioning 354 using UNPRECONDITIONED norm type for convergence test 355 PC Object: 1 MPI processes 356 type: mg 357 type is MULTIPLICATIVE, levels=2 cycles=v 358 Cycles per PCApply=1 359 Not using Galerkin computed coarse grid matrices 360 Coarse grid solver -- level ------------------------------- 361 KSP Object: (mg_coarse_) 1 MPI processes 362 type: preonly 363 maximum iterations=10000, initial guess is zero 364 tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 365 left preconditioning 366 using NONE norm type for convergence test 367 PC Object: (mg_coarse_) 1 MPI processes 368 type: lu 369 out-of-place factorization 370 tolerance for zero pivot 2.22045e-14 371 using diagonal shift on blocks to prevent zero pivot [INBLOCKS] 372 matrix ordering: nd 373 factor fill ratio given 5., needed 1. 374 Factored matrix follows: 375 Mat Object: 1 MPI processes 376 type: seqaij 377 rows=1, cols=1 378 package used to perform factorization: petsc 379 total: nonzeros=1, allocated nonzeros=1 380 not using I-node routines 381 linear system matrix = precond matrix: 382 Mat Object: 1 MPI processes 383 type: seqaij 384 rows=1, cols=1 385 total: nonzeros=1, allocated nonzeros=1 386 total number of mallocs used during MatSetValues calls=0 387 not using I-node routines 388 Down solver (pre-smoother) on level 1 ------------------------------- 389 KSP Object: (mg_levels_1_) 1 MPI processes 390 type: chebyshev 391 eigenvalue estimates used: min = 0.151841, max = 1.67025 392 eigenvalues estimate via gmres min 0.994613, max 1.51841 393 eigenvalues estimated using gmres with translations [0. 0.1; 0. 1.1] 394 KSP Object: (mg_levels_1_esteig_) 1 MPI processes 395 type: gmres 396 restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement 397 happy breakdown tolerance 1e-30 398 maximum iterations=10, initial guess is zero 399 tolerances: relative=1e-12, absolute=1e-50, divergence=10000. 400 left preconditioning 401 using PRECONDITIONED norm type for convergence test 402 estimating eigenvalues using noisy right hand side 403 maximum iterations=2, nonzero initial guess 404 tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 405 left preconditioning 406 using NONE norm type for convergence test 407 PC Object: (mg_levels_1_) 1 MPI processes 408 type: sor 409 type = local_symmetric, iterations = 1, local iterations = 1, omega = 1. 410 linear system matrix = precond matrix: 411 Mat Object: 1 MPI processes 412 type: seqaij 413 rows=9, cols=9 414 total: nonzeros=41, allocated nonzeros=41 415 total number of mallocs used during MatSetValues calls=0 416 not using I-node routines 417 Up solver (post-smoother) same as down solver (pre-smoother) 418 linear system matrix = precond matrix: 419 Mat Object: 1 MPI processes 420 type: seqaij 421 rows=9, cols=9 422 total: nonzeros=41, allocated nonzeros=41 423 total number of mallocs used during MatSetValues calls=0 424 not using I-node routines 425Timestep: 0004 time = 0.4 L_2 Error: 0.16664 4264 TS dt 0.1 time 0.4 427 0 SNES Function norm 0.37449 428 0 KSP unpreconditioned resid norm 3.744902213650e-01 true resid norm 3.744902213650e-01 ||r(i)||/||b|| 1.000000000000e+00 429 1 KSP unpreconditioned resid norm 5.982974938131e-02 true resid norm 5.982974938131e-02 ||r(i)||/||b|| 1.597631819684e-01 430 2 KSP unpreconditioned resid norm 5.840834302842e-03 true resid norm 5.840834302842e-03 ||r(i)||/||b|| 1.559676052836e-02 431 3 KSP unpreconditioned resid norm 4.612964201615e-04 true resid norm 4.612964201614e-04 ||r(i)||/||b|| 1.231798305654e-03 432 4 KSP unpreconditioned resid norm 8.073771822390e-05 true resid norm 8.073771822391e-05 ||r(i)||/||b|| 2.155936620444e-04 433 5 KSP unpreconditioned resid norm 5.730193368363e-06 true resid norm 5.730193368351e-06 ||r(i)||/||b|| 1.530131640678e-05 434 6 KSP unpreconditioned resid norm 6.010616249637e-07 true resid norm 6.010616249715e-07 ||r(i)||/||b|| 1.605012869978e-06 435 7 KSP unpreconditioned resid norm 4.858498871957e-09 true resid norm 4.858498872870e-09 ||r(i)||/||b|| 1.297363347743e-08 436 8 KSP unpreconditioned resid norm 1.751406135113e-11 true resid norm 1.751403567040e-11 ||r(i)||/||b|| 4.676767154710e-11 437 1 SNES Function norm 1.751e-11 438 Nonlinear solve converged due to CONVERGED_FNORM_RELATIVE iterations 1 439SNES Object: 1 MPI processes 440 type: newtonls 441 maximum iterations=50, maximum function evaluations=10000 442 tolerances: relative=1e-08, absolute=1e-50, solution=1e-08 443 total number of linear solver iterations=8 444 total number of function evaluations=2 445 norm schedule ALWAYS 446 SNESLineSearch Object: 1 MPI processes 447 type: bt 448 interpolation: cubic 449 alpha=1.000000e-04 450 maxstep=1.000000e+08, minlambda=1.000000e-12 451 tolerances: relative=1.000000e-08, absolute=1.000000e-15, lambda=1.000000e-08 452 maximum iterations=40 453 KSP Object: 1 MPI processes 454 type: fgmres 455 restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement 456 happy breakdown tolerance 1e-30 457 maximum iterations=10000, initial guess is zero 458 tolerances: relative=1e-09, absolute=1e-50, divergence=10000. 459 right preconditioning 460 using UNPRECONDITIONED norm type for convergence test 461 PC Object: 1 MPI processes 462 type: mg 463 type is MULTIPLICATIVE, levels=2 cycles=v 464 Cycles per PCApply=1 465 Not using Galerkin computed coarse grid matrices 466 Coarse grid solver -- level ------------------------------- 467 KSP Object: (mg_coarse_) 1 MPI processes 468 type: preonly 469 maximum iterations=10000, initial guess is zero 470 tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 471 left preconditioning 472 using NONE norm type for convergence test 473 PC Object: (mg_coarse_) 1 MPI processes 474 type: lu 475 out-of-place factorization 476 tolerance for zero pivot 2.22045e-14 477 using diagonal shift on blocks to prevent zero pivot [INBLOCKS] 478 matrix ordering: nd 479 factor fill ratio given 5., needed 1. 480 Factored matrix follows: 481 Mat Object: 1 MPI processes 482 type: seqaij 483 rows=1, cols=1 484 package used to perform factorization: petsc 485 total: nonzeros=1, allocated nonzeros=1 486 not using I-node routines 487 linear system matrix = precond matrix: 488 Mat Object: 1 MPI processes 489 type: seqaij 490 rows=1, cols=1 491 total: nonzeros=1, allocated nonzeros=1 492 total number of mallocs used during MatSetValues calls=0 493 not using I-node routines 494 Down solver (pre-smoother) on level 1 ------------------------------- 495 KSP Object: (mg_levels_1_) 1 MPI processes 496 type: chebyshev 497 eigenvalue estimates used: min = 0.151841, max = 1.67025 498 eigenvalues estimate via gmres min 0.994613, max 1.51841 499 eigenvalues estimated using gmres with translations [0. 0.1; 0. 1.1] 500 KSP Object: (mg_levels_1_esteig_) 1 MPI processes 501 type: gmres 502 restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement 503 happy breakdown tolerance 1e-30 504 maximum iterations=10, initial guess is zero 505 tolerances: relative=1e-12, absolute=1e-50, divergence=10000. 506 left preconditioning 507 using PRECONDITIONED norm type for convergence test 508 estimating eigenvalues using noisy right hand side 509 maximum iterations=2, nonzero initial guess 510 tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 511 left preconditioning 512 using NONE norm type for convergence test 513 PC Object: (mg_levels_1_) 1 MPI processes 514 type: sor 515 type = local_symmetric, iterations = 1, local iterations = 1, omega = 1. 516 linear system matrix = precond matrix: 517 Mat Object: 1 MPI processes 518 type: seqaij 519 rows=9, cols=9 520 total: nonzeros=41, allocated nonzeros=41 521 total number of mallocs used during MatSetValues calls=0 522 not using I-node routines 523 Up solver (post-smoother) same as down solver (pre-smoother) 524 linear system matrix = precond matrix: 525 Mat Object: 1 MPI processes 526 type: seqaij 527 rows=9, cols=9 528 total: nonzeros=41, allocated nonzeros=41 529 total number of mallocs used during MatSetValues calls=0 530 not using I-node routines 531Timestep: 0005 time = 0.5 L_2 Error: 0.2076 5325 TS dt 0.1 time 0.5 533 0 SNES Function norm 0.372928 534 0 KSP unpreconditioned resid norm 3.729281667988e-01 true resid norm 3.729281667988e-01 ||r(i)||/||b|| 1.000000000000e+00 535 1 KSP unpreconditioned resid norm 6.002541949214e-02 true resid norm 6.002541949214e-02 ||r(i)||/||b|| 1.609570551010e-01 536 2 KSP unpreconditioned resid norm 5.775316419172e-03 true resid norm 5.775316419172e-03 ||r(i)||/||b|| 1.548640444284e-02 537 3 KSP unpreconditioned resid norm 4.576993940278e-04 true resid norm 4.576993940278e-04 ||r(i)||/||b|| 1.227312482070e-03 538 4 KSP unpreconditioned resid norm 8.039211062475e-05 true resid norm 8.039211062473e-05 ||r(i)||/||b|| 2.155699616760e-04 539 5 KSP unpreconditioned resid norm 5.708171590220e-06 true resid norm 5.708171590238e-06 ||r(i)||/||b|| 1.530635682265e-05 540 6 KSP unpreconditioned resid norm 5.991553964591e-07 true resid norm 5.991553964691e-07 ||r(i)||/||b|| 1.606624143229e-06 541 7 KSP unpreconditioned resid norm 4.791537579858e-09 true resid norm 4.791537591102e-09 ||r(i)||/||b|| 1.284841966278e-08 542 8 KSP unpreconditioned resid norm 1.772048749492e-11 true resid norm 1.772051502236e-11 ||r(i)||/||b|| 4.751723414853e-11 543 1 SNES Function norm 1.772e-11 544 Nonlinear solve converged due to CONVERGED_FNORM_RELATIVE iterations 1 545SNES Object: 1 MPI processes 546 type: newtonls 547 maximum iterations=50, maximum function evaluations=10000 548 tolerances: relative=1e-08, absolute=1e-50, solution=1e-08 549 total number of linear solver iterations=8 550 total number of function evaluations=2 551 norm schedule ALWAYS 552 SNESLineSearch Object: 1 MPI processes 553 type: bt 554 interpolation: cubic 555 alpha=1.000000e-04 556 maxstep=1.000000e+08, minlambda=1.000000e-12 557 tolerances: relative=1.000000e-08, absolute=1.000000e-15, lambda=1.000000e-08 558 maximum iterations=40 559 KSP Object: 1 MPI processes 560 type: fgmres 561 restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement 562 happy breakdown tolerance 1e-30 563 maximum iterations=10000, initial guess is zero 564 tolerances: relative=1e-09, absolute=1e-50, divergence=10000. 565 right preconditioning 566 using UNPRECONDITIONED norm type for convergence test 567 PC Object: 1 MPI processes 568 type: mg 569 type is MULTIPLICATIVE, levels=2 cycles=v 570 Cycles per PCApply=1 571 Not using Galerkin computed coarse grid matrices 572 Coarse grid solver -- level ------------------------------- 573 KSP Object: (mg_coarse_) 1 MPI processes 574 type: preonly 575 maximum iterations=10000, initial guess is zero 576 tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 577 left preconditioning 578 using NONE norm type for convergence test 579 PC Object: (mg_coarse_) 1 MPI processes 580 type: lu 581 out-of-place factorization 582 tolerance for zero pivot 2.22045e-14 583 using diagonal shift on blocks to prevent zero pivot [INBLOCKS] 584 matrix ordering: nd 585 factor fill ratio given 5., needed 1. 586 Factored matrix follows: 587 Mat Object: 1 MPI processes 588 type: seqaij 589 rows=1, cols=1 590 package used to perform factorization: petsc 591 total: nonzeros=1, allocated nonzeros=1 592 not using I-node routines 593 linear system matrix = precond matrix: 594 Mat Object: 1 MPI processes 595 type: seqaij 596 rows=1, cols=1 597 total: nonzeros=1, allocated nonzeros=1 598 total number of mallocs used during MatSetValues calls=0 599 not using I-node routines 600 Down solver (pre-smoother) on level 1 ------------------------------- 601 KSP Object: (mg_levels_1_) 1 MPI processes 602 type: chebyshev 603 eigenvalue estimates used: min = 0.151841, max = 1.67025 604 eigenvalues estimate via gmres min 0.994613, max 1.51841 605 eigenvalues estimated using gmres with translations [0. 0.1; 0. 1.1] 606 KSP Object: (mg_levels_1_esteig_) 1 MPI processes 607 type: gmres 608 restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement 609 happy breakdown tolerance 1e-30 610 maximum iterations=10, initial guess is zero 611 tolerances: relative=1e-12, absolute=1e-50, divergence=10000. 612 left preconditioning 613 using PRECONDITIONED norm type for convergence test 614 estimating eigenvalues using noisy right hand side 615 maximum iterations=2, nonzero initial guess 616 tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 617 left preconditioning 618 using NONE norm type for convergence test 619 PC Object: (mg_levels_1_) 1 MPI processes 620 type: sor 621 type = local_symmetric, iterations = 1, local iterations = 1, omega = 1. 622 linear system matrix = precond matrix: 623 Mat Object: 1 MPI processes 624 type: seqaij 625 rows=9, cols=9 626 total: nonzeros=41, allocated nonzeros=41 627 total number of mallocs used during MatSetValues calls=0 628 not using I-node routines 629 Up solver (post-smoother) same as down solver (pre-smoother) 630 linear system matrix = precond matrix: 631 Mat Object: 1 MPI processes 632 type: seqaij 633 rows=9, cols=9 634 total: nonzeros=41, allocated nonzeros=41 635 total number of mallocs used during MatSetValues calls=0 636 not using I-node routines 637Timestep: 0006 time = 0.6 L_2 Error: 0.24859 6386 TS dt 0.1 time 0.6 639 0 SNES Function norm 0.372986 640 0 KSP unpreconditioned resid norm 3.729860798745e-01 true resid norm 3.729860798745e-01 ||r(i)||/||b|| 1.000000000000e+00 641 1 KSP unpreconditioned resid norm 5.998286066860e-02 true resid norm 5.998286066860e-02 ||r(i)||/||b|| 1.608179605222e-01 642 2 KSP unpreconditioned resid norm 5.687596001641e-03 true resid norm 5.687596001641e-03 ||r(i)||/||b|| 1.524881572941e-02 643 3 KSP unpreconditioned resid norm 4.534382695579e-04 true resid norm 4.534382695579e-04 ||r(i)||/||b|| 1.215697566275e-03 644 4 KSP unpreconditioned resid norm 7.976964452116e-05 true resid norm 7.976964452115e-05 ||r(i)||/||b|| 2.138676181910e-04 645 5 KSP unpreconditioned resid norm 5.651631937849e-06 true resid norm 5.651631937845e-06 ||r(i)||/||b|| 1.515239372940e-05 646 6 KSP unpreconditioned resid norm 5.948808411325e-07 true resid norm 5.948808411183e-07 ||r(i)||/||b|| 1.594914323126e-06 647 7 KSP unpreconditioned resid norm 4.702487381111e-09 true resid norm 4.702487348039e-09 ||r(i)||/||b|| 1.260767519695e-08 648 8 KSP unpreconditioned resid norm 1.779657220875e-11 true resid norm 1.779659688198e-11 ||r(i)||/||b|| 4.771383663424e-11 649 1 SNES Function norm 1.780e-11 650 Nonlinear solve converged due to CONVERGED_FNORM_RELATIVE iterations 1 651SNES Object: 1 MPI processes 652 type: newtonls 653 maximum iterations=50, maximum function evaluations=10000 654 tolerances: relative=1e-08, absolute=1e-50, solution=1e-08 655 total number of linear solver iterations=8 656 total number of function evaluations=2 657 norm schedule ALWAYS 658 SNESLineSearch Object: 1 MPI processes 659 type: bt 660 interpolation: cubic 661 alpha=1.000000e-04 662 maxstep=1.000000e+08, minlambda=1.000000e-12 663 tolerances: relative=1.000000e-08, absolute=1.000000e-15, lambda=1.000000e-08 664 maximum iterations=40 665 KSP Object: 1 MPI processes 666 type: fgmres 667 restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement 668 happy breakdown tolerance 1e-30 669 maximum iterations=10000, initial guess is zero 670 tolerances: relative=1e-09, absolute=1e-50, divergence=10000. 671 right preconditioning 672 using UNPRECONDITIONED norm type for convergence test 673 PC Object: 1 MPI processes 674 type: mg 675 type is MULTIPLICATIVE, levels=2 cycles=v 676 Cycles per PCApply=1 677 Not using Galerkin computed coarse grid matrices 678 Coarse grid solver -- level ------------------------------- 679 KSP Object: (mg_coarse_) 1 MPI processes 680 type: preonly 681 maximum iterations=10000, initial guess is zero 682 tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 683 left preconditioning 684 using NONE norm type for convergence test 685 PC Object: (mg_coarse_) 1 MPI processes 686 type: lu 687 out-of-place factorization 688 tolerance for zero pivot 2.22045e-14 689 using diagonal shift on blocks to prevent zero pivot [INBLOCKS] 690 matrix ordering: nd 691 factor fill ratio given 5., needed 1. 692 Factored matrix follows: 693 Mat Object: 1 MPI processes 694 type: seqaij 695 rows=1, cols=1 696 package used to perform factorization: petsc 697 total: nonzeros=1, allocated nonzeros=1 698 not using I-node routines 699 linear system matrix = precond matrix: 700 Mat Object: 1 MPI processes 701 type: seqaij 702 rows=1, cols=1 703 total: nonzeros=1, allocated nonzeros=1 704 total number of mallocs used during MatSetValues calls=0 705 not using I-node routines 706 Down solver (pre-smoother) on level 1 ------------------------------- 707 KSP Object: (mg_levels_1_) 1 MPI processes 708 type: chebyshev 709 eigenvalue estimates used: min = 0.151841, max = 1.67025 710 eigenvalues estimate via gmres min 0.994613, max 1.51841 711 eigenvalues estimated using gmres with translations [0. 0.1; 0. 1.1] 712 KSP Object: (mg_levels_1_esteig_) 1 MPI processes 713 type: gmres 714 restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement 715 happy breakdown tolerance 1e-30 716 maximum iterations=10, initial guess is zero 717 tolerances: relative=1e-12, absolute=1e-50, divergence=10000. 718 left preconditioning 719 using PRECONDITIONED norm type for convergence test 720 estimating eigenvalues using noisy right hand side 721 maximum iterations=2, nonzero initial guess 722 tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 723 left preconditioning 724 using NONE norm type for convergence test 725 PC Object: (mg_levels_1_) 1 MPI processes 726 type: sor 727 type = local_symmetric, iterations = 1, local iterations = 1, omega = 1. 728 linear system matrix = precond matrix: 729 Mat Object: 1 MPI processes 730 type: seqaij 731 rows=9, cols=9 732 total: nonzeros=41, allocated nonzeros=41 733 total number of mallocs used during MatSetValues calls=0 734 not using I-node routines 735 Up solver (post-smoother) same as down solver (pre-smoother) 736 linear system matrix = precond matrix: 737 Mat Object: 1 MPI processes 738 type: seqaij 739 rows=9, cols=9 740 total: nonzeros=41, allocated nonzeros=41 741 total number of mallocs used during MatSetValues calls=0 742 not using I-node routines 743Timestep: 0007 time = 0.7 L_2 Error: 0.28962 7447 TS dt 0.1 time 0.7 745 0 SNES Function norm 0.374276 746 0 KSP unpreconditioned resid norm 3.742759385298e-01 true resid norm 3.742759385298e-01 ||r(i)||/||b|| 1.000000000000e+00 747 1 KSP unpreconditioned resid norm 5.979030900298e-02 true resid norm 5.979030900298e-02 ||r(i)||/||b|| 1.597492727901e-01 748 2 KSP unpreconditioned resid norm 5.617806615202e-03 true resid norm 5.617806615202e-03 ||r(i)||/||b|| 1.500979901959e-02 749 3 KSP unpreconditioned resid norm 4.495138691066e-04 true resid norm 4.495138691066e-04 ||r(i)||/||b|| 1.201022622166e-03 750 4 KSP unpreconditioned resid norm 7.920536127941e-05 true resid norm 7.920536127941e-05 ||r(i)||/||b|| 2.116229047225e-04 751 5 KSP unpreconditioned resid norm 5.595911671793e-06 true resid norm 5.595911671801e-06 ||r(i)||/||b|| 1.495129955129e-05 752 6 KSP unpreconditioned resid norm 5.892156061787e-07 true resid norm 5.892156061538e-07 ||r(i)||/||b|| 1.574281286872e-06 753 7 KSP unpreconditioned resid norm 4.617371734791e-09 true resid norm 4.617371702498e-09 ||r(i)||/||b|| 1.233681149965e-08 754 8 KSP unpreconditioned resid norm 1.776419283029e-11 true resid norm 1.776421370272e-11 ||r(i)||/||b|| 4.746287931973e-11 755 1 SNES Function norm 1.776e-11 756 Nonlinear solve converged due to CONVERGED_FNORM_RELATIVE iterations 1 757SNES Object: 1 MPI processes 758 type: newtonls 759 maximum iterations=50, maximum function evaluations=10000 760 tolerances: relative=1e-08, absolute=1e-50, solution=1e-08 761 total number of linear solver iterations=8 762 total number of function evaluations=2 763 norm schedule ALWAYS 764 SNESLineSearch Object: 1 MPI processes 765 type: bt 766 interpolation: cubic 767 alpha=1.000000e-04 768 maxstep=1.000000e+08, minlambda=1.000000e-12 769 tolerances: relative=1.000000e-08, absolute=1.000000e-15, lambda=1.000000e-08 770 maximum iterations=40 771 KSP Object: 1 MPI processes 772 type: fgmres 773 restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement 774 happy breakdown tolerance 1e-30 775 maximum iterations=10000, initial guess is zero 776 tolerances: relative=1e-09, absolute=1e-50, divergence=10000. 777 right preconditioning 778 using UNPRECONDITIONED norm type for convergence test 779 PC Object: 1 MPI processes 780 type: mg 781 type is MULTIPLICATIVE, levels=2 cycles=v 782 Cycles per PCApply=1 783 Not using Galerkin computed coarse grid matrices 784 Coarse grid solver -- level ------------------------------- 785 KSP Object: (mg_coarse_) 1 MPI processes 786 type: preonly 787 maximum iterations=10000, initial guess is zero 788 tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 789 left preconditioning 790 using NONE norm type for convergence test 791 PC Object: (mg_coarse_) 1 MPI processes 792 type: lu 793 out-of-place factorization 794 tolerance for zero pivot 2.22045e-14 795 using diagonal shift on blocks to prevent zero pivot [INBLOCKS] 796 matrix ordering: nd 797 factor fill ratio given 5., needed 1. 798 Factored matrix follows: 799 Mat Object: 1 MPI processes 800 type: seqaij 801 rows=1, cols=1 802 package used to perform factorization: petsc 803 total: nonzeros=1, allocated nonzeros=1 804 not using I-node routines 805 linear system matrix = precond matrix: 806 Mat Object: 1 MPI processes 807 type: seqaij 808 rows=1, cols=1 809 total: nonzeros=1, allocated nonzeros=1 810 total number of mallocs used during MatSetValues calls=0 811 not using I-node routines 812 Down solver (pre-smoother) on level 1 ------------------------------- 813 KSP Object: (mg_levels_1_) 1 MPI processes 814 type: chebyshev 815 eigenvalue estimates used: min = 0.151841, max = 1.67025 816 eigenvalues estimate via gmres min 0.994613, max 1.51841 817 eigenvalues estimated using gmres with translations [0. 0.1; 0. 1.1] 818 KSP Object: (mg_levels_1_esteig_) 1 MPI processes 819 type: gmres 820 restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement 821 happy breakdown tolerance 1e-30 822 maximum iterations=10, initial guess is zero 823 tolerances: relative=1e-12, absolute=1e-50, divergence=10000. 824 left preconditioning 825 using PRECONDITIONED norm type for convergence test 826 estimating eigenvalues using noisy right hand side 827 maximum iterations=2, nonzero initial guess 828 tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 829 left preconditioning 830 using NONE norm type for convergence test 831 PC Object: (mg_levels_1_) 1 MPI processes 832 type: sor 833 type = local_symmetric, iterations = 1, local iterations = 1, omega = 1. 834 linear system matrix = precond matrix: 835 Mat Object: 1 MPI processes 836 type: seqaij 837 rows=9, cols=9 838 total: nonzeros=41, allocated nonzeros=41 839 total number of mallocs used during MatSetValues calls=0 840 not using I-node routines 841 Up solver (post-smoother) same as down solver (pre-smoother) 842 linear system matrix = precond matrix: 843 Mat Object: 1 MPI processes 844 type: seqaij 845 rows=9, cols=9 846 total: nonzeros=41, allocated nonzeros=41 847 total number of mallocs used during MatSetValues calls=0 848 not using I-node routines 849Timestep: 0008 time = 0.8 L_2 Error: 0.33068 8508 TS dt 0.1 time 0.8 851 0 SNES Function norm 0.376133 852 0 KSP unpreconditioned resid norm 3.761327824878e-01 true resid norm 3.761327824878e-01 ||r(i)||/||b|| 1.000000000000e+00 853 1 KSP unpreconditioned resid norm 5.960286328784e-02 true resid norm 5.960286328784e-02 ||r(i)||/||b|| 1.584622932721e-01 854 2 KSP unpreconditioned resid norm 5.587587079708e-03 true resid norm 5.587587079708e-03 ||r(i)||/||b|| 1.485535784132e-02 855 3 KSP unpreconditioned resid norm 4.467633715519e-04 true resid norm 4.467633715519e-04 ||r(i)||/||b|| 1.187780997437e-03 856 4 KSP unpreconditioned resid norm 7.887710489546e-05 true resid norm 7.887710489544e-05 ||r(i)||/||b|| 2.097054778734e-04 857 5 KSP unpreconditioned resid norm 5.566404901123e-06 true resid norm 5.566404901117e-06 ||r(i)||/||b|| 1.479904214756e-05 858 6 KSP unpreconditioned resid norm 5.843952130956e-07 true resid norm 5.843952130851e-07 ||r(i)||/||b|| 1.553693908890e-06 859 7 KSP unpreconditioned resid norm 4.557400918821e-09 true resid norm 4.557400934472e-09 ||r(i)||/||b|| 1.211646829699e-08 860 8 KSP unpreconditioned resid norm 1.769382827441e-11 true resid norm 1.769385166385e-11 ||r(i)||/||b|| 4.704150365949e-11 861 1 SNES Function norm 1.769e-11 862 Nonlinear solve converged due to CONVERGED_FNORM_RELATIVE iterations 1 863SNES Object: 1 MPI processes 864 type: newtonls 865 maximum iterations=50, maximum function evaluations=10000 866 tolerances: relative=1e-08, absolute=1e-50, solution=1e-08 867 total number of linear solver iterations=8 868 total number of function evaluations=2 869 norm schedule ALWAYS 870 SNESLineSearch Object: 1 MPI processes 871 type: bt 872 interpolation: cubic 873 alpha=1.000000e-04 874 maxstep=1.000000e+08, minlambda=1.000000e-12 875 tolerances: relative=1.000000e-08, absolute=1.000000e-15, lambda=1.000000e-08 876 maximum iterations=40 877 KSP Object: 1 MPI processes 878 type: fgmres 879 restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement 880 happy breakdown tolerance 1e-30 881 maximum iterations=10000, initial guess is zero 882 tolerances: relative=1e-09, absolute=1e-50, divergence=10000. 883 right preconditioning 884 using UNPRECONDITIONED norm type for convergence test 885 PC Object: 1 MPI processes 886 type: mg 887 type is MULTIPLICATIVE, levels=2 cycles=v 888 Cycles per PCApply=1 889 Not using Galerkin computed coarse grid matrices 890 Coarse grid solver -- level ------------------------------- 891 KSP Object: (mg_coarse_) 1 MPI processes 892 type: preonly 893 maximum iterations=10000, initial guess is zero 894 tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 895 left preconditioning 896 using NONE norm type for convergence test 897 PC Object: (mg_coarse_) 1 MPI processes 898 type: lu 899 out-of-place factorization 900 tolerance for zero pivot 2.22045e-14 901 using diagonal shift on blocks to prevent zero pivot [INBLOCKS] 902 matrix ordering: nd 903 factor fill ratio given 5., needed 1. 904 Factored matrix follows: 905 Mat Object: 1 MPI processes 906 type: seqaij 907 rows=1, cols=1 908 package used to perform factorization: petsc 909 total: nonzeros=1, allocated nonzeros=1 910 not using I-node routines 911 linear system matrix = precond matrix: 912 Mat Object: 1 MPI processes 913 type: seqaij 914 rows=1, cols=1 915 total: nonzeros=1, allocated nonzeros=1 916 total number of mallocs used during MatSetValues calls=0 917 not using I-node routines 918 Down solver (pre-smoother) on level 1 ------------------------------- 919 KSP Object: (mg_levels_1_) 1 MPI processes 920 type: chebyshev 921 eigenvalue estimates used: min = 0.151841, max = 1.67025 922 eigenvalues estimate via gmres min 0.994613, max 1.51841 923 eigenvalues estimated using gmres with translations [0. 0.1; 0. 1.1] 924 KSP Object: (mg_levels_1_esteig_) 1 MPI processes 925 type: gmres 926 restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement 927 happy breakdown tolerance 1e-30 928 maximum iterations=10, initial guess is zero 929 tolerances: relative=1e-12, absolute=1e-50, divergence=10000. 930 left preconditioning 931 using PRECONDITIONED norm type for convergence test 932 estimating eigenvalues using noisy right hand side 933 maximum iterations=2, nonzero initial guess 934 tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 935 left preconditioning 936 using NONE norm type for convergence test 937 PC Object: (mg_levels_1_) 1 MPI processes 938 type: sor 939 type = local_symmetric, iterations = 1, local iterations = 1, omega = 1. 940 linear system matrix = precond matrix: 941 Mat Object: 1 MPI processes 942 type: seqaij 943 rows=9, cols=9 944 total: nonzeros=41, allocated nonzeros=41 945 total number of mallocs used during MatSetValues calls=0 946 not using I-node routines 947 Up solver (post-smoother) same as down solver (pre-smoother) 948 linear system matrix = precond matrix: 949 Mat Object: 1 MPI processes 950 type: seqaij 951 rows=9, cols=9 952 total: nonzeros=41, allocated nonzeros=41 953 total number of mallocs used during MatSetValues calls=0 954 not using I-node routines 955Timestep: 0009 time = 0.9 L_2 Error: 0.37179 9569 TS dt 0.1 time 0.9 957 0 SNES Function norm 0.377927 958 0 KSP unpreconditioned resid norm 3.779274613460e-01 true resid norm 3.779274613460e-01 ||r(i)||/||b|| 1.000000000000e+00 959 1 KSP unpreconditioned resid norm 5.952775794612e-02 true resid norm 5.952775794612e-02 ||r(i)||/||b|| 1.575110676904e-01 960 2 KSP unpreconditioned resid norm 5.596783441033e-03 true resid norm 5.596783441033e-03 ||r(i)||/||b|| 1.480914729271e-02 961 3 KSP unpreconditioned resid norm 4.456392946867e-04 true resid norm 4.456392946866e-04 ||r(i)||/||b|| 1.179166216447e-03 962 4 KSP unpreconditioned resid norm 7.879630242934e-05 true resid norm 7.879630242930e-05 ||r(i)||/||b|| 2.084958371341e-04 963 5 KSP unpreconditioned resid norm 5.569022502593e-06 true resid norm 5.569022502581e-06 ||r(i)||/||b|| 1.473569156035e-05 964 6 KSP unpreconditioned resid norm 5.821644196935e-07 true resid norm 5.821644196750e-07 ||r(i)||/||b|| 1.540413119496e-06 965 7 KSP unpreconditioned resid norm 4.532618828776e-09 true resid norm 4.532618817180e-09 ||r(i)||/||b|| 1.199335661144e-08 966 8 KSP unpreconditioned resid norm 1.763889785026e-11 true resid norm 1.763888741317e-11 ||r(i)||/||b|| 4.667267985859e-11 967 1 SNES Function norm 1.764e-11 968 Nonlinear solve converged due to CONVERGED_FNORM_RELATIVE iterations 1 969SNES Object: 1 MPI processes 970 type: newtonls 971 maximum iterations=50, maximum function evaluations=10000 972 tolerances: relative=1e-08, absolute=1e-50, solution=1e-08 973 total number of linear solver iterations=8 974 total number of function evaluations=2 975 norm schedule ALWAYS 976 SNESLineSearch Object: 1 MPI processes 977 type: bt 978 interpolation: cubic 979 alpha=1.000000e-04 980 maxstep=1.000000e+08, minlambda=1.000000e-12 981 tolerances: relative=1.000000e-08, absolute=1.000000e-15, lambda=1.000000e-08 982 maximum iterations=40 983 KSP Object: 1 MPI processes 984 type: fgmres 985 restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement 986 happy breakdown tolerance 1e-30 987 maximum iterations=10000, initial guess is zero 988 tolerances: relative=1e-09, absolute=1e-50, divergence=10000. 989 right preconditioning 990 using UNPRECONDITIONED norm type for convergence test 991 PC Object: 1 MPI processes 992 type: mg 993 type is MULTIPLICATIVE, levels=2 cycles=v 994 Cycles per PCApply=1 995 Not using Galerkin computed coarse grid matrices 996 Coarse grid solver -- level ------------------------------- 997 KSP Object: (mg_coarse_) 1 MPI processes 998 type: preonly 999 maximum iterations=10000, initial guess is zero 1000 tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 1001 left preconditioning 1002 using NONE norm type for convergence test 1003 PC Object: (mg_coarse_) 1 MPI processes 1004 type: lu 1005 out-of-place factorization 1006 tolerance for zero pivot 2.22045e-14 1007 using diagonal shift on blocks to prevent zero pivot [INBLOCKS] 1008 matrix ordering: nd 1009 factor fill ratio given 5., needed 1. 1010 Factored matrix follows: 1011 Mat Object: 1 MPI processes 1012 type: seqaij 1013 rows=1, cols=1 1014 package used to perform factorization: petsc 1015 total: nonzeros=1, allocated nonzeros=1 1016 not using I-node routines 1017 linear system matrix = precond matrix: 1018 Mat Object: 1 MPI processes 1019 type: seqaij 1020 rows=1, cols=1 1021 total: nonzeros=1, allocated nonzeros=1 1022 total number of mallocs used during MatSetValues calls=0 1023 not using I-node routines 1024 Down solver (pre-smoother) on level 1 ------------------------------- 1025 KSP Object: (mg_levels_1_) 1 MPI processes 1026 type: chebyshev 1027 eigenvalue estimates used: min = 0.151841, max = 1.67025 1028 eigenvalues estimate via gmres min 0.994613, max 1.51841 1029 eigenvalues estimated using gmres with translations [0. 0.1; 0. 1.1] 1030 KSP Object: (mg_levels_1_esteig_) 1 MPI processes 1031 type: gmres 1032 restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement 1033 happy breakdown tolerance 1e-30 1034 maximum iterations=10, initial guess is zero 1035 tolerances: relative=1e-12, absolute=1e-50, divergence=10000. 1036 left preconditioning 1037 using PRECONDITIONED norm type for convergence test 1038 estimating eigenvalues using noisy right hand side 1039 maximum iterations=2, nonzero initial guess 1040 tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 1041 left preconditioning 1042 using NONE norm type for convergence test 1043 PC Object: (mg_levels_1_) 1 MPI processes 1044 type: sor 1045 type = local_symmetric, iterations = 1, local iterations = 1, omega = 1. 1046 linear system matrix = precond matrix: 1047 Mat Object: 1 MPI processes 1048 type: seqaij 1049 rows=9, cols=9 1050 total: nonzeros=41, allocated nonzeros=41 1051 total number of mallocs used during MatSetValues calls=0 1052 not using I-node routines 1053 Up solver (post-smoother) same as down solver (pre-smoother) 1054 linear system matrix = precond matrix: 1055 Mat Object: 1 MPI processes 1056 type: seqaij 1057 rows=9, cols=9 1058 total: nonzeros=41, allocated nonzeros=41 1059 total number of mallocs used during MatSetValues calls=0 1060 not using I-node routines 1061Timestep: 0010 time = 1 L_2 Error: 0.41293 106210 TS dt 0.1 time 1. 1063