Bachelor's Thesis "Synthesizing Ranking Functions for LASW programs with Goblint"

src/analyses/apron/relationAnalysis.apron.ml

@@ -613,7 +613,7 @@ struct
         if (one_var || GobApron.Lincons1.num_vars lincons1 >= 2) && (exact || Apron.Lincons1.get_typ lincons1 <> EQ) then
           RD.cil_exp_of_lincons1 lincons1
           |> Option.map e_inv
-          |> Option.filter (InvariantCil.exp_is_suitable ~scope)
+          (*|> Option.filter (InvariantCil.exp_is_suitable ~scope)*)
         else
           None
       )

src/analyses/loopTermination.ml

@@ -5,10 +5,317 @@
 open Analyses
 open GoblintCil
 open TerminationPreprocessing
+open ListMatrix
+open SparseVector
 
 (** Contains all loop counter variables (varinfo) and maps them to their corresponding loop statement. *)
 let loop_counters : stmt VarToStmt.t ref = ref VarToStmt.empty
 
+module V = struct
+  type t = string
+  let is_int _ = true
+  let compare = String.compare
+  let print fmt v = Format.fprintf fmt "%s" v
+end
+
+module Rat = struct
+  include Mpqf
+  let zero = of_int 0
+  let one = of_int 1
+  let m_one = of_int (-1)
+  let is_zero x = cmp_int x 0 = 0
+  let is_one x = cmp_int x 1 = 0
+  let is_m_one x = cmp_int x (-1) = 0
+  let is_int x = is_one (of_mpz (get_den x)) (* I assume the numbers to be normalized *)
+  let min x y = if cmp x y <= 0 then x else y
+  let floor x =
+    let n = get_num x in
+    let d = get_den x in
+    of_mpz (Mpzf.fdiv_q n d)
+  let ceiling x =
+    let n = get_num x in
+    let d = get_den x in
+    of_mpz (Mpzf.cdiv_q n d)
+  let minus = neg
+  let sign = sgn
+  let compare = cmp
+  let mult = mul
+  let hash _ = 0
+  let get_den x = Z.of_string (Mpzf.to_string (get_den x))
+  let get_num x = Z.of_string (Mpzf.to_string (get_num x))
+end
+
+module Ex = struct
+  include Set.Make(String)
+
+  let print fmt s = match elements s with
+    | [] -> Format.fprintf fmt "()"
+    | e::l ->
+      Format.fprintf fmt "%s" e;
+      List.iter (Format.fprintf fmt ", %s") l
+end
+
+module Sim = OcplibSimplex.Basic.Make (V) (Rat) (Ex)
+
+module Matrix = ListMatrix (Rat) (SparseVector)
+
+module SparseVec = SparseVector (Rat)
+
+module VarSet = Set.Make(String)
+
+let string_of_constraints (constraints: Invariant.t) =
+  "CONSTRAINTS: " ^
+  match constraints with
+  | `Top -> "Top"
+  | `Bot -> "Bot"
+  | `Lifted constraints -> constraints |> Invariant.Exp.process |> List.map Invariant.Exp.show |> String.concat "; "
+
+let exp_list_of_constraints (constraints : Invariant.t) =
+  match constraints with
+  | `Top -> []
+  | `Bot -> failwith "Constraints are Bot"
+  | `Lifted constraints ->
+    (* Split Eq into two Le's and mirror Ge into Le *)
+    let to_leq acc = function
+      | BinOp (Le, _, _, _) as e -> e :: acc
+      | BinOp (Ge, e1, e2, t) -> (BinOp (Le, e2, e1, t)) :: acc
+      | BinOp (Eq, e1, e2, t) -> (BinOp (Le, e1, e2, t)) :: (BinOp (Le, e2, e1, t)) :: acc
+      | BinOp (Ne, _, _, _) -> acc
+      | _ -> failwith "Found something else, help me" in
+    let cs =
+      constraints |> Invariant.Exp.process |> List.map Invariant.Exp.to_cil |> List.fold_left to_leq []
+    in
+    Printf.printf "Number of constraints after conversion to Leq list: %i\n" (List.length cs);
+    cs
+
+let coeff_in_exp vname e =
+  let rec inner = function
+    | BinOp (PlusA, e1, e2, _) ->
+      (match inner e1, inner e2 with
+       | Some c1, Some c2 -> failwith "This shouldn't happen"
+       | None, None -> None
+       | c, None | None, c -> c)
+    | BinOp (MinusA, e1, e2, _) ->
+      failwith "found a Minus"
+    (*(match find_coeff v e1, find_coeff v e2 with
+      | Some c1, Some c2 -> failwith "This shouldn't happen"
+      | c, None -> c
+      | None, Some c -> Some (Rat.mult Rat.m_one c))*)
+    | BinOp (Mult, Const (CInt (c, _, _)), Lval (Var v, _), _) -> if v.vname = vname then Some (c |> cilint_to_int |> Rat.of_int) else None
+    | BinOp (Mult, _, _, _) -> failwith "found another Mult"
+    | Const (CInt (c, _, _)) -> None
+    | Lval (Var v, _) -> if v.vname = vname then Some Rat.one else None
+    | _ -> failwith "found something else"
+  in
+  match e with
+  | BinOp (Le, e1, e2, _) ->
+    (match inner e1, inner e2 with
+     | Some c1, Some c2 -> failwith "This shouldn't happen"
+     | None, None -> Rat.zero
+     | Some c, None -> c
+     | None, Some c -> Rat.neg c)
+  | _ -> failwith "This shouldn't happen"
+
+let const_in_exp e =
+  let rec inner = function
+    | BinOp (PlusA, e1, e2, _) ->
+      (match inner e1, inner e2 with
+       | Some c1, Some c2 -> failwith "This shouldn't happen"
+       | None, None -> None
+       | c, None | None, c -> c)
+    | Const (CInt (c, _, _)) -> if Z.compare c Z.zero = 0 then None else Some (c |> cilint_to_int |> Rat.of_int)
+    | BinOp (Mult, _, _, _) -> None
+    | Lval (Var v, _) -> None
+    | _ -> failwith "found something else"
+  in
+  match e with
+  | BinOp (Le, e1, e2, _) ->
+    (match inner e1, inner e2 with
+     | Some c1, Some c2 -> failwith "This shouldn't happen"
+     | None, None -> Rat.zero
+     | Some c, None -> Rat.neg c
+     | None, Some c -> c)
+  | _ -> failwith "This shouldn't happen"
+
+let rec transposed_matrices_of_exp_list (cs : exp list) =
+  let rec vars_from_exp acc = function
+    | BinOp (_, e1, e2, _) -> VarSet.union (vars_from_exp acc e1) (vars_from_exp acc e2)
+    | UnOp (_, e, _) -> vars_from_exp acc e
+    | Lval (Var v, _) ->
+      let len = String.length v.vname in
+      if String.starts_with ~prefix:"old_" v.vname then
+        VarSet.add (String.sub v.vname 4 (len - 4)) acc
+      else
+        VarSet.add v.vname acc
+    | CastE _ -> failwith "Encountered cast"
+    | Question _ -> failwith "Encountered question"
+    | _ -> acc
+  in
+  let vars = cs |> List.fold_left vars_from_exp VarSet.empty |> VarSet.elements in
+  let old_vars = List.map (String.cat "old_") vars in
+  let atr = List.fold_left (fun m v -> cs |> List.map (coeff_in_exp v) |> SparseVec.of_list |> Matrix.append_row m) (Matrix.empty ()) old_vars in
+  let a'tr = List.fold_left (fun m v -> cs |> List.map (coeff_in_exp v) |> SparseVec.of_list |> Matrix.append_row m) (Matrix.empty ()) vars in
+  let b = cs |> List.map const_in_exp |> SparseVec.of_list in
+  atr, a'tr, b
+
+
+let termination_provable a_transposed a'_transposed b =
+  let num_constraints = Matrix.num_cols a_transposed in
+  let zero, one, m_one = Sim.Core.R2.zero, Sim.Core.R2.of_r Rat.one, Sim.Core.R2.of_r Rat.m_one in
+  let ass_l1_a'tr_eq_zero sim =
+    let rec aux i sim =
+      if i >= Matrix.num_rows a'_transposed then sim else
+        let l1_ri = Matrix.get_row a'_transposed i
+          |> SparseVec.to_sparse_list
+          |> List.map (fun (idx, v) -> "x" ^ string_of_int idx, v)
+
+        in
+        l1_ri |> List.map (fun (var, value) -> (Printf.sprintf "(%s * %s)" (Rat.to_string value) var)) |> String.concat " + " |> Printf.printf "l1_r%i: %s\n" i;
+
+        let sim = fst @@
+          match l1_ri with
+          | [] -> sim, false
+          | [(var, _)] -> Sim.Assert.var sim var
+                            ~min:{Sim.Core.bvalue = zero; explanation = Ex.singleton (var ^ ">=0")}
+                            ~max:{Sim.Core.bvalue = zero; explanation = Ex.singleton (var ^ "<=0")}
+          | _ ->
+            let l1_ri = l1_ri |> Sim.Core.P.from_list in
+            Sim.Assert.poly sim l1_ri ("l1_r" ^ string_of_int i)
+              ~min:{Sim.Core.bvalue = zero; explanation = Ex.singleton ("l1_r" ^ string_of_int i ^ ">=0")}
+              ~max:{Sim.Core.bvalue = zero; explanation = Ex.singleton ("l1_r" ^ string_of_int i ^ "<=0")}
+        in aux (i + 1) sim
+    in aux 0 sim
+  in
+  let rec ass_l1_minus_l2_atr_eq_zero sim =
+    let rec aux i sim =
+      if i >= Matrix.num_rows a_transposed then sim else
+        let l1_m_l2_ri = Matrix.get_row a_transposed i
+          |> SparseVec.to_sparse_list
+          (* Example: (2, 0, -3) -> x0 * 2 + y0 * (-2) + x2 * (-3) + y2 * (-(-3))
+             map doesn't work here because we double the number of list elements *)
+          |> List.fold_left (fun l (idx, v) -> ("x" ^ string_of_int idx, v) :: ("y" ^ string_of_int idx, Rat.mul Rat.m_one v) :: l) []
+        in
+        l1_m_l2_ri |> List.map (fun (var, value) -> (Printf.sprintf "(%s * %s)" (Rat.to_string value) var)) |> String.concat " + " |> Printf.printf "l1_m_l2_r%i: %s\n" i;
+
+        let sim = fst @@
+          match l1_m_l2_ri with
+          | [] -> sim, false
+          | _ ->
+            let l1_m_l2_ri = Sim.Core.P.from_list l1_m_l2_ri in
+            Sim.Assert.poly sim l1_m_l2_ri ("l1_m_l2_r" ^ string_of_int i)
+              ~min:{Sim.Core.bvalue = zero; explanation = Ex.singleton ("l1_m_l2_r" ^ string_of_int i ^ ">=0")}
+              ~max:{Sim.Core.bvalue = zero; explanation = Ex.singleton ("l1_m_l2_r" ^ string_of_int i ^ "<=0")}
+        in aux (i + 1) sim
+    in aux 0 sim
+  in    
+  let rec ass_l2_atr_plus_a'tr_eq_zero sim =
+    let m = Matrix.add a_transposed a'_transposed in
+    let rec aux i sim =
+      if i >= Matrix.num_rows m then sim else
+        let l2_ri = Matrix.get_row m i
+          |> SparseVec.to_sparse_list
+          |> List.map (fun (idx, v) -> "y" ^ string_of_int idx, v)
+        in
+        l2_ri |> List.map (fun (var, value) -> (Printf.sprintf "(%s * %s)" (Rat.to_string value) var)) |> String.concat " + " |> Printf.printf "l2_r%i: %s\n" i;
+
+        let sim = fst @@
+          match l2_ri with
+          | [] -> sim, false
+          | [(var, _)] -> Sim.Assert.var sim var
+                            ~min:{Sim.Core.bvalue = zero; explanation = Ex.singleton (var ^ ">=0")}
+                            ~max:{Sim.Core.bvalue = zero; explanation = Ex.singleton (var ^ "<=0")}
+          | _ -> let l2_ri = l2_ri |> Sim.Core.P.from_list in
+            Sim.Assert.poly sim l2_ri ("l2_r" ^ string_of_int i)
+              ~min:{Sim.Core.bvalue = zero; explanation = Ex.singleton ("l2_r" ^ string_of_int i ^ ">=0")}
+              ~max:{Sim.Core.bvalue = zero; explanation = Ex.singleton ("l2_r" ^ string_of_int i ^ "<=0")}
+        in aux (i + 1) sim
+    in aux 0 sim
+  in
+  let ass_l2_b_lt_zero sim =
+    let l2_b = b
+      |> SparseVec.to_sparse_list
+      |> List.map (fun (i, v) -> "y" ^ string_of_int i, v)
+    in
+    l2_b |> List.map (fun (var, value) -> (Printf.sprintf "(%s * %s)" (Rat.to_string value) var)) |> String.concat " + " |> Printf.printf "l2_b: %s\n";
+
+    fst @@
+    match l2_b with
+    | [] -> sim, false
+    | [(var, value)] when Rat.compare value Rat.zero > 0 -> Sim.Assert.var sim var
+                                                              ~max:{Sim.Core.bvalue = m_one; explanation = Ex.singleton (var ^ "<0")} (* creates a contradiction *)
+    | [(var, value)] when Rat.compare value Rat.zero < 0 -> Sim.Assert.var sim var
+                                                              ~min:{Sim.Core.bvalue = one; explanation = Ex.singleton (var ^ ">0")}
+    | [(var, _)] -> failwith "This shouldn't happen"
+    | _ -> let l2_b = Sim.Core.P.from_list l2_b in
+      Sim.Assert.poly sim l2_b "l2_b"
+        ~max:{Sim.Core.bvalue = m_one; explanation = Ex.singleton "l2_b<0"}
+  in
+  let ass_li_ge_zero sim =
+    let f sim i =
+      let sim', _ =
+        Sim.Assert.var sim ("x" ^ string_of_int i)
+          ~min:{Sim.Core.bvalue = zero; explanation = Ex.singleton ("x" ^ string_of_int i ^ ">=0")}
+      in
+      fst @@ Sim.Assert.var sim' ("y" ^ string_of_int i)
+        ~min:{Sim.Core.bvalue = zero; explanation = Ex.singleton ("y" ^ string_of_int i ^ ">=0")}
+    in
+    List.fold_left f sim (List.init num_constraints Fun.id)
+  in
+  let sim = Sim.Core.empty ~is_int:true ~check_invs:true
+    |> ass_l1_a'tr_eq_zero |> ass_l1_minus_l2_atr_eq_zero |> ass_l2_atr_plus_a'tr_eq_zero |> ass_l2_b_lt_zero |> ass_li_ge_zero
+    |> Sim.Solve.solve in
+  match sim.status with
+  | UNSAT _ -> false
+  | SAT -> true
+  | UNK -> failwith "Simplex returned UNK, this shouldn't happen"
+
+let test_termination_provable () =
+  let zero = Rat.zero in
+  let one = Rat.one in
+  let two = Rat.of_int 2 in
+  let m_one = Rat.m_one in
+  let m_two = Rat.of_int (-2) in
+
+  (* From regression test 52: *)
+  (*let a_transposed = Matrix.append_row (Matrix.append_row (Matrix.append_row (Matrix.append_row (Matrix.empty ())
+    (SparseVec.of_list [m_one; zero; zero; zero; m_one; zero; zero; zero; zero; zero; one; m_one; zero; zero])) (*i*)
+    (SparseVec.of_list [zero; zero; zero; one; one; zero; zero; zero; one; one; zero; zero; one; m_one])) (*j*)
+    (SparseVec.of_list [zero; m_one; zero; zero; zero; zero; zero; zero; zero; zero; zero; zero; zero; zero])) (*nat*)
+    (SparseVec.of_list [zero; zero; m_one; zero; zero; zero; zero; zero; zero; zero; zero; zero; zero; zero]) (*pos*)
+    in
+    let a'_transposed = Matrix.append_row (Matrix.append_row (Matrix.append_row (Matrix.append_row (Matrix.empty ())
+    (SparseVec.of_list [zero; zero; zero; zero; zero; zero; zero; zero; zero; m_one; m_one; one; zero; zero])) (*i'*)
+    (SparseVec.of_list [zero; zero; zero; zero; zero; m_one; zero; zero; m_one; zero; zero; zero; m_one; one])) (*j'*)
+    (SparseVec.of_list [zero; zero; zero; zero; zero; zero; m_one; zero; zero; m_one; m_one; one; zero; zero])) (*nat'*)
+    (SparseVec.of_list [zero; zero; zero; zero; zero; zero; zero; m_one; zero; zero; zero; zero; one; m_one]) (*pos'*)
+    in
+    let b = SparseVec.of_list [Rat.of_int 2147483647; zero; m_one; Rat.of_int 2147483646; m_one;
+    Rat.of_int 2147483647; zero; m_one; m_one; m_one;
+    zero; zero; zero; zero]*)
+
+  (* From regression test 53: *)
+  (*let a_transposed = Matrix.append_row (Matrix.append_row (Matrix.empty ())
+    (SparseVec.of_list [zero; m_one; zero; zero; zero; m_one; one; zero; zero])) (*x1*)
+    (SparseVec.of_list [m_one; zero; zero; zero; zero; zero; zero; m_one; one]) (*x2*)
+    in
+    let a'_transposed = Matrix.append_row (Matrix.append_row (Matrix.empty ())
+    (SparseVec.of_list [zero; zero; m_one; zero; one; two; m_two; zero; zero])) (*x1'*)
+    (SparseVec.of_list [zero; zero; zero; m_one; zero; zero; zero; one; m_one]) (*x2'*)
+    in
+    let b = SparseVec.of_list [zero; m_two; m_one; m_one; Rat.of_int 1073741823;
+    zero; two; one; m_one]*)
+
+  (* From regression test 53, simplified: *)
+  let a_transposed = Matrix.append_row (Matrix.empty ())
+      (SparseVec.of_list [m_one; zero; m_one; one]) (*x1*)
+  in
+  let a'_transposed = Matrix.append_row (Matrix.empty ())
+      (SparseVec.of_list [zero; one; two; m_two]) (*x1'*)
+  in
+  let b = SparseVec.of_list [m_one; Rat.of_int 1073741823; zero; two]
+  in
+  termination_provable a_transposed a'_transposed b
+
 (** Checks whether a variable can be bounded. *)
 let ask_bound man varinfo =
   let open IntDomain.IntDomTuple in
@@ -62,7 +369,13 @@
                   M.success ~category:Termination ~loc "Loop terminates: bounded by %a iteration(s)" IntDomain.IntDomTuple.pretty bound;
             end
           | None ->
-            failwith "Encountered a call to __goblint_bounded with an unknown loop counter variable."
+            (*failwith "Encountered a call to __goblint_bounded with an unknown loop counter variable."*)
+            Printf.printf "%s\n" (string_of_constraints (man.ask (Queries.Invariant Invariant.default_context)));
+            let atr, a'tr, b = man.ask (Queries.Invariant Invariant.default_context) |> exp_list_of_constraints |> transposed_matrices_of_exp_list in
+            if Matrix.is_empty atr then Printf.printf "nothing to see here" else
+              let term_provable = termination_provable atr a'tr b in
+              Printf.printf "\n%s\n%s\n%s\nTermination%s provable\n" (Matrix.show atr) (Matrix.show a'tr) (SparseVec.show b) (if term_provable then "" else " not");
+
         end
       | "__goblint_bounded", _ ->
         failwith "__goblint_bounded call unexpected arguments"

src/cdomains/affineEquality/sparseImplementation/listMatrix.ml

@@ -11,6 +11,7 @@ let timing_wrap = Vector.timing_wrap
 module type SparseMatrix =
 sig
   include Matrix
+  val add: t -> t -> t
   val get_col_upper_triangular: t -> int -> vec
 
   val swap_rows: t -> int -> int -> t
@@ -508,4 +509,6 @@ module ListMatrix: SparseMatrixFunctor =
       if M.tracing then M.tracel "linear_disjunct" "linear_disjunct between \n%s and \n%s =>\n%s" (show m1)  (show m2) (show res);
       res
 
+  let add = List.map2 V.add
+
   end

src/cdomains/affineEquality/sparseImplementation/sparseVector.ml

@@ -8,6 +8,7 @@ open Batteries
 module type SparseVector =
 sig
   include Vector
+  val add: t -> t -> t
   val push_first: t -> int -> num -> t
 
   val is_zero_vec: t -> bool
@@ -367,4 +368,18 @@ module SparseVector: SparseVectorFunctor =
     let rev v =
       let entries = List.rev_map (fun (idx, value) -> (v.len - 1 - idx, value)) v.entries in
       {entries; len = v.len}
+
+    let add a b =
+      if a.len <> b.len then raise (Invalid_argument "Unequal lengths") else
+        let rec aux a b =
+          match a, b with
+          | [], b -> b
+          | a, [] -> a
+          | (ix, vx)::xs, (iy, vy)::ys when ix = iy ->
+            let sum = A.add vx vy in
+            if A.compare sum A.zero = 0 then aux xs ys
+            else (ix, A.add vx vy) :: aux xs ys
+          | (ix, vx)::xs, (iy, _)::_ when ix < iy -> (ix, vx) :: aux xs b
+          | _, y::ys -> y :: aux a ys
+        in {entries = aux a.entries b.entries; len = a.len}
   end