CristalCaveGem » LustreC

(********************************************************************)

(*                                                                  *)

(*  The LustreC compiler toolset   /  The LustreC Development Team  *)

(*  Copyright 2012 -    --   ONERA - CNRS - INPT                    *)

(*                                                                  *)

(*  LustreC is free software, distributed WITHOUT ANY WARRANTY      *)

(*  under the terms of the GNU Lesser General Public License        *)

(*  version 2.1.                                                    *)

(*                                                                  *)

(********************************************************************)

open Utils.Format

open Lustre_types

open Corelang

open Machine_code_types

(*open Machine_code_common*)

module Mpfr = Lustrec_mpfr

let pp_print_version fmt () =

  fprintf fmt

    "/* @[<v>\

     C code generated by %s@,\

     Version number %s@,\

     Code is %s compliant@,\

     Using %s numbers */@,\

     @]"

    (Filename.basename Sys.executable_name) 

    Version.number 

    (if !Options.ansi then "ANSI C90" else "C99")

    (if !Options.mpfr then "MPFR multi-precision" else "(double) floating-point")

let protect_filename s =

  Str.global_replace (Str.regexp "\\.\\|\\ ") "_" s

let file_to_module_name basename =

  let baseNAME = Ocaml_utils.uppercase basename in

  let baseNAME = protect_filename baseNAME in

  baseNAME

let var_is name v =

  v.var_id = name

let mk_local n m =

  let used name =

    let open List in

    exists (var_is name) m.mstep.step_inputs

    || exists (var_is name) m.mstep.step_outputs

    || exists (var_is name) m.mstep.step_locals

    || exists (var_is name) m.mmemory in

  mk_new_name used n

(* Generation of a non-clashing name for the self memory variable (for step and reset functions) *)

let mk_self = mk_local "self"

let mk_mem = mk_local "mem"

let mk_mem_in = mk_local "mem_in"

let mk_mem_out = mk_local "mem_out"

(* Generation of a non-clashing name for the instance variable of static allocation macro *)

let mk_instance m =

  let used name =

    let open List in

    exists (var_is name) m.mstep.step_inputs

    || exists (var_is name) m.mmemory in

  mk_new_name used "inst"

(* Generation of a non-clashing name for the attribute variable of static allocation macro *)

let mk_attribute m =

  let used name =

    let open List in

    exists (var_is name) m.mstep.step_inputs

    || exists (var_is name) m.mmemory in

  mk_new_name used "attr"

let mk_call_var_decl loc id =

  { var_id = id;

    var_orig = false;

    var_dec_type = mktyp Location.dummy_loc Tydec_any;

    var_dec_clock = mkclock Location.dummy_loc Ckdec_any;

    var_dec_const = false;

    var_dec_value = None;

    var_parent_nodeid = None;

    var_type = Type_predef.type_arrow (Types.new_var ()) (Types.new_var ());

    var_clock = Clocks.new_var true;

    var_loc = loc }

(* counter for loop variable creation *)

let loop_cpt = ref (-1)

let reset_loop_counter () =

 loop_cpt := -1

let mk_loop_var m () =

  let vars = m.mstep.step_inputs

             @ m.mstep.step_outputs

             @ m.mstep.step_locals

             @ m.mmemory in

  let rec aux () =

    incr loop_cpt;

    let s = sprintf "__%s_%d" "i" !loop_cpt in

    if List.exists (var_is s) vars then aux () else s

  in aux ()

(*

let addr_cpt = ref (-1)

let reset_addr_counter () =

 addr_cpt := -1

let mk_addr_var m var =

  let vars = m.mmemory in

  let rec aux () =

    incr addr_cpt;

    let s = Printf.sprintf "%s_%s_%d" var "addr" !addr_cpt in

    if List.exists (fun v -> v.var_id = s) vars then aux () else s

  in aux ()

*)

let pp_global_init_name fmt id = fprintf fmt "%s_INIT" id

let pp_global_clear_name fmt id = fprintf fmt "%s_CLEAR" id

let pp_machine_memtype_name ?(ghost=false) fmt id =

  fprintf fmt "struct %s_mem%s" id (if ghost then "_ghost" else "")

let pp_machine_regtype_name fmt id = fprintf fmt "struct %s_reg" id

let pp_machine_alloc_name fmt id = fprintf fmt "%s_alloc" id

let pp_machine_dealloc_name fmt id = fprintf fmt "%s_dealloc" id

let pp_machine_static_declare_name fmt id = fprintf fmt "%s_DECLARE" id

let pp_machine_static_link_name fmt id = fprintf fmt "%s_LINK" id

let pp_machine_static_alloc_name fmt id = fprintf fmt "%s_ALLOC" id

let pp_machine_reset_name fmt id = fprintf fmt "%s_reset" id

let pp_machine_init_name fmt id = fprintf fmt "%s_init" id

let pp_machine_clear_name fmt id = fprintf fmt "%s_clear" id

let pp_machine_step_name fmt id = fprintf fmt "%s_step" id

let pp_mod pp_val v1 v2 fmt =

  if !Options.integer_div_euclidean then

    (* (a mod_C b) + (a mod_C b < 0 ? abs(b) : 0) *)

    fprintf fmt "((%a %% %a) + ((%a %% %a) < 0?(abs(%a)):0))"

      pp_val v1 pp_val v2

      pp_val v1 pp_val v2

      pp_val v2

  else (* Regular behavior: printing a % *)

    fprintf fmt "(%a %% %a)" pp_val v1 pp_val v2

let pp_div pp_val v1 v2 fmt =

  if !Options.integer_div_euclidean then

    (* (a - ((a mod_C b) + (a mod_C b < 0 ? abs(b) : 0))) div_C b *)

    fprintf fmt "(%a - %t) / %a"

      pp_val v1

      (pp_mod pp_val v1 v2)

      pp_val v2

  else (* Regular behavior: printing a / *)

    fprintf fmt "(%a / %a)" pp_val v1 pp_val v2

let pp_basic_lib_fun is_int i pp_val fmt vl =

  match i, vl with

  (*  | "ite", [v1; v2; v3] -> fprintf fmt "(%a?(%a):(%a))" pp_val v1 pp_val v2 pp_val v3 *)

  | "uminus", [v] ->

    fprintf fmt "(- %a)" pp_val v

  | "not", [v] ->

    fprintf fmt "(!%a)" pp_val v

  | "impl", [v1; v2] ->

    fprintf fmt "(!%a || %a)" pp_val v1 pp_val v2

  | "=", [v1; v2] ->

    fprintf fmt "(%a == %a)" pp_val v1 pp_val v2

  | "mod", [v1; v2] ->

     if is_int then

       pp_mod pp_val v1 v2 fmt 

     else

       fprintf fmt "(%a %% %a)" pp_val v1 pp_val v2

  | "equi", [v1; v2] ->

    fprintf fmt "(!%a == !%a)" pp_val v1 pp_val v2

  | "xor", [v1; v2] ->

    fprintf fmt "(!%a != !%a)" pp_val v1 pp_val v2

  | "/", [v1; v2] ->

     if is_int then

       pp_div pp_val v1 v2 fmt

     else

       fprintf fmt "(%a / %a)" pp_val v1 pp_val v2

  | _, [v1; v2] ->

    fprintf fmt "(%a %s %a)" pp_val v1 i pp_val v2

  | _ ->

    (* TODO: raise proper error *)

    eprintf "internal error: Basic_library.pp_c %s@." i;

    assert false

let rec pp_c_dimension fmt dim =

  let open Dimension in

  match dim.dim_desc with

  | Dident id ->

    fprintf fmt "%s" id

  | Dint i ->

    fprintf fmt "%d" i

  | Dbool b ->

    fprintf fmt "%B" b

  | Dite (i, t, e) ->

    fprintf fmt "((%a)?%a:%a)"

      pp_c_dimension i pp_c_dimension t pp_c_dimension e

  | Dappl (f, args) ->

    fprintf fmt "%a"

      (pp_basic_lib_fun (Basic_library.is_numeric_operator f) f pp_c_dimension)

      args

  | Dlink dim' ->

    fprintf fmt "%a" pp_c_dimension dim'

  | Dvar ->

    fprintf fmt "_%s" (Utils.name_of_dimension dim.dim_id)

  | Dunivar ->

    fprintf fmt "'%s" (Utils.name_of_dimension dim.dim_id)

let is_basic_c_type t =

  Types.(is_int_type t || is_real_type t || is_bool_type t)

let pp_c_basic_type_desc t_desc =

  if Types.is_bool_type t_desc then

    if !Options.cpp then "bool" else "_Bool"

  else if Types.is_int_type t_desc then !Options.int_type

  else if Types.is_real_type t_desc then

    if !Options.mpfr then Mpfr.mpfr_t else !Options.real_type

  else

    assert false (* Not a basic C type. Do not handle arrays or pointers *)

let pp_basic_c_type ?(var_opt=None) fmt t =

  match var_opt with

  | Some v when Machine_types.is_exportable v ->

     Machine_types.pp_c_var_type fmt v

  | _ ->

     fprintf fmt "%s" (pp_c_basic_type_desc t)

let pp_c_type ?var_opt var_id fmt t =

  let rec aux t pp_suffix =

    if is_basic_c_type  t then

       fprintf fmt "%a %s%a"

         (pp_basic_c_type ~var_opt) t

         var_id

         pp_suffix ()

    else

      let open Types in

      match (repr t).tdesc with

      | Tclock t' ->

        aux t' pp_suffix

      | Tarray (d, t')  ->

        let pp_suffix' fmt () =

          fprintf fmt "%a[%a]" pp_suffix () pp_c_dimension d in

        aux t' pp_suffix'

      | Tstatic (_, t') ->

        fprintf fmt "const "; aux t' pp_suffix

      | Tconst ty ->

        fprintf fmt "%s %s" ty var_id

      | Tarrow (_, _) ->

        fprintf fmt "void (*%s)()" var_id

      | _ ->

        (* TODO: raise proper error *)

        eprintf "internal error: C_backend_common.pp_c_type %a@." print_ty t;

        assert false

  in aux t (fun _ () -> ())

(*

let rec pp_c_initialize fmt t = 

  match (Types.repr t).Types.tdesc with

  | Types.Tint -> pp_print_string fmt "0"

  | Types.Tclock t' -> pp_c_initialize fmt t'

  | Types.Tbool -> pp_print_string fmt "0" 

  | Types.Treal when not !Options.mpfr -> pp_print_string fmt "0."

  | Types.Tarray (d, t') when Dimension.is_dimension_const d ->

    fprintf fmt "{%a}"

      (Utils.fprintf_list ~sep:"," (fun fmt _ -> pp_c_initialize fmt t'))

      (Utils.duplicate 0 (Dimension.size_const_dimension d))

  | _ -> assert false

 *)

let pp_c_tag fmt t =

  pp_print_string fmt

    (if t = tag_true then "1" else if t = tag_false then "0" else t)

(* Prints a constant value *)

let rec pp_c_const fmt c =

  match c with

  | Const_int i ->

    pp_print_int fmt i

  | Const_real r ->

    Real.pp fmt r

  (* | Const_float r   -> pp_print_float fmt r *)

  | Const_tag t ->

    pp_c_tag fmt t

  | Const_array ca ->

    pp_print_braced pp_c_const fmt ca

  | Const_struct fl ->

    pp_print_braced (fun fmt (_, c) -> pp_c_const fmt c) fmt fl

  | Const_string _

  | Const_modeid _ -> assert false (* string occurs in annotations not in C *)

(* Prints a value expression [v], with internal function calls only.

   [pp_var] is a printer for variables (typically [pp_c_var_read]),

   but an offset suffix may be added for array variables

*)

let rec pp_c_val m self pp_var fmt v =

  let pp_c_val = pp_c_val m self pp_var in

  match v.value_desc with

  | Cst c ->

    pp_c_const fmt c

  | Array vl ->

    pp_print_braced pp_c_val fmt vl

  | Access (t, i) ->

    fprintf fmt "%a[%a]" pp_c_val t pp_c_val i

  | Power (v, _) ->

    (* TODO: raise proper error *)

    eprintf "internal error: C_backend_common.pp_c_val %a@."

      (Machine_code_common.pp_val m) v;

    assert false

  | Var v ->

     if Machine_code_common.is_memory m v then

       (* array memory vars are represented by an indirection to a local var

        *  with the right type, in order to avoid casting everywhere. *)

       if Types.is_array_type v.var_type

       && not (Types.is_real_type v.var_type && !Options.mpfr)

       then fprintf fmt "%a" pp_var v

       else fprintf fmt "%s->_reg.%a" self pp_var v

     else

       pp_var fmt v

  | Fun (n, vl) ->

    pp_basic_lib_fun (Types.is_int_type v.value_type) n pp_c_val fmt vl

(* Access to the value of a variable:

   - if it's not a scalar output, then its name is enough

   - otherwise, dereference it (it has been declared as a pointer,

     despite its scalar Lustre type)

   - moreover, dereference memory array variables.

*)

let pp_c_var_read m fmt id =

  (* mpfr_t is a static array, not treated as general arrays *)

  if Types.is_address_type id.var_type

  then

    if Machine_code_common.is_memory m id

    && not (Types.is_real_type id.var_type && !Options.mpfr)

    then fprintf fmt "(*%s)" id.var_id

    else fprintf fmt "%s" id.var_id

  else

    if Machine_code_common.is_output m id

    then fprintf fmt "*%s" id.var_id

    else fprintf fmt "%s" id.var_id

(* Addressable value of a variable, the one that is passed around in calls:

   - if it's not a scalar non-output, then its name is enough

   - otherwise, reference it (it must be passed as a pointer,

     despite its scalar Lustre type)

*)

let pp_c_var_write m fmt id =

  if Types.is_address_type id.var_type

  then

    fprintf fmt "%s" id.var_id

  else

    if Machine_code_common.is_output m id

    then

      fprintf fmt "%s" id.var_id

    else

      fprintf fmt "&%s" id.var_id

(* Declaration of an input variable:

   - if its type is array/matrix/etc, then declare it as a mere pointer,

     in order to cope with unknown/parametric array dimensions, 

     as it is the case for generics

*)

let pp_c_decl_input_var fmt id =

  if !Options.ansi && Types.is_address_type id.var_type

  then

    pp_c_type ~var_opt:id (sprintf "(*%s)" id.var_id) fmt

      (Types.array_base_type id.var_type)

  else

    pp_c_type ~var_opt:id id.var_id fmt id.var_type

(* Declaration of an output variable:

   - if its type is scalar, then pass its address

   - if its type is array/matrix/struct/etc, then declare it as a mere pointer,

     in order to cope with unknown/parametric array dimensions, 

     as it is the case for generics

*)

let pp_c_decl_output_var fmt id =

  if (not !Options.ansi) && Types.is_address_type id.var_type

  then

    pp_c_type ~var_opt:id id.var_id fmt id.var_type

  else

    pp_c_type ~var_opt:id (sprintf "(*%s)" id.var_id) fmt

      (Types.array_base_type id.var_type)

(* Declaration of a local/mem variable:

   - if it's an array/matrix/etc, its size(s) should be

     known in order to statically allocate memory, 

     so we print the full type

*)

let pp_c_decl_local_var m fmt id =

  if id.var_dec_const

  then

    fprintf fmt "%a = %a"

      (pp_c_type ~var_opt:id id.var_id)

      id.var_type

      (pp_c_val m "" (pp_c_var_read m))

      (Machine_code_common.get_const_assign m id)

  else

    fprintf fmt "%a"

      (pp_c_type ~var_opt:id id.var_id) id.var_type

(* Declaration of a struct variable:

   - if it's an array/matrix/etc, we declare it as a pointer

*)

let pp_c_decl_struct_var fmt id =

  if Types.is_array_type id.var_type

  then

    pp_c_type (sprintf "(*%s)" id.var_id) fmt

      (Types.array_base_type id.var_type)

  else

    pp_c_type id.var_id  fmt id.var_type

let pp_c_decl_instance_var ?(ghost=false) fmt (name, (node, _)) =

  fprintf fmt "%a %s%s"

    (pp_machine_memtype_name ~ghost) (node_name node)

    (if ghost then "" else "*")

    name

(* let pp_c_checks self fmt m =

 *   pp_print_list

 *     (fun fmt (loc, check) ->

 *        fprintf fmt

 *          "@[<v>%a@,assert (%a);@]"

 *          Location.pp_c_loc loc

 *          (pp_c_val m self (pp_c_var_read m)) check)

 *     fmt

 *     m.mstep.step_checks *)

let has_c_prototype funname dependencies =

  (* We select the last imported node with the name funname.

     The order of evaluation of dependencies should be

     compatible with overloading. (Not checked yet) *)

  let imported_node_opt =

    List.fold_left

      (fun res dep ->

         match res with

         | Some _ -> res

         | None ->

           let decls = dep.content in

           let matched = fun t -> match t.top_decl_desc with

             | ImportedNode nd -> nd.nodei_id = funname

             | _ -> false

           in

           if List.exists matched decls then

             match (List.find matched decls).top_decl_desc with

             | ImportedNode nd -> Some nd

             | _ -> assert false

           else

             None) None dependencies in

  match imported_node_opt with

  | None -> false

  | Some nd -> (match nd.nodei_prototype with Some "C" -> true | _ -> false)

(********************************************************************************************)

(*                       Struct Printing functions                                          *)

(********************************************************************************************)

(* let pp_registers_struct fmt m =

 *   pp_print_braced

 *     ~pp_prologue:(fun fmt () ->

 *         fprintf fmt "@[%a " pp_machine_regtype_name m.mname.node_id)

 *     ~pp_open_box:pp_open_vbox0

 *     ~pp_sep:pp_print_semicolon

 *     ~pp_eol:pp_print_semicolon

 *     ~pp_epilogue:(fun fmt () -> pp_print_string fmt "@] _reg;")

 *     pp_c_decl_struct_var

 *     fmt m.mmemory *)

let print_machine_struct ?(ghost=false) fmt m =

  if not (fst (Machine_code_common.get_stateless_status m)) then

    (* Define struct *)

    fprintf fmt "@[<v 2>%a {%a%a@]@,};"

      (pp_machine_memtype_name ~ghost) m.mname.node_id

      (if ghost then

         (fun fmt -> function

            | [] -> pp_print_nothing fmt ()

            | _ -> fprintf fmt "@,%a _reg;"

                     pp_machine_regtype_name m.mname.node_id)

       else

         pp_print_list

           ~pp_open_box:pp_open_vbox0

           ~pp_prologue:(fun fmt () ->

               fprintf fmt "@,@[%a {" pp_machine_regtype_name m.mname.node_id)

           ~pp_sep:pp_print_semicolon

           ~pp_eol:pp_print_semicolon'

           ~pp_epilogue:(fun fmt () -> fprintf fmt "}@] _reg;")

           pp_c_decl_struct_var)

      m.mmemory

      (pp_print_list

         ~pp_open_box:pp_open_vbox0

         ~pp_prologue:pp_print_cut

         ~pp_sep:pp_print_semicolon

         ~pp_eol:pp_print_semicolon'

         (pp_c_decl_instance_var ~ghost))

      m.minstances

(********************************************************************************************)

(*                      Prototype Printing functions                                        *)

(********************************************************************************************)

let print_global_init_prototype fmt baseNAME =

  fprintf fmt "void %a ()"

    pp_global_init_name baseNAME

let print_global_clear_prototype fmt baseNAME =

  fprintf fmt "void %a ()"

    pp_global_clear_name baseNAME

let print_alloc_prototype fmt (name, static) =

  fprintf fmt "%a * %a %a"

    (pp_machine_memtype_name ~ghost:false) name

    pp_machine_alloc_name name

    (pp_print_parenthesized pp_c_decl_input_var) static

let print_dealloc_prototype fmt name =

  fprintf fmt "void %a (%a * _alloc)"

    pp_machine_dealloc_name name

    (pp_machine_memtype_name ~ghost:false) name

let print_reset_prototype self fmt (name, static) =

  fprintf fmt "void %a (%a%a *%s)"

    pp_machine_reset_name name

    (pp_print_list ~pp_sep:pp_print_comma ~pp_eol:pp_print_comma

       pp_c_decl_input_var) static

    (pp_machine_memtype_name ~ghost:false) name

    self

let print_init_prototype self fmt (name, static) =

  fprintf fmt "void %a (%a%a *%s)"

    pp_machine_init_name name

    (pp_print_list ~pp_sep:pp_print_comma ~pp_eol:pp_print_comma

       pp_c_decl_input_var) static

    (pp_machine_memtype_name ~ghost:false) name

    self

let print_clear_prototype self fmt (name, static) =

  fprintf fmt "void %a (%a%a *%s)"

    pp_machine_clear_name name

    (pp_print_list ~pp_sep:pp_print_comma ~pp_eol:pp_print_comma

       pp_c_decl_input_var) static

    (pp_machine_memtype_name ~ghost:false) name

    self

let print_stateless_prototype fmt (name, inputs, outputs) =

  fprintf fmt "void %a (@[<v>%a%a@])"

    pp_machine_step_name name

    (pp_print_list ~pp_sep:pp_print_comma ~pp_eol:pp_print_comma

       ~pp_epilogue:pp_print_cut pp_c_decl_input_var) inputs

    (pp_print_list ~pp_sep:pp_print_comma pp_c_decl_output_var) outputs

let print_step_prototype self fmt (name, inputs, outputs) =

  fprintf fmt "void %a (@[<v>%a%a%a *%s@])"

    pp_machine_step_name name

    (pp_print_list ~pp_sep:pp_print_comma ~pp_eol:pp_print_comma

       ~pp_epilogue:pp_print_cut pp_c_decl_input_var) inputs

    (pp_print_list ~pp_sep:pp_print_comma ~pp_eol:pp_print_comma

       ~pp_epilogue:pp_print_cut pp_c_decl_output_var) outputs

    (pp_machine_memtype_name ~ghost:false) name

    self

let print_import_prototype fmt dep =

  fprintf fmt "#include \"%s.h\"" dep.name

let print_import_alloc_prototype fmt dep =

  if dep.is_stateful then

    fprintf fmt "#include \"%s_alloc.h\"" dep.name

let pp_c_var m self pp_var fmt var =

    pp_c_val m self pp_var fmt (Machine_code_common.mk_val (Var var) var.var_type)

let pp_array_suffix =

  pp_print_list ~pp_sep:pp_print_nothing (fun fmt v -> fprintf fmt "[%s]" v)

let mpfr_vars vars =

  if !Options.mpfr then

    List.filter (fun v -> Types.(is_real_type (array_base_type v.var_type))) vars

  else []

let mpfr_consts consts =

  if !Options.mpfr then

    List.filter (fun c -> Types.(is_real_type (array_base_type c.const_type))) consts

  else []

(* type directed initialization: useless wrt the lustre compilation model,

   except for MPFR injection, where values are dynamically allocated

*)

let pp_initialize m self pp_var fmt var =

  let rec aux indices fmt typ =

    if Types.is_array_type typ

    then

      let dim = Types.array_type_dimension typ in

      let idx = mk_loop_var m () in

      fprintf fmt "@[<v 2>{@,int %s;@,for(%s=0;%s<%a;%s++)@,%a @]@,}"

        idx idx idx pp_c_dimension dim idx

        (aux (idx::indices)) (Types.array_element_type typ)

    else

      let indices = List.rev indices in

      let pp_var_suffix fmt var =

        fprintf fmt "%a%a" (pp_c_var m self pp_var) var pp_array_suffix indices in

      Mpfr.pp_inject_init pp_var_suffix fmt var

  in

  reset_loop_counter ();

  aux [] fmt var.var_type

(* type directed clear: useless wrt the lustre compilation model,

   except for MPFR injection, where values are dynamically allocated

*)

let pp_clear m self pp_var fmt var =

  let rec aux indices fmt typ =

    if Types.is_array_type typ

    then

      let dim = Types.array_type_dimension typ in

      let idx = mk_loop_var m () in

      fprintf fmt "@[<v 2>{@,int %s;@,for(%s=0;%s<%a;%s++)@,%a @]@,}"

        idx idx idx pp_c_dimension dim idx

        (aux (idx::indices)) (Types.array_element_type typ)

    else

      let indices = List.rev indices in

      let pp_var_suffix fmt var =

        fprintf fmt "%a%a" (pp_c_var m self pp_var) var pp_array_suffix indices in

      Mpfr.pp_inject_clear pp_var_suffix fmt var

  in

  reset_loop_counter ();

  aux [] fmt var.var_type

  (*** Common functions for main ***)

let pp_print_file file_suffix fmt (typ, arg) =

  fprintf fmt

    "@[<v 2>if (traces) {@,\

     fprintf(f_%s, \"%%%s\\n\", %s);@,\

     fflush(f_%s);@]@,\

     }"

    file_suffix typ arg

    file_suffix

let print_put_var fmt file_suffix name var_type var_id =

  let pp_file = pp_print_file ("out" ^ file_suffix) in

  let unclocked_t = Types.unclock_type var_type in

  fprintf fmt "@[<v>%a@]"

    (fun fmt () ->

       if Types.is_int_type unclocked_t then

         fprintf fmt "_put_int(\"%s\", %s);@,%a"

           name var_id

           pp_file ("d", var_id)

       else if Types.is_bool_type unclocked_t then

         fprintf fmt "_put_bool(\"%s\", %s);@,%a"

           name var_id

           pp_file ("i", var_id)

       else if Types.is_real_type unclocked_t then

         if !Options.mpfr then

           fprintf fmt "_put_double(\"%s\", mpfr_get_d(%s, %s), %i);@,%a"

             name var_id (Mpfr.mpfr_rnd ()) !Options.print_prec_double

             pp_file (".*f",

                      string_of_int !Options.print_prec_double

                      ^ ", mpfr_get_d(" ^ var_id ^ ", MPFR_RNDN)")

         else

           fprintf fmt "_put_double(\"%s\", %s, %i);@,%a"

             name var_id !Options.print_prec_double

             pp_file (".*f",

                      string_of_int !Options.print_prec_double ^ ", " ^ var_id)

       else begin

         eprintf "Impossible to print the _put_xx for type %a@.@?"

           Types.print_ty var_type;

         assert false

       end) ()

let pp_file_decl fmt inout idx =

  let idx = idx + 1 in (* we start from 1: in1, in2, ... *)

  fprintf fmt "FILE *f_%s%i;" inout idx

let pp_file_open fmt inout idx =

  let idx = idx + 1 in (* we start from 1: in1, in2, ... *)

  fprintf fmt

    "@[<v>const char* cst_char_suffix_%s%i = \"_simu.%s%i\";@,\

     size_t l%s%i = strlen(dir) + strlen(prefix) + strlen(cst_char_suffix_%s%i);@,\

     char* f_%s%i_name = malloc((l%s%i+2) * sizeof(char));@,\

     strcpy (f_%s%i_name, dir);@,\

     strcat(f_%s%i_name, \"/\");@,\

     strcat(f_%s%i_name, prefix);@,\

     strcat(f_%s%i_name, cst_char_suffix_%s%i);@,\

     f_%s%i = fopen(f_%s%i_name, \"w\");@,\

     free(f_%s%i_name);\

     @]"

    inout idx inout idx

    inout idx inout idx

    inout idx inout idx

    inout idx

    inout idx

    inout idx

    inout idx inout idx

    inout idx inout idx

    inout idx;

  "f_" ^ inout ^ string_of_int idx

(* Local Variables: *)

(* compile-command:"make -C ../../.." *)

(* End: *)