Commit df87c2b0 authored by Heiko Becker's avatar Heiko Becker
Browse files

Remove temporary files

parent 77e7f380
open Interval
open Opt_func
let start_interval = Array.init 3 (function
| 0 -> {low = 3.99999999999999966693e-01; high = 9.00000000000000022204e-01}
| 1 -> {low = 3.79999999999999982236e+00; high = 7.80000000000000071054e+00}
| 2 -> {low = -4.50000000000000000000e+00; high = -2.99999999999999988898e-01}
| _ -> failwith "Out of boundaries"
)
let f_X input_array =
let var_w0 = input_array.(0) in
let var_r0 = input_array.(1) in
let var_v0 = input_array.(2) in
((({low = 3.00000000000000000000e+00; high = 3.00000000000000000000e+00} +$ ({low = 2.00000000000000000000e+00; high = 2.00000000000000000000e+00} *$ ({low = 1.00000000000000000000e+00; high = 1.00000000000000000000e+00} /$ (var_r0 *$ var_r0)))) -$ (({low = 1.25000000000000000000e-01; high = 1.25000000000000000000e-01} *$ (({low = 3.00000000000000000000e+00; high = 3.00000000000000000000e+00} -$ ({low = 2.00000000000000000000e+00; high = 2.00000000000000000000e+00} *$ var_v0)) *$ (((var_w0 *$ var_w0) *$ var_r0) *$ var_r0))) *$ ({low = 1.00000000000000000000e+00; high = 1.00000000000000000000e+00} /$ ({low = 1.00000000000000000000e+00; high = 1.00000000000000000000e+00} -$ var_v0)))) -$ {low = 4.50000000000000000000e+00; high = 4.50000000000000000000e+00})
let _ =
let x_tol = size_max_X start_interval *. 0.000000e+00 +. 1.000000e-02 in
let upper_bound, lower_bound, c = Opt0.opt f_X start_interval x_tol (1.000000e-02) (1.000000e-02) (1000000) in
let () = Printf.printf "iter_max = %d\n" c in
let () = Printf.printf "max = %0.20e\n" upper_bound in
let () = Printf.printf "lower_max = %0.20e\n" lower_bound in
let upper_bound, lower_bound, c = Opt0.opt (fun x -> ~-$ (f_X x)) start_interval x_tol (1.000000e-02) (1.000000e-02) (1000000) in
let () = Printf.printf "iter_min = %d\n" c in
let () = Printf.printf "min = %0.20e\n" (-. upper_bound) in
let () = Printf.printf "lower_min = %0.20e\n" (-. lower_bound) in
flush stdout
open Interval
open Opt_func
let start_interval = Array.init 0 (function
| _ -> failwith "Out of boundaries"
)
let f_X input_array =
{low = 5.32817243572231868460e-13; high = 5.32817243572231868460e-13}
let _ =
let x_tol = size_max_X start_interval *. 0.000000e+00 +. 1.000000e-02 in
let upper_bound, lower_bound, c = Opt0.opt f_X start_interval x_tol (1.000000e-02) (1.000000e-02) (1000000) in
let () = Printf.printf "iter_max = %d\n" c in
let () = Printf.printf "max = %0.20e\n" upper_bound in
let () = Printf.printf "lower_max = %0.20e\n" lower_bound in
let upper_bound, lower_bound, c = Opt0.opt (fun x -> ~-$ (f_X x)) start_interval x_tol (1.000000e-02) (1.000000e-02) (1000000) in
let () = Printf.printf "iter_min = %d\n" c in
let () = Printf.printf "min = %0.20e\n" (-. upper_bound) in
let () = Printf.printf "lower_min = %0.20e\n" (-. lower_bound) in
flush stdout
open Interval
open Opt_func
let start_interval = Array.init 3 (function
| 0 -> {low = 3.99999999999999966693e-01; high = 9.00000000000000022204e-01}
| 1 -> {low = 3.79999999999999982236e+00; high = 7.80000000000000071054e+00}
| 2 -> {low = -4.50000000000000000000e+00; high = -2.99999999999999988898e-01}
| _ -> failwith "Out of boundaries"
)
let f_X input_array =
let var_w0 = input_array.(0) in
let var_r0 = input_array.(1) in
let var_v0 = input_array.(2) in
(abs_I((((({low = 2.00000000000000000000e+00; high = 2.00000000000000000000e+00} *$ (~-$(((var_r0 *$ floor_power2_I(var_r0)) /$ ((var_r0 *$ var_r0) *$ (var_r0 *$ var_r0)))))) +$ ({low = 2.00000000000000000000e+00; high = 2.00000000000000000000e+00} *$ (~-$(((var_r0 *$ floor_power2_I(var_r0)) /$ ((var_r0 *$ var_r0) *$ (var_r0 *$ var_r0))))))) +$ (~-$((({low = 1.00000000000000000000e+00; high = 1.00000000000000000000e+00} /$ ({low = 1.00000000000000000000e+00; high = 1.00000000000000000000e+00} -$ var_v0)) *$ ({low = 1.25000000000000000000e-01; high = 1.25000000000000000000e-01} *$ (({low = 3.00000000000000000000e+00; high = 3.00000000000000000000e+00} -$ ({low = 2.00000000000000000000e+00; high = 2.00000000000000000000e+00} *$ var_v0)) *$ (((var_w0 *$ var_w0) *$ var_r0) *$ floor_power2_I(var_r0)))))))) +$ (~-$((({low = 1.00000000000000000000e+00; high = 1.00000000000000000000e+00} /$ ({low = 1.00000000000000000000e+00; high = 1.00000000000000000000e+00} -$ var_v0)) *$ ({low = 1.25000000000000000000e-01; high = 1.25000000000000000000e-01} *$ (({low = 3.00000000000000000000e+00; high = 3.00000000000000000000e+00} -$ ({low = 2.00000000000000000000e+00; high = 2.00000000000000000000e+00} *$ var_v0)) *$ (var_r0 *$ ((var_w0 *$ var_w0) *$ floor_power2_I(var_r0)))))))))) +$ (abs_I(({low = 2.00000000000000000000e+00; high = 2.00000000000000000000e+00} *$ (~-$((floor_power2_I(((var_r0 *$ var_r0) +$ {low = -6.92779167366097823219e-15; high = 6.92779167366097823219e-15})) /$ ((var_r0 *$ var_r0) *$ (var_r0 *$ var_r0))))))) +$ (abs_I(floor_power2_I((({low = 2.00000000000000000000e+00; high = 2.00000000000000000000e+00} *$ ({low = 1.00000000000000000000e+00; high = 1.00000000000000000000e+00} /$ (var_r0 *$ var_r0))) +$ {low = -1.00525868360255573858e-16; high = 1.00525868360255573858e-16}))) +$ (abs_I(floor_power2_I((({low = 3.00000000000000000000e+00; high = 3.00000000000000000000e+00} +$ ({low = 2.00000000000000000000e+00; high = 2.00000000000000000000e+00} *$ ({low = 1.00000000000000000000e+00; high = 1.00000000000000000000e+00} /$ (var_r0 *$ var_r0)))) +$ {low = -1.14403656168070042939e-16; high = 1.14403656168070042939e-16}))) +$ (abs_I(((~-$((({low = 1.25000000000000000000e-01; high = 1.25000000000000000000e-01} *$ (({low = 3.00000000000000000000e+00; high = 3.00000000000000000000e+00} -$ ({low = 2.00000000000000000000e+00; high = 2.00000000000000000000e+00} *$ var_v0)) *$ (((var_w0 *$ var_w0) *$ var_r0) *$ var_r0))) *$ (~-$(((~-$(floor_power2_I(var_v0))) /$ (({low = 1.00000000000000000000e+00; high = 1.00000000000000000000e+00} -$ var_v0) *$ ({low = 1.00000000000000000000e+00; high = 1.00000000000000000000e+00} -$ var_v0)))))))) +$ (~-$((({low = 1.00000000000000000000e+00; high = 1.00000000000000000000e+00} /$ ({low = 1.00000000000000000000e+00; high = 1.00000000000000000000e+00} -$ var_v0)) *$ ({low = 1.25000000000000000000e-01; high = 1.25000000000000000000e-01} *$ ((((var_w0 *$ var_w0) *$ var_r0) *$ var_r0) *$ (~-$(({low = 2.00000000000000000000e+00; high = 2.00000000000000000000e+00} *$ floor_power2_I(var_v0))))))))))) +$ (abs_I((~-$((({low = 1.00000000000000000000e+00; high = 1.00000000000000000000e+00} /$ ({low = 1.00000000000000000000e+00; high = 1.00000000000000000000e+00} -$ var_v0)) *$ ({low = 1.25000000000000000000e-01; high = 1.25000000000000000000e-01} *$ ((((var_w0 *$ var_w0) *$ var_r0) *$ var_r0) *$ floor_power2_I((({low = 3.00000000000000000000e+00; high = 3.00000000000000000000e+00} -$ ({low = 2.00000000000000000000e+00; high = 2.00000000000000000000e+00} *$ var_v0)) +$ {low = -8.88178419700125232339e-16; high = 8.88178419700125232339e-16})))))))) +$ (abs_I(((~-$((({low = 1.00000000000000000000e+00; high = 1.00000000000000000000e+00} /$ ({low = 1.00000000000000000000e+00; high = 1.00000000000000000000e+00} -$ var_v0)) *$ ({low = 1.25000000000000000000e-01; high = 1.25000000000000000000e-01} *$ (({low = 3.00000000000000000000e+00; high = 3.00000000000000000000e+00} -$ ({low = 2.00000000000000000000e+00; high = 2.00000000000000000000e+00} *$ var_v0)) *$ (var_r0 *$ (var_r0 *$ (var_w0 *$ floor_power2_I(var_w0))))))))) +$ (~-$((({low = 1.00000000000000000000e+00; high = 1.00000000000000000000e+00} /$ ({low = 1.00000000000000000000e+00; high = 1.00000000000000000000e+00} -$ var_v0)) *$ ({low = 1.25000000000000000000e-01; high = 1.25000000000000000000e-01} *$ (({low = 3.00000000000000000000e+00; high = 3.00000000000000000000e+00} -$ ({low = 2.00000000000000000000e+00; high = 2.00000000000000000000e+00} *$ var_v0)) *$ (var_r0 *$ (var_r0 *$ (var_w0 *$ floor_power2_I(var_w0))))))))))) +$ (abs_I((~-$((({low = 1.00000000000000000000e+00; high = 1.00000000000000000000e+00} /$ ({low = 1.00000000000000000000e+00; high = 1.00000000000000000000e+00} -$ var_v0)) *$ ({low = 1.25000000000000000000e-01; high = 1.25000000000000000000e-01} *$ (({low = 3.00000000000000000000e+00; high = 3.00000000000000000000e+00} -$ ({low = 2.00000000000000000000e+00; high = 2.00000000000000000000e+00} *$ var_v0)) *$ (var_r0 *$ (var_r0 *$ floor_power2_I(((var_w0 *$ var_w0) +$ {low = -9.99200722162641034293e-17; high = 9.99200722162641034293e-17})))))))))) +$ (abs_I((~-$((({low = 1.00000000000000000000e+00; high = 1.00000000000000000000e+00} /$ ({low = 1.00000000000000000000e+00; high = 1.00000000000000000000e+00} -$ var_v0)) *$ ({low = 1.25000000000000000000e-01; high = 1.25000000000000000000e-01} *$ (({low = 3.00000000000000000000e+00; high = 3.00000000000000000000e+00} -$ ({low = 2.00000000000000000000e+00; high = 2.00000000000000000000e+00} *$ var_v0)) *$ (var_r0 *$ floor_power2_I((((var_w0 *$ var_w0) *$ var_r0) +$ {low = -1.57207580286922207145e-15; high = 1.57207580286922207145e-15}))))))))) +$ (abs_I((~-$((({low = 1.00000000000000000000e+00; high = 1.00000000000000000000e+00} /$ ({low = 1.00000000000000000000e+00; high = 1.00000000000000000000e+00} -$ var_v0)) *$ ({low = 1.25000000000000000000e-01; high = 1.25000000000000000000e-01} *$ (({low = 3.00000000000000000000e+00; high = 3.00000000000000000000e+00} -$ ({low = 2.00000000000000000000e+00; high = 2.00000000000000000000e+00} *$ var_v0)) *$ floor_power2_I(((((var_w0 *$ var_w0) *$ var_r0) *$ var_r0) +$ {low = -1.85318427270431242061e-14; high = 1.85318427270431242061e-14})))))))) +$ (abs_I((~-$((({low = 1.00000000000000000000e+00; high = 1.00000000000000000000e+00} /$ ({low = 1.00000000000000000000e+00; high = 1.00000000000000000000e+00} -$ var_v0)) *$ ({low = 1.25000000000000000000e-01; high = 1.25000000000000000000e-01} *$ floor_power2_I(((({low = 3.00000000000000000000e+00; high = 3.00000000000000000000e+00} -$ ({low = 2.00000000000000000000e+00; high = 2.00000000000000000000e+00} *$ var_v0)) *$ (((var_w0 *$ var_w0) *$ var_r0) *$ var_r0)) +$ {low = -3.52554252458503743828e-13; high = 3.52554252458503743828e-13}))))))) +$ (abs_I((~-$((({low = 1.25000000000000000000e-01; high = 1.25000000000000000000e-01} *$ (({low = 3.00000000000000000000e+00; high = 3.00000000000000000000e+00} -$ ({low = 2.00000000000000000000e+00; high = 2.00000000000000000000e+00} *$ var_v0)) *$ (((var_w0 *$ var_w0) *$ var_r0) *$ var_r0))) *$ (~-$((floor_power2_I((({low = 1.00000000000000000000e+00; high = 1.00000000000000000000e+00} -$ var_v0) +$ {low = -4.44089209850062616169e-16; high = 4.44089209850062616169e-16})) /$ (({low = 1.00000000000000000000e+00; high = 1.00000000000000000000e+00} -$ var_v0) *$ ({low = 1.00000000000000000000e+00; high = 1.00000000000000000000e+00} -$ var_v0))))))))) +$ (abs_I((~-$(floor_power2_I(((({low = 1.25000000000000000000e-01; high = 1.25000000000000000000e-01} *$ (({low = 3.00000000000000000000e+00; high = 3.00000000000000000000e+00} -$ ({low = 2.00000000000000000000e+00; high = 2.00000000000000000000e+00} *$ var_v0)) *$ (((var_w0 *$ var_w0) *$ var_r0) *$ var_r0))) *$ ({low = 1.00000000000000000000e+00; high = 1.00000000000000000000e+00} /$ ({low = 1.00000000000000000000e+00; high = 1.00000000000000000000e+00} -$ var_v0))) +$ {low = -7.82140847814635974942e-14; high = 7.82140847814635974942e-14}))))) +$ (abs_I(floor_power2_I(((({low = 3.00000000000000000000e+00; high = 3.00000000000000000000e+00} +$ ({low = 2.00000000000000000000e+00; high = 2.00000000000000000000e+00} *$ ({low = 1.00000000000000000000e+00; high = 1.00000000000000000000e+00} /$ (var_r0 *$ var_r0)))) -$ (({low = 1.25000000000000000000e-01; high = 1.25000000000000000000e-01} *$ (({low = 3.00000000000000000000e+00; high = 3.00000000000000000000e+00} -$ ({low = 2.00000000000000000000e+00; high = 2.00000000000000000000e+00} *$ var_v0)) *$ (((var_w0 *$ var_w0) *$ var_r0) *$ var_r0))) *$ ({low = 1.00000000000000000000e+00; high = 1.00000000000000000000e+00} /$ ({low = 1.00000000000000000000e+00; high = 1.00000000000000000000e+00} -$ var_v0)))) +$ {low = -8.21032467213572644859e-14; high = 8.21032467213572644859e-14}))) +$ abs_I(floor_power2_I((((({low = 3.00000000000000000000e+00; high = 3.00000000000000000000e+00} +$ ({low = 2.00000000000000000000e+00; high = 2.00000000000000000000e+00} *$ ({low = 1.00000000000000000000e+00; high = 1.00000000000000000000e+00} /$ (var_r0 *$ var_r0)))) -$ (({low = 1.25000000000000000000e-01; high = 1.25000000000000000000e-01} *$ (({low = 3.00000000000000000000e+00; high = 3.00000000000000000000e+00} -$ ({low = 2.00000000000000000000e+00; high = 2.00000000000000000000e+00} *$ var_v0)) *$ (((var_w0 *$ var_w0) *$ var_r0) *$ var_r0))) *$ ({low = 1.00000000000000000000e+00; high = 1.00000000000000000000e+00} /$ ({low = 1.00000000000000000000e+00; high = 1.00000000000000000000e+00} -$ var_v0)))) -$ {low = 4.50000000000000000000e+00; high = 4.50000000000000000000e+00}) +$ {low = -8.56559604001577654152e-14; high = 8.56559604001577654152e-14})))))))))))))))))
let _ =
let x_tol = size_max_X start_interval *. 0.000000e+00 +. 1.000000e-02 in
let upper_bound, lower_bound, c = Opt0.opt f_X start_interval x_tol (1.000000e-02) (1.000000e-02) (1000000) in
let () = Printf.printf "iter_max = %d\n" c in
let () = Printf.printf "max = %0.20e\n" upper_bound in
let () = Printf.printf "lower_max = %0.20e\n" lower_bound in
let () = Printf.printf "iter_min = 0\n" in
let () = Printf.printf "min = 0\n" in
let () = Printf.printf "lower_min = 0\n" in
flush stdout
open Interval
open Opt_func
let start_interval = Array.init 3 (function
| 0 -> {low = 3.99999999999999966693e-01; high = 9.00000000000000022204e-01}
| 1 -> {low = 3.79999999999999982236e+00; high = 7.80000000000000071054e+00}
| 2 -> {low = -4.50000000000000000000e+00; high = -2.99999999999999988898e-01}
| _ -> failwith "Out of boundaries"
)
let f_X input_array =
let var_w0 = input_array.(0) in
let var_r0 = input_array.(1) in
let var_v0 = input_array.(2) in
((({low = 6.00000000000000000000e+00; high = 6.00000000000000000000e+00} *$ var_v0) -$ (({low = 5.00000000000000000000e-01; high = 5.00000000000000000000e-01} *$ (var_v0 *$ (((var_w0 *$ var_w0) *$ var_r0) *$ var_r0))) *$ ({low = 1.00000000000000000000e+00; high = 1.00000000000000000000e+00} /$ ({low = 1.00000000000000000000e+00; high = 1.00000000000000000000e+00} -$ var_v0)))) -$ {low = 2.50000000000000000000e+00; high = 2.50000000000000000000e+00})
let _ =
let x_tol = size_max_X start_interval *. 0.000000e+00 +. 1.000000e-02 in
let upper_bound, lower_bound, c = Opt0.opt f_X start_interval x_tol (1.000000e-02) (1.000000e-02) (1000000) in
let () = Printf.printf "iter_max = %d\n" c in
let () = Printf.printf "max = %0.20e\n" upper_bound in
let () = Printf.printf "lower_max = %0.20e\n" lower_bound in
let upper_bound, lower_bound, c = Opt0.opt (fun x -> ~-$ (f_X x)) start_interval x_tol (1.000000e-02) (1.000000e-02) (1000000) in
let () = Printf.printf "iter_min = %d\n" c in
let () = Printf.printf "min = %0.20e\n" (-. upper_bound) in
let () = Printf.printf "lower_min = %0.20e\n" (-. lower_bound) in
flush stdout
open Interval
open Opt_func
let start_interval = Array.init 0 (function
| _ -> failwith "Out of boundaries"
)
let f_X input_array =
{low = 7.38509559225976221917e-13; high = 7.38509559225976221917e-13}
let _ =
let x_tol = size_max_X start_interval *. 0.000000e+00 +. 1.000000e-02 in
let upper_bound, lower_bound, c = Opt0.opt f_X start_interval x_tol (1.000000e-02) (1.000000e-02) (1000000) in
let () = Printf.printf "iter_max = %d\n" c in
let () = Printf.printf "max = %0.20e\n" upper_bound in
let () = Printf.printf "lower_max = %0.20e\n" lower_bound in
let upper_bound, lower_bound, c = Opt0.opt (fun x -> ~-$ (f_X x)) start_interval x_tol (1.000000e-02) (1.000000e-02) (1000000) in
let () = Printf.printf "iter_min = %d\n" c in
let () = Printf.printf "min = %0.20e\n" (-. upper_bound) in
let () = Printf.printf "lower_min = %0.20e\n" (-. lower_bound) in
flush stdout
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment