DNDSR 0.2.1
Distributed Numeric Data Structure for CFV
Loading...
Searching...
No Matches
Classes | Macros | Typedefs | Enumerations | Functions | Variables
test_Reconstruction.cpp File Reference

Phase-0 regression tests for CFV Variational Reconstruction. More...

#include "doctest.h"
#include "CFV/VariationalReconstruction.hpp"
#include <array>
#include <cmath>
#include <iostream>
#include <iomanip>
#include <functional>
#include <nlohmann/json.hpp>
#include <map>
#include <string>
#include <vector>
#include "CFV/VariationalReconstruction_LimiterProcedure.hxx"
Include dependency graph for test_Reconstruction.cpp:

Go to the source code of this file.

Classes

struct  LimiterTestCase
 

Macros

#define DOCTEST_CONFIG_IMPLEMENT
 
#define POLY_TEST(testName, method, polyFunc, polyDeg)
 

Typedefs

using tVR = CFV::VariationalReconstruction< g_dim >
 
using ScalarFunc = std::function< DNDS::real(const tPoint &)>
 

Enumerations

enum class  RecMethod {
  GaussGreen , VFV_P1_HQM , VFV_P2_HQM , VFV_P3_HQM ,
  VFV_P1_Default , VFV_P2_Default , VFV_P3_Default
}
 

Functions

int main (int argc, char **argv)
 
 POLY_TEST ("const", GaussGreen, [](const tPoint &) { return 1.0;}, 0) POLY_TEST(" const "
 
const tPoint POLY_TEST ("const", VFV_P2_HQM, [](const tPoint &) { return 1.0;}, 0) POLY_TEST(" const "
 
const tPoint const tPoint POLY_TEST ("const", VFV_P1_Default, [](const tPoint &) { return 1.0;}, 0) POLY_TEST("linear"
 
 POLY_TEST ("linear", VFV_P1_HQM, [](const tPoint &p) { return p[0]+2 *p[1];}, 1) POLY_TEST("linear"
 
 POLY_TEST ("linear", VFV_P3_HQM, [](const tPoint &p) { return p[0]+2 *p[1];}, 1) POLY_TEST("linear"
 
 POLY_TEST ("quad", VFV_P2_HQM, [](const tPoint &p) { return p[0] *p[0]+p[1] *p[1];}, 2) POLY_TEST("quad"
 
 POLY_TEST ("quad", VFV_P2_Default, [](const tPoint &p) { return p[0] *p[0]+p[1] *p[1];}, 2) POLY_TEST("cubic"
 
 POLY_TEST ("cubic", VFV_P3_Default, [](const tPoint &p) { return p[0] *p[0] *p[0]+p[0] *p[1] *p[1];}, 3) struct PeriodicTestCase
 
 TEST_CASE ("Periodic reconstruction convergence series")
 
 TEST_CASE ("VFV P2 HQM converges on IV10 base mesh")
 
 TEST_CASE ("DEBUG compare InterpolateFace vs Legacy face2cell")
 
 TEST_CASE ("Limiter procedure: reconstruction + smooth indicator + WBAP limiter")
 

Variables

 VFV_P1_HQM
 
const tPoint VFV_P3_HQM
 
const tPoint const tPoint GaussGreen
 
const tPoint const tPoint const tPoint & p { return p[0] + 2 * p[1]
 
 VFV_P2_HQM
 
 VFV_P1_Default
 

Detailed Description

Phase-0 regression tests for CFV Variational Reconstruction.

Parameterized over [mesh+function, reconstruction_method].

All iterative VR uses Jacobi iteration (SORInstead=false) so that golden values are deterministic across MPI partitionings (np=1,2,4).

Error metric: L1 pointwise error at 6th-degree quadrature points, divided by domain volume (so the dimension matches the field function).

Meshes:

  • Uniform_3x3_wall (9 quads, wall BC) – polynomial exactness tests
  • IV10_10 (100 quads, periodic) + bisections -> 400, 1600 cells
  • IV10U_10 (322 tris, periodic) + bisections -> ~1288, ~5152 cells

Golden values captured from commit c774b89 on dev/harry_refac1.

Definition in file test_Reconstruction.cpp.

Macro Definition Documentation

◆ DOCTEST_CONFIG_IMPLEMENT

#define DOCTEST_CONFIG_IMPLEMENT

Definition at line 21 of file test_Reconstruction.cpp.

◆ POLY_TEST

#define POLY_TEST (   testName,
  method,
  polyFunc,
  polyDeg 
)
Value:
TEST_CASE("Wall/" testName "/" #method) \
{ \
ScalarFunc f = polyFunc; \
auto vr = buildVR(g_wall_mesh, RecMethod::method); \
auto bc = makeDirichletBC(f); \
DNDS::real err = runTest(vr, RecMethod::method, f, bc, 100, 1e-15, false); \
if (g_mpi.rank == 0) \
std::cout << "[Wall/" testName "/" #method "] err = " \
<< std::scientific << err << std::endl; \
if (polyDeg == 0) \
CHECK(err < 1e-12); \
else \
CHECK(err < 10.0); /* golden-value regression checked separately */ \
}
DNDS::real
double real
Canonical floating-point scalar used throughout DNDSR (double precision).
Definition Defines.hpp:110
DNDS::MPIInfo::rank
int rank
This rank's 0-based index within comm (-1 until initialised).
Definition MPI.hpp:235
err
real err
Definition test_Direct.cpp:315
TEST_CASE
TEST_CASE("3D: VFV P2 HQM error < P1 on sinCos3D")
Definition test_Reconstruction3D.cpp:509
ScalarFunc
std::function< DNDS::real(const tPoint &)> ScalarFunc
Definition test_Reconstruction.cpp:232

Definition at line 454 of file test_Reconstruction.cpp.

Typedef Documentation

◆ ScalarFunc

using ScalarFunc = std::function<DNDS::real(const tPoint &)>

Definition at line 232 of file test_Reconstruction.cpp.

◆ tVR

using tVR = CFV::VariationalReconstruction<g_dim>

Definition at line 41 of file test_Reconstruction.cpp.

Enumeration Type Documentation

◆ RecMethod

enum class RecMethod
strong
Enumerator
GaussGreen 
VFV_P1_HQM 
VFV_P2_HQM 
VFV_P3_HQM 
VFV_P1_Default 
VFV_P2_Default 
VFV_P3_Default 

Definition at line 128 of file test_Reconstruction.cpp.

Function Documentation

◆ main()

int main ( int  argc,
char **  argv 
)

Definition at line 399 of file test_Reconstruction.cpp.

Here is the call graph for this function:

◆ POLY_TEST() [1/8]

POLY_TEST ( "const"  ,
GaussGreen  ,
[] (const tPoint &) { return 1.0;}  ,
0   
) const

◆ POLY_TEST() [2/8]

const tPoint const tPoint POLY_TEST ( "const"  ,
VFV_P1_Default  ,
[] (const tPoint &) { return 1.0;}  ,
0   
)

◆ POLY_TEST() [3/8]

const tPoint POLY_TEST ( "const"  ,
VFV_P2_HQM  ,
[] (const tPoint &) { return 1.0;}  ,
0   
) const

◆ POLY_TEST() [4/8]

POLY_TEST ( "cubic"  ,
VFV_P3_Default  ,
[] (const tPoint &p) { return p[0] *p[0] *p[0]+p[0] *p[1] *p[1];}  ,
3   
)

Definition at line 505 of file test_Reconstruction.cpp.

◆ POLY_TEST() [5/8]

POLY_TEST ( "linear"  ,
VFV_P1_HQM  ,
[] (const tPoint &p) { return p[0]+2 *p[1];}  ,
1   
)

◆ POLY_TEST() [6/8]

POLY_TEST ( "linear"  ,
VFV_P3_HQM  ,
[] (const tPoint &p) { return p[0]+2 *p[1];}  ,
1   
)

◆ POLY_TEST() [7/8]

POLY_TEST ( "quad"  ,
VFV_P2_Default  ,
[] (const tPoint &p) { return p[0] *p[0]+p[1] *p[1];}  ,
2   
)

◆ POLY_TEST() [8/8]

POLY_TEST ( "quad"  ,
VFV_P2_HQM  ,
[] (const tPoint &p) { return p[0] *p[0]+p[1] *p[1];}  ,
2   
)

◆ TEST_CASE() [1/4]

TEST_CASE ( "DEBUG compare InterpolateFace vs Legacy face2cell"  )

Definition at line 653 of file test_Reconstruction.cpp.

Here is the call graph for this function:

◆ TEST_CASE() [2/4]

TEST_CASE ( "Limiter procedure: reconstruction + smooth indicator + WBAP limiter"  )

Definition at line 961 of file test_Reconstruction.cpp.

Here is the call graph for this function:

◆ TEST_CASE() [3/4]

TEST_CASE ( "Periodic reconstruction convergence series"  )

Definition at line 552 of file test_Reconstruction.cpp.

Here is the call graph for this function:

◆ TEST_CASE() [4/4]

TEST_CASE ( "VFV P2 HQM converges on IV10 base mesh"  )

Definition at line 589 of file test_Reconstruction.cpp.

Here is the call graph for this function:

Variable Documentation

◆ GaussGreen

const tPoint const tPoint GaussGreen

Definition at line 483 of file test_Reconstruction.cpp.

◆ p

const tPoint & p { return p[0] + 2 * p[1]

Definition at line 483 of file test_Reconstruction.cpp.

◆ VFV_P1_Default

VFV_P1_Default

Definition at line 491 of file test_Reconstruction.cpp.

◆ VFV_P1_HQM

VFV_P1_HQM

Definition at line 473 of file test_Reconstruction.cpp.

◆ VFV_P2_HQM

VFV_P2_HQM

Definition at line 487 of file test_Reconstruction.cpp.

◆ VFV_P3_HQM

VFV_P3_HQM

Definition at line 477 of file test_Reconstruction.cpp.