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

Split Form Class #3

Closed
wants to merge 99 commits into from
Closed

Split Form Class #3

wants to merge 99 commits into from

Conversation

abtin98
Copy link

@abtin98 abtin98 commented Jul 9, 2019

Still need to figure out what variables I need and how/where to store them, so this is very rough and subject to change.

dougshidong and others added 23 commits June 26, 2019 21:18
Use FEValues as done in deal.II tutorials
Goal is to derive DGStrong and DGWeak from DGBase.

The loop over the cells will be the same. The volume, boundary and internal face integrals will differ.
No need to manually destruct the pointer anymore
To see the difference, try using

vim -d dg/weak_dg.cpp dg/strong_dg.cpp
Using a ChartManifold to create the bump that decays.
Added test case.

TODO:

Code up wall and inlet/outlet boundary conditions
Non-dimensionalize initial conditions to use Mach number

Check entropy convergence.
Faster compilation times.

Used to recompile everything anytime physics.h was touched.
dougshidong and others added 21 commits August 20, 2019 10:13
Change some of the longer test to use MPI. All the tests pass

Use cmake ../ -DMPIMAX to set the default number of processor in those tests. Not the same as ctest -jX.

Test project /home/ddong/Codes/PHiLiP_temp/PHiLiP/build_release
      Start  1: 1D_numerical_flux_conservation
 1/39 Test  dougshidong#1: 1D_numerical_flux_conservation ..................................   Passed    0.26 sec
      Start  2: 2D_numerical_flux_conservation
 2/39 Test  dougshidong#2: 2D_numerical_flux_conservation ..................................   Passed    0.22 sec
      Start  3: 3D_numerical_flux_conservation
 3/39 Test  dougshidong#3: 3D_numerical_flux_conservation ..................................   Passed    0.27 sec
      Start  4: 1D_jacobian_matrix_regression
 4/39 Test  dougshidong#4: 1D_jacobian_matrix_regression ...................................   Passed    0.69 sec
      Start  5: 2D_jacobian_matrix_regression
 5/39 Test  dougshidong#5: 2D_jacobian_matrix_regression ...................................   Passed    0.74 sec
      Start  6: 3D_jacobian_matrix_regression
 6/39 Test  dougshidong#6: 3D_jacobian_matrix_regression ...................................   Passed    3.27 sec
      Start  7: 1D_euler_convert_primitive_conservative
 7/39 Test  dougshidong#7: 1D_euler_convert_primitive_conservative .........................   Passed    0.19 sec
      Start  8: 2D_euler_convert_primitive_conservative
 8/39 Test  dougshidong#8: 2D_euler_convert_primitive_conservative .........................   Passed    0.22 sec
      Start  9: 3D_euler_convert_primitive_conservative
 9/39 Test  dougshidong#9: 3D_euler_convert_primitive_conservative .........................   Passed    0.57 sec
      Start 10: 1D_euler_manufactured_solution_source
10/39 Test dougshidong#10: 1D_euler_manufactured_solution_source ...........................   Passed    0.27 sec
      Start 11: 2D_euler_manufactured_solution_source
11/39 Test dougshidong#11: 2D_euler_manufactured_solution_source ...........................   Passed    0.24 sec
      Start 12: 3D_euler_manufactured_solution_source
12/39 Test dougshidong#12: 3D_euler_manufactured_solution_source ...........................   Passed    0.61 sec
      Start 13: 1D_euler_convective_jacobian
13/39 Test dougshidong#13: 1D_euler_convective_jacobian ....................................   Passed    0.20 sec
      Start 14: 2D_euler_convective_jacobian
14/39 Test dougshidong#14: 2D_euler_convective_jacobian ....................................   Passed    0.24 sec
      Start 15: 3D_euler_convective_jacobian
15/39 Test dougshidong#15: 3D_euler_convective_jacobian ....................................   Passed    1.04 sec
      Start 16: 1D_ADVECTION_IMPLICIT_MANUFACTURED_SOLUTION
16/39 Test dougshidong#16: 1D_ADVECTION_IMPLICIT_MANUFACTURED_SOLUTION .....................   Passed    0.79 sec
      Start 17: 2D_ADVECTION_IMPLICIT_MANUFACTURED_SOLUTION
17/39 Test dougshidong#17: 2D_ADVECTION_IMPLICIT_MANUFACTURED_SOLUTION .....................   Passed    7.34 sec
      Start 18: MPI_2D_ADVECTION_IMPLICIT_MANUFACTURED_SOLUTION
18/39 Test dougshidong#18: MPI_2D_ADVECTION_IMPLICIT_MANUFACTURED_SOLUTION .................   Passed    3.57 sec
      Start 19: MPI_3D_ADVECTION_IMPLICIT_MANUFACTURED_SOLUTION
19/39 Test dougshidong#19: MPI_3D_ADVECTION_IMPLICIT_MANUFACTURED_SOLUTION .................   Passed    7.67 sec
      Start 20: 1D_ADVECTION_IMPLICIT_MANUFACTURED_SOLUTION_STRONG
20/39 Test dougshidong#20: 1D_ADVECTION_IMPLICIT_MANUFACTURED_SOLUTION_STRONG ..............   Passed    0.82 sec
      Start 21: 2D_ADVECTION_IMPLICIT_MANUFACTURED_SOLUTION_STRONG
21/39 Test dougshidong#21: 2D_ADVECTION_IMPLICIT_MANUFACTURED_SOLUTION_STRONG ..............   Passed   10.93 sec
      Start 22: MPI_3D_ADVECTION_IMPLICIT_MANUFACTURED_SOLUTION_STRONG
22/39 Test dougshidong#22: MPI_3D_ADVECTION_IMPLICIT_MANUFACTURED_SOLUTION_STRONG ..........   Passed   37.31 sec
      Start 23: 1D_DIFFUSION_IMPLICIT_MANUFACTURED_SOLUTION
23/39 Test dougshidong#23: 1D_DIFFUSION_IMPLICIT_MANUFACTURED_SOLUTION .....................   Passed    0.80 sec
      Start 24: 2D_DIFFUSION_IMPLICIT_MANUFACTURED_SOLUTION
24/39 Test dougshidong#24: 2D_DIFFUSION_IMPLICIT_MANUFACTURED_SOLUTION .....................   Passed    2.56 sec
      Start 25: MPI_3D_DIFFUSION_IMPLICIT_MANUFACTURED_SOLUTION
25/39 Test dougshidong#25: MPI_3D_DIFFUSION_IMPLICIT_MANUFACTURED_SOLUTION .................   Passed   21.02 sec
      Start 26: 1D_CONVECTION_DIFFUSION_IMPLICIT_MANUFACTURED_SOLUTION
26/39 Test dougshidong#26: 1D_CONVECTION_DIFFUSION_IMPLICIT_MANUFACTURED_SOLUTION ..........   Passed    0.91 sec
      Start 27: 2D_CONVECTION_DIFFUSION_IMPLICIT_MANUFACTURED_SOLUTION
27/39 Test dougshidong#27: 2D_CONVECTION_DIFFUSION_IMPLICIT_MANUFACTURED_SOLUTION ..........   Passed    2.62 sec
      Start 28: MPI_3D_CONVECTION_DIFFUSION_IMPLICIT_MANUFACTURED_SOLUTION
28/39 Test dougshidong#28: MPI_3D_CONVECTION_DIFFUSION_IMPLICIT_MANUFACTURED_SOLUTION ......   Passed   19.84 sec
      Start 29: 1D_ADVECTION_VECTOR_VALUED_IMPLICIT_MANUFACTURED_SOLUTION
29/39 Test dougshidong#29: 1D_ADVECTION_VECTOR_VALUED_IMPLICIT_MANUFACTURED_SOLUTION .......   Passed    0.83 sec
      Start 30: 2D_ADVECTION_VECTOR_VALUED_IMPLICIT_MANUFACTURED_SOLUTION
30/39 Test dougshidong#30: 2D_ADVECTION_VECTOR_VALUED_IMPLICIT_MANUFACTURED_SOLUTION .......   Passed   12.23 sec
      Start 31: MPI_3D_ADVECTION_VECTOR_VALUED_IMPLICIT_MANUFACTURED_SOLUTION
31/39 Test dougshidong#31: MPI_3D_ADVECTION_VECTOR_VALUED_IMPLICIT_MANUFACTURED_SOLUTION ...   Passed   12.63 sec
      Start 32: 1D_BURGERS_INVISCID_IMPLICIT_MANUFACTURED_SOLUTION
32/39 Test dougshidong#32: 1D_BURGERS_INVISCID_IMPLICIT_MANUFACTURED_SOLUTION ..............   Passed    0.84 sec
      Start 33: 1D_EULER_LAXFRIEDRICHS_MANUFACTURED_SOLUTION
33/39 Test dougshidong#33: 1D_EULER_LAXFRIEDRICHS_MANUFACTURED_SOLUTION ....................   Passed    1.64 sec
      Start 34: 2D_EULER_LAXFRIEDRICHS_MANUFACTURED_SOLUTION
34/39 Test dougshidong#34: 2D_EULER_LAXFRIEDRICHS_MANUFACTURED_SOLUTION ....................   Passed   38.39 sec
      Start 35: MPI_2D_EULER_LAXFRIEDRICHS_MANUFACTURED_SOLUTION
35/39 Test dougshidong#35: MPI_2D_EULER_LAXFRIEDRICHS_MANUFACTURED_SOLUTION ................   Passed   15.22 sec
      Start 36: 1D_EULER_ROE_MANUFACTURED_SOLUTION
36/39 Test dougshidong#36: 1D_EULER_ROE_MANUFACTURED_SOLUTION ..............................   Passed    1.69 sec
      Start 37: MPI_2D_EULER_ROE_MANUFACTURED_SOLUTION
37/39 Test dougshidong#37: MPI_2D_EULER_ROE_MANUFACTURED_SOLUTION ..........................   Passed   10.56 sec
      Start 38: MPI_2D_EULER_INTEGRATION_CYLINDER
38/39 Test dougshidong#38: MPI_2D_EULER_INTEGRATION_CYLINDER ...............................   Passed  120.18 sec
      Start 39: MPI_2D_EULER_INTEGRATION_GAUSSIAN_BUMP
39/39 Test dougshidong#39: MPI_2D_EULER_INTEGRATION_GAUSSIAN_BUMP ..........................   Passed  312.47 sec

100% tests passed, 0 tests failed out of 39

Total Test time (real) = 652.82 sec
Very slow... 600 seconds instead of 3 seconds.
I would suggest searching the word 'Passed' to see the orders retrieved.
Test project /home/abtinameri/Desktop/thesis/PHiLiP/build
      Start  1: 1D_numerical_flux_conservation
 1/42 Test  dougshidong#1: 1D_numerical_flux_conservation ..................................   Passed    0.42 sec
      Start  2: 2D_numerical_flux_conservation
 2/42 Test  dougshidong#2: 2D_numerical_flux_conservation ..................................   Passed    0.45 sec
      Start  3: 3D_numerical_flux_conservation
 3/42 Test  dougshidong#3: 3D_numerical_flux_conservation ..................................   Passed    0.40 sec
      Start  4: 1D_jacobian_matrix_regression
 4/42 Test  dougshidong#4: 1D_jacobian_matrix_regression ...................................   Passed    0.42 sec
      Start  5: 2D_jacobian_matrix_regression
 5/42 Test  dougshidong#5: 2D_jacobian_matrix_regression ...................................   Passed    0.47 sec
      Start  6: 3D_jacobian_matrix_regression
 6/42 Test  dougshidong#6: 3D_jacobian_matrix_regression ...................................   Passed    2.48 sec
      Start  7: 1D_euler_convert_primitive_conservative
 7/42 Test  dougshidong#7: 1D_euler_convert_primitive_conservative .........................   Passed    0.22 sec
      Start  8: 2D_euler_convert_primitive_conservative
 8/42 Test  dougshidong#8: 2D_euler_convert_primitive_conservative .........................   Passed    0.21 sec
      Start  9: 3D_euler_convert_primitive_conservative
 9/42 Test  dougshidong#9: 3D_euler_convert_primitive_conservative .........................   Passed    0.48 sec
      Start 10: 1D_euler_manufactured_solution_source
10/42 Test dougshidong#10: 1D_euler_manufactured_solution_source ...........................   Passed    0.21 sec
      Start 11: 2D_euler_manufactured_solution_source
11/42 Test dougshidong#11: 2D_euler_manufactured_solution_source ...........................   Passed    0.25 sec
      Start 12: 3D_euler_manufactured_solution_source
12/42 Test dougshidong#12: 3D_euler_manufactured_solution_source ...........................   Passed    0.51 sec
      Start 13: 1D_euler_convective_jacobian
13/42 Test dougshidong#13: 1D_euler_convective_jacobian ....................................   Passed    0.21 sec
      Start 14: 2D_euler_convective_jacobian
14/42 Test dougshidong#14: 2D_euler_convective_jacobian ....................................   Passed    0.23 sec
      Start 15: 3D_euler_convective_jacobian
15/42 Test dougshidong#15: 3D_euler_convective_jacobian ....................................   Passed    0.87 sec
      Start 16: 1D_ADVECTION_IMPLICIT_MANUFACTURED_SOLUTION
16/42 Test dougshidong#16: 1D_ADVECTION_IMPLICIT_MANUFACTURED_SOLUTION .....................   Passed    0.52 sec
      Start 17: 2D_ADVECTION_IMPLICIT_MANUFACTURED_SOLUTION
17/42 Test dougshidong#17: 2D_ADVECTION_IMPLICIT_MANUFACTURED_SOLUTION .....................   Passed    4.98 sec
      Start 18: MPI_2D_ADVECTION_EXPLICIT_MANUFACTURED_SOLUTION
18/42 Test dougshidong#18: MPI_2D_ADVECTION_EXPLICIT_MANUFACTURED_SOLUTION .................   Passed  284.29 sec
      Start 19: MPI_2D_ADVECTION_IMPLICIT_MANUFACTURED_SOLUTION
19/42 Test dougshidong#19: MPI_2D_ADVECTION_IMPLICIT_MANUFACTURED_SOLUTION .................   Passed    3.93 sec
      Start 20: MPI_3D_ADVECTION_IMPLICIT_MANUFACTURED_SOLUTION
20/42 Test dougshidong#20: MPI_3D_ADVECTION_IMPLICIT_MANUFACTURED_SOLUTION .................   Passed    9.86 sec
      Start 21: 1D_ADVECTION_IMPLICIT_MANUFACTURED_SOLUTION_STRONG
21/42 Test dougshidong#21: 1D_ADVECTION_IMPLICIT_MANUFACTURED_SOLUTION_STRONG ..............   Passed    0.59 sec
      Start 22: 2D_ADVECTION_IMPLICIT_MANUFACTURED_SOLUTION_STRONG
22/42 Test dougshidong#22: 2D_ADVECTION_IMPLICIT_MANUFACTURED_SOLUTION_STRONG ..............   Passed    7.40 sec
      Start 23: MPI_3D_ADVECTION_IMPLICIT_MANUFACTURED_SOLUTION_STRONG
23/42 Test dougshidong#23: MPI_3D_ADVECTION_IMPLICIT_MANUFACTURED_SOLUTION_STRONG ..........   Passed   41.23 sec
      Start 24: 1D_DIFFUSION_IMPLICIT_MANUFACTURED_SOLUTION
24/42 Test dougshidong#24: 1D_DIFFUSION_IMPLICIT_MANUFACTURED_SOLUTION .....................   Passed    0.53 sec
      Start 25: 2D_DIFFUSION_IMPLICIT_MANUFACTURED_SOLUTION
25/42 Test dougshidong#25: 2D_DIFFUSION_IMPLICIT_MANUFACTURED_SOLUTION .....................   Passed    1.86 sec
      Start 26: MPI_3D_DIFFUSION_IMPLICIT_MANUFACTURED_SOLUTION
26/42 Test dougshidong#26: MPI_3D_DIFFUSION_IMPLICIT_MANUFACTURED_SOLUTION .................   Passed   29.47 sec
      Start 27: 1D_CONVECTION_DIFFUSION_IMPLICIT_MANUFACTURED_SOLUTION
27/42 Test dougshidong#27: 1D_CONVECTION_DIFFUSION_IMPLICIT_MANUFACTURED_SOLUTION ..........   Passed    0.58 sec
      Start 28: 2D_CONVECTION_DIFFUSION_IMPLICIT_MANUFACTURED_SOLUTION
28/42 Test dougshidong#28: 2D_CONVECTION_DIFFUSION_IMPLICIT_MANUFACTURED_SOLUTION ..........   Passed    1.94 sec
      Start 29: MPI_3D_CONVECTION_DIFFUSION_IMPLICIT_MANUFACTURED_SOLUTION
29/42 Test dougshidong#29: MPI_3D_CONVECTION_DIFFUSION_IMPLICIT_MANUFACTURED_SOLUTION ......   Passed   28.02 sec
      Start 30: 1D_ADVECTION_VECTOR_VALUED_IMPLICIT_MANUFACTURED_SOLUTION
30/42 Test dougshidong#30: 1D_ADVECTION_VECTOR_VALUED_IMPLICIT_MANUFACTURED_SOLUTION .......   Passed    0.57 sec
      Start 31: 2D_ADVECTION_VECTOR_VALUED_IMPLICIT_MANUFACTURED_SOLUTION
31/42 Test dougshidong#31: 2D_ADVECTION_VECTOR_VALUED_IMPLICIT_MANUFACTURED_SOLUTION .......   Passed    9.53 sec
      Start 32: MPI_3D_ADVECTION_VECTOR_VALUED_IMPLICIT_MANUFACTURED_SOLUTION
32/42 Test dougshidong#32: MPI_3D_ADVECTION_VECTOR_VALUED_IMPLICIT_MANUFACTURED_SOLUTION ...   Passed   17.23 sec
      Start 33: 1D_BURGERS_INVISCID_IMPLICIT_MANUFACTURED_SOLUTION
33/42 Test dougshidong#33: 1D_BURGERS_INVISCID_IMPLICIT_MANUFACTURED_SOLUTION ..............   Passed    0.54 sec
      Start 34: 1D_burgers_energy_stability
34/42 Test dougshidong#34: 1D_burgers_energy_stability .....................................   Passed  211.12 sec
      Start 35: 2D_advection_explicit_periodic
35/42 Test dougshidong#35: 2D_advection_explicit_periodic ..................................   Passed  170.45 sec
      Start 36: 1D_EULER_LAXFRIEDRICHS_MANUFACTURED_SOLUTION
36/42 Test dougshidong#36: 1D_EULER_LAXFRIEDRICHS_MANUFACTURED_SOLUTION ....................   Passed    1.13 sec
      Start 37: 2D_EULER_LAXFRIEDRICHS_MANUFACTURED_SOLUTION
37/42 Test dougshidong#37: 2D_EULER_LAXFRIEDRICHS_MANUFACTURED_SOLUTION ....................   Passed   29.69 sec
      Start 38: MPI_2D_EULER_LAXFRIEDRICHS_MANUFACTURED_SOLUTION
38/42 Test dougshidong#38: MPI_2D_EULER_LAXFRIEDRICHS_MANUFACTURED_SOLUTION ................   Passed   22.17 sec
      Start 39: 1D_EULER_ROE_MANUFACTURED_SOLUTION
39/42 Test dougshidong#39: 1D_EULER_ROE_MANUFACTURED_SOLUTION ..............................   Passed    1.13 sec
      Start 40: MPI_2D_EULER_ROE_MANUFACTURED_SOLUTION
40/42 Test dougshidong#40: MPI_2D_EULER_ROE_MANUFACTURED_SOLUTION ..........................   Passed   14.58 sec
      Start 41: MPI_2D_EULER_INTEGRATION_CYLINDER
41/42 Test dougshidong#41: MPI_2D_EULER_INTEGRATION_CYLINDER ...............................   Passed  167.23 sec
      Start 42: MPI_2D_EULER_INTEGRATION_GAUSSIAN_BUMP
42/42 Test dougshidong#42: MPI_2D_EULER_INTEGRATION_GAUSSIAN_BUMP ..........................   Passed  245.55 sec

100% tests passed, 0 tests failed out of 42
dougshidong added a commit that referenced this pull request May 19, 2022
Physics and numerical fluxes are now computed outside of the DG discretization.
Partial implementation of vector-valued problems by adding the std::array structures

Diffusive terms require a function that will evaluate the solution numerical flux and the viscous numerical flux.
Current implementation of SIPG is working to recover optimal convergence orders of the solution but not of the output functional

Suboptimal output error is perceived for diffusion (maybe, I get 2p) and convection-diffusion (definitely)

TODO: Find the adjoint consistent boundary condition for the diffusion problem and convection-diffusion problem
TODO: Find out why the linear solver is soooo much slower for the diffusion problems

      Start  1: 1D_ADVECTION_IMPLICIT_MANUFACTURED_SOLUTION
 1/12 Test  #1: 1D_ADVECTION_IMPLICIT_MANUFACTURED_SOLUTION ..............   Passed    0.53 sec
      Start  2: 2D_ADVECTION_IMPLICIT_MANUFACTURED_SOLUTION
 2/12 Test  #2: 2D_ADVECTION_IMPLICIT_MANUFACTURED_SOLUTION ..............   Passed    5.84 sec
      Start  3: 3D_ADVECTION_IMPLICIT_MANUFACTURED_SOLUTION
 3/12 Test  #3: 3D_ADVECTION_IMPLICIT_MANUFACTURED_SOLUTION ..............   Passed   47.98 sec
      Start  4: 1D_DIFFUSION_IMPLICIT_MANUFACTURED_SOLUTION
 4/12 Test  #4: 1D_DIFFUSION_IMPLICIT_MANUFACTURED_SOLUTION ..............   Passed    0.64 sec
      Start  5: 2D_DIFFUSION_IMPLICIT_MANUFACTURED_SOLUTION
 5/12 Test  #5: 2D_DIFFUSION_IMPLICIT_MANUFACTURED_SOLUTION ..............   Passed    3.23 sec
      Start  6: 3D_DIFFUSION_IMPLICIT_MANUFACTURED_SOLUTION
 6/12 Test  #6: 3D_DIFFUSION_IMPLICIT_MANUFACTURED_SOLUTION ..............   Passed  237.84 sec
      Start  7: 1D_CONVECTION_DIFFUSION_IMPLICIT_MANUFACTURED_SOLUTION
 7/12 Test  #7: 1D_CONVECTION_DIFFUSION_IMPLICIT_MANUFACTURED_SOLUTION ...   Passed    0.75 sec
      Start  8: 2D_CONVECTION_DIFFUSION_IMPLICIT_MANUFACTURED_SOLUTION
 8/12 Test  #8: 2D_CONVECTION_DIFFUSION_IMPLICIT_MANUFACTURED_SOLUTION ...   Passed    3.45 sec
      Start  9: 3D_CONVECTION_DIFFUSION_IMPLICIT_MANUFACTURED_SOLUTION
 9/12 Test  #9: 3D_CONVECTION_DIFFUSION_IMPLICIT_MANUFACTURED_SOLUTION ...   Passed  266.84 sec
      Start 10: 1D_numerical_flux_conservation
10/12 Test #10: 1D_numerical_flux_conservation ...........................   Passed    0.27 sec
      Start 11: 2D_numerical_flux_conservation
11/12 Test #11: 2D_numerical_flux_conservation ...........................   Passed    0.25 sec
      Start 12: 3D_numerical_flux_conservation
12/12 Test #12: 3D_numerical_flux_conservation ...........................   Passed    0.24 sec

100% tests passed, 0 tests failed out of 12

Total Test time (real) = 567.87 sec
dougshidong added a commit that referenced this pull request May 19, 2022
Checks source term against a finite-difference of the convective fluxes.

Since the source term is analytically derived for the manufactured solution,
it also checks for the physical flux in a way.

Also, new test for grid convegence of Euler. Not supposed to pass

Test project /home/ddong/Codes/PHiLiP/build_release
      Start  1: 1D_numerical_flux_conservation
 1/27 Test  #1: 1D_numerical_flux_conservation ..............................   Passed    0.17 sec
      Start  2: 2D_numerical_flux_conservation
 2/27 Test  #2: 2D_numerical_flux_conservation ..............................   Passed    0.17 sec
      Start  3: 3D_numerical_flux_conservation
 3/27 Test  #3: 3D_numerical_flux_conservation ..............................   Passed    0.16 sec
      Start  4: 1D_jacobian_matrix_regression
 4/27 Test  #4: 1D_jacobian_matrix_regression ...............................   Passed    0.45 sec
      Start  5: 2D_jacobian_matrix_regression
 5/27 Test  #5: 2D_jacobian_matrix_regression ...............................   Passed    0.50 sec
      Start  6: 3D_jacobian_matrix_regression
 6/27 Test  #6: 3D_jacobian_matrix_regression ...............................   Passed    1.57 sec
      Start  7: 1D_euler_convert_primitive_conservative
 7/27 Test  #7: 1D_euler_convert_primitive_conservative .....................   Passed    0.16 sec
      Start  8: 2D_euler_convert_primitive_conservative
 8/27 Test  #8: 2D_euler_convert_primitive_conservative .....................   Passed    0.20 sec
      Start  9: 3D_euler_convert_primitive_conservative
 9/27 Test  #9: 3D_euler_convert_primitive_conservative .....................   Passed    0.44 sec
      Start 10: 1D_euler_manufactured_solution_source
10/27 Test #10: 1D_euler_manufactured_solution_source .......................   Passed    0.16 sec
      Start 11: 2D_euler_manufactured_solution_source
11/27 Test #11: 2D_euler_manufactured_solution_source .......................   Passed    0.21 sec
      Start 12: 3D_euler_manufactured_solution_source
12/27 Test #12: 3D_euler_manufactured_solution_source .......................   Passed    0.54 sec
      Start 13: 1D_ADVECTION_IMPLICIT_MANUFACTURED_SOLUTION
13/27 Test #13: 1D_ADVECTION_IMPLICIT_MANUFACTURED_SOLUTION .................   Passed    0.49 sec
      Start 14: 2D_ADVECTION_IMPLICIT_MANUFACTURED_SOLUTION
14/27 Test #14: 2D_ADVECTION_IMPLICIT_MANUFACTURED_SOLUTION .................   Passed    2.38 sec
      Start 15: 3D_ADVECTION_IMPLICIT_MANUFACTURED_SOLUTION
15/27 Test #15: 3D_ADVECTION_IMPLICIT_MANUFACTURED_SOLUTION .................   Passed   25.05 sec
      Start 16: 1D_DIFFUSION_IMPLICIT_MANUFACTURED_SOLUTION
16/27 Test #16: 1D_DIFFUSION_IMPLICIT_MANUFACTURED_SOLUTION .................   Passed    0.58 sec
      Start 17: 2D_DIFFUSION_IMPLICIT_MANUFACTURED_SOLUTION
17/27 Test #17: 2D_DIFFUSION_IMPLICIT_MANUFACTURED_SOLUTION .................   Passed    1.49 sec
      Start 18: 3D_DIFFUSION_IMPLICIT_MANUFACTURED_SOLUTION
18/27 Test #18: 3D_DIFFUSION_IMPLICIT_MANUFACTURED_SOLUTION .................   Passed   48.04 sec
      Start 19: 1D_CONVECTION_DIFFUSION_IMPLICIT_MANUFACTURED_SOLUTION
19/27 Test #19: 1D_CONVECTION_DIFFUSION_IMPLICIT_MANUFACTURED_SOLUTION ......   Passed    0.51 sec
      Start 20: 2D_CONVECTION_DIFFUSION_IMPLICIT_MANUFACTURED_SOLUTION
20/27 Test #20: 2D_CONVECTION_DIFFUSION_IMPLICIT_MANUFACTURED_SOLUTION ......   Passed    1.50 sec
      Start 21: 3D_CONVECTION_DIFFUSION_IMPLICIT_MANUFACTURED_SOLUTION
21/27 Test #21: 3D_CONVECTION_DIFFUSION_IMPLICIT_MANUFACTURED_SOLUTION ......   Passed   43.90 sec
      Start 22: 1D_ADVECTION_VECTOR_VALUED_IMPLICIT_MANUFACTURED_SOLUTION
22/27 Test #22: 1D_ADVECTION_VECTOR_VALUED_IMPLICIT_MANUFACTURED_SOLUTION ...   Passed    0.49 sec
      Start 23: 2D_ADVECTION_VECTOR_VALUED_IMPLICIT_MANUFACTURED_SOLUTION
23/27 Test #23: 2D_ADVECTION_VECTOR_VALUED_IMPLICIT_MANUFACTURED_SOLUTION ...   Passed    6.13 sec
      Start 24: 3D_ADVECTION_VECTOR_VALUED_IMPLICIT_MANUFACTURED_SOLUTION
24/27 Test #24: 3D_ADVECTION_VECTOR_VALUED_IMPLICIT_MANUFACTURED_SOLUTION ...   Passed   20.49 sec
      Start 25: 1D_EULER_IMPLICIT_MANUFACTURED_SOLUTION
25/27 Test #25: 1D_EULER_IMPLICIT_MANUFACTURED_SOLUTION .....................***Failed    0.42 sec
      Start 26: 2D_EULER_IMPLICIT_MANUFACTURED_SOLUTION
26/27 Test #26: 2D_EULER_IMPLICIT_MANUFACTURED_SOLUTION .....................***Failed    0.41 sec
      Start 27: 3D_EULER_IMPLICIT_MANUFACTURED_SOLUTION
27/27 Test #27: 3D_EULER_IMPLICIT_MANUFACTURED_SOLUTION .....................***Failed    0.41 sec

89% tests passed, 3 tests failed out of 27

Total Test time (real) = 157.22 sec

The following tests FAILED:
	 25 - 1D_EULER_IMPLICIT_MANUFACTURED_SOLUTION (Failed)
	 26 - 2D_EULER_IMPLICIT_MANUFACTURED_SOLUTION (Failed)
	 27 - 3D_EULER_IMPLICIT_MANUFACTURED_SOLUTION (Failed)
Errors while running CTest
dougshidong added a commit that referenced this pull request May 19, 2022
Change some of the longer test to use MPI. All the tests pass

Use cmake ../ -DMPIMAX to set the default number of processor in those tests. Not the same as ctest -jX.

Test project /home/ddong/Codes/PHiLiP_temp/PHiLiP/build_release
      Start  1: 1D_numerical_flux_conservation
 1/39 Test  #1: 1D_numerical_flux_conservation ..................................   Passed    0.26 sec
      Start  2: 2D_numerical_flux_conservation
 2/39 Test  #2: 2D_numerical_flux_conservation ..................................   Passed    0.22 sec
      Start  3: 3D_numerical_flux_conservation
 3/39 Test  #3: 3D_numerical_flux_conservation ..................................   Passed    0.27 sec
      Start  4: 1D_jacobian_matrix_regression
 4/39 Test  #4: 1D_jacobian_matrix_regression ...................................   Passed    0.69 sec
      Start  5: 2D_jacobian_matrix_regression
 5/39 Test  #5: 2D_jacobian_matrix_regression ...................................   Passed    0.74 sec
      Start  6: 3D_jacobian_matrix_regression
 6/39 Test  #6: 3D_jacobian_matrix_regression ...................................   Passed    3.27 sec
      Start  7: 1D_euler_convert_primitive_conservative
 7/39 Test  #7: 1D_euler_convert_primitive_conservative .........................   Passed    0.19 sec
      Start  8: 2D_euler_convert_primitive_conservative
 8/39 Test  #8: 2D_euler_convert_primitive_conservative .........................   Passed    0.22 sec
      Start  9: 3D_euler_convert_primitive_conservative
 9/39 Test  #9: 3D_euler_convert_primitive_conservative .........................   Passed    0.57 sec
      Start 10: 1D_euler_manufactured_solution_source
10/39 Test #10: 1D_euler_manufactured_solution_source ...........................   Passed    0.27 sec
      Start 11: 2D_euler_manufactured_solution_source
11/39 Test #11: 2D_euler_manufactured_solution_source ...........................   Passed    0.24 sec
      Start 12: 3D_euler_manufactured_solution_source
12/39 Test #12: 3D_euler_manufactured_solution_source ...........................   Passed    0.61 sec
      Start 13: 1D_euler_convective_jacobian
13/39 Test #13: 1D_euler_convective_jacobian ....................................   Passed    0.20 sec
      Start 14: 2D_euler_convective_jacobian
14/39 Test #14: 2D_euler_convective_jacobian ....................................   Passed    0.24 sec
      Start 15: 3D_euler_convective_jacobian
15/39 Test #15: 3D_euler_convective_jacobian ....................................   Passed    1.04 sec
      Start 16: 1D_ADVECTION_IMPLICIT_MANUFACTURED_SOLUTION
16/39 Test #16: 1D_ADVECTION_IMPLICIT_MANUFACTURED_SOLUTION .....................   Passed    0.79 sec
      Start 17: 2D_ADVECTION_IMPLICIT_MANUFACTURED_SOLUTION
17/39 Test #17: 2D_ADVECTION_IMPLICIT_MANUFACTURED_SOLUTION .....................   Passed    7.34 sec
      Start 18: MPI_2D_ADVECTION_IMPLICIT_MANUFACTURED_SOLUTION
18/39 Test #18: MPI_2D_ADVECTION_IMPLICIT_MANUFACTURED_SOLUTION .................   Passed    3.57 sec
      Start 19: MPI_3D_ADVECTION_IMPLICIT_MANUFACTURED_SOLUTION
19/39 Test #19: MPI_3D_ADVECTION_IMPLICIT_MANUFACTURED_SOLUTION .................   Passed    7.67 sec
      Start 20: 1D_ADVECTION_IMPLICIT_MANUFACTURED_SOLUTION_STRONG
20/39 Test #20: 1D_ADVECTION_IMPLICIT_MANUFACTURED_SOLUTION_STRONG ..............   Passed    0.82 sec
      Start 21: 2D_ADVECTION_IMPLICIT_MANUFACTURED_SOLUTION_STRONG
21/39 Test #21: 2D_ADVECTION_IMPLICIT_MANUFACTURED_SOLUTION_STRONG ..............   Passed   10.93 sec
      Start 22: MPI_3D_ADVECTION_IMPLICIT_MANUFACTURED_SOLUTION_STRONG
22/39 Test #22: MPI_3D_ADVECTION_IMPLICIT_MANUFACTURED_SOLUTION_STRONG ..........   Passed   37.31 sec
      Start 23: 1D_DIFFUSION_IMPLICIT_MANUFACTURED_SOLUTION
23/39 Test #23: 1D_DIFFUSION_IMPLICIT_MANUFACTURED_SOLUTION .....................   Passed    0.80 sec
      Start 24: 2D_DIFFUSION_IMPLICIT_MANUFACTURED_SOLUTION
24/39 Test #24: 2D_DIFFUSION_IMPLICIT_MANUFACTURED_SOLUTION .....................   Passed    2.56 sec
      Start 25: MPI_3D_DIFFUSION_IMPLICIT_MANUFACTURED_SOLUTION
25/39 Test #25: MPI_3D_DIFFUSION_IMPLICIT_MANUFACTURED_SOLUTION .................   Passed   21.02 sec
      Start 26: 1D_CONVECTION_DIFFUSION_IMPLICIT_MANUFACTURED_SOLUTION
26/39 Test #26: 1D_CONVECTION_DIFFUSION_IMPLICIT_MANUFACTURED_SOLUTION ..........   Passed    0.91 sec
      Start 27: 2D_CONVECTION_DIFFUSION_IMPLICIT_MANUFACTURED_SOLUTION
27/39 Test #27: 2D_CONVECTION_DIFFUSION_IMPLICIT_MANUFACTURED_SOLUTION ..........   Passed    2.62 sec
      Start 28: MPI_3D_CONVECTION_DIFFUSION_IMPLICIT_MANUFACTURED_SOLUTION
28/39 Test #28: MPI_3D_CONVECTION_DIFFUSION_IMPLICIT_MANUFACTURED_SOLUTION ......   Passed   19.84 sec
      Start 29: 1D_ADVECTION_VECTOR_VALUED_IMPLICIT_MANUFACTURED_SOLUTION
29/39 Test #29: 1D_ADVECTION_VECTOR_VALUED_IMPLICIT_MANUFACTURED_SOLUTION .......   Passed    0.83 sec
      Start 30: 2D_ADVECTION_VECTOR_VALUED_IMPLICIT_MANUFACTURED_SOLUTION
30/39 Test #30: 2D_ADVECTION_VECTOR_VALUED_IMPLICIT_MANUFACTURED_SOLUTION .......   Passed   12.23 sec
      Start 31: MPI_3D_ADVECTION_VECTOR_VALUED_IMPLICIT_MANUFACTURED_SOLUTION
31/39 Test #31: MPI_3D_ADVECTION_VECTOR_VALUED_IMPLICIT_MANUFACTURED_SOLUTION ...   Passed   12.63 sec
      Start 32: 1D_BURGERS_INVISCID_IMPLICIT_MANUFACTURED_SOLUTION
32/39 Test #32: 1D_BURGERS_INVISCID_IMPLICIT_MANUFACTURED_SOLUTION ..............   Passed    0.84 sec
      Start 33: 1D_EULER_LAXFRIEDRICHS_MANUFACTURED_SOLUTION
33/39 Test #33: 1D_EULER_LAXFRIEDRICHS_MANUFACTURED_SOLUTION ....................   Passed    1.64 sec
      Start 34: 2D_EULER_LAXFRIEDRICHS_MANUFACTURED_SOLUTION
34/39 Test #34: 2D_EULER_LAXFRIEDRICHS_MANUFACTURED_SOLUTION ....................   Passed   38.39 sec
      Start 35: MPI_2D_EULER_LAXFRIEDRICHS_MANUFACTURED_SOLUTION
35/39 Test #35: MPI_2D_EULER_LAXFRIEDRICHS_MANUFACTURED_SOLUTION ................   Passed   15.22 sec
      Start 36: 1D_EULER_ROE_MANUFACTURED_SOLUTION
36/39 Test #36: 1D_EULER_ROE_MANUFACTURED_SOLUTION ..............................   Passed    1.69 sec
      Start 37: MPI_2D_EULER_ROE_MANUFACTURED_SOLUTION
37/39 Test #37: MPI_2D_EULER_ROE_MANUFACTURED_SOLUTION ..........................   Passed   10.56 sec
      Start 38: MPI_2D_EULER_INTEGRATION_CYLINDER
38/39 Test #38: MPI_2D_EULER_INTEGRATION_CYLINDER ...............................   Passed  120.18 sec
      Start 39: MPI_2D_EULER_INTEGRATION_GAUSSIAN_BUMP
39/39 Test #39: MPI_2D_EULER_INTEGRATION_GAUSSIAN_BUMP ..........................   Passed  312.47 sec

100% tests passed, 0 tests failed out of 39

Total Test time (real) = 652.82 sec
dougshidong pushed a commit that referenced this pull request May 19, 2022
Test project /home/abtinameri/Desktop/thesis/PHiLiP/build
      Start  1: 1D_numerical_flux_conservation
 1/42 Test  #1: 1D_numerical_flux_conservation ..................................   Passed    0.42 sec
      Start  2: 2D_numerical_flux_conservation
 2/42 Test  #2: 2D_numerical_flux_conservation ..................................   Passed    0.45 sec
      Start  3: 3D_numerical_flux_conservation
 3/42 Test  #3: 3D_numerical_flux_conservation ..................................   Passed    0.40 sec
      Start  4: 1D_jacobian_matrix_regression
 4/42 Test  #4: 1D_jacobian_matrix_regression ...................................   Passed    0.42 sec
      Start  5: 2D_jacobian_matrix_regression
 5/42 Test  #5: 2D_jacobian_matrix_regression ...................................   Passed    0.47 sec
      Start  6: 3D_jacobian_matrix_regression
 6/42 Test  #6: 3D_jacobian_matrix_regression ...................................   Passed    2.48 sec
      Start  7: 1D_euler_convert_primitive_conservative
 7/42 Test  #7: 1D_euler_convert_primitive_conservative .........................   Passed    0.22 sec
      Start  8: 2D_euler_convert_primitive_conservative
 8/42 Test  #8: 2D_euler_convert_primitive_conservative .........................   Passed    0.21 sec
      Start  9: 3D_euler_convert_primitive_conservative
 9/42 Test  #9: 3D_euler_convert_primitive_conservative .........................   Passed    0.48 sec
      Start 10: 1D_euler_manufactured_solution_source
10/42 Test #10: 1D_euler_manufactured_solution_source ...........................   Passed    0.21 sec
      Start 11: 2D_euler_manufactured_solution_source
11/42 Test #11: 2D_euler_manufactured_solution_source ...........................   Passed    0.25 sec
      Start 12: 3D_euler_manufactured_solution_source
12/42 Test #12: 3D_euler_manufactured_solution_source ...........................   Passed    0.51 sec
      Start 13: 1D_euler_convective_jacobian
13/42 Test #13: 1D_euler_convective_jacobian ....................................   Passed    0.21 sec
      Start 14: 2D_euler_convective_jacobian
14/42 Test #14: 2D_euler_convective_jacobian ....................................   Passed    0.23 sec
      Start 15: 3D_euler_convective_jacobian
15/42 Test #15: 3D_euler_convective_jacobian ....................................   Passed    0.87 sec
      Start 16: 1D_ADVECTION_IMPLICIT_MANUFACTURED_SOLUTION
16/42 Test #16: 1D_ADVECTION_IMPLICIT_MANUFACTURED_SOLUTION .....................   Passed    0.52 sec
      Start 17: 2D_ADVECTION_IMPLICIT_MANUFACTURED_SOLUTION
17/42 Test #17: 2D_ADVECTION_IMPLICIT_MANUFACTURED_SOLUTION .....................   Passed    4.98 sec
      Start 18: MPI_2D_ADVECTION_EXPLICIT_MANUFACTURED_SOLUTION
18/42 Test #18: MPI_2D_ADVECTION_EXPLICIT_MANUFACTURED_SOLUTION .................   Passed  284.29 sec
      Start 19: MPI_2D_ADVECTION_IMPLICIT_MANUFACTURED_SOLUTION
19/42 Test #19: MPI_2D_ADVECTION_IMPLICIT_MANUFACTURED_SOLUTION .................   Passed    3.93 sec
      Start 20: MPI_3D_ADVECTION_IMPLICIT_MANUFACTURED_SOLUTION
20/42 Test #20: MPI_3D_ADVECTION_IMPLICIT_MANUFACTURED_SOLUTION .................   Passed    9.86 sec
      Start 21: 1D_ADVECTION_IMPLICIT_MANUFACTURED_SOLUTION_STRONG
21/42 Test #21: 1D_ADVECTION_IMPLICIT_MANUFACTURED_SOLUTION_STRONG ..............   Passed    0.59 sec
      Start 22: 2D_ADVECTION_IMPLICIT_MANUFACTURED_SOLUTION_STRONG
22/42 Test #22: 2D_ADVECTION_IMPLICIT_MANUFACTURED_SOLUTION_STRONG ..............   Passed    7.40 sec
      Start 23: MPI_3D_ADVECTION_IMPLICIT_MANUFACTURED_SOLUTION_STRONG
23/42 Test #23: MPI_3D_ADVECTION_IMPLICIT_MANUFACTURED_SOLUTION_STRONG ..........   Passed   41.23 sec
      Start 24: 1D_DIFFUSION_IMPLICIT_MANUFACTURED_SOLUTION
24/42 Test #24: 1D_DIFFUSION_IMPLICIT_MANUFACTURED_SOLUTION .....................   Passed    0.53 sec
      Start 25: 2D_DIFFUSION_IMPLICIT_MANUFACTURED_SOLUTION
25/42 Test #25: 2D_DIFFUSION_IMPLICIT_MANUFACTURED_SOLUTION .....................   Passed    1.86 sec
      Start 26: MPI_3D_DIFFUSION_IMPLICIT_MANUFACTURED_SOLUTION
26/42 Test #26: MPI_3D_DIFFUSION_IMPLICIT_MANUFACTURED_SOLUTION .................   Passed   29.47 sec
      Start 27: 1D_CONVECTION_DIFFUSION_IMPLICIT_MANUFACTURED_SOLUTION
27/42 Test #27: 1D_CONVECTION_DIFFUSION_IMPLICIT_MANUFACTURED_SOLUTION ..........   Passed    0.58 sec
      Start 28: 2D_CONVECTION_DIFFUSION_IMPLICIT_MANUFACTURED_SOLUTION
28/42 Test #28: 2D_CONVECTION_DIFFUSION_IMPLICIT_MANUFACTURED_SOLUTION ..........   Passed    1.94 sec
      Start 29: MPI_3D_CONVECTION_DIFFUSION_IMPLICIT_MANUFACTURED_SOLUTION
29/42 Test #29: MPI_3D_CONVECTION_DIFFUSION_IMPLICIT_MANUFACTURED_SOLUTION ......   Passed   28.02 sec
      Start 30: 1D_ADVECTION_VECTOR_VALUED_IMPLICIT_MANUFACTURED_SOLUTION
30/42 Test #30: 1D_ADVECTION_VECTOR_VALUED_IMPLICIT_MANUFACTURED_SOLUTION .......   Passed    0.57 sec
      Start 31: 2D_ADVECTION_VECTOR_VALUED_IMPLICIT_MANUFACTURED_SOLUTION
31/42 Test #31: 2D_ADVECTION_VECTOR_VALUED_IMPLICIT_MANUFACTURED_SOLUTION .......   Passed    9.53 sec
      Start 32: MPI_3D_ADVECTION_VECTOR_VALUED_IMPLICIT_MANUFACTURED_SOLUTION
32/42 Test #32: MPI_3D_ADVECTION_VECTOR_VALUED_IMPLICIT_MANUFACTURED_SOLUTION ...   Passed   17.23 sec
      Start 33: 1D_BURGERS_INVISCID_IMPLICIT_MANUFACTURED_SOLUTION
33/42 Test #33: 1D_BURGERS_INVISCID_IMPLICIT_MANUFACTURED_SOLUTION ..............   Passed    0.54 sec
      Start 34: 1D_burgers_energy_stability
34/42 Test #34: 1D_burgers_energy_stability .....................................   Passed  211.12 sec
      Start 35: 2D_advection_explicit_periodic
35/42 Test #35: 2D_advection_explicit_periodic ..................................   Passed  170.45 sec
      Start 36: 1D_EULER_LAXFRIEDRICHS_MANUFACTURED_SOLUTION
36/42 Test #36: 1D_EULER_LAXFRIEDRICHS_MANUFACTURED_SOLUTION ....................   Passed    1.13 sec
      Start 37: 2D_EULER_LAXFRIEDRICHS_MANUFACTURED_SOLUTION
37/42 Test #37: 2D_EULER_LAXFRIEDRICHS_MANUFACTURED_SOLUTION ....................   Passed   29.69 sec
      Start 38: MPI_2D_EULER_LAXFRIEDRICHS_MANUFACTURED_SOLUTION
38/42 Test #38: MPI_2D_EULER_LAXFRIEDRICHS_MANUFACTURED_SOLUTION ................   Passed   22.17 sec
      Start 39: 1D_EULER_ROE_MANUFACTURED_SOLUTION
39/42 Test #39: 1D_EULER_ROE_MANUFACTURED_SOLUTION ..............................   Passed    1.13 sec
      Start 40: MPI_2D_EULER_ROE_MANUFACTURED_SOLUTION
40/42 Test #40: MPI_2D_EULER_ROE_MANUFACTURED_SOLUTION ..........................   Passed   14.58 sec
      Start 41: MPI_2D_EULER_INTEGRATION_CYLINDER
41/42 Test #41: MPI_2D_EULER_INTEGRATION_CYLINDER ...............................   Passed  167.23 sec
      Start 42: MPI_2D_EULER_INTEGRATION_GAUSSIAN_BUMP
42/42 Test #42: MPI_2D_EULER_INTEGRATION_GAUSSIAN_BUMP ..........................   Passed  245.55 sec

100% tests passed, 0 tests failed out of 42
jbrillon added a commit that referenced this pull request Feb 11, 2023
* ESFR Split forms. Contains all unsteady unit and integration tests. Makes use of operators in DG strong. Provides weight-adjusted curvilinear Mass matrices. Includes free-stream preserving metric terms.

* Removed the accidentally pushed mesh file. Will rebase last 4 commits into 1.

* Cleaning up the accidental added line at end of some files and removed #if 0 #endif commented blocks.

* And this file deleted last line.

* Finished removing #if 0 #endif blocks, added some comments and included new numerical fluxes in the numerical_flux_conservation test.

* Apply suggestions from code review

dg.cpp minor edits

* Implemented sum factorization in operators class. All operations can be written as a combination of matrix-vector products and inner products. A unit test verifies by testing the cputime for A*u and also M*u where M is a mass matrix--uses both matrix-vector product and inner product forms. All that is left is to clean up the unit test a bit to fail at wrong orders.

* Apply suggestions from code review

* Apply suggestions from code review

* Apply suggestions from code review

* Update src/parameters/all_parameters.cpp

* Update with all changes requested from pull request. Only missing merge conflict with ode solver base then pull request changes complete.

* Small cleanup changes before I merge upstream master in my branch.

* Made corrections to explicit ode solver with Butcher tableau c values. Updated unsteady tests. Continued to clean up and removed deprecated Burgers Rewienski lines.

* Cleaned up last few comments on pull request wrt Flux reconstruction operators and flags.

* Removed old references to Runge-Kutta order.

* Cleaned up re-used code in operators, and other comments from PR.

* Removed collection tuple from dg. Removed finite element and quadrature collections from DG. Made consistent throughout the code to call operators finite element and quadrature collections from the operators collection tuple.

* Fixed typos in error messages.

* Created initial condition base class with interpolate and projection capabilities for different initial condition functions. This generalizes the initial condition class.

* Turned axuiliary mass matrix into a single rather than array of dim because we use the same mass matrix for each dim direction. This reduces memory consumption.

* Templated initial condition base by nstate and the functions therein.

* Added a strong form Navier-Stokes integration test.

* Made mapping support points arrays instead of vectors. Added flag in dg allocate system for the implicit matrix/vectors to allocate to reduce memory costs.

* Fix minor changes from Pranshul.

* Small update. Switching branch for other work.

* Auxiliary variable correctly computed for NAvier-Stokes equation now. Verified through viscous TGV test strong.

* Fixed small compiling error.

* Fix TGV viscous strong control file

* Strong viscous TGV energy check test confirmed pass.

* Added definitions for projection.

* Changed initial condition base to just initial condition since nothing is derived from it.

* Current status of new operators backup

* Made entropy conserving flux work for all physics types. Will return to my other branch with restructuring operators class for efficiency now.

* Recent version, new operators quarter done.

* For energy tests, only initializes ode solver on first iteration now.

* For Julien to debug his flow solver case. I will revert these changes after so that the test case is mroe reasonable when using ctests.

* Fixed the Burgers' energy inviscid initial condition.

* Fixed stash error

* Apply suggestions from code review

* cleaning up InitialCondition; changed to SetInitialCondition

* removing commented lines from previous commit

* fixing bug that would run two tgv energy check quick tests

* cleaning up src/parameter/ files; fixing comments marked as fixed that were not actually fixed

* minor editorial changes

* correcting all_parameters.cpp from last commit; indentation and doxygen fixes

* indentation clean up in unit tests and euler split tgv

* cleaning up src/testing/ files

* removing unused CurvManifold class declaration in advection_explicit_periodic.h

* changing tabs to spaces in burgers_stability.cpp

* cleaning up ConvectionDiffusionPeriodic class files in src/testing/

* reducing code repetition in CMakeList for flow_variable_tests

* reviewing src/dg, src/numerical_flux, and src/ode_solver/ files; moving all convective num fluxes to convective_numerical_flux.<cpp/hpp>; converting tabs to spaces in convective_numerical_flux.cpp; general clean up in numerical flux directory; indentation fixes

* New operators class completed with passing unit tests. Only unit tests to change atm are sum factorization, metric splitting and weight adjusted inverse. Next step is to deprecate the previous operators class and make respective changes throughout dg, tests etc.

* reviewing src/mesh/grids files

* making repeated constants class members in NonsymmetricCurvedGridManifold

* correcting "set test_type = flow_solver" to "set run_type = flow_simulation" in prm files

* The sum factorization test.

* Changed mapping shape functions class so metric operators less variables pass to construct.

* fixing time refinement study parameter file

* New operators class complete with top-nothc efficiency :)

* adding print statements for operators allocation status

* adding forgotten print statements

* comment fix

* resolving a PR comment

* Implemented new operators class throughout DG/DG strong. All that is left is having it pass the tests. It currently compiles.

* I have DG strong conservative working, now need to tweak split forms and pass all tests.

* Conservative DG works 2D curvilinear. Now need cuvilinear split form to work properly, then will do auxtesting/tgv and finito.

* Got curvilinear stability and weight-adjusted curvilinear working perfectly. Now only needs to pass auxiliary/conv_diff explicit and it's all done!

* 3D TGV works. Now just auxiliary.

* bug fix for writing pvtu files in subdirectory

* setting default value for WALE model constant to that recommended by literature

* Auxiliary and conv diffusion passes.

* Inverse mass on the fly and mass evaluated on the fly in application per cell done. Only doesnt work for curvilinear at the moment. Has flag to use. For curvilinear, do not use inverse on the fly yet. Works perfectly for linear grids.

* 2pt flux Hadamard sum after the Hadamard product now O n^d+1. All is finished.

* Made mass matrix computation more efficient. FINALLY PASSES NACA 0012 OPTIMIZATION I'm so happy!

* compiles. Have not yet tested since the merge.

* Fixed TGV Euler split form test that got stuff deleted during the merge. Also, in mass matrix DG, changed condition of factoring out det of Jac by Cartesian element rather than linear. This is because unstructured linear meshes do NOT have a constant metric Jacobian determinant, only Cartesian meshes do. Will run ctest and clean up documentation now.

* Moved the entropy consevring, entropy stable split, and central numerical fluxes into convective numerical flux.

* Changed tabs to spaces in burgers epxlicit test.

* Added documentation equations to DG strong.

* Modified TGV split test to use initial condition TGV function.

* Fixed Polynomial degree ramping issue.

* Run curvilinear mesh on Narval. Also, made 3D periodic nonsymm curv grid.

* Added energy file declaration from control file. Also increased warping of grid.

* small clean up

* fixing pde_type flag in strong dg

* correcting initial condition setting in split tgv

* Finished with curv TGV work for conference

* 1D_BURGERS_STABILITY_ORDERS_OF_ACCURACY_LONG passing

* 1D_TIME_REFINEMENT_STUDY_ADVECTION_EXPLICIT now passes and behaves as supposed to.

* 1D periodic now CONFIRMED to work with dealii matched faces. Thus, the hardcoded 1D periodic was removed. All tests that are affected by this pass.

* Viscous TGV energy check quick STRONG passing.

* Smoothed out auxiliary test to reach asymptotic convergence in L infinity norm in a reasonable time/grid level.

* Unit test verifying DG Strong Auxiliary equations' right-hand-side for d\in[1,3], for Navier-Stokes equations.

* Fixed CMake files and cleaned up the include header files in unit test.

* Resolved conflicts from merge that caused the tests to fail. Tests are passing now.

* Removed use_strong_form flag in dg.cpp. Instead uses polymorphism to build operators needed for strong or fevalues needed for weak and assembles the residual functions therein. Also, corrected parameter files for advection unsteady weight adjusted curvilinear tests.

* Made a few small minor changes to operators.

* Small bug fix I saw.

* fixing Gaussian bump flow case merge commit; related tests now passing

* adding physics_model to use_aux_equation flag

* Refactoring NumericalFluxConvective

* introducing has_nonzero_diffusion bool in physics

* assigning has_nonzero_diffusion to use_auxiliary_eqn

* bug fix

* Added viscous pseudotime to converge the diffusion with explicit steady state (both weak and strong needed this). Added the max_dt_cell for strong form calculation. Verified the boundary integrals in strong form. Fully verified viscous strong form, with energy tests recovering the exact provable stability bound to machine precision. Added a central viscous flux to achieve the exact provable stability condition.

* Fixed 3D periodic nonsymmetric grid. Setup TGV for curvilinear grid. Fixed bug in weight-adjusted on-the-fly for curvilinear and now passes both 2D and 3D tests. Added condition for sum-factorization tests for cpu orders of convergence.

* improving MPI_VISCOUS_TAYLOR_GREEN_VORTEX_RESTART_FROM_PARAMETER_FILE_CHECK

* Update. Switching branches.

* Added check for convergence in Hadamard product sum-factorized test. The test is still the same but it checks if the computed error was off because it didn't consider enough polynomial points.

* This always passes.

* Speed up TGV run.

* Fixed bug with diffusion implicit convergence times.

* Cleaned up TGV test for curvilinear grid.

* Fix nonsymmetric grid for 2D advection, the 1/5 was too much warping, changed back to 1/10.

* Adding Ismail-Roe split form convective numerical flux; introducing parameter 'two_point_num_flux_type'.

* FlowSolver using set_initial_condition; can be used for curvilinear grids now

* adding safeguard for parameter: number_of_grid_elements_per_dimension

* Fixed projection IC for multi-state.

* The test parameter file for multi-state curvilinear projection.

* renaming "entropy_conserving_flux" to "two_point_flux"

* using pcout in ODE factory

* adding output_high_order_grid parameter

* adding constant_time_step parameter for TGV DNS verification

* minor merge fix

* bug fix in parameters for constant_time_step

* updating ode_solver_type in missed prm files from explicit to runge_kutta; all tests expected to pass are passing

* allowing for negative grid bounds; needed for verifying viscous TGV

* Fixing bug which adds dissipation to Euler

* Test to confirm conservation of entropy

* PR comments re: entropy check test

* Changing physics per PR comments; using isothermal TGV initialization

* Adding comment re: MPI behaviour of dot product.

* Move entropy calculation into physics

* Fixing an accidental change

* Addressing final PR comments

* fixing ctest command for TGV restart tests when MPIMAX is not a power of 2

* setting output_high_order_grid to false by default

* removing commented code in tgv shell script

* Bug fixed for auxiliary rhs test to work with general processors.

* Small fix from last commit.

* Sum-factorized Hadamard product works with general diagonal weights for other tensor directions. This enables the use of general matrices for 2pt fluxes (stiffness operator rather than derivative).

* Setup volume skew-symmetric form for DG solver. Setup surface Hadamard products using sum-factorization type algorithms with corresponding unit tests. Hadamard products handle basis of variable sizes with arbitrary diagonal products.

* Put surface two point flux in dg strong. Uncollocated updated version works perfectly for Burgers' and curvilinear. Exploring why not yet machine precision entropy conservation for Euler TGV uncollocated IR.

* PeriodicTurbulence Adaptive Time Step MPI Bug Fix (#3)

* flipping sign on viscous flux

* reverting sign flip on viscous flux

* renaming mean_specific_energy and trying equivalent pressure gradient method

* unused var fix

* removing misleading compute_dimensional_temperature

* const member initialization for sutherlands law temperatures

* bug fix for adaptive time step; using MPI max

* making temperature_inf an input parameter that must be consistent with the Prandtl number

* removing commented code in header

* doxygen fix

* Added Chandrashekar flux. General uncollocated two-point fluxes works perfectly. Verified with 3D TGV curvilinear FR c+ uncollocated entropy conservation. Found a bug with weight-adjusted curvilinear mass inverse with uncollocated and overintegration. To be fixed in future only overintegration doesn't work as expected with uncollocated curvilinear weight-adjusted.

* Clean up. Giong to merge.

* Deleted commented out debugging code.

* Cleaned up euler/physics_model as per comments. Added in strong DG to interpolate the conservative variables to the face for the numerical flux if conservative to pass the TGV quick tets. In the future, I think that it should be the interpolation of the entropy variables to the face even for conservative DG.

* Added Ranocha entropy conserving, KEP and PEP numerical flux. Added mapping support points extraction that scales on order (grid_degre+1)^dim to alleviate previous bottleneck of computing the mapping support nodes that prevnted total scaling from (p+1)^dim.

* Added test verifying that the global inverse mass matrix and global mass matrix on the fly recover same vector for multi state and any FR scheme fully veryfing the mass matrix application.

* Integrated quantities for periodic turbulence now uses sum-factorization algorithms to reduce round-off error build up and provide scaling speed-up.

* improvements and addressing PR comments

* using constexpr where possible in operators.cpp

* editorial changes

* adding abort condition for aux equation

* Replacing exit(1) by std::abort(), editorial changes to weak_dg.hpp

* correcting abort for non-diffusive problems

* cleaning up added polymorphism to DGBase (functions that are not called in DGBase and only in the derived DGStrong class); 1D, 2D, 3D_AUXILIARY_RHS_TEST are passing

* Unifying the two versions of compute_entropy_variables in Euler physics

* moving allocate_auxiliary_equation() definition to DGBase

* Changing auxiliary explicit timestepping check to NOT-implicit. This will allow ODE solvers such as rrk_explicit to solve diffusive terms.

---------

Co-authored-by: Alexander <alexander@localhost.localdomain>
Co-authored-by: Julien Brillon <43619715+jbrillon@users.noreply.github.com>
Co-authored-by: Carolyn <c.pethrick@gmail.com>
PranshulThakur pushed a commit that referenced this pull request Feb 24, 2023
* ESFR Split forms. Contains all unsteady unit and integration tests. Makes use of operators in DG strong. Provides weight-adjusted curvilinear Mass matrices. Includes free-stream preserving metric terms.

* Removed the accidentally pushed mesh file. Will rebase last 4 commits into 1.

* Cleaning up the accidental added line at end of some files and removed #if 0 #endif commented blocks.

* And this file deleted last line.

* Finished removing #if 0 #endif blocks, added some comments and included new numerical fluxes in the numerical_flux_conservation test.

* Apply suggestions from code review

dg.cpp minor edits

* Implemented sum factorization in operators class. All operations can be written as a combination of matrix-vector products and inner products. A unit test verifies by testing the cputime for A*u and also M*u where M is a mass matrix--uses both matrix-vector product and inner product forms. All that is left is to clean up the unit test a bit to fail at wrong orders.

* Apply suggestions from code review

* Apply suggestions from code review

* Apply suggestions from code review

* Update src/parameters/all_parameters.cpp

* Update with all changes requested from pull request. Only missing merge conflict with ode solver base then pull request changes complete.

* Small cleanup changes before I merge upstream master in my branch.

* Made corrections to explicit ode solver with Butcher tableau c values. Updated unsteady tests. Continued to clean up and removed deprecated Burgers Rewienski lines.

* Cleaned up last few comments on pull request wrt Flux reconstruction operators and flags.

* Removed old references to Runge-Kutta order.

* Cleaned up re-used code in operators, and other comments from PR.

* Removed collection tuple from dg. Removed finite element and quadrature collections from DG. Made consistent throughout the code to call operators finite element and quadrature collections from the operators collection tuple.

* Fixed typos in error messages.

* Created initial condition base class with interpolate and projection capabilities for different initial condition functions. This generalizes the initial condition class.

* Turned axuiliary mass matrix into a single rather than array of dim because we use the same mass matrix for each dim direction. This reduces memory consumption.

* Templated initial condition base by nstate and the functions therein.

* Added a strong form Navier-Stokes integration test.

* Made mapping support points arrays instead of vectors. Added flag in dg allocate system for the implicit matrix/vectors to allocate to reduce memory costs.

* Fix minor changes from Pranshul.

* Small update. Switching branch for other work.

* Auxiliary variable correctly computed for NAvier-Stokes equation now. Verified through viscous TGV test strong.

* Fixed small compiling error.

* Fix TGV viscous strong control file

* Strong viscous TGV energy check test confirmed pass.

* Added definitions for projection.

* Changed initial condition base to just initial condition since nothing is derived from it.

* Current status of new operators backup

* Made entropy conserving flux work for all physics types. Will return to my other branch with restructuring operators class for efficiency now.

* Recent version, new operators quarter done.

* For energy tests, only initializes ode solver on first iteration now.

* For Julien to debug his flow solver case. I will revert these changes after so that the test case is mroe reasonable when using ctests.

* Fixed the Burgers' energy inviscid initial condition.

* Fixed stash error

* Apply suggestions from code review

* cleaning up InitialCondition; changed to SetInitialCondition

* removing commented lines from previous commit

* fixing bug that would run two tgv energy check quick tests

* cleaning up src/parameter/ files; fixing comments marked as fixed that were not actually fixed

* minor editorial changes

* correcting all_parameters.cpp from last commit; indentation and doxygen fixes

* indentation clean up in unit tests and euler split tgv

* cleaning up src/testing/ files

* removing unused CurvManifold class declaration in advection_explicit_periodic.h

* changing tabs to spaces in burgers_stability.cpp

* cleaning up ConvectionDiffusionPeriodic class files in src/testing/

* reducing code repetition in CMakeList for flow_variable_tests

* reviewing src/dg, src/numerical_flux, and src/ode_solver/ files; moving all convective num fluxes to convective_numerical_flux.<cpp/hpp>; converting tabs to spaces in convective_numerical_flux.cpp; general clean up in numerical flux directory; indentation fixes

* New operators class completed with passing unit tests. Only unit tests to change atm are sum factorization, metric splitting and weight adjusted inverse. Next step is to deprecate the previous operators class and make respective changes throughout dg, tests etc.

* reviewing src/mesh/grids files

* making repeated constants class members in NonsymmetricCurvedGridManifold

* correcting "set test_type = flow_solver" to "set run_type = flow_simulation" in prm files

* The sum factorization test.

* Changed mapping shape functions class so metric operators less variables pass to construct.

* fixing time refinement study parameter file

* New operators class complete with top-nothc efficiency :)

* adding print statements for operators allocation status

* adding forgotten print statements

* comment fix

* resolving a PR comment

* Implemented new operators class throughout DG/DG strong. All that is left is having it pass the tests. It currently compiles.

* I have DG strong conservative working, now need to tweak split forms and pass all tests.

* Conservative DG works 2D curvilinear. Now need cuvilinear split form to work properly, then will do auxtesting/tgv and finito.

* Got curvilinear stability and weight-adjusted curvilinear working perfectly. Now only needs to pass auxiliary/conv_diff explicit and it's all done!

* 3D TGV works. Now just auxiliary.

* bug fix for writing pvtu files in subdirectory

* setting default value for WALE model constant to that recommended by literature

* Auxiliary and conv diffusion passes.

* Inverse mass on the fly and mass evaluated on the fly in application per cell done. Only doesnt work for curvilinear at the moment. Has flag to use. For curvilinear, do not use inverse on the fly yet. Works perfectly for linear grids.

* 2pt flux Hadamard sum after the Hadamard product now O n^d+1. All is finished.

* Made mass matrix computation more efficient. FINALLY PASSES NACA 0012 OPTIMIZATION I'm so happy!

* compiles. Have not yet tested since the merge.

* Fixed TGV Euler split form test that got stuff deleted during the merge. Also, in mass matrix DG, changed condition of factoring out det of Jac by Cartesian element rather than linear. This is because unstructured linear meshes do NOT have a constant metric Jacobian determinant, only Cartesian meshes do. Will run ctest and clean up documentation now.

* Moved the entropy consevring, entropy stable split, and central numerical fluxes into convective numerical flux.

* Changed tabs to spaces in burgers epxlicit test.

* Added documentation equations to DG strong.

* Modified TGV split test to use initial condition TGV function.

* Fixed Polynomial degree ramping issue.

* Run curvilinear mesh on Narval. Also, made 3D periodic nonsymm curv grid.

* Added energy file declaration from control file. Also increased warping of grid.

* small clean up

* fixing pde_type flag in strong dg

* correcting initial condition setting in split tgv

* Finished with curv TGV work for conference

* 1D_BURGERS_STABILITY_ORDERS_OF_ACCURACY_LONG passing

* 1D_TIME_REFINEMENT_STUDY_ADVECTION_EXPLICIT now passes and behaves as supposed to.

* 1D periodic now CONFIRMED to work with dealii matched faces. Thus, the hardcoded 1D periodic was removed. All tests that are affected by this pass.

* Viscous TGV energy check quick STRONG passing.

* Smoothed out auxiliary test to reach asymptotic convergence in L infinity norm in a reasonable time/grid level.

* Unit test verifying DG Strong Auxiliary equations' right-hand-side for d\in[1,3], for Navier-Stokes equations.

* Fixed CMake files and cleaned up the include header files in unit test.

* Resolved conflicts from merge that caused the tests to fail. Tests are passing now.

* Removed use_strong_form flag in dg.cpp. Instead uses polymorphism to build operators needed for strong or fevalues needed for weak and assembles the residual functions therein. Also, corrected parameter files for advection unsteady weight adjusted curvilinear tests.

* Made a few small minor changes to operators.

* Small bug fix I saw.

* fixing Gaussian bump flow case merge commit; related tests now passing

* adding physics_model to use_aux_equation flag

* Refactoring NumericalFluxConvective

* introducing has_nonzero_diffusion bool in physics

* assigning has_nonzero_diffusion to use_auxiliary_eqn

* bug fix

* Added viscous pseudotime to converge the diffusion with explicit steady state (both weak and strong needed this). Added the max_dt_cell for strong form calculation. Verified the boundary integrals in strong form. Fully verified viscous strong form, with energy tests recovering the exact provable stability bound to machine precision. Added a central viscous flux to achieve the exact provable stability condition.

* Fixed 3D periodic nonsymmetric grid. Setup TGV for curvilinear grid. Fixed bug in weight-adjusted on-the-fly for curvilinear and now passes both 2D and 3D tests. Added condition for sum-factorization tests for cpu orders of convergence.

* improving MPI_VISCOUS_TAYLOR_GREEN_VORTEX_RESTART_FROM_PARAMETER_FILE_CHECK

* Update. Switching branches.

* Added check for convergence in Hadamard product sum-factorized test. The test is still the same but it checks if the computed error was off because it didn't consider enough polynomial points.

* This always passes.

* Speed up TGV run.

* Fixed bug with diffusion implicit convergence times.

* Cleaned up TGV test for curvilinear grid.

* Fix nonsymmetric grid for 2D advection, the 1/5 was too much warping, changed back to 1/10.

* Adding Ismail-Roe split form convective numerical flux; introducing parameter 'two_point_num_flux_type'.

* FlowSolver using set_initial_condition; can be used for curvilinear grids now

* adding safeguard for parameter: number_of_grid_elements_per_dimension

* Fixed projection IC for multi-state.

* The test parameter file for multi-state curvilinear projection.

* renaming "entropy_conserving_flux" to "two_point_flux"

* using pcout in ODE factory

* adding set_initial_condition_method parameter

* progress on reading the file

* adding read_values_from_file_and_project to SetInitialCondition class; code compiles -- ready to be tested for reading in serial

* dhit test files

* successfully reading DHIT flow setup file and projecting values to initialize solution; need to clean up this commit

* reading files in parallel; MPI test passing

* cleaning up set_initial_condition constructor, adding parameters for reading files and DHIT flow case, cleaning up test file, adding gdown script to get files

* adding output_velocity_field()

* bug fixes

* setter/getter for time_step in FlowSolverCaseBase, outputting velocity field at hard-coded times for DHIT case

* Outputting velocity field at a general set of times via input parameters

* fixing seg fault

* successfully preparing velocity field output files

* adding output_high_order_grid parameter

* interpolating velocity field to equidistant nodes for output

* introducing parameter output_velocity_at_equidistant_nodes

* adding constant_time_step parameter for TGV DNS verification

* minor merge fix

* bug fix in parameters for constant_time_step

* updating ode_solver_type in missed prm files from explicit to runge_kutta; all tests expected to pass are passing

* allowing for negative grid bounds; needed for verifying viscous TGV

* Fixing bug which adds dissipation to Euler

* Test to confirm conservation of entropy

* PR comments re: entropy check test

* Changing physics per PR comments; using isothermal TGV initialization

* Adding comment re: MPI behaviour of dot product.

* Move entropy calculation into physics

* Fixing an accidental change

* Addressing final PR comments

* fixing ctest command for TGV restart tests when MPIMAX is not a power of 2

* setting output_high_order_grid to false by default

* removing commented code in tgv shell script

* Bug fixed for auxiliary rhs test to work with general processors.

* Small fix from last commit.

* Modify vtk output to write higher-order on straight grids

* Add time data in unsteady flow solver

* PeriodicTurbulence Adaptive Time Step MPI Bug Fix (#3)

* flipping sign on viscous flux

* reverting sign flip on viscous flux

* renaming mean_specific_energy and trying equivalent pressure gradient method

* unused var fix

* removing misleading compute_dimensional_temperature

* const member initialization for sutherlands law temperatures

* bug fix for adaptive time step; using MPI max

* making temperature_inf an input parameter that must be consistent with the Prandtl number

* removing commented code in header

* doxygen fix

* renaming parameter set_initial_condition_method to apply_initial_condition_method

* adding strain rate based dissipation rate to the integrated quantities

* writing vorticity field

* updating job_compile_PHiLiP for faster queue times

* adding script for getting flow initialization files on cluster

* adding documentation for getting the flow initialization files

* fixing initialization of adaptive time step

* error fix in last commit

* adding data table for storing the exact output times for the velocity field files

* adding output_flow_field_files_directory_name parameter

* adding option to use constant viscosity

* Printing message when negative pressure, density, or entropy is encountered

* bug fix for outputting at equidistant nodes

* adding option to output the solution files (.vtu) at the velocity field output times

* changing int to unsigned int for file indexing in flow solver

* setting proper index of current desired time for output velocity field when restarting the computation

* adding timer to flow solver unsteady component

* correcting LES filter width

* closing issue 72

* additional merge fix

* minor edit

* fixing DHIT prm files; all tests passing that should be

* adding timer start/stop message

* addressing comment on src/dg/artificial_dissipation.h

* adding check_positive functions in euler.cpp

* merge fix

* correcting comment header

* removing mpi comment for table handler

* addressing PR comments

* adding DHIT test case files input files to gitignore

* Correction to `output_flow_field` (#7)

* Initial progress towards fixing the way soln_grad was computed at equidistant nodes

* cleaning up code, removing output_velocity_field_at_equidistant_nodes parameter since always will be true

* writing vorticity for mpi test as well

---------

Co-authored-by: Alexander <alexander@localhost.localdomain>
Co-authored-by: AlexanderCicchino <alexander.cicchino@mail.mcgill.ca>
Co-authored-by: A. Cicchino <39313421+AlexanderCicchino@users.noreply.github.com>
Co-authored-by: Carolyn <c.pethrick@gmail.com>
DominicRoy2 pushed a commit to DominicRoy2/PHiLiP that referenced this pull request Mar 11, 2024
* flipping sign on viscous flux

* reverting sign flip on viscous flux

* renaming mean_specific_energy and trying equivalent pressure gradient method

* unused var fix

* removing misleading compute_dimensional_temperature

* const member initialization for sutherlands law temperatures

* bug fix for adaptive time step; using MPI max

* making temperature_inf an input parameter that must be consistent with the Prandtl number

* removing commented code in header

* doxygen fix
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

None yet

2 participants