summaryrefslogtreecommitdiffstats
path: root/Help/guide/tutorial/Step8/MathFunctions/mysqrt.cxx
diff options
context:
space:
mode:
Diffstat (limited to 'Help/guide/tutorial/Step8/MathFunctions/mysqrt.cxx')
-rw-r--r--Help/guide/tutorial/Step8/MathFunctions/mysqrt.cxx25
1 files changed, 14 insertions, 11 deletions
diff --git a/Help/guide/tutorial/Step8/MathFunctions/mysqrt.cxx b/Help/guide/tutorial/Step8/MathFunctions/mysqrt.cxx
index 7d80ee9..28ab042 100644
--- a/Help/guide/tutorial/Step8/MathFunctions/mysqrt.cxx
+++ b/Help/guide/tutorial/Step8/MathFunctions/mysqrt.cxx
@@ -1,10 +1,10 @@
-#include <iostream>
-
-#include "MathFunctions.h"
+#include "mysqrt.h"
-// include the generated table
-#include "Table.h"
+#include <cmath>
+#include <iostream>
+namespace mathfunctions {
+namespace detail {
// a hack square root calculation using simple operations
double mysqrt(double x)
{
@@ -12,12 +12,13 @@ double mysqrt(double x)
return 0;
}
- // use the table to help find an initial value
+ // if we have both log and exp then use them
+#if defined(HAVE_LOG) && defined(HAVE_EXP)
+ double result = std::exp(std::log(x) * 0.5);
+ std::cout << "Computing sqrt of " << x << " to be " << result
+ << " using log and exp" << std::endl;
+#else
double result = x;
- if (x >= 1 && x < 10) {
- std::cout << "Use the table to help find an initial value " << std::endl;
- result = sqrtTable[static_cast<int>(x)];
- }
// do ten iterations
for (int i = 0; i < 10; ++i) {
@@ -28,6 +29,8 @@ double mysqrt(double x)
result = result + 0.5 * delta / result;
std::cout << "Computing sqrt of " << x << " to be " << result << std::endl;
}
-
+#endif
return result;
}
+}
+}