Licence change... man that was boring work...
[olsrd.git] / src / mantissa.c
1 /*
2  * The olsr.org Optimized Link-State Routing daemon(olsrd)
3  * Copyright (c) 2004, Andreas T√łnnesen(andreto@olsr.org)
4  * All rights reserved.
5  *
6  * Redistribution and use in source and binary forms, with or without 
7  * modification, are permitted provided that the following conditions 
8  * are met:
9  *
10  * * Redistributions of source code must retain the above copyright 
11  *   notice, this list of conditions and the following disclaimer.
12  * * Redistributions in binary form must reproduce the above copyright 
13  *   notice, this list of conditions and the following disclaimer in 
14  *   the documentation and/or other materials provided with the 
15  *   distribution.
16  * * Neither the name of olsr.org, olsrd nor the names of its 
17  *   contributors may be used to endorse or promote products derived 
18  *   from this software without specific prior written permission.
19  *
20  * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS 
21  * "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT 
22  * LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS 
23  * FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE 
24  * COPYRIGHT OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, 
25  * INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, 
26  * BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; 
27  * LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER 
28  * CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT 
29  * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN 
30  * ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE 
31  * POSSIBILITY OF SUCH DAMAGE.
32  *
33  * Visit http://www.olsr.org for more information.
34  *
35  * If you find this software useful feel free to make a donation
36  * to the project. For more information see the website or contact
37  * the copyright holders.
38  *
39  * $Id: mantissa.c,v 1.6 2004/11/21 11:28:56 kattemat Exp $
40  */
41
42
43 #include "mantissa.h"
44 #include "math.h"
45
46 /**
47  *Function that converts a double to a mantissa/exponent
48  *product as described in RFC3626:
49  *
50  * value = C*(1+a/16)*2^b [in seconds]
51  *
52  *  where a is the integer represented by the four highest bits of the
53  *  field and b the integer represented by the four lowest bits of the
54  *  field.
55  *
56  *@param interval the time interval to process
57  *
58  *@return a 8-bit mantissa/exponent product
59  */
60
61 olsr_u8_t
62 double_to_me(double interval)
63 {
64   int a, b;
65
66   b = 0;
67
68   while(interval / VTIME_SCALE_FACTOR >= pow((double)2, (double)b))
69     b++;
70
71   b--;
72   if(b < 0)
73     {
74       a = 1;
75       b = 0;
76     } 
77   else 
78     if (b > 15)
79       {
80         a = 15;
81         b = 15;
82       } 
83     else 
84       { 
85         a = (int)(16*((double)interval/(VTIME_SCALE_FACTOR*(double)pow(2,b))-1));
86         while(a >= 16)
87           {
88             a -= 16;
89             b++;
90           }
91       }
92
93   //printf("Generated mantissa/exponent(%d/%d): %d from %f\n", a, b, (olsr_u8_t) (a*16+b), interval);
94   //printf("Resolves back to: %f\n", me_to_double((olsr_u8_t) (a*16+b)));
95   return (olsr_u8_t) (a*16+b);
96 }
97
98
99
100
101 /**
102  *Function that converts a mantissa/exponent 8bit value back
103  *to double as described in RFC3626:
104  *
105  * value = C*(1+a/16)*2^b [in seconds]
106  *
107  *  where a is the integer represented by the four highest bits of the
108  *  field and b the integer represented by the four lowest bits of the
109  *  field.
110  *
111  *@param me the 8 bit mantissa/exponen value
112  *
113  *@return a double value
114  */
115 double
116 me_to_double(olsr_u8_t me)
117 {
118   int a, b;
119
120   a = me>>4;
121   b = me - a*16;
122
123   return (double)(VTIME_SCALE_FACTOR*(1+(double)a/16)*(double)pow(2,b));
124 }