logarithme & uint64_t

Version imprimable

25/11/2012, 22h54
Array

logarithme & uint64_t
Bonjour,

Je souhaite avoir une fonction log() qui puisse calculer le logarithme d'un uin64_t.
Est-ce possible? Math.h n'a pas l'air à supporter cela donc...
Quelqu'un a une solution?

Une fonction déclarée du genre
Code:

1 2 double log(uint64_t n) { ... };
Merci!
25/11/2012, 23h44
kwariz

Bonsoir,

un cast ? double x=log( (double) n );
27/11/2012, 09h50
mith06

Salut,

Attention en utilisant les double tu ne peux plus représenter les nombre entiers au delà de 2^53.

Tu ne pourras donc pas calculer le log de la plage [(2^53)+1 : (2^64)-1]

Citation:

Envoyé par mith06

Salut,

Attention en utilisant les double tu ne peux plus représenter les nombre entiers au delà de 2^53.

Tu ne pourras donc pas calculer le log de la plage [(2^53)+1 : (2^64)-1]

Bonjour,
je me permets un petite rectification : «en utilisant les double tu ne peux plus représenter exactement les nombre entiers au delà de 2^53.»
Cela ne signifie pas qu'on ne peut plus calculer le log avec une précision suffisante (celles des double) en utilisant un simple cast.
On peut soit calculer l'erreur commise, soit se convaincre qu'un cast est largement suffisant avec un petit programme de test comme :

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
#include <stdio.h>
#include <math.h>
#include <stdint.h>
#include <inttypes.h>
 
int main()                                                                          
{                                                                                   
  uint64_t max_u64=UINT64_MAX;                                                      
  double max_d=18446744073709551615.0;                                              
 
  printf("avec cast  : log2(%"PRIu64")=%g\n", max_u64, log( (double) max_u64 ) / lo\
g(2) );                                                                             
  printf("sans cast  : log2(%g)=%g\n", max_d, log(max_d) / log(2) );                
  printf("différence : %g\n", (log(max_d)-log((double)max_u64)));                   
  return 0;                                                                         
}

qui donne :

Code:

1
2
3
4
> ./log_test 
avec cast  : log2(18446744073709551615)=64
sans cast  : log2(1.84467e+19)=64
différence : 0

27/11/2012, 14h56
diogene

Le cast ne peut rien apporter si le code compile sans lui : Dans log(x), l'argument de log() doit être un double. Si le type de x admet un cast implicite en double, il sera effectué et ce sera équivalent à log((double)x). Sinon, on aura une erreur de compilation.

Pas étonnant donc que les résultats du test soient identiques