ampl/tp057.mod


If Cut/Copy and Paste fails, then click here for download.


# Model tp057
# Source version 1
# Simple APM canonicalizer version 1.3
# AMPL back end version 1.1
param a_1 := 8;
param a_2 := 8;
param a_3 := 10;
param a_4 := 10;
param a_5 := 10;
param a_6 := 10;
param a_7 := 12;
param a_8 := 12;
param a_9 := 12;
param a_10 := 12;
param a_11 := 14;
param a_12 := 14;
param a_13 := 14;
param a_14 := 16;
param a_15 := 16;
param a_16 := 16;
param a_17 := 18;
param a_18 := 18;
param a_19 := 20;
param a_20 := 20;
param a_21 := 20;
param a_22 := 22;
param a_23 := 22;
param a_24 := 22;
param a_25 := 24;
param a_26 := 24;
param a_27 := 24;
param a_28 := 26;
param a_29 := 26;
param a_30 := 26;
param a_31 := 28;
param a_32 := 28;
param a_33 := 30;
param a_34 := 30;
param a_35 := 30;
param a_36 := 32;
param a_37 := 32;
param a_38 := 34;
param a_39 := 36;
param a_40 := 36;
param a_41 := 38;
param a_42 := 38;
param a_43 := 40;
param a_44 := 42;
param b_1 := 0.49;
param b_2 := 0.49;
param b_3 := 0.48;
param b_4 := 0.47;
param b_5 := 0.48;
param b_6 := 0.47;
param b_7 := 0.46;
param b_8 := 0.46;
param b_9 := 0.45;
param b_10 := 0.43;
param b_11 := 0.45;
param b_12 := 0.43;
param b_13 := 0.43;
param b_14 := 0.44;
param b_15 := 0.43;
param b_16 := 0.43;
param b_17 := 0.46;
param b_18 := 0.45;
param b_19 := 0.42;
param b_20 := 0.42;
param b_21 := 0.43;
param b_22 := 0.41;
param b_23 := 0.41;
param b_24 := 0.40;
param b_25 := 0.42;
param b_26 := 0.40;
param b_27 := 0.40;
param b_28 := 0.41;
param b_29 := 0.40;
param b_30 := 0.41;
param b_31 := 0.41;
param b_32 := 0.40;
param b_33 := 0.40;
param b_34 := 0.40;
param b_35 := 0.38;
param b_36 := 0.41;
param b_37 := 0.40;
param b_38 := 0.40;
param b_39 := 0.41;
param b_40 := 0.38;
param b_41 := 0.40;
param b_42 := 0.40;
param b_43 := 0.39;
param b_44 := 0.39;
var x_1 >= 0.4;
let x_1 := 0.42;
var x_2 >= -4;
let x_2 := 5;
var aux_1 = b_1 - x_1 - (0.49 - x_1)*exp(-x_2*(a_1 - 8));
var aux_2 = b_2 - x_1 - (0.49 - x_1)*exp(-x_2*(a_2 - 8));
var aux_3 = b_3 - x_1 - (0.49 - x_1)*exp(-x_2*(a_3 - 8));
var aux_4 = b_4 - x_1 - (0.49 - x_1)*exp(-x_2*(a_4 - 8));
var aux_5 = b_5 - x_1 - (0.49 - x_1)*exp(-x_2*(a_5 - 8));
var aux_6 = b_6 - x_1 - (0.49 - x_1)*exp(-x_2*(a_6 - 8));
var aux_7 = b_7 - x_1 - (0.49 - x_1)*exp(-x_2*(a_7 - 8));
var aux_8 = b_8 - x_1 - (0.49 - x_1)*exp(-x_2*(a_8 - 8));
var aux_9 = b_9 - x_1 - (0.49 - x_1)*exp(-x_2*(a_9 - 8));
var aux_10 = b_10 - x_1 - (0.49 - x_1)*exp(-x_2*(a_10 - 8));
var aux_11 = b_11 - x_1 - (0.49 - x_1)*exp(-x_2*(a_11 - 8));
var aux_12 = b_12 - x_1 - (0.49 - x_1)*exp(-x_2*(a_12 - 8));
var aux_13 = b_13 - x_1 - (0.49 - x_1)*exp(-x_2*(a_13 - 8));
var aux_14 = b_14 - x_1 - (0.49 - x_1)*exp(-x_2*(a_14 - 8));
var aux_15 = b_15 - x_1 - (0.49 - x_1)*exp(-x_2*(a_15 - 8));
var aux_16 = b_16 - x_1 - (0.49 - x_1)*exp(-x_2*(a_16 - 8));
var aux_17 = b_17 - x_1 - (0.49 - x_1)*exp(-x_2*(a_17 - 8));
var aux_18 = b_18 - x_1 - (0.49 - x_1)*exp(-x_2*(a_18 - 8));
var aux_19 = b_19 - x_1 - (0.49 - x_1)*exp(-x_2*(a_19 - 8));
var aux_20 = b_20 - x_1 - (0.49 - x_1)*exp(-x_2*(a_20 - 8));
var aux_21 = b_21 - x_1 - (0.49 - x_1)*exp(-x_2*(a_21 - 8));
var aux_22 = b_22 - x_1 - (0.49 - x_1)*exp(-x_2*(a_22 - 8));
var aux_23 = b_23 - x_1 - (0.49 - x_1)*exp(-x_2*(a_23 - 8));
var aux_24 = b_24 - x_1 - (0.49 - x_1)*exp(-x_2*(a_24 - 8));
var aux_25 = b_25 - x_1 - (0.49 - x_1)*exp(-x_2*(a_25 - 8));
var aux_26 = b_26 - x_1 - (0.49 - x_1)*exp(-x_2*(a_26 - 8));
var aux_27 = b_27 - x_1 - (0.49 - x_1)*exp(-x_2*(a_27 - 8));
var aux_28 = b_28 - x_1 - (0.49 - x_1)*exp(-x_2*(a_28 - 8));
var aux_29 = b_29 - x_1 - (0.49 - x_1)*exp(-x_2*(a_29 - 8));
var aux_30 = b_30 - x_1 - (0.49 - x_1)*exp(-x_2*(a_30 - 8));
var aux_31 = b_31 - x_1 - (0.49 - x_1)*exp(-x_2*(a_31 - 8));
var aux_32 = b_32 - x_1 - (0.49 - x_1)*exp(-x_2*(a_32 - 8));
var aux_33 = b_33 - x_1 - (0.49 - x_1)*exp(-x_2*(a_33 - 8));
var aux_34 = b_34 - x_1 - (0.49 - x_1)*exp(-x_2*(a_34 - 8));
var aux_35 = b_35 - x_1 - (0.49 - x_1)*exp(-x_2*(a_35 - 8));
var aux_36 = b_36 - x_1 - (0.49 - x_1)*exp(-x_2*(a_36 - 8));
var aux_37 = b_37 - x_1 - (0.49 - x_1)*exp(-x_2*(a_37 - 8));
var aux_38 = b_38 - x_1 - (0.49 - x_1)*exp(-x_2*(a_38 - 8));
var aux_39 = b_39 - x_1 - (0.49 - x_1)*exp(-x_2*(a_39 - 8));
var aux_40 = b_40 - x_1 - (0.49 - x_1)*exp(-x_2*(a_40 - 8));
var aux_41 = b_41 - x_1 - (0.49 - x_1)*exp(-x_2*(a_41 - 8));
var aux_42 = b_42 - x_1 - (0.49 - x_1)*exp(-x_2*(a_42 - 8));
var aux_43 = b_43 - x_1 - (0.49 - x_1)*exp(-x_2*(a_43 - 8));
var aux_44 = b_44 - x_1 - (0.49 - x_1)*exp(-x_2*(a_44 - 8));
var s_1 = (aux_1)^2;
var s_2 = s_1 + (aux_2)^2;
var s_3 = s_2 + (aux_3)^2;
var s_4 = s_3 + (aux_4)^2;
var s_5 = s_4 + (aux_5)^2;
var s_6 = s_5 + (aux_6)^2;
var s_7 = s_6 + (aux_7)^2;
var s_8 = s_7 + (aux_8)^2;
var s_9 = s_8 + (aux_9)^2;
var s_10 = s_9 + (aux_10)^2;
var s_11 = s_10 + (aux_11)^2;
var s_12 = s_11 + (aux_12)^2;
var s_13 = s_12 + (aux_13)^2;
var s_14 = s_13 + (aux_14)^2;
var s_15 = s_14 + (aux_15)^2;
var s_16 = s_15 + (aux_16)^2;
var s_17 = s_16 + (aux_17)^2;
var s_18 = s_17 + (aux_18)^2;
var s_19 = s_18 + (aux_19)^2;
var s_20 = s_19 + (aux_20)^2;
var s_21 = s_20 + (aux_21)^2;
var s_22 = s_21 + (aux_22)^2;
var s_23 = s_22 + (aux_23)^2;
var s_24 = s_23 + (aux_24)^2;
var s_25 = s_24 + (aux_25)^2;
var s_26 = s_25 + (aux_26)^2;
var s_27 = s_26 + (aux_27)^2;
var s_28 = s_27 + (aux_28)^2;
var s_29 = s_28 + (aux_29)^2;
var s_30 = s_29 + (aux_30)^2;
var s_31 = s_30 + (aux_31)^2;
var s_32 = s_31 + (aux_32)^2;
var s_33 = s_32 + (aux_33)^2;
var s_34 = s_33 + (aux_34)^2;
var s_35 = s_34 + (aux_35)^2;
var s_36 = s_35 + (aux_36)^2;
var s_37 = s_36 + (aux_37)^2;
var s_38 = s_37 + (aux_38)^2;
var s_39 = s_38 + (aux_39)^2;
var s_40 = s_39 + (aux_40)^2;
var s_41 = s_40 + (aux_41)^2;
var s_42 = s_41 + (aux_42)^2;
var s_43 = s_42 + (aux_43)^2;
var s_44 = s_43 + (aux_44)^2;
var mf = s_44;
var myminfun = mf;
var myabsdevnod0_0 = abs(myminfun - 0.02845966972298671);
var myreldevnod0_0 = abs(myminfun/0.02845966972298671 - 1);
var myabsdevnod0_1 = abs(x_1 - 0.4199526507578012);
var myreldevnod0_1 = abs(x_1/0.4199526507578012 - 1);
var myabsdevnod0_2 = abs(x_2 - 1.284845193624845);
var myreldevnod0_2 = abs(x_2/1.284845193624845 - 1);
var myabsdevnod1_0 = (1/2)*(abs(myabsdevnod0_0 - myabsdevnod0_1) + (myabsdevnod0_0 + myabsdevnod0_1));
var myreldevnod1_0 = (1/2)*(abs(myreldevnod0_0 - myreldevnod0_1) + (myreldevnod0_0 + myreldevnod0_1));
var myabsdevnod1_2 = myabsdevnod0_2;
var myreldevnod1_2 = myreldevnod0_2;
var myabsdevnod2_0 = (1/2)*(abs(myabsdevnod1_0 - myabsdevnod1_2) + (myabsdevnod1_0 + myabsdevnod1_2));
var myreldevnod2_0 = (1/2)*(abs(myreldevnod1_0 - myreldevnod1_2) + (myreldevnod1_0 + myreldevnod1_2));
var zmyabsdevmax = myabsdevnod2_0;
var zmyreldevmax = myreldevnod2_0;
minimize obj: myminfun;
subject to ci_1: 0.49*x_2 - x_1*x_2 - 0.09 >= 0;
solve;
display x_1;
display x_2;
display myminfun;
display zmyabsdevmax;
display zmyreldevmax;
# End Model

Stephan K.H. Seidl