e r l a n g : c o o k b o o k

/ Erlang.NumberComparing

This Web


WebHome 
WebChanges 
TOC
NewRecipe 
WebTopicList 
WebStatistics 

All Webs


Chicken
Cookbook
Erlang
Know
Main
Plugins
Sandbox
Scm
TWiki  

Erlang Links


Erlang.org
Erlang Wiki
ErlMan
Erlang Wiki
The Jungerl
Erlang-fr.org
Joe Armstrong
Lambda the Ultimate

Erlang Web Ring


[Prev]: Joe Armstrong's Page
[Next]: Joe Armstrong's Page

Comparing Floating-Point Numbers

Problem

You want to compare two floating-point numbers and know if they are equal. Unfortunately, floating-point arithmetic is not precise so very few results will match exactly. Consequently, we usually want to compare floating point values up to a certain number of decimal places.

Solution

Erlang (unforunately) does not implement the full tower of numerics, which would provide built-in support for access to exact and inexact integers, bignum support, and a raft of other tools making this particular recipe somewhat unnecessary.

Consequently, we will implement a "fuzzy match" on two real numbers where the difference is below some epsilon threshhold.

In these cases, you can use floating-point byte strings to represent and compare numbers:

1> Aval = 8.001e-3 * 9.001e5.
7201.70
2> Bval = 8.0011e-3 * 9.001e5.
7201.79
3> Aval == Bval.
false
4> ABin = << Aval/float >>.
<<64,188,33,179,57,192,235,237>>
5> BBin = << Bval/float >>.
<<64,188,33,202,68,166,34,63>>
6> << ABinTest:3/binary, ARest:5/binary >> = ABin.
<<64,188,33,179,57,192,235,237>>
7> ABinTest.
<<64,188,33>>
8> << BBinTest:3/binary, BRest:5/binary >> = BBin.
<<64,188,33,202,68,166,34,63>>
9> BBinTest.
<<64,188,33>>
10> BinTest == BinTest2.
true
11> << ABinTest2:4/binary, ARest2:4/binary >> = ABin.
<<64,188,33,179,57,192,235,237>>
12> << BBinTest2:4/binary, BRest2:4/binary >> = BBin. 
<<64,188,33,202,68,166,34,63>>
13> ABinTest2 == BBinTest2.
false

Another option is to convert the numbers into strings and then compare the portions of the numbers of interest:

14> [Nsa]=io_lib:format("~.12f", [Aval]).
["7201.700099999999"]
15> [NSB]=io_lib:format("~.12f", [Bval]).  
["7201.790110000001"]
16> string:substr(NSAa, 1, 6) == string:substr(NSB, 1, 6).
true

Discussion

It would obviously be better to make use of some (imaginary) Erlang bignum and exact integer facility to carry calculations through your programs with the highest level of precision. However, Erlang does not consistently provide the level of detail necessary for certain mathematical operations. Consequently, it is frequently useful to provide "fuzzy" matching on numerical terms. Hopefully this recipe can help deal with those situations.

A final thought:

20> Equal_to_digit = fun(Digit, Numa, Numb) ->                     
20> [N_a] = io_lib:format("~.12f", [Numa]),                   
20> [N_b] = io_lib:format("~.12f", [Numb]),   
20> string:substr(N_a, 1, Digit) == string:substr(N_b,1,Digit) end.
#Fun<erl_eval.18.83214871>
21> Equal_to_digit(6, Aval, Bval).
true

Note: Some error handling would obviously be necessary to handle cases where the digits are insufficient for the match.

See Also

Volume 2, Section 4.2.2 of The Art of Computer Programming

Contributors

Based on work by BrentAFulgham and NoelWelsh.

-- BrentAFulgham - 23 Aug 2004

CookbookForm
TopicType: Recipe
ParentTopic: NumberRecipes
TopicOrder: 020

 
 
Copyright © 2004 by the contributing authors. All material on the Erlang Cookbook web site is the property of the contributing authors.
This material can be redistributed and/or modified under the terms of the GNU Lesser General Public License (LGPL), version 2.1, as published by the Free Software Foundation.
Ideas, requests, problems regarding Schematics Cookbook? Send feedback.
/ You are Main.guest