main.cpp

// #include "Tests/Test_1D.h"
// #include "Tests/Test_2D.h"
// #include "Tests/Test_3D.h"
// #include "Tests/Test_3DV.h"
#include "Tests/Test_Vortex_Rings.h"

static const bool Export_VTK = true;

int main()
{

    // int nx = 8*100;
    // int ny = 8*16;
    // int nz = 8*16;
    int nx = 32;
    int ny = 32;
    int nz = 32;

    // --- 1D solver tests
    // Test_Dirichlet_1D(nx);
    // Test_Dirichlet_1D_IHBC(nx);          // Not working VKFFT->
    // Test_Neumann_1D(nx);
    // Test_Periodic_1D(nx);
    // Test_Unbounded_1D(nx);

    //--- 2D scalar solver tests
    // Test_Dirichlet_2D(nx,ny);
    // Test_Dirichlet_2D_IHBC(nx,ny);       // Not working VKFFT->
    // Test_Neumann_2D(nx,ny);
    // Test_Neumann_2D_IHBC(nx,ny);         // Not working VKFFT->
    // Test_Periodic_2D(nx,ny);
    // Test_Unbounded_2D(nx,ny);

    //--- 3D scalar solver tests
    // Test_Dirichlet_3D(nx,ny,nz);
    // Test_Dirichlet_3D_IHBC(nx,ny,nz);    // Not working VKFFT->
    // Test_Neumann_3D(nx,ny,nz);
    // Test_Neumann_3D_IHBC(nx,ny,nz);      // Not working VKFFT->
    // Test_Periodic_3D(nx,ny,nz);
    // Test_Unbounded_3D(nx,ny,nz);

    //--- 3D vector solver tests
    // Test_Dirichlet_3DV(nx,ny,nz);
    // Test_Neumann_3DV(nx,ny,nz);
    // Test_Periodic_3DV(nx,ny,nz);        // Not working VKFFT-> Convolution?
    // Test_Unbounded_3DV(nx,ny,nz);
    // Test_Unbounded_3DV_Curl(nx,ny,nz);

    // --- VPM tests
    Vortex_Ring_Evolution_Test();

    return 0;
}