Discussion :

[AlexiZ]

Bonjour,

Voilà j'ai un petit souci avec ce type de données (plus d'info ici : http://caml.inria.fr/pub/docs/manual...ibref/Num.html), j'aimerais simplement faire une division de 2 entiers et en sortir un nombre décimal très grand.
Mon exemple : 103993 / 33102, qui doit donner une approximation de pi.
En caml, ça me donne ça :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
 
#load "nums.cma";;
open Num;;
 
(*Version 1 :*)
Int 103993 // Int 33102;;
(*Résultat*)
- : Num.num = Ratio <abstr>
 
(*Version 2 :*)
let printnum n = print_string (string_of_num n);;
#install_printer printnum;;
let n = Int 103993 // Int 33102;;
printnum n;;
(*Résultat*)
103993/33102- : unit = ()

Si je remplace mon "//" par une addition (+/), une multiplication (*/) ou une soustraction (-/), j'ai un résultat qui s'affiche correctement et tout va bien dans les 2 cas.
Mais pour la division, rien à faire.
Quelqu'un a une solution pour ma simple division ?

[SpiceGuid]

j'aimerais simplement faire une division de 2 entiers et en sortir un nombre décimal très grand.

Le résultat de Num.(//) est un nombre rationnel c'est-à-dire ∈ ℚ.

Si tu veux calculer dans ℕ il faut utiliser le module Big_int.

PI en goutte-à-goutte (exact, temps ∞) :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

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
#load "nums.cma";;
let add  = Big_int.add_big_int
and sub  = Big_int.sub_big_int
and succ = Big_int.succ_big_int
and pred = Big_int.pred_big_int
and mult = Big_int.mult_big_int
and div  = Big_int.div_big_int
and add_int  = Big_int.add_int_big_int
and mult_int = Big_int.mult_int_big_int
and big_int  = Big_int.big_int_of_int
and int_of   = Big_int.int_of_big_int
;;
let pi () =
  let rec g q r t i =
    let i3 = mult_int 3 i in
    let u = mult_int 3 (mult (add_int 1 i3) (add_int 2 i3))
    and y = int_of (div (add (mult q (add_int (-12) (mult_int 27 i))) (mult_int 5 r)) (mult_int 5 t)) 
    in begin
      print_int y;
      flush stdout;
      g
      (mult_int 10 (mult q (mult i (add_int (-1) (mult_int 2 i)))))
      (mult_int 10 (mult u (sub (add (mult q (add_int (-2) (mult_int 5 i))) r) (mult_int y t))))
      (mult t u)
      (add_int 1 i);
      ()
    end
  in g (big_int 1) (big_int 180) (big_int 60) (big_int 2);; 

pi ();;

[AlexiZ]

Merci d'avoir répondu.

L'objectif n'est pas de cacluler dans ℕ mais bien dans ℚ. Est-ce que c'est possible ?
La méthode du goutte à goutte donne bien un nombre important de décimales mais je ne comprends pas vraiment comment cela fonctionne, une précision sur les paramètres de g ?

[SpiceGuid]

L'objectif n'est pas de cacluler dans ℕ mais bien dans ℚ. Est-ce que c'est possible ?

Oui, c'est possible en utilisant le module Num.

une précision sur les paramètres de g ?

http://web.comlab.ox.ac.uk/oucl/work...ons/spigot.pdf

[AlexiZ]

Oui, c'est possible en utilisant Num.(//)

ça va pas être simple
Mais en tout cas merci, je vais voir avec tout ça pour trouver quelque chose de probant.
Juste une dernière question en passant, je ne peux pas utiliser "power_num" sans avoir une "Exception : invalid argument power_num", une idée ?

[SpiceGuid]

L'objectif n'est pas de cacluler dans ℕ mais bien dans ℚ. Est-ce que c'est possible ?

Ne reviens pas demain pour nous dire qu'en fait tu veux calculer ∏ dans ℝ