public
Last active

obsolete: see 293169

  • Download Gist
bignum_sqrt.erl
Erlang
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64
%% square root function for integers in Erlang
%% by Kenji Rikitake <kenji.rikitake@acm.org> 26-JAN-2010
%% Distributed under MIT license at the end of the source code.
 
-module(bignum_sqrt).
 
-export([sqrt/1]).
 
%% significand digits for float estimation
-define(DIGITS, 10).
 
%% computing square root of N
%% with Newton-Raphson method
 
sqrt(N) when (is_integer(N) and (N > 0)) ->
% estimation of the root by the float calc
L = integer_to_list(N),
Len = length(L),
case Len > ?DIGITS of
true ->
Exp = Len - ?DIGITS,
M = list_to_integer(lists:sublist(L, ?DIGITS)),
XM = math:exp(math:log(M) / 2) *
math:pow(10, (Exp rem 2) / 2),
X = list_to_integer(integer_to_list(trunc(XM)) ++
lists:duplicate(Exp div 2, $0));
false ->
X = trunc(math:sqrt(N))
end,
sqrt(N, X).
 
sqrt(N, X) ->
% Newton-Raphson method of solving square root
X2 = ((N div X) + X) div 2,
case X2 =/= X of
true ->
sqrt(N, X2);
false ->
X
end.
 
%% MIT License:
%% Copyright (c) 2010 by Kenji Rikitake.
%%
%% Permission is hereby granted, free of charge, to any person
%% obtaining a copy of this software and associated documentation
%% files (the "Software"), to deal in the Software without
%% restriction, including without limitation the rights to use,
%% copy, modify, merge, publish, distribute, sublicense, and/or sell
%% copies of the Software, and to permit persons to whom the
%% Software is furnished to do so, subject to the following
%% conditions:
%%
%% The above copyright notice and this permission notice shall be
%% included in all copies or substantial portions of the Software.
%%
%% THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND,
%% EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES
%% OF MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND
%% NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT
%% HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY,
%% WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING
%% FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR
%% OTHER DEALINGS IN THE SOFTWARE.

Please sign in to comment on this gist.

Something went wrong with that request. Please try again.