reconstructing units

Rezept.mod

@@ -1,11 +1,11 @@
 set ZUTATEN;
 set STOFFE;
 set INDEX; #Indexmenge
-param preis {z in ZUTATEN} >= 0;
-param werte {z in ZUTATEN, s in STOFFE} >= 0;
-param minwert {s in STOFFE} >= 0 default 0;
-param maxwert {s in STOFFE} >= 0 default Infinity;
-var menge {z in ZUTATEN} >= 0;
+param preis {z in ZUTATEN} >= 0; #in €/kg oder €/10Stück
+param werte {z in ZUTATEN, s in STOFFE} >= 0; #je 100g oder Stück
+param minwert {s in STOFFE} >= 0 default 0;   #ebenso
+param maxwert {s in STOFFE} >= 0 default Infinity; #ebenso
+var menge {z in ZUTATEN} >= 0; # in g oder Stück
 param ref {z in ZUTATEN} >= 0 default 0;
 
 param bezugsmenge {z in ZUTATEN} default 100; #nicht immer ist menge in 100g, teilweise in Stück
@@ -13,15 +13,15 @@ param weight_penalty >= 0 <= 1 default 0.05;       #Transportkosten je 100g
 param pill_penalty {z in ZUTATEN, i in INDEX} default 0; #für einzelne Zutaten können Strafkosten bei Überschreiten einer Höchstmenge verhängt werden
 #Die Höchstmenge wird in pill_penalty[z,1], die Strafkosten für jede darüberliegende Einheit in pill_penalty[z,2] gespeichert.
 
-var kaufpreis = sum {z in ZUTATEN} preis[z]/10*menge[z]; #Alnatura-Preise sind in 1kg statt 100g…
+var kaufpreis = sum {z in ZUTATEN} preis[z]/(10*bezugsmenge[z])*menge[z]; #Alnatura-Preise sind in €/1kg statt 100g…
 var omega_ratio = sum{z in ZUTATEN} menge[z]*werte[z,"Omega-6"]/sum{y in ZUTATEN} menge[y]*werte[z,"Omega-3"];
-var rezeptwert {s in STOFFE} = sum {z in ZUTATEN} menge[z]*werte[z,s];
-var total_weight = sum {z in ZUTATEN} menge[z]*bezugsmenge[z];
+var rezeptwert {s in STOFFE} = sum {z in ZUTATEN} menge[z]*werte[z,s]/bezugsmenge[z];
+var total_weight = sum {z in ZUTATEN} menge[z];
 var total_pill_penalty = sum {z in ZUTATEN} max(menge[z]-pill_penalty[z,1],0)*pill_penalty[z,2];
 #Wenn die zutatenweise Obergrenze in pill_penalty[z,1] überschritten ist, werden Strafkosten für die Überschussmenge berechnet.
 
 minimize gesamtpreis: kaufpreis*(1-weight_penalty) + total_weight/100*weight_penalty + total_pill_penalty; 
 subject to ugrenze {s in STOFFE}: rezeptwert[s] >= minwert[s]; #undefined renders 0 => no problem here
 subject to ogrenze {s in STOFFE}: rezeptwert[s] <= maxwert[s]; #undefined should render Infinity here.
-subject to omega_ratio_constraint: omega_ratio <= 6;
+#subject to omega_ratio_constraint: omega_ratio <= 6;
 

config.dat

@@ -1,4 +1,4 @@
-param weight_penalty 0.5455;       #Transportkosten je 100g
+param weight_penalty 0;       #Transportkosten je 100g
 #should be below 1. If equal 1, prices will not be taken into account.
 
 #Um einzelne Zutaten nach Stück statt 100g zu berechnen:

grenzen-orig.dat

@@ -7,11 +7,13 @@ Omega-3 1.67
 Omega-6 17
 Fiber 28
 Cholesterol 0
-#VitaminA 3000
+#VitaminA 3000 #orig min
+VitaminA 2132 #orig Recipe
 VitaminB6 1.3
 VitaminB12 2.4 #nnot hold
 VitaminC  90 #cannot hold
 #VitaminD 600 #cannot hold
+VitaminD 300 #orig Recipe
 VitaminE 20.00
 VitaminK 120
 Thiamin 1.2
@@ -20,7 +22,9 @@ Niacin 16
 Folate 400
 PantothenicAcid 5
 #Choline 550  #not in swiss data
-Calcium 1
+Choline 255 #orig Recipe
+#Calcium 1
+Calcium 0.61 #orig Recipe
 Chloride 2.3
 Chromium 35 #not in swiss data
 Copper 0.90 #not in swiss data
@@ -31,6 +35,7 @@ Manganese 2.3 #not in swiss data
 Molybdenum 45 #not in swiss data
 Phosphorus 0.7
 #Potassium 3.5
+Potassium 3.4 #orig Recipe
 Selenium 55 #cannot hold
 Sodium 1.5
 Sulfur 2  #not in swiss data

main.ampl

@@ -3,4 +3,4 @@ model Rezept.mod;
 data Rezept.dat;
 option solver "ampl.linux64/minos";
 solve;
-display {i in ZUTATEN} menge[i]*bezugsmenge[i], ref, total_weight, kaufpreis ;
+display {i in ZUTATEN} menge[i], ref, total_weight, kaufpreis ;

preise.dat

@@ -1,7 +1,8 @@
+#alle Preise in €/(bezugsmenge*10)
 feingemahleneHaferflockenOR  00.78
 SojaflockenOR  04.3
 JodsalzOR  00.22
-MultivitamineOR  24.1258741258741
+MultivitamineOR  0.241258741258741
 SonnenblumenölOR  03.38
 PuderzuckerOR  00.65
 ErbsenproteinOR  09.676

zutaten.dat

@@ -5,7 +5,7 @@ MultivitamineOR
 SonnenblumenölOR
 PuderzuckerOR
 ErbsenproteinOR
-##WeizenmehlOR
+WeizenmehlOR
 ##OlivenölOR
 ##MaltodextrinOR
 ##Sojaflocken