Source file matrix.ml

(* File generated from matrix.idl *)

type t
and equation = {
  var : int;
  expr : Vector.t
}

external make : int -> int -> t
	= "camlidl_matrix_matrix_make"

let dummy = make 0 0
external copy : t -> t
	= "camlidl_matrix_matrix_copy"

external _print : t -> unit
	= "camlidl_matrix_matrix__print"

external compare : t -> t -> int
	= "camlidl_matrix_matrix_compare"

external hash : t -> int
	= "camlidl_matrix_matrix_hash"

external compare_sort : t -> t -> int
	= "camlidl_matrix_matrix_compare_sort"

external hash_sort : t -> int
	= "camlidl_matrix_matrix_hash_sort"

external get : t -> int -> int -> int
	= "camlidl_matrix_matrix_get"

external get_str10 : t -> int -> int -> string
	= "camlidl_matrix_matrix_get_str10"

external set : t -> int -> int -> int -> unit
	= "camlidl_matrix_matrix_set"

external set_str10 : t -> int -> int -> string -> unit
	= "camlidl_matrix_matrix_set_str10"

external get_row : t -> int -> Vector.t
	= "camlidl_matrix_matrix_get_row"

external set_row : t -> int -> Vector.t -> unit
	= "camlidl_matrix_matrix_set_row"

external nbrows : t -> int
	= "camlidl_matrix_matrix_nbrows"

external nbcolumns : t -> int
	= "camlidl_matrix_matrix_nbcolumns"

external sort_rows : t -> unit
	= "camlidl_matrix_matrix_sort_rows"

external merge_sort : t -> t -> t
	= "camlidl_matrix_matrix_merge_sort"

external assign_variable : t -> int -> Vector.t -> t
	= "camlidl_matrix_matrix_assign_variable"

external substitute_variable : t -> int -> Vector.t -> t
	= "camlidl_matrix_matrix_substitute_variable"

external assign_variables : t -> equation array -> t
	= "camlidl_matrix_matrix_assign_variables"

external substitute_variables : t -> equation array -> t
	= "camlidl_matrix_matrix_substitute_variables"

external add_dims : t -> int -> t
	= "camlidl_matrix_matrix_add_dims"

external add_dims_multi : t -> Vector.dimsup array -> t
	= "camlidl_matrix_matrix_add_dims_multi"

external del_dims_multi : t -> Vector.dimsup array -> t
	= "camlidl_matrix_matrix_del_dims_multi"

external add_permute_dims : t -> int -> int array -> t
	= "camlidl_matrix_matrix_add_permute_dims"

external permute_del_dims : t -> int -> int array -> t
	= "camlidl_matrix_matrix_permute_del_dims"

external is_row_dummy_constraint : t -> int -> bool
	= "camlidl_matrix_matrix_is_row_dummy_constraint"

 
let get_big_int mat i j = 
try Big_int.big_int_of_int (get mat i j) 
with Polka.Overflow(s) -> Big_int.big_int_of_string s 
let set_big_int mat i j v = 
if Big_int.is_int_big_int v 
then set mat i j (Big_int.int_of_big_int v) 
else set_str10 mat i j (Big_int.string_of_big_int v) 
let print formatter mat = 
let nbrows = nbrows mat in 
Format.fprintf formatter "@[<v>"; 
for i=0 to nbrows-1 do 
Format.fprintf formatter "%a@ " Vector.print (get_row mat i); 
done; 
Format.fprintf formatter "@]"; 
() 
let print_constraints assoc formatter mat =  
let nbrows = nbrows mat 
and nbcolumns = nbcolumns mat 
in 
if nbrows > !Polka.print_limit then 
Format.fprintf formatter "matrix with %i constraints" nbrows 
else begin 
let lchaines = ref [] in 
for i = nbrows-1 downto 0 do 
lchaines := (Polka.to_constraint assoc 
(fun j -> get_big_int mat i j)  
nbcolumns) :: !lchaines 
done; 
Polka.print_list formatter "{@[<hov>" ",@," "@]}"  
Format.pp_print_string 
!lchaines 
end 

let print_frames assoc formatter mat =  
let nbrows = nbrows mat                
and nbcolumns = nbcolumns mat          
in 
if nbrows > !Polka.print_limit then 
Format.fprintf formatter "matrix with %i frames" nbrows 
else begin 
let lchaines = ref [] in 
for i = nbrows-1 downto 0 do 
lchaines := (Polka.to_frame assoc 
(fun j -> get_big_int mat i j)  
nbcolumns) :: !lchaines 
done; 
Polka.print_list formatter "{@[<hov>" ",@," "@]}" 
Format.pp_print_string 
!lchaines 
end 

let of_lconstraints assoc dim lstr = 
let mat = make (List.length lstr) (!Polka.dec + dim) in 
let i = ref 0 and l = ref lstr in 
while !l<>[] do 
let vec = Vector.of_constraint assoc dim (List.hd !l) in 
set_row mat !i vec; 
incr i;  
l := List.tl !l 
done; 
mat 

let of_lframes assoc dim lstr = 
let mat = make (List.length lstr) (!Polka.dec + dim) in 
let i = ref 0 and l = ref lstr in 
while !l<>[] do 
let vec = Vector.of_frame assoc dim (List.hd !l) in 
set_row mat !i vec; 
incr i;  
l := List.tl !l 
done; 
mat