X-Git-Url: https://scm.cri.mines-paristech.fr/git/Faustine.git/blobdiff_plain/06159b51a934937f647ec7119b47cb466d8e50b1..73e86244c539dd6ff8af5d8527b11fde18c8d070:/interpretor/faustexp.ml?ds=sidebyside

diff --git a/interpretor/faustexp.ml b/interpretor/faustexp.ml
index 5904569..d0b5ec0 100644
--- a/interpretor/faustexp.ml
+++ b/interpretor/faustexp.ml
@@ -1,20 +1,361 @@
 (**
  	Module: Faustexp	
-	Description: dimension estimation and delay estimation of faust expressions.
+	Description: Faust expression evaluation
 	@author WANG Haisheng	
-	Created: 03/06/2013	Modified: 04/06/2013
+	Created: 03/06/2013	Modified: 04/08/2013
 *)
 
 open Types;;
+open Basic;;
+open Symbol;;
 open Value;;
+open Signal;;
+open Beam;;
+
+exception NotYetDone;;
+exception Dimension_error of string;;
+exception Process_error of string;;
+
+class dimension : int * int -> dimension_type = 
+  fun (init : int * int) -> 
+    object (self)
+      val _input = fst init
+      val _output = snd init
+
+      method input = _input
+      method output = _output
+
+      method par : dimension_type -> dimension_type = 
+	fun dim -> 
+	  new dimension 
+	    ((self#input + dim#input), (self#output + dim#output))
+	    
+      method seq : dimension_type -> dimension_type = 
+	fun dim -> 
+	  if self#output = dim#input then
+	    new dimension (self#input, dim#output)
+	  else raise (Dimension_error "seq dimension not matched.")
+	      
+      method split : dimension_type -> dimension_type =
+	fun dim ->
+	  if dim#input mod self#output = 0 then
+	    new dimension (self#input, dim#output)
+	  else raise (Dimension_error "split dimension not matched.")
+
+      method merge : dimension_type -> dimension_type =
+	fun dim ->
+	  if self#output mod dim#input = 0 then
+	    new dimension (self#input, dim#output)
+	  else raise (Dimension_error "merge dimension not matched.")
+
+      method _rec : dimension_type -> dimension_type = 
+	fun dim ->
+	  if self#output >= dim#input && self#input >= dim#output then
+	    new dimension (self#input - dim#output, self#output)
+	  else raise (Dimension_error "rec dimension not matched.")
+    end;;
+
+(*
+class process : faust_exp -> process_type = 
+  fun (exp_init : faust_exp) ->
+  object (self)
+    val exp = exp_init
+    val left = 
+      match exp_init with
+      |	Const b -> exp_init
+      |	Ident s -> exp_init
+      |	Par (e1, e2) -> e1
+      |	Seq (e1, e2) -> e1
+      |	Split (e1, e2) -> e1
+      |	Merge (e1, e2) -> e1
+      |	Rec (e1, e2) -> e1
+
+    val right = 
+      match exp_init with
+      |	Const b -> exp_init
+      |	Ident s -> exp_init
+      |	Par (e1, e2) -> e2
+      |	Seq (e1, e2) -> e2
+      |	Split (e1, e2) -> e2
+      |	Merge (e1, e2) -> e2
+      |	Rec (e1, e2) -> e2
+
+    val proc_left = 
+
+      val dim = new dimension
+      val delay = 0
+      method get_exp = exp
+      method get_dim = dim
+      method get_delay = delay
+      method to_string = "NotYetDone"
+      method virtual evaluate : beam_type -> beam_type
+  end;;
+*)
+
+class proc_const : faust_exp -> process_type = 
+  fun (exp_init : faust_exp) ->
+    object (self)
+      val _exp = exp_init
+      val _dim = new dimension (0,1)
+      val _delay = 0
+      val _const = 
+	match exp_init with
+	| Const b -> b
+	| _ -> raise (Process_error "const process constructor.")
+
+      method exp = _exp
+      method dim = _dim
+      method delay = _delay
+      method private const = _const
+
+      method eval : beam_type -> beam_type = 
+	fun (input : beam_type) ->
+	  if input#get = [||] then
+	    new beam [| new signal 0 (fun t -> new value self#const)|]
+	  else
+	    raise (Process_error "proc_const accepts no input.")
+      end;;
+
+
+class proc_ident : faust_exp -> process_type = 
+  fun (exp_init : faust_exp) ->
+    object (self)
+      val _exp = exp_init
+      val _symbol = 
+	match exp_init with
+	| Ident s -> s
+	| _ -> raise (Process_error "ident process constructor.")
+
+      val _dim = new dimension (dimension_of_symbol self#symb)
+      val _delay = delay_of_symbol _symbol      
+
+      method exp = _exp
+      method dim = _dim
+      method delay = _delay
+      method symb = _symbol
+
+      method private beam_of_ident : int -> signal_type -> beam_type = 
+	fun (n : int) ->
+	  fun (s, signal_type) ->
+	    if n = (self#dim)#input then 
+	      new beam [|s|]
+	    else raise (Process_error ("Ident " ^ string_of_symbol self#symb))
+
+      method eval : beam_type -> beam_type = 
+	fun (input : beam_type) ->
+	  let n = Array.length input#get in
+	  match self#symb with
+	  | Pass -> self#beam_of_ident n input#get.(0)
+	  | Stop -> if n = 1 then new beam [||] 
+	            else raise (Process_error "Ident !")
+	  | Add -> self#beam_of_ident n ((input#get.(0))#add input#get.(1))
+	  | Sub -> self#beam_of_ident n ((input#get.(0))#sub input#get.(1))
+	  | Mul -> self#beam_of_ident n ((input#get.(0))#mul input#get.(1))
+	  | Div -> self#beam_of_ident n ((input#get.(0))#div input#get.(1))
+	  | Mem -> self#beam_of_ident n ((input#get.(0))#mem)
+	  | Delay -> self#beam_of_ident n ((input#get.(0))#delay input#get.(1))
+	  | Floor -> self#beam_of_ident n ((input#get.(0))#floor)
+	  | Int -> self#beam_of_ident n ((input#get.(0))#int)
+	  | Sin -> self#beam_of_ident n ((input#get.(0))#sin)
+	  | Cos -> self#beam_of_ident n ((input#get.(0))#cos)
+	  | Atan -> self#beam_of_ident n ((input#get.(0))#atan)
+	  | Atan2 -> self#beam_of_ident n ((input#get.(0))#atan2 input#get.(1))
+	  | Sqrt -> self#beam_of_ident n ((input#get.(0))#sqrt)
+	  | Rdtable -> self#beam_of_ident n 
+		((input#get.(1))#rdtable input#get.(0) input#get.(2))
+	  | Mod -> self#beam_of_ident n 
+		((input#get.(0))#_mod input#get.(1))
+	  | Vectorize -> self#beam_of_ident n 
+		((input#get.(0))#vectorzie input#get.(1))
+	  | Vconcat -> self#beam_of_ident n 
+		((input#get.(0))#vconcat input#get.(1))
+	  | Vpick -> self#beam_of_ident n 
+		((input#get.(0))#vpick input#get.(1))
+	  | Serialize -> self#beam_of_ident n 
+		(input#get.(0))#serialize
+	  | Larger -> self#beam_of_ident n 
+		((input#get.(0))#larger input#get.(1))
+	  | Smaller -> self#beam_of_ident n 
+		((input#get.(0))#smaller input#get.(1))
+	  | Prefix -> self#beam_of_ident n 
+		((input#get.(1))#prefix input#get.(0))
+	  | Selec2 -> self#beam_of_ident n 
+		((input#get.(0))#select2 input#get.(1) input#get.(2))
+	  | Select3 -> self#beam_of_ident n 
+		((input#get.(0))#select3 input#get.(1) 
+		   input#get.(2) input#get.(3))
+      end;;
+
+
+class proc_par : faust_exp -> process_type = 
+  fun (exp_init : faust_exp) ->
+    object (self)
+      val _exp = exp_init
+      val _exp_left = 
+	match exp_init with
+	| Par (e1, e2) -> e1
+	| _ -> raise (Process_error "par process constructor.")
+      val _exp_right = 
+	match exp_init with
+	| Par (e1, e2) -> e2
+	| _ -> raise (Process_error "par process constructor.")
+
+      val proc_left = (new proc_factory)#make _exp_left
+      val proc_right = (new proc_factory)#make _exp_right
+
+      val _dim = (proc_left#dim)#par proc_right#dim
+      val _delay = max proc_left#delay proc_right#delay
+
+      method exp = _exp
+      method dim = _dim
+      method delay = _delay
+
+      method eval : beam_type -> beam_type = 
+	fun (input : beam_type) ->
+	  let (sub_input1, sub_input2) = input#cut proc_left#dim#input in
+	  let sub_output1 = proc_left#eval sub_input1 in
+	  let sub_output2 = proc_right#eval sub_input2 in
+	  sub_output1#append sub_output2
+      end
+
+and proc_split : faust_exp -> process_type =
+  fun (exp_init : faust_exp) ->
+    object (self)
+      val _exp = exp_init
+      val _exp_left = 
+	match exp_init with
+	| Split (e1, e2) -> e1
+	| _ -> raise (Process_error "par process constructor.")
+      val _exp_right = 
+	match exp_init with
+	| Split (e1, e2) -> e2
+	| _ -> raise (Process_error "par process constructor.")
+
+      val proc_left = (new proc_factory)#make _exp_left
+      val proc_right = (new proc_factory)#make _exp_right
+
+      val _dim = (proc_left#dim)#split proc_right#dim
+      val _delay = proc_left#delay + proc_right#delay
+
+      method exp = _exp
+      method dim = _dim
+      method delay = _delay
+
+      method eval : beam_type -> beam_type = 
+	fun (input : beam_type) ->
+	  let mid_output = proc_left#eval input in
+	  let mid_input = mid_output#matching proc_right#dim#input in
+	  proc_right#eval mid_input
+      end
+
+
+and proc_merge : faust_exp -> process_type =
+  fun (exp_init : faust_exp) -> 
+    object (self)
+      val _exp = exp_init
+      val _exp_left = 
+	match exp_init with
+	| Merge (e1, e2) -> e1
+	| _ -> raise (Process_error "merge process constructor.")
+      val _exp_right = 
+	match exp_init with
+	| Merge (e1, e2) -> e2
+	| _ -> raise (Process_error "merge process constructor.")
+
+      val proc_left = (new proc_factory)#make _exp_left
+      val proc_right = (new proc_factory)#make _exp_right
+
+      val _dim = (proc_left#dim)#merge proc_right#dim
+      val _delay = proc_left#delay + proc_right#delay
+
+      method exp = _exp
+      method dim = _dim
+      method delay = _delay
+
+      method eval : beam_type -> beam_type = 
+	fun (input : beam_type) ->
+	  let mid_output = proc_left#eval input in
+	  let mid_input = mid_output#matching proc_right#dim#input in
+	  proc_right#eval mid_input
+
+      end
+
+and proc_seq : faust_exp -> process_type =
+  fun (exp_init : faust_exp) -> 
+    object (self)
+      val _exp = exp_init
+      val _exp_left = 
+	match exp_init with
+	| Seq (e1, e2) -> e1
+	| _ -> raise (Process_error "seq process constructor.")
+      val _exp_right = 
+	match exp_init with
+	| Seq (e1, e2) -> e2
+	| _ -> raise (Process_error "seq process constructor.")
+
+      val proc_left = (new proc_factory)#make _exp_left
+      val proc_right = (new proc_factory)#make _exp_right
+
+      val _dim = (proc_left#dim)#seq proc_right#dim
+      val _delay = proc_left#delay + proc_right#delay
+
+      method exp = _exp
+      method dim = _dim
+      method delay = _delay
+
+      method eval : beam_type -> beam_type = 
+	fun (input : beam_type) ->
+	  let mid_output = proc_left#eval input in
+	  proc_right#eval mid_output
+      end
+
+and proc_rec : faust_exp -> process_type =
+  fun (exp_init : faust_exp) -> 
+  object (self)
+    val _exp = exp_init
+    val _exp_left = 
+      match exp_init with
+      | Rec (e1, e2) -> e1
+      | _ -> raise (Process_error "rec process constructor.")
+    val _exp_right = 
+      match exp_init with
+      | Rec (e1, e2) -> e2
+      | _ -> raise (Process_error "rec process constructor.")
+	    
+    val proc_left = (new proc_factory)#make _exp_left
+    val proc_right = (new proc_factory)#make _exp_right
+
+    val _dim = (proc_left#dim)#_rec proc_right#dim
+    val _delay = proc_left#delay
+	
+    method exp = _exp
+    method dim = _dim
+    method delay = _delay
+	
+    method eval : beam_type -> beam_type = 
+      fun (input : beam_type) ->
+	let mid_output = proc_left#eval input in
+	proc_right#eval mid_output
+	  
+
+  end
+
+and proc_factory = 
+  object
+      method make : faust_exp -> process_type = 
+	fun (exp : faust_exp) ->
+	  match exp with
+	  | Const b -> new proc_const exp
+	  | Ident s -> new proc_ident exp
+	  | Par (e1, e2) -> new proc_par exp
+	  | Seq (e1, e2) -> new proc_seq exp
+	  | Split (e1, e2) -> new proc_split exp
+	  | Merge (e1, e2) -> new proc_merge exp
+	  | Rec (e1, e2) -> new proc_rec exp
+  end;;
 
-(* EXCEPTIONS *)
 
-(** Exception raised in beam matching of faust expressions.*)
-exception Beam_Matching_Error of string;;
 
-(** Exception raised in case that the branch under call hasn't yet been programed.*)
-exception NotYetDone;;
 
 
 (* PROCESS DELAY ESTIMATION *)
@@ -35,19 +376,15 @@ let rec delay exp_faust = match exp_faust with
 			|Pass			->	0
 			|Stop			->	0
 			|Mem			->	1
-			|Delay			->	100000 (* danger! *)
+			|Delay			->	100000 
 			|Floor		        ->	0
 			|Int			->	0
 			|Sin			->	0
-			|Cos			->	0
-			|Atan			->	0
-			|Atantwo                ->      0
-			|Sqrt			->	0
-			|Rdtable		->	100000 (* danger! *)
+			|Rdtable		->	100000
 			|Mod			->	0
 			|Larger			->	0
 			|Smaller		->	0
-			|Vectorize		->	100 (* danger! *)
+			|Vectorize		->	100
 			|Concat			->	0
 			|Nth			->	0
 			|Serialize		->	0
@@ -94,46 +431,6 @@ let subtree_left = fun d_tree -> subtree d_tree 0;;
 (** val subtree_right : dimension -> dimension, returns the right subtree of dimension tree.*)
 let subtree_right = fun d_tree -> subtree d_tree 1;;
 
-
-(** val d_par : int * int -> int * int -> int * int, process dimension for constructor "par(,)", 
-which is the addition of two dimensions.*)
-let d_par a b = (((fst a) + (fst b)), ((snd a) + (snd b)));;
-
-
-(** val d_seq : int * int -> int * int -> int * int, process dimension for constructor "seq(:)", 
-which is (size of input beam of first exp, size of output beam of second exp) 
-along with beam matching.*)
-let d_seq a b = if (snd a) = (fst b) then (fst a, snd b) else raise (Beam_Matching_Error "seq");;
-
-
-(** val d_split : int * int -> int * int -> int * int, process dimension for constructor "split(<:)", 
-which is (size of input beam of first exp, size of output beam of second exp) 
-along with beam matching.*)
-let d_split a b = 
-  if ((fst b) mod (snd a)) = 0 then 
-    (fst a, snd b) 
-  else raise (Beam_Matching_Error "split");;
-
-
-(** val d_merge : int * int -> int * int -> int * int, process dimension for constructor "merge(:>)", 
-which is (size of input beam of first exp, size of output beam of second exp) 
-along with beam matching. *)
-let d_merge a b = 
-  if ((snd a) mod (fst b)) = 0 then 
-    (fst a, snd b) 
-  else raise (Beam_Matching_Error "merge");;
-
-
-(** val d_rec : int * int -> int * int -> int * int, process dimension for constructor "rec(~)", 
-which is (size of input beam of first exp - size of output beam of second exp, 
-size of output beam of first exp) 
-along with beam matching.*)
-let d_rec a b = 
-  if (fst a) >= (snd b) && (snd a) >= (fst b) then 
-    ((fst a) - (snd b), snd a) 
-  else raise (Beam_Matching_Error "rec");;
-
-
 (** val dim : faust_exp -> int * int, returns dimension for faust expression, 
 along with beam matching.*)
 let rec dim exp_faust = 
@@ -163,10 +460,6 @@ returns the dimension tree of constructor(e1, e2).*)
 			|Floor	    	        ->      End (1, 1)
 			|Int			->	End (1, 1)
 			|Sin			->	End (1, 1)
-			|Cos			->	End (1, 1)
-			|Atan			->	End (1, 1)
-			|Atantwo                ->      End (2, 1)
-			|Sqrt			->	End (1, 1)
 			|Rdtable		->	End (3, 1)
 			|Mod  		        ->	End (2, 1)
 			|Vectorize		->	End (2, 1)