Commit 90e8aaa8 authored by Pietro Incardona's avatar Pietro Incardona

Compiling 3d lid driven cavity

parent 4d432843
......@@ -327,7 +327,7 @@ public:
// get the position
auto key = it.get();
grid_key_dx<2> gkey = g_map.getGKey(key);
grid_key_dx<Sys_eqs::dims> gkey = g_map.getGKey(key);
// Calculate the non-zero colums
T::value(g_map,key,gs,spacing,cols,1.0);
......
......@@ -17,6 +17,7 @@
constexpr unsigned int x = 0;
constexpr unsigned int y = 1;
constexpr unsigned int z = 2;
constexpr unsigned int V = 0;
struct sys_nn
......@@ -148,11 +149,11 @@ BOOST_AUTO_TEST_CASE( der_central_non_periodic)
BOOST_REQUIRE_EQUAL(cols_x.size(),2ul);
BOOST_REQUIRE_EQUAL(cols_y.size(),2ul);
BOOST_REQUIRE_EQUAL(cols_x[17+1],1/spacing[0]);
BOOST_REQUIRE_EQUAL(cols_x[17-1],-1/spacing[0]);
BOOST_REQUIRE_EQUAL(cols_x[17+1],1/spacing[0]/2.0);
BOOST_REQUIRE_EQUAL(cols_x[17-1],-1/spacing[0]/2.0);
BOOST_REQUIRE_EQUAL(cols_y[17+16],1/spacing[1]);
BOOST_REQUIRE_EQUAL(cols_y[17-16],-1/spacing[1]);
BOOST_REQUIRE_EQUAL(cols_y[17+16],1/spacing[1]/2.0);
BOOST_REQUIRE_EQUAL(cols_y[17-16],-1/spacing[1]/2.0);
// filled colums
......@@ -173,9 +174,9 @@ BOOST_AUTO_TEST_CASE( der_central_non_periodic)
BOOST_REQUIRE_EQUAL(cols_yx.size(),4ul);
BOOST_REQUIRE_EQUAL(cols_yy.size(),3ul);
BOOST_REQUIRE_EQUAL(cols_xx[32],1/spacing[0]/spacing[0]);
BOOST_REQUIRE_EQUAL(cols_xx[34],-2/spacing[0]/spacing[0]);
BOOST_REQUIRE_EQUAL(cols_xx[36],1/spacing[0]/spacing[0]);
BOOST_REQUIRE_EQUAL(cols_xx[32],1/spacing[0]/spacing[0]/2.0/2.0);
BOOST_REQUIRE_EQUAL(cols_xx[34],-2/spacing[0]/spacing[0]/2.0/2.0);
BOOST_REQUIRE_EQUAL(cols_xx[36],1/spacing[0]/spacing[0]/2.0/2.0);
BOOST_REQUIRE_EQUAL(cols_xy[17],1/spacing[0]/spacing[1]);
BOOST_REQUIRE_EQUAL(cols_xy[19],-1/spacing[0]/spacing[1]);
......
......@@ -19,6 +19,8 @@
#include "Solvers/umfpack_solver.hpp"
#include "data_type/aggregate.hpp"
BOOST_AUTO_TEST_SUITE( eq_test_suite )
// Stokes flow
struct lid_nn
......@@ -103,8 +105,6 @@ typedef Avg<y,v_x,lid_nn,FORWARD> avg_vx_f;
#define EQ_2 1
#define EQ_3 2
BOOST_AUTO_TEST_SUITE( eq_test_suite )
// Lid driven cavity, uncompressible fluid
BOOST_AUTO_TEST_CASE(lid_driven_cavity)
......
......@@ -6,4 +6,5 @@
#include "FiniteDifference/FDScheme_unit_tests.hpp"
#include "FiniteDifference/util/common_test.hpp"
#include "FiniteDifference/eq_unit_test.hpp"
#include "FiniteDifference/eq_unit_test_3d.hpp"
#include "util/util_num_unit_tests.hpp"
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment