summaryrefslogtreecommitdiffstats
path: root/Help/guide/tutorial/Step11/MathFunctions/mysqrt.cxx
diff options
context:
space:
mode:
authorBrad King <brad.king@kitware.com>2019-06-25 11:05:18 (GMT)
committerKitware Robot <kwrobot@kitware.com>2019-06-25 11:05:35 (GMT)
commit2e8b8eb857037a4c0f08a1c7c0ce5a778ebfbc54 (patch)
treec3c8880699a5237170dd122d8987b8a99a487a34 /Help/guide/tutorial/Step11/MathFunctions/mysqrt.cxx
parentedf232dbe487e36cec2843c9495df9914b41ae28 (diff)
parent8ab8563533e57f4f4af769daf4590c6d262e5c41 (diff)
downloadCMake-2e8b8eb857037a4c0f08a1c7c0ce5a778ebfbc54.zip
CMake-2e8b8eb857037a4c0f08a1c7c0ce5a778ebfbc54.tar.gz
CMake-2e8b8eb857037a4c0f08a1c7c0ce5a778ebfbc54.tar.bz2
Merge topic 'help-guide-tutorial'
8ab8563533 Help: Add tutorial guide table of contents 0a6c9c417a Help: Update tutorial to include text from CMake book eef3e020c2 Help: Populate tutorial guide text 862cfc0e6c Help/guide/tutorial: Adopt tutorial code d2fde94809 Help: Add infrastructure for guide-level documentation f067af8e38 Tests/Tutorial: Remove trailing blank lines d5a81ef747 Tests/Tutorial: Remove unused file Acked-by: Kitware Robot <kwrobot@kitware.com> Merge-request: !3439
Diffstat (limited to 'Help/guide/tutorial/Step11/MathFunctions/mysqrt.cxx')
-rw-r--r--Help/guide/tutorial/Step11/MathFunctions/mysqrt.cxx45
1 files changed, 45 insertions, 0 deletions
diff --git a/Help/guide/tutorial/Step11/MathFunctions/mysqrt.cxx b/Help/guide/tutorial/Step11/MathFunctions/mysqrt.cxx
new file mode 100644
index 0000000..96d9421
--- /dev/null
+++ b/Help/guide/tutorial/Step11/MathFunctions/mysqrt.cxx
@@ -0,0 +1,45 @@
+#include "MathFunctions.h"
+#include <iostream>
+
+// include the generated table
+#include "Table.h"
+
+#include <cmath>
+
+namespace mathfunctions {
+namespace detail {
+// a hack square root calculation using simple operations
+double mysqrt(double x)
+{
+ if (x <= 0) {
+ return 0;
+ }
+
+ // if we have both log and exp then use them
+#if defined(HAVE_LOG) && defined(HAVE_EXP)
+ double result = exp(log(x) * 0.5);
+ std::cout << "Computing sqrt of " << x << " to be " << result << " using log"
+ << std::endl;
+#else
+ // use the table to help find an initial value
+ double result = x;
+ if (x >= 1 && x < 10) {
+ result = sqrtTable[static_cast<int>(x)];
+ }
+
+ // if we have both log and exp then use them
+
+ // do ten iterations
+ for (int i = 0; i < 10; ++i) {
+ if (result <= 0) {
+ result = 0.1;
+ }
+ double delta = x - (result * result);
+ result = result + 0.5 * delta / result;
+ std::cout << "Computing sqrt of " << x << " to be " << result << std::endl;
+ }
+#endif
+ return result;
+}
+}
+}