summaryrefslogtreecommitdiff
path: root/tests/equilibrium.cpp
diff options
context:
space:
mode:
Diffstat (limited to 'tests/equilibrium.cpp')
-rw-r--r--tests/equilibrium.cpp40
1 files changed, 0 insertions, 40 deletions
diff --git a/tests/equilibrium.cpp b/tests/equilibrium.cpp
deleted file mode 100644
index 9201e55..0000000
--- a/tests/equilibrium.cpp
+++ /dev/null
@@ -1,40 +0,0 @@
-#include <forstio/test/suite.hpp>
-
-#include "../c++/equilibrium.hpp"
-
-
-namespace {
-
-template<typename Descriptor>
-void check_equilibrium(){
- using namespace kel;
-
- using dfi = lbm::df_info<lbm::sch::Float64,Descriptor>;
-
- saw::data<lbm::sch::Float64> rho{1.0};
- saw::data<lbm::sch::FixedArray<lbm::sch::Float64,Descriptor::D>> vel;
- for(saw::data<lbm::sch::UInt64> i{0u}; i.get() < Descriptor::D; ++i){
- vel.at(i) = {0.0};
- }
- auto eq = lbm::equilibrium<lbm::sch::Float64,Descriptor>(rho,vel);
-
- for(saw::data<lbm::sch::UInt64> i{0u}; i.get() < Descriptor::Q; ++i){
- SAW_EXPECT(eq.at(i).get() == dfi::weights[i.get()], std::string{"No velocity and normalized rho should be exactly the weights: "} + std::to_string(eq.at(i).get()) + std::string{" "} + std::to_string(dfi::weights[i.get()]));
- }
-}
-
-SAW_TEST("Equilibrium at rest D1Q3"){
- using namespace kel;
- check_equilibrium<lbm::sch::Descriptor<1,3>>();
-}
-
-SAW_TEST("Equilibrium at rest D2Q5"){
- using namespace kel;
- check_equilibrium<lbm::sch::Descriptor<2,5>>();
-}
-
-SAW_TEST("Equilibrium at rest D2Q9"){
- using namespace kel;
- check_equilibrium<lbm::sch::Descriptor<2,9>>();
-}
-}