92cd935d4bbaff0f3fa3061a90f35a3eb4c1a043
[olsrd.git] / src / mantissa.c
1 /*
2  * OLSR ad-hoc routing table management protocol
3  * Copyright (C) 2004 Andreas T√łnnesen (andreto@ifi.uio.no)
4  *
5  * This file is part of the olsr.org OLSR daemon.
6  *
7  * olsr.org is free software; you can redistribute it and/or modify
8  * it under the terms of the GNU General Public License as published by
9  * the Free Software Foundation; either version 2 of the License, or
10  * (at your option) any later version.
11  *
12  * olsr.org is distributed in the hope that it will be useful,
13  * but WITHOUT ANY WARRANTY; without even the implied warranty of
14  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
15  * GNU General Public License for more details.
16  *
17  * You should have received a copy of the GNU General Public License
18  * along with olsr.org; if not, write to the Free Software
19  * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
20  * 
21  * 
22  * $ Id $
23  *
24  */
25
26
27 #include "mantissa.h"
28 #include "math.h"
29
30 /**
31  *Function that converts a double to a mantissa/exponent
32  *product as described in RFC3626:
33  *
34  * value = C*(1+a/16)*2^b [in seconds]
35  *
36  *  where a is the integer represented by the four highest bits of the
37  *  field and b the integer represented by the four lowest bits of the
38  *  field.
39  *
40  *@param interval the time interval to process
41  *
42  *@return a 8-bit mantissa/exponent product
43  */
44
45 olsr_u8_t
46 double_to_me(double interval)
47 {
48   int a, b;
49
50   b = 0;
51
52   while(interval / VTIME_SCALE_FACTOR >= pow((double)2, (double)b))
53     b++;
54
55   b--;
56   if(b < 0)
57     {
58       a = 1;
59       b = 0;
60     } 
61   else 
62     if (b > 15)
63       {
64         a = 15;
65         b = 15;
66       } 
67     else 
68       { 
69         a = (int)(16*((double)interval/(VTIME_SCALE_FACTOR*(double)pow(2,b))-1));
70         while(a >= 16)
71           {
72             a -= 16;
73             b++;
74           }
75       }
76
77   //printf("Generated mantissa/exponent(%d/%d): %d from %f\n", a, b, (olsr_u8_t) (a*16+b), interval);
78   //printf("Resolves back to: %f\n", me_to_double((olsr_u8_t) (a*16+b)));
79   return (olsr_u8_t) (a*16+b);
80 }
81
82
83
84
85 /**
86  *Function that converts a mantissa/exponent 8bit value back
87  *to double as described in RFC3626:
88  *
89  * value = C*(1+a/16)*2^b [in seconds]
90  *
91  *  where a is the integer represented by the four highest bits of the
92  *  field and b the integer represented by the four lowest bits of the
93  *  field.
94  *
95  *@param me the 8 bit mantissa/exponen value
96  *
97  *@return a double value
98  */
99 double
100 me_to_double(olsr_u8_t me)
101 {
102   int a, b;
103
104   a = me>>4;
105   b = me - a*16;
106
107   return (double)(VTIME_SCALE_FACTOR*(1+(double)a/16)*(double)pow(2,b));
108 }