diff options
Diffstat (limited to 'Help/guide/tutorial/Step6/MathFunctions/mysqrt.cxx')
| -rw-r--r-- | Help/guide/tutorial/Step6/MathFunctions/mysqrt.cxx | 33 | 
1 files changed, 33 insertions, 0 deletions
diff --git a/Help/guide/tutorial/Step6/MathFunctions/mysqrt.cxx b/Help/guide/tutorial/Step6/MathFunctions/mysqrt.cxx new file mode 100644 index 0000000..b9ad20a --- /dev/null +++ b/Help/guide/tutorial/Step6/MathFunctions/mysqrt.cxx @@ -0,0 +1,33 @@ +#include "MathFunctions.h" +#include "TutorialConfig.h" +#include <iostream> + +#include <cmath> + +// 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 +  double result = x; + +  // 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; +}  | 
