| /libCEED/examples/fluids/qfunctions/ |
| H A D | newtonian_types.h | e334ad8fca0b75b5d678f79cfd10f4127eb21e08 Fri May 13 03:44:12 UTC 2022 Jed Brown <jed@jedbrown.org> examples/fluids: add IJacobian for SUPG
I think this is exact for Galerkin on channel.yaml and blasius.yaml.
Specifically, SNES converges almost perfectly in this test than when
-snes_mf_operator is added.
-options_file examples/fluids/channel.yaml -snes_monitor -ksp_converged_reason -stab none -ts_max_steps 3 -newtonian_unit_tests -ksp_rtol 1e-10 -order 2 -ts_dt 3e-6 -pc_type lu
I say "think" because something is fishy with the viscous energy term
velocity*stress. If I disable that term, then convergence is clearly
better with the analytic Jacobian (makes sense because it's accuracy is
machine epsilon) than with -snes_mf_operator. This is unexpected, but
the difference is near the limit of numerical stability so I'm not
confident it's wrong. Also, the term itself is so simple.
The Jacobian is inexact for SUPG in that Tau_d is frozen.
WIP: examples/fluids: support MatShell and Blasius_Inflow_Jacobian
|
| H A D | blasius.h | e334ad8fca0b75b5d678f79cfd10f4127eb21e08 Fri May 13 03:44:12 UTC 2022 Jed Brown <jed@jedbrown.org> examples/fluids: add IJacobian for SUPG
I think this is exact for Galerkin on channel.yaml and blasius.yaml.
Specifically, SNES converges almost perfectly in this test than when
-snes_mf_operator is added.
-options_file examples/fluids/channel.yaml -snes_monitor -ksp_converged_reason -stab none -ts_max_steps 3 -newtonian_unit_tests -ksp_rtol 1e-10 -order 2 -ts_dt 3e-6 -pc_type lu
I say "think" because something is fishy with the viscous energy term
velocity*stress. If I disable that term, then convergence is clearly
better with the analytic Jacobian (makes sense because it's accuracy is
machine epsilon) than with -snes_mf_operator. This is unexpected, but
the difference is near the limit of numerical stability so I'm not
confident it's wrong. Also, the term itself is so simple.
The Jacobian is inexact for SUPG in that Tau_d is frozen.
WIP: examples/fluids: support MatShell and Blasius_Inflow_Jacobian
|
| H A D | newtonian.h | e334ad8fca0b75b5d678f79cfd10f4127eb21e08 Fri May 13 03:44:12 UTC 2022 Jed Brown <jed@jedbrown.org> examples/fluids: add IJacobian for SUPG
I think this is exact for Galerkin on channel.yaml and blasius.yaml.
Specifically, SNES converges almost perfectly in this test than when
-snes_mf_operator is added.
-options_file examples/fluids/channel.yaml -snes_monitor -ksp_converged_reason -stab none -ts_max_steps 3 -newtonian_unit_tests -ksp_rtol 1e-10 -order 2 -ts_dt 3e-6 -pc_type lu
I say "think" because something is fishy with the viscous energy term
velocity*stress. If I disable that term, then convergence is clearly
better with the analytic Jacobian (makes sense because it's accuracy is
machine epsilon) than with -snes_mf_operator. This is unexpected, but
the difference is near the limit of numerical stability so I'm not
confident it's wrong. Also, the term itself is so simple.
The Jacobian is inexact for SUPG in that Tau_d is frozen.
WIP: examples/fluids: support MatShell and Blasius_Inflow_Jacobian
|
| /libCEED/examples/fluids/ |
| H A D | Makefile | e334ad8fca0b75b5d678f79cfd10f4127eb21e08 Fri May 13 03:44:12 UTC 2022 Jed Brown <jed@jedbrown.org> examples/fluids: add IJacobian for SUPG
I think this is exact for Galerkin on channel.yaml and blasius.yaml.
Specifically, SNES converges almost perfectly in this test than when
-snes_mf_operator is added.
-options_file examples/fluids/channel.yaml -snes_monitor -ksp_converged_reason -stab none -ts_max_steps 3 -newtonian_unit_tests -ksp_rtol 1e-10 -order 2 -ts_dt 3e-6 -pc_type lu
I say "think" because something is fishy with the viscous energy term
velocity*stress. If I disable that term, then convergence is clearly
better with the analytic Jacobian (makes sense because it's accuracy is
machine epsilon) than with -snes_mf_operator. This is unexpected, but
the difference is near the limit of numerical stability so I'm not
confident it's wrong. Also, the term itself is so simple.
The Jacobian is inexact for SUPG in that Tau_d is frozen.
WIP: examples/fluids: support MatShell and Blasius_Inflow_Jacobian
|
| H A D | navierstokes.h | e334ad8fca0b75b5d678f79cfd10f4127eb21e08 Fri May 13 03:44:12 UTC 2022 Jed Brown <jed@jedbrown.org> examples/fluids: add IJacobian for SUPG
I think this is exact for Galerkin on channel.yaml and blasius.yaml.
Specifically, SNES converges almost perfectly in this test than when
-snes_mf_operator is added.
-options_file examples/fluids/channel.yaml -snes_monitor -ksp_converged_reason -stab none -ts_max_steps 3 -newtonian_unit_tests -ksp_rtol 1e-10 -order 2 -ts_dt 3e-6 -pc_type lu
I say "think" because something is fishy with the viscous energy term
velocity*stress. If I disable that term, then convergence is clearly
better with the analytic Jacobian (makes sense because it's accuracy is
machine epsilon) than with -snes_mf_operator. This is unexpected, but
the difference is near the limit of numerical stability so I'm not
confident it's wrong. Also, the term itself is so simple.
The Jacobian is inexact for SUPG in that Tau_d is frozen.
WIP: examples/fluids: support MatShell and Blasius_Inflow_Jacobian
|
| H A D | navierstokes.c | e334ad8fca0b75b5d678f79cfd10f4127eb21e08 Fri May 13 03:44:12 UTC 2022 Jed Brown <jed@jedbrown.org> examples/fluids: add IJacobian for SUPG
I think this is exact for Galerkin on channel.yaml and blasius.yaml.
Specifically, SNES converges almost perfectly in this test than when
-snes_mf_operator is added.
-options_file examples/fluids/channel.yaml -snes_monitor -ksp_converged_reason -stab none -ts_max_steps 3 -newtonian_unit_tests -ksp_rtol 1e-10 -order 2 -ts_dt 3e-6 -pc_type lu
I say "think" because something is fishy with the viscous energy term
velocity*stress. If I disable that term, then convergence is clearly
better with the analytic Jacobian (makes sense because it's accuracy is
machine epsilon) than with -snes_mf_operator. This is unexpected, but
the difference is near the limit of numerical stability so I'm not
confident it's wrong. Also, the term itself is so simple.
The Jacobian is inexact for SUPG in that Tau_d is frozen.
WIP: examples/fluids: support MatShell and Blasius_Inflow_Jacobian
|
| /libCEED/examples/fluids/problems/ |
| H A D | blasius.c | e334ad8fca0b75b5d678f79cfd10f4127eb21e08 Fri May 13 03:44:12 UTC 2022 Jed Brown <jed@jedbrown.org> examples/fluids: add IJacobian for SUPG
I think this is exact for Galerkin on channel.yaml and blasius.yaml.
Specifically, SNES converges almost perfectly in this test than when
-snes_mf_operator is added.
-options_file examples/fluids/channel.yaml -snes_monitor -ksp_converged_reason -stab none -ts_max_steps 3 -newtonian_unit_tests -ksp_rtol 1e-10 -order 2 -ts_dt 3e-6 -pc_type lu
I say "think" because something is fishy with the viscous energy term
velocity*stress. If I disable that term, then convergence is clearly
better with the analytic Jacobian (makes sense because it's accuracy is
machine epsilon) than with -snes_mf_operator. This is unexpected, but
the difference is near the limit of numerical stability so I'm not
confident it's wrong. Also, the term itself is so simple.
The Jacobian is inexact for SUPG in that Tau_d is frozen.
WIP: examples/fluids: support MatShell and Blasius_Inflow_Jacobian
|
| H A D | newtonian.c | e334ad8fca0b75b5d678f79cfd10f4127eb21e08 Fri May 13 03:44:12 UTC 2022 Jed Brown <jed@jedbrown.org> examples/fluids: add IJacobian for SUPG
I think this is exact for Galerkin on channel.yaml and blasius.yaml.
Specifically, SNES converges almost perfectly in this test than when
-snes_mf_operator is added.
-options_file examples/fluids/channel.yaml -snes_monitor -ksp_converged_reason -stab none -ts_max_steps 3 -newtonian_unit_tests -ksp_rtol 1e-10 -order 2 -ts_dt 3e-6 -pc_type lu
I say "think" because something is fishy with the viscous energy term
velocity*stress. If I disable that term, then convergence is clearly
better with the analytic Jacobian (makes sense because it's accuracy is
machine epsilon) than with -snes_mf_operator. This is unexpected, but
the difference is near the limit of numerical stability so I'm not
confident it's wrong. Also, the term itself is so simple.
The Jacobian is inexact for SUPG in that Tau_d is frozen.
WIP: examples/fluids: support MatShell and Blasius_Inflow_Jacobian
|
| /libCEED/examples/fluids/src/ |
| H A D | setupts.c | e334ad8fca0b75b5d678f79cfd10f4127eb21e08 Fri May 13 03:44:12 UTC 2022 Jed Brown <jed@jedbrown.org> examples/fluids: add IJacobian for SUPG
I think this is exact for Galerkin on channel.yaml and blasius.yaml.
Specifically, SNES converges almost perfectly in this test than when
-snes_mf_operator is added.
-options_file examples/fluids/channel.yaml -snes_monitor -ksp_converged_reason -stab none -ts_max_steps 3 -newtonian_unit_tests -ksp_rtol 1e-10 -order 2 -ts_dt 3e-6 -pc_type lu
I say "think" because something is fishy with the viscous energy term
velocity*stress. If I disable that term, then convergence is clearly
better with the analytic Jacobian (makes sense because it's accuracy is
machine epsilon) than with -snes_mf_operator. This is unexpected, but
the difference is near the limit of numerical stability so I'm not
confident it's wrong. Also, the term itself is so simple.
The Jacobian is inexact for SUPG in that Tau_d is frozen.
WIP: examples/fluids: support MatShell and Blasius_Inflow_Jacobian
|
| H A D | setuplibceed.c | e334ad8fca0b75b5d678f79cfd10f4127eb21e08 Fri May 13 03:44:12 UTC 2022 Jed Brown <jed@jedbrown.org> examples/fluids: add IJacobian for SUPG
I think this is exact for Galerkin on channel.yaml and blasius.yaml.
Specifically, SNES converges almost perfectly in this test than when
-snes_mf_operator is added.
-options_file examples/fluids/channel.yaml -snes_monitor -ksp_converged_reason -stab none -ts_max_steps 3 -newtonian_unit_tests -ksp_rtol 1e-10 -order 2 -ts_dt 3e-6 -pc_type lu
I say "think" because something is fishy with the viscous energy term
velocity*stress. If I disable that term, then convergence is clearly
better with the analytic Jacobian (makes sense because it's accuracy is
machine epsilon) than with -snes_mf_operator. This is unexpected, but
the difference is near the limit of numerical stability so I'm not
confident it's wrong. Also, the term itself is so simple.
The Jacobian is inexact for SUPG in that Tau_d is frozen.
WIP: examples/fluids: support MatShell and Blasius_Inflow_Jacobian
|