From 17e4dbbd46c80e463cef4453c1267e8c9041e4b4 Mon Sep 17 00:00:00 2001 From: Klaus Basan Date: Fri, 26 Jan 2018 04:28:33 +0100 Subject: [PATCH] Ref T231, disable ARTISTIC forammting for some areas --- src/blackmisc/simulation/interpolatorspline.cpp | 5 +++++ 1 file changed, 5 insertions(+) diff --git a/src/blackmisc/simulation/interpolatorspline.cpp b/src/blackmisc/simulation/interpolatorspline.cpp index f46813327..1aa957b47 100644 --- a/src/blackmisc/simulation/interpolatorspline.cpp +++ b/src/blackmisc/simulation/interpolatorspline.cpp @@ -27,6 +27,7 @@ namespace BlackMisc template std::array solveTridiagonal(std::array, N> &matrix, std::array &d) { + // *INDENT-OFF* const auto a = [&matrix](auto i) -> double& { return matrix[i][i-1]; }; // subdiagonal const auto b = [&matrix](auto i) -> double& { return matrix[i][i ]; }; // main diagonal const auto c = [&matrix](auto i) -> double& { return matrix[i][i+1]; }; // superdiagonal @@ -47,6 +48,7 @@ namespace BlackMisc d[i] -= c(i) * d[i+1]; } return d; + // *INDENT-ON* } //! \private Linear equation expressed as tridiagonal matrix. @@ -58,6 +60,7 @@ namespace BlackMisc std::array, N> a {{}}; std::array b {{}}; + // *INDENT-OFF* a[0][0] = 2.0 / (x[1] - x[0]); a[0][1] = 1.0 / (x[1] - x[0]); b[0] = 3.0 * (y[1] - y[0]) / ((x[1] - x[0]) * (x[1] - x[0])); @@ -74,6 +77,8 @@ namespace BlackMisc b[i] = 3.0 * (y[i] - y[i-1]) / ((x[i] - x[i-1]) * (x[i] - x[i-1])) + 3.0 * (y[i+1] - y[i]) / ((x[i+1] - x[i]) * (x[i+1] - x[i])); } + // *INDENT-ON* + solveTridiagonal(a, b); return b; }