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

Use weak form by default

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.

Instead of p+1 for entropy

Getting vortex to work. I can see it move and it diffuses super quickly

Most of the EULER_INTEGRATION test do not pass. Need to tweak some things.

Mainly to merge back into master so that other students call get the latest versions

2D_EULER_IMPLICIT_MANUFACTURED_SOLUTION now has a tolerance of 0.3 instead of 0.1

It's possible to get optimal orders on even polynomial degrees, we just need big grids

Post processor now allows the physics to specify what it wants to output

Also, the explicit ode solver will only assemble the residual without assembling the Jacobian.

Got the code to compile and run successfully for 1D linear advection weak/strong

Will solve p0, p1, ..., pn successively.
Each time interpolating from coarser to finer solution space.

Also, linear solver now uses RCM re-ordering

 ********************************************
 Convergence summary
 ********************************************
p cells     dx     L2_entropy_error   area_error
1    18 3.6459e+01 1.0503e+00     - 6.3992e+00    -
1    72 2.0509e+01 4.8038e-01  1.13 4.0986e-01 3.96
1   288 1.0740e+01 1.0588e-01  2.18 2.5774e-02 3.99
1  1152 5.4795e+00 2.1981e-02  2.27 1.6133e-03 4.00
 ********************************************
p cells     dx     L2_entropy_error   area_error
2    18 3.6459e+01 3.9484e-01     - 4.4870e-02    -
2    72 2.0509e+01 5.0201e-02  2.98 7.1469e-04 5.97
2   288 1.0740e+01 6.8890e-03  2.87 1.1221e-05 5.99
2  1152 5.4795e+00 9.4486e-04  2.87 1.7576e-07 6.00
 ********************************************
p cells     dx     L2_entropy_error   area_error
3    18 3.6459e+01 2.3019e-02     - 1.8383e-04    -
3    72 2.0509e+01 5.6555e-03  2.03 7.2973e-07 7.98
3   288 1.0740e+01 3.4905e-04  4.02 2.8499e-09 8.00
3  1152 5.4795e+00 1.9742e-05  4.14 9.0949e-11 4.97
 ********************************************
p cells     dx     L2_entropy_error    area_error
4    18 3.6459e+01 1.9260e-02     - 4.9017e-07     -
4    72 2.0509e+01 7.9555e-04  4.60 4.9749e-10  9.94
4   288 1.0740e+01 3.0078e-05  4.73 5.6843e-11  3.13
4  1152 5.4795e+00 1.5808e-06  4.25 6.9122e-11 -0.28
 ********************************************

Took too long otherwise

37:  ********************************************
37:  Convergence summary
37:  ********************************************
37: p cells DoFs     dx     L2_entropy_error
37: 1    27  432 4.8113e-02 6.4831e-02     -
37: 1   108 1728 2.4056e-02 7.2022e-03  3.17
37: 1   432 6912 1.2028e-02 6.8091e-04  3.40
37:  ********************************************
37: p cells DoFs      dx     L2_entropy_error
37: 2    27   972 3.2075e-02 4.6651e-03     -
37: 2   108  3888 1.6038e-02 6.7912e-04  2.78
37: 2   432 15552 8.0188e-03 5.1894e-05  3.71
37:  ********************************************
37: p cells DoFs      dx     L2_entropy_error
37: 3    27  1728 2.4056e-02 2.7273e-03     -
37: 3   108  6912 1.2028e-02 7.4464e-05  5.19
37: 3   432 27648 6.0141e-03 2.9142e-06  4.68
37:  ********************************************

Also add 2019_08_12_tests_outputs.log which lists all the passing tests with the orders.

Suggest opening it with vim, and searching for 'Passed'. The orders will be printed right above of the corresponding passed test.

Had to use Philip's version of the vortex where density is constant

 ********************************************
 Convergence summary
 ********************************************
p cells DoFs      dx      soln_L2_error
1   400  6400 1.2500e-02 2.4245e-01    -
1   900 14400 8.3333e-03 1.1197e-01 1.91
1  1600 25600 6.2500e-03 6.4467e-02 1.92
1  2500 40000 5.0000e-03 4.1902e-02 1.93
 ********************************************
p cells DoFs      dx      soln_L2_error
2   400 14400 8.3333e-03 2.5350e-02    -
2   900 32400 5.5556e-03 7.3640e-03 3.05
2  1600 57600 4.1667e-03 3.0701e-03 3.04
2  2500 90000 3.3333e-03 1.5792e-03 2.98
 ********************************************
p cells  DoFs      dx      soln_L2_error
3   400  25600 6.2500e-03 2.1770e-03    -
3   900  57600 4.1667e-03 4.3163e-04 3.99
3  1600 102400 3.1250e-03 1.3717e-04 3.98
3  2500 160000 2.5000e-03 5.9569e-05 3.74
 ********************************************

Takes forever. Do not run if you don't have to. Run overnight if you do.

See discussion in the following thread: https://groups.google.com/forum/#!topic/dealii/T_Ykn8KXxME

deal.II replaced it with compute_n_locally_owned_dofs_per_processor()

Mac don't have this feature. Won't be able to catch NaNs and overflows on Mac.

1D can't be run with multiple processors and therefore must use a different object.

Earlier, when using the base Triangulation, n_dofs() would thing that the mesh is just local and return the wrong number of dofs and active cells.

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