Skip to content
Snippets Groups Projects
Commit db986e47 authored by Jean-Baptiste Keck's avatar Jean-Baptiste Keck
Browse files

modified poissonSolver main

parent f50a9732
No related branches found
No related tags found
No related merge requests found
...@@ -39,6 +39,18 @@ std::function<T(T)> func(std::size_t k) { ...@@ -39,6 +39,18 @@ std::function<T(T)> func(std::size_t k) {
} }
} }
std::string bdsToStr(domain::Boundary bd) {
switch(bd) {
case(NONE): return "None ";
case(PERIODIC) : return "Periodic ";
case(HOMOGENEOUS_NEUMANN): return "Hom_Neum. ";
case(HOMOGENEOUS_DIRICHLET): return "Hom_Diric.";
case(NEUMANN): return "Neumann ";
case(DIRICHLET): return "Dirichlet ";
}
return "";
}
template <typename T, std::size_t Dim, bool verbose=false> template <typename T, std::size_t Dim, bool verbose=false>
void test(bool includePeriodicBds=false) { void test(bool includePeriodicBds=false) {
Shape<Dim> shape; Shape<Dim> shape;
...@@ -114,8 +126,8 @@ void test(bool includePeriodicBds=false) { ...@@ -114,8 +126,8 @@ void test(bool includePeriodicBds=false) {
std::stringstream ss; std::stringstream ss;
ss << "["; ss << "[";
for (std::size_t k=0; k<Dim-1; k++) for (std::size_t k=0; k<Dim-1; k++)
ss << bdsConfig[k].first << "/" << bdsConfig[k].second << ","; ss << bdsToStr(bdsConfig[k].first) << "/" << bdsToStr(bdsConfig[k].second) << ",";
ss << bdsConfig[Dim-1].first << "/" << bdsConfig[Dim-1].second; ss << bdsToStr(bdsConfig[Dim-1].first) << "/" << bdsToStr(bdsConfig[Dim-1].second);
ss << "]"; ss << "]";
const auto dist = out.distance(ref); const auto dist = out.distance(ref);
...@@ -155,11 +167,11 @@ void test(bool includePeriodicBds=false) { ...@@ -155,11 +167,11 @@ void test(bool includePeriodicBds=false) {
int main(int argc, const char *argv[]) { int main(int argc, const char *argv[]) {
std::cout << "== TEST 1D - float =="; std::cout << "== TEST 1D - float ==";
test<float,1,true>(true); test<float,1,true>();
std::cout << "== TEST 2D - float =="; std::cout << "== TEST 2D - float ==";
test<float,2,false>(); test<float,2,true>();
std::cout << "== TEST 3D - float =="; std::cout << "== TEST 3D - float ==";
test<float,3,false>(); test<float,3,true>();
std::cout << std::endl; std::cout << std::endl;
std::cout << "== TEST 1D - double =="; std::cout << "== TEST 1D - double ==";
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment