From 095fa728b6277006840c9bae12d2bc44b6d688b3 Mon Sep 17 00:00:00 2001
From: Gennady Pospelov <g.pospelov@fz-juelich.de>
Date: Thu, 14 Dec 2017 18:06:37 +0100
Subject: [PATCH] Prove of principles for qLoopedTest replacement
---
Tests/UnitTests/Core/CMakeLists.txt | 1 +
.../Core/Numeric2/SpecializationTest.h | 91 +++++++++++++++++++
Tests/UnitTests/Core/Numeric2/testlist.h | 3 +
3 files changed, 95 insertions(+)
create mode 100644 Tests/UnitTests/Core/Numeric2/SpecializationTest.h
create mode 100644 Tests/UnitTests/Core/Numeric2/testlist.h
diff --git a/Tests/UnitTests/Core/CMakeLists.txt b/Tests/UnitTests/Core/CMakeLists.txt
index f2b228a4dd2..9aa0602db98 100644
--- a/Tests/UnitTests/Core/CMakeLists.txt
+++ b/Tests/UnitTests/Core/CMakeLists.txt
@@ -32,3 +32,4 @@ ADD_GTEST(Core "DataStructure" ${libs} 0)
ADD_GTEST(Core "Other" ${libs} 0)
ADD_GTEST(Core "Numeric0" ${libs} 1)
ADD_GTEST(Core "Numeric1" ${libs} 1)
+ADD_GTEST(Core "Numeric2" ${libs} 0)
diff --git a/Tests/UnitTests/Core/Numeric2/SpecializationTest.h b/Tests/UnitTests/Core/Numeric2/SpecializationTest.h
new file mode 100644
index 00000000000..b8b0d1a5ad6
--- /dev/null
+++ b/Tests/UnitTests/Core/Numeric2/SpecializationTest.h
@@ -0,0 +1,91 @@
+#include "google_test.h"
+#include "Vectors3D.h"
+#include "Complex.h"
+#include "gtest/internal/gtest-param-util.h"
+#include "HardParticles.h"
+
+const complex_t I(0,1);
+double eps_polyh = 7.5e-13;
+
+using ::testing::Values;
+using ::testing::internal::ParamGenerator;
+using ::testing::Combine;
+
+auto qlist = testing::Combine(
+ testing::Values(
+ cvector_t({ 1, 0, 0 }),
+ cvector_t({ 0, 1, 0 }),
+ cvector_t({ 0, 0, 1 }),
+ cvector_t({ 1, 1, 0 }),
+ cvector_t({ 1, 0, 1 }),
+ cvector_t({ 1, 0, 1 }),
+ cvector_t({ 1, 1, 1 })
+ ),
+ testing::Values(
+ cvector_t({ 1, 0, 0 }),
+ cvector_t({ 0, 1, 0 }),
+ cvector_t({ 0, 0, 1 }),
+ cvector_t({ 1, 1, 0 }),
+ cvector_t({ 1, 0, 1 }),
+ cvector_t({ 1, 0, 1 }),
+ cvector_t({ 1, 1, 1 })
+ ),
+ testing::Values(
+ 1e-19, 1e-17, 1e-15, 1e-13, 1e-11, 1e-9, 1e-7, 1e-5, 1e-4, 1e-3, 1e-2, .1,
+ 1., 1e1, 1e2, 1e3, 1e4 ),
+ testing::Values(
+ -1e-15, +1e-14, -1e-13*I, +1e-12*I,
+ -1e-11, +1e-10, -1e-9*I, +1e-8*I,
+ -1e-7, +1e-6, -1e-5*I, +1e-4*I,
+ -1e-3, +1e-2, -1e-1*I, +1e-1*I,
+ .9, -.99, .999, -.9999 )
+ );
+
+
+class SpecializationTest : public ::testing::Test
+{
+protected:
+ SpecializationTest() : gen(qlist) {}
+ ~SpecializationTest();
+
+ void test_all(IFormFactorBorn* p0, IFormFactorBorn* p1, double eps=1e-12) {
+ for (auto it : gen) {
+ cvector_t qdir = std::get<0>(it);
+ cvector_t qdev = std::get<1>(it);
+ double qmag = std::get<2>(it);
+ complex_t qeps = std::get<3>(it);
+ m_q = qmag * (qdir + qeps*qdev).unit();
+
+ if (skip_q(1e-99, 5e2))
+ continue;
+
+ test_ff_eq(p0, p1, eps_polyh);
+ }
+ }
+
+ void test_ff_eq(IFormFactorBorn* p0, IFormFactorBorn* p1, double eps=1e-12) {
+ complex_t f0 = p0->evaluate_for_q(m_q);
+ complex_t f1 = p1->evaluate_for_q(m_q);
+ double avge = (std::abs(f0) + std::abs(f1))/2;
+ //std::cout<<"q="<<q<<" -> "<<std::setprecision(16)<<" f0="<<f0<<", f1="<<f1<<"\n";
+ EXPECT_NEAR( real(f0), real(f1), eps*avge );
+ EXPECT_NEAR( imag(f0), imag(f1), eps*avge );
+ }
+
+ bool skip_q( double qmag_begin=1e-99, double qmag_end=1e99 ) {
+ return m_q.mag()<=qmag_begin || m_q.mag()>=qmag_end;
+ }
+
+ cvector_t m_q;
+ ParamGenerator<std::tuple<cvector_t, cvector_t, double, complex_t>> gen;
+};
+
+SpecializationTest::~SpecializationTest() = default;
+
+TEST_F(SpecializationTest, Start)
+{
+ double L = 1.5, H = .24, alpha = .6;
+ FormFactorAnisoPyramid p0(L, L, H, alpha);
+ FormFactorPyramid p1(L, H, alpha);
+ test_all(&p0, &p1, eps_polyh);
+}
diff --git a/Tests/UnitTests/Core/Numeric2/testlist.h b/Tests/UnitTests/Core/Numeric2/testlist.h
new file mode 100644
index 00000000000..7e4055c3006
--- /dev/null
+++ b/Tests/UnitTests/Core/Numeric2/testlist.h
@@ -0,0 +1,3 @@
+// To renew this file, run /G/ba/dev-tools/code-tools/update-gtestlist.py <directory>
+
+#include "SpecializationTest.h"
--
GitLab