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
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
|
(**************************************************************************)
(* *)
(* Ocamlgraph: a generic graph library for OCaml *)
(* Copyright (C) 2004-2010 *)
(* Sylvain Conchon, Jean-Christophe Filliatre and Julien Signoles *)
(* *)
(* This software is free software; you can redistribute it and/or *)
(* modify it under the terms of the GNU Library General Public *)
(* License version 2.1, with the special exception on linking *)
(* described in file LICENSE. *)
(* *)
(* This software is distributed in the hope that it will be useful, *)
(* but WITHOUT ANY WARRANTY; without even the implied warranty of *)
(* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. *)
(* *)
(**************************************************************************)
open Outils_math
let debug_outil_tort = ref false
let (w,h)= (600.,600.)
(*** Tortue Hyperbolique ***)
type coord = float * float
type turtle =
{
pos : coord ; (* with |pos| < 1 *)
dir : coord (* with |dir| = 1 *)
}
let make_turtle pos angle =
{
pos = pos ;
dir = expi angle
}
let make_turtle_dir pos dir =
{
pos = pos ;
dir = dir
}
let advance turt step =
{ pos = gamma turt.pos turt.dir step ;
dir = delta turt.pos turt.dir step }
let turn turtle u =
{ turtle with dir = turtle.dir *& u }
let turn_left turtle angle =
turn turtle (expi angle) (*** a comprendre pourquoi je dois inverser + et - de l'angle ***)
let turn_right turtle angle =
turn turtle (expi (-.angle)) (*** a comprendre pourquoi je dois inverser + et - de l'angle ***)
let to_tortue(x,y)=
((float x*.(2./.w) -. 1.),(1. -. float y *.(2./.h)))
(* ((float x*.(2./.w) ),(float y *.(2./.h) ))*)
let from_tortue (x,y) =
let xzoom = (w/.2.)
and yzoom = (h/.2.) in
(truncate (x*.xzoom +. xzoom), truncate(yzoom -. y*.yzoom))
let depart = to_tortue (truncate(w/.2.), truncate(h/.2.))
let origine =ref depart
(* GTK *)
let point_courant = ref (0,0)
(*let canvas = graphEdGTK.root *)
let moveto_gtk x y = point_courant := (x,y)
let tmoveto_gtk tor =
let (x,y)= from_tortue tor.pos in
point_courant := (x,y)
let tlineto_gtk tor line =
let (x',y')= from_tortue tor.pos in
point_courant := (x',y');
List.append line [(float x'); (float y') ]
let tdraw_string_gtk tor (ellipse : GnoCanvas.ellipse) =
let (x,y) = from_tortue tor.pos in
(* debug *)
if !debug_outil_tort then Format.eprintf "tdraw_string_gtk x=%d y=%d@." x y;
(* /debug *)
moveto_gtk x y;
ellipse#parent#move ~x:(float x) ~y:(float y);
ellipse#parent#set [`X (float x); `Y (float y)]
(* avance la tortue en tra�ant, d'une distance d, en un certain nombre d'etapes,
et retourne la nouvelle position de la tortue *)
let tdraw_edge_gtk tor d etapes line =
let d = d /. (float etapes) in
let rec list_points t liste = function
| 0 -> (t,liste)
| n ->let t = advance t d in
list_points t (tlineto_gtk t liste) (n-1)
in
let l = let (x,y) =from_tortue tor.pos in [(float x); (float y)] in
let t,lpoints = list_points tor l etapes in
(* debug *)
if (!debug_outil_tort)
then
(let ltext=
let rec chaine = function
|[]->""
|e::l->(string_of_float e)^" "^chaine l
in chaine lpoints in
Format.eprintf "taille %d %s @." (List.length lpoints) ltext);
(* /debug *)
let p = Array.of_list lpoints in
line#set [`POINTS p];
t
|
