X-Git-Url: https://scm.cri.ensmp.fr/git/Faustine.git/blobdiff_plain/440ca0ba93966e89b68dc54207c461afc0d56264..73e86244c539dd6ff8af5d8527b11fde18c8d070:/interpretor/faustexp.ml diff --git a/interpretor/faustexp.ml b/interpretor/faustexp.ml index 6b324a0..d0b5ec0 100644 --- a/interpretor/faustexp.ml +++ b/interpretor/faustexp.ml @@ -1,22 +1,87 @@ (** 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;; - - - -class virtual expression : expression_type = +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 - val exp = exp_init + 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 @@ -24,56 +89,270 @@ class virtual expression : expression_type = method get_delay = delay method to_string = "NotYetDone" method virtual evaluate : beam_type -> beam_type -end;; - - -class exp_const = - object - inherit expression - method evaluate = fun b1 -> - - end;; - - -class exp_ident = - object - inherit expression - - end;; - - -class exp_par = - object - inherit expression + 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 exp_split = - object - inherit expression - +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 exp_merge = - object - inherit expression - - 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 + -class exp_seq = - object - inherit expression + end - end;; - -class exp_rec = +and proc_factory = object - inherit expression - - end;; - + 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;; @@ -97,15 +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 - |Rdtable -> 100000 (* danger! *) + |Rdtable -> 100000 |Mod -> 0 |Larger -> 0 |Smaller -> 0 - |Vectorize -> 100 (* danger! *) + |Vectorize -> 100 |Concat -> 0 |Nth -> 0 |Serialize -> 0 @@ -152,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 =