From 81cc0094a2c7478144702e1532bcb067faaebf26 Mon Sep 17 00:00:00 2001 From: Jaron Kent-Dobias Date: Wed, 17 Feb 2021 12:08:10 +0100 Subject: Got the variation calculation fixed and working. --- stokes_test.cpp | 25 +++++++++++++++++++++++++ 1 file changed, 25 insertions(+) create mode 100644 stokes_test.cpp (limited to 'stokes_test.cpp') diff --git a/stokes_test.cpp b/stokes_test.cpp new file mode 100644 index 0000000..7d96c8c --- /dev/null +++ b/stokes_test.cpp @@ -0,0 +1,25 @@ +#include "stokes.hpp" + +using namespace std::complex_literals; + +int main() { + + Vector> z(2), dz(2), ddz(2), dH(2); + Matrix> ddH(2, 2); + + z << -0.75067 - 0.190132 * 1i, -0.625994 + 0.665987 * 1i; + dz << -0.0636149 - 0.469166 * 1i, 0.820037 - 0.449064 * 1i; + ddz << 0.55777 + 0.730164 * 1i, 0.361959 - 0.463245 * 1i; + dH << 0.967613 - 0.907519 * 1i, 0.712336 - 0.649056 * 1i; + + ddH << -0.371925 - 0.280788 * 1i, -0.163888 - 0.141297 * 1i, + -0.163888 - 0.141297 * 1i, 0.230969 + 0.942449 * 1i; + + std::cout << zDot(z, dH) << std::endl << std::endl; + + std::cout << segmentCost(z, dz, dH) << std::endl << std::endl; + + std::cout << variation(z, dz, ddz, dH, ddH) << std::endl; + + return 0; +} -- cgit v1.2.3-70-g09d2