diff options
| author | Brad King <brad.king@kitware.com> | 2019-06-25 11:05:18 (GMT) | 
|---|---|---|
| committer | Kitware Robot <kwrobot@kitware.com> | 2019-06-25 11:05:35 (GMT) | 
| commit | 2e8b8eb857037a4c0f08a1c7c0ce5a778ebfbc54 (patch) | |
| tree | c3c8880699a5237170dd122d8987b8a99a487a34 /Help/guide/tutorial/Step8/MathFunctions/mysqrt.cxx | |
| parent | edf232dbe487e36cec2843c9495df9914b41ae28 (diff) | |
| parent | 8ab8563533e57f4f4af769daf4590c6d262e5c41 (diff) | |
| download | CMake-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/Step8/MathFunctions/mysqrt.cxx')
| -rw-r--r-- | Help/guide/tutorial/Step8/MathFunctions/mysqrt.cxx | 42 | 
1 files changed, 42 insertions, 0 deletions
diff --git a/Help/guide/tutorial/Step8/MathFunctions/mysqrt.cxx b/Help/guide/tutorial/Step8/MathFunctions/mysqrt.cxx new file mode 100644 index 0000000..5b862fb --- /dev/null +++ b/Help/guide/tutorial/Step8/MathFunctions/mysqrt.cxx @@ -0,0 +1,42 @@ +#include "MathFunctions.h" +#include "TutorialConfig.h" +#include <iostream> + +// include the generated table +#include "Table.h" + +#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 +  // 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; +}  | 
