Homec4science
Diffusion GooseFEM (rstcheck) Browse

rGOOSEFEM/tests/basic/Iterate.cpp
f762324f0857rstcheck

Last Change
Raw File

Iterate.cpp
View Options

#include <catch2/catch.hpp>
#include <xtensor/xrandom.hpp>
#include <xtensor/xmath.hpp>
#include <Eigen/Eigen>
#include <GooseFEM/GooseFEM.h>
#define ISCLOSE(a,b) REQUIRE_THAT((a), Catch::WithinAbs((b), 1.e-12));
TEST_CASE("GooseFEM::Iterate", "Iterate.h")
{
SECTION("StopList::all_less - sorted input")
{
GooseFEM::Iterate::StopList residuals(5);
std::vector<double> res = {5e+0, 5e+1, 5e-1, 5e-2, 5e-3, 5e-4, 4e-4, 3e-4, 2e-4, 1e-4};
std::vector<bool> conv = {false, false, false, false, false, false, false, false, false, true};
for (size_t i = 0; i < res.size(); ++i) {
residuals.roll_insert(res[i]);
REQUIRE(residuals.all_less(1e-3) == conv[i]);
}
}
SECTION("StopList::all_less - unsorted input")
{
GooseFEM::Iterate::StopList residuals(5);
std::vector<double> res = {5e+0, 5e+1, 5e-1, 5e-2, 5e-3, 5e-4, 4e-4, 5e-4, 4e-4, 4e-4};
std::vector<bool> conv = {false, false, false, false, false, false, false, false, false, true};
for (size_t i = 0; i < res.size(); ++i) {
residuals.roll_insert(res[i]);
REQUIRE(residuals.all_less(1e-3) == conv[i]);
}
}
SECTION("StopList::descending && all_less - sorted input")
{
GooseFEM::Iterate::StopList residuals(5);
std::vector<double> res = {5e+0, 5e+1, 5e-1, 5e-2, 5e-3, 5e-4, 4e-4, 3e-4, 2e-4, 1e-4};
std::vector<bool> conv = {false, false, false, false, false, false, false, false, false, true};
for (size_t i = 0; i < res.size(); ++i) {
residuals.roll_insert(res[i]);
REQUIRE((residuals.descending() && residuals.all_less(1e-3)) == conv[i]);
}
}
SECTION("StopList::descending && all_less - unsorted input")
{
GooseFEM::Iterate::StopList residuals(5);
std::vector<double> res = {5e+0, 5e+1, 5e-1, 5e-2, 5e-3, 5e-4, 4e-4, 3e-4, 2e-4, 3e-4, 2e-4, 1e-4, 9e-5, 8e-5};
std::vector<bool> conv = {false, false, false, false, false, false, false, false, false, false, false, false, false, true};
for (size_t i = 0; i < res.size(); ++i) {
residuals.roll_insert(res[i]);
REQUIRE((residuals.descending() && residuals.all_less(1e-3)) == conv[i]);
}
}
SECTION("StopList::stop_simple - sorted input")
{
GooseFEM::Iterate::StopList stop(5);
std::vector<double> res = {5e+0, 5e+1, 5e-1, 5e-2, 5e-3, 5e-4, 4e-4, 3e-4, 2e-4, 1e-4};
std::vector<bool> conv = {false, false, false, false, false, false, false, false, false, true};
for (size_t i = 0; i < res.size(); ++i) {
REQUIRE(stop.stop_simple(res[i], 1e-3) == conv[i]);
}
}
SECTION("StopList::stop_simple - unsorted input")
{
GooseFEM::Iterate::StopList a(5);
GooseFEM::Iterate::StopList b(5);
std::vector<double> res = {5e+0, 5e+1, 5e-1, 5e-2, 5e-3, 5e-4, 4e-4, 5e-4, 4e-4, 4e-4};
std::vector<bool> conv = {false, false, false, false, false, false, false, false, false, true};
for (size_t i = 0; i < res.size(); ++i) {
REQUIRE(a.stop_simple(res[i], 1e-3) == conv[i]);
REQUIRE(b.stop(res[i], 1e-3) == false);
}
}
SECTION("StopList::stop - sorted input")
{
GooseFEM::Iterate::StopList stop(5);
std::vector<double> res = {5e+0, 5e+1, 5e-1, 5e-2, 5e-3, 5e-4, 4e-4, 3e-4, 2e-4, 1e-4};
std::vector<bool> conv = {false, false, false, false, false, false, false, false, false, true};
for (size_t i = 0; i < res.size(); ++i) {
REQUIRE(stop.stop(res[i], 1e-3) == conv[i]);
}
}
SECTION("StopList::stop - unsorted input")
{
GooseFEM::Iterate::StopList stop(5);
std::vector<double> res = {5e+0, 5e+1, 5e-1, 5e-2, 5e-3, 5e-4, 4e-4, 3e-4, 2e-4, 3e-4, 2e-4, 1e-4, 9e-5, 8e-5};
std::vector<bool> conv = {false, false, false, false, false, false, false, false, false, false, false, false, false, true};
for (size_t i = 0; i < res.size(); ++i) {
REQUIRE(stop.stop(res[i], 1e-3) == conv[i]);
}
}
}
c4science · Help