How is the square root function implemented? [closed]

by

Simple implementation using Binary Search with C++

double root(double n){
  // Max and min are used to take into account numbers less than 1
  double lo = min(1, n), hi = max(1, n), mid;

  // Update the bounds to be off the target by a factor of 10
  while(100 * lo * lo < n) lo *= 10;
  while(0.01 * hi * hi > n) hi *= 0.1;

  for(int i = 0 ; i < 100 ; i++){
      mid = (lo+hi)/2;
      if(mid*mid == n) return mid;
      if(mid*mid > n) hi = mid;
      else lo = mid;
  }
  return mid;
}

Note that while loop is most common with the binary search but Personally I prefer using for when dealing with decimal numbers, it saves some special cases handling and gets pretty accurate result from small loops like that 1000 or even 500 (Both will give the same result for almost all numbers but just to be safe).

Edit: Check out this Wikipedia article for various -special purpose- methods specialized in computing the square root.

Edit 2: Apply updates suggested by @jorgbrown to fix the function in case of input less than 1. Also, apply optimization to make the bounds off the target root by a factor of 10

More Related Contents:

  1. Writing your own square root function
  2. how can i increase the speed of my C++ calculation [closed]
  3. How to create a custom EL function to invoke a static method?
  4. Understanding how recursive functions work
  5. How could I create a function with a completion handler in Swift?
  6. Difference between . and : in Lua
  7. Value receiver vs. pointer receiver
  8. Correct use of modules, subroutines and functions in Fortran
  9. Passing external function of multiple variables as a function of one variable in Fortran
  10. Jump to function definition
  11. How to write a Rust function that takes an iterator?
  12. Does Fortran preserve the value of internal variables through function and subroutine calls?
  13. Is finding the equivalence of two functions undecidable?
  14. How to deal with name/value pairs of function arguments in MATLAB
  15. How to compose `not` with a function of arbitrary arity?
  16. Uses for Haskell id function
  17. Return value from a VBScript function
  18. Is there a way to get a vector with the name of all functions that one could use in R?
  19. F# generics / function overloading syntax
  20. How to correct “Function definitions are not permitted at the prompt or in scripts”
  21. What are the differences between functions and methods in Swift?
  22. How to avoid SQL injection in CodeIgniter?
  23. No function matches the given name and argument types
  24. Return multiple values from a function in swift
  25. Why can’t ‘kotlin.Result’ be used as a return type?
  26. How many lines of code should a function/procedure/method have? [duplicate]
  27. How to evaluate functions in GDB?
  28. Proper use of modules in Fortran
  29. Matlab ode45. How to change a parameter inside it while calling it?
  30. How do you strip the unit from any number in SASS?

Leave a Comment