Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Newton solver not converging, np=2, 3D_Skew_Steady #34

Open
eyoung55 opened this issue May 12, 2021 · 2 comments
Open

Newton solver not converging, np=2, 3D_Skew_Steady #34

eyoung55 opened this issue May 12, 2021 · 2 comments
Labels
bug Something isn't working

Comments

@eyoung55
Copy link
Collaborator

eyoung55 commented May 12, 2021

The regression test 3D_Skew_Steady fails when np=2 as the Newton solver blows up. Possibly an error in marking boundaries or parent-child facets being unreliable in parallel? Output is below.

Solving
|    
Process 0: No Jacobian form specified for nonlinear variational problem.
Process 0: Differentiating residual form F to obtain Jacobian J = F'.
Process 0: Solving nonlinear variational problem.
  Process 0: Newton iteration 0: r (abs) = 1.120e+06 (tol = 1.000e-06) r (rel) = 1.000e+00 (tol = 1.000e-05)
  Process 0: Newton iteration 1: r (abs) = 6.355e+05 (tol = 1.000e-06) r (rel) = 5.674e-01 (tol = 1.000e-05)
  Process 0: Newton iteration 2: r (abs) = 1.274e+06 (tol = 1.000e-06) r (rel) = 1.137e+00 (tol = 1.000e-05)
  Process 0: Newton iteration 3: r (abs) = 8.578e+05 (tol = 1.000e-06) r (rel) = 7.658e-01 (tol = 1.000e-05)
  Process 0: Newton iteration 4: r (abs) = 1.393e+07 (tol = 1.000e-06) r (rel) = 1.244e+01 (tol = 1.000e-05)
  Process 0: Newton iteration 5: r (abs) = 2.407e+08 (tol = 1.000e-06) r (rel) = 2.149e+02 (tol = 1.000e-05)
  Process 0: Newton iteration 6: r (abs) = 8.855e+07 (tol = 1.000e-06) r (rel) = 7.906e+01 (tol = 1.000e-05)
  Process 0: Newton iteration 7: r (abs) = 4.285e+07 (tol = 1.000e-06) r (rel) = 3.825e+01 (tol = 1.000e-05)
  Process 0: Newton iteration 8: r (abs) = 9.166e+07 (tol = 1.000e-06) r (rel) = 8.183e+01 (tol = 1.000e-05)
  Process 0: Newton iteration 9: r (abs) = 2.379e+13 (tol = 1.000e-06) r (rel) = 2.124e+07 (tol = 1.000e-05)
  Process 0: Newton iteration 10: r (abs) = 6.146e+12 (tol = 1.000e-06) r (rel) = 5.487e+06 (tol = 1.000e-05)
  Process 0: Newton iteration 11: r (abs) = 1.653e+12 (tol = 1.000e-06) r (rel) = 1.476e+06 (tol = 1.000e-05)
  Process 0: Newton iteration 12: r (abs) = 2.040e+12 (tol = 1.000e-06) r (rel) = 1.821e+06 (tol = 1.000e-05)
  Process 0: Newton iteration 13: r (abs) = 1.371e+14 (tol = 1.000e-06) r (rel) = 1.224e+08 (tol = 1.000e-05)
  Process 0: Newton iteration 14: r (abs) = 3.583e+13 (tol = 1.000e-06) r (rel) = 3.198e+07 (tol = 1.000e-05)
  Process 0: Newton iteration 15: r (abs) = 1.344e+15 (tol = 1.000e-06) r (rel) = 1.200e+09 (tol = 1.000e-05)
  Process 0: Newton iteration 16: r (abs) = 3.470e+14 (tol = 1.000e-06) r (rel) = 3.098e+08 (tol = 1.000e-05)
  Process 0: Newton iteration 17: r (abs) = 1.248e+14 (tol = 1.000e-06) r (rel) = 1.114e+08 (tol = 1.000e-05)
  Process 0: Newton iteration 18: r (abs) = 8.591e+14 (tol = 1.000e-06) r (rel) = 7.670e+08 (tol = 1.000e-05)
  Process 0: Newton iteration 19: r (abs) = 2.622e+15 (tol = 1.000e-06) r (rel) = 2.341e+09 (tol = 1.000e-05)
  Process 0: Newton iteration 20: r (abs) = 7.180e+14 (tol = 1.000e-06) r (rel) = 6.410e+08 (tol = 1.000e-05)
  Process 0: Newton iteration 21: r (abs) = 7.329e+14 (tol = 1.000e-06) r (rel) = 6.543e+08 (tol = 1.000e-05)
  Process 0: Newton iteration 22: r (abs) = 2.961e+14 (tol = 1.000e-06) r (rel) = 2.643e+08 (tol = 1.000e-05)
  Process 0: Newton iteration 23: r (abs) = 7.600e+13 (tol = 1.000e-06) r (rel) = 6.785e+07 (tol = 1.000e-05)
  Process 0: Newton iteration 24: r (abs) = 2.846e+14 (tol = 1.000e-06) r (rel) = 2.541e+08 (tol = 1.000e-05)
  Process 0: Newton iteration 25: r (abs) = 3.824e+14 (tol = 1.000e-06) r (rel) = 3.414e+08 (tol = 1.000e-05)
  Process 0: Newton iteration 26: r (abs) = 1.306e+14 (tol = 1.000e-06) r (rel) = 1.166e+08 (tol = 1.000e-05)
  Process 0: Newton iteration 27: r (abs) = 4.862e+13 (tol = 1.000e-06) r (rel) = 4.340e+07 (tol = 1.000e-05)
  Process 0: Newton iteration 28: r (abs) = 2.073e+13 (tol = 1.000e-06) r (rel) = 1.851e+07 (tol = 1.000e-05)
  Process 0: Newton iteration 29: r (abs) = 3.471e+14 (tol = 1.000e-06) r (rel) = 3.099e+08 (tol = 1.000e-05)
  Process 0: Newton iteration 30: r (abs) = 1.154e+14 (tol = 1.000e-06) r (rel) = 1.030e+08 (tol = 1.000e-05)
  Process 0: Newton iteration 31: r (abs) = 3.055e+13 (tol = 1.000e-06) r (rel) = 2.728e+07 (tol = 1.000e-05)
  Process 0: Newton iteration 32: r (abs) = 3.287e+13 (tol = 1.000e-06) r (rel) = 2.934e+07 (tol = 1.000e-05)
  Process 0: Newton iteration 33: r (abs) = 1.092e+13 (tol = 1.000e-06) r (rel) = 9.745e+06 (tol = 1.000e-05)
  Process 0: Newton iteration 34: r (abs) = 5.564e+12 (tol = 1.000e-06) r (rel) = 4.967e+06 (tol = 1.000e-05)
  Process 0: Newton iteration 35: r (abs) = 9.925e+12 (tol = 1.000e-06) r (rel) = 8.861e+06 (tol = 1.000e-05)
  Process 0: Newton iteration 36: r (abs) = 1.294e+13 (tol = 1.000e-06) r (rel) = 1.155e+07 (tol = 1.000e-05)
  Process 0: Newton iteration 37: r (abs) = 9.416e+13 (tol = 1.000e-06) r (rel) = 8.407e+07 (tol = 1.000e-05)
  Process 0: Newton iteration 38: r (abs) = 2.460e+13 (tol = 1.000e-06) r (rel) = 2.196e+07 (tol = 1.000e-05)
  Process 0: Newton iteration 39: r (abs) = 6.238e+14 (tol = 1.000e-06) r (rel) = 5.569e+08 (tol = 1.000e-05)
  Process 0: Newton iteration 40: r (abs) = 1.583e+14 (tol = 1.000e-06) r (rel) = 1.414e+08 (tol = 1.000e-05)
@eyoung55 eyoung55 added the bug Something isn't working label May 12, 2021
@eyoung55 eyoung55 added this to To do in Parallelization May 12, 2021
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
bug Something isn't working
Projects
Development

No branches or pull requests

2 participants
@eyoung55 and others