Changeset 421 for experiments/frams/foraminifera/data/scripts/foraminifera.inc

Timestamp:

08/03/15 12:06:54 (9 years ago)

Author:

oriona

Message:

refactoring, removed unused params, creatures and food chanfe to forams and nutirents

File:

• experiments/frams/foraminifera/data/scripts/foraminifera.inc (modified) (7 diffs)

experiments/frams/foraminifera/data/scripts/foraminifera.inc

@@ -6 +6 @@
 function reproduce_haploid(parent, parent2)
 {
-        var number = parent.user2["Va"] / ExpParams.ofnumd + parent2.user2["Va"] / ExpParams.ofnumd;
+        var number = parent.user2["Va"] / ExpParams.energy0d + parent2.user2["Va"] / ExpParams.energy0d;
 
         for (var j = 0; j < number; j++)
         {
-                createOffspring(ExpParams.gend, ExpParams.ofnumd, [parent.user1, parent2.user1], {"Va" : ExpParams.ofnumd, "gen" : 1 , "growth_step" : ExpParams.growth_step, "rsize": ExpParams.rads[1], "vinit": ExpParams.ofnumd}); //TODO genes from both generations in user1
+                createOffspring(ExpParams.gend, ExpParams.energy0d, [parent.user1, parent2.user1], {"Va" : ExpParams.energy0d, "gen" : 1 , "growth_step" : ExpParams.growth_step, "rsize": ExpParams.rads[1], "vinit": ExpParams.energy0d}); //TODO genes from both generations in user1
         }
 }
@@ -16 +16 @@
 function reproduce_diploid(parent)
 {
-        var number = parent.user2["Va"] / ExpParams.ofnumh;
+        var number = parent.user2["Va"] / ExpParams.energy0h;
         for (var j = 0; j < number / 2; j++)
         {
@@ -23 +23 @@
                 if (Math.rnd01 < ExpParams.crossprob)
                 {
-                        crossover(parent, "vamin");
+                        crossover(parent, "minenergy");
                         crossed = 1;
                 }
@@ -29 +29 @@
                 for (var k = 0; k < 2; k++)
                 {
-                        createOffspring(ExpParams.genh, ExpParams.ofnumh, {"vamin" : ExpParams.v_min_h, "amin": ExpParams.age_min_h}, {"Va" : ExpParams.ofnumh, "gen" : 0 , "growth_step" : ExpParams.growth_step, "rsize" : ExpParams.rads[0], "vinit": ExpParams.ofnumh});
+                        createOffspring(ExpParams.genh, ExpParams.energy0h, {"minenergy" : ExpParams.minenerg[0], "minage": ExpParams.minage[0]}, {"Va" : ExpParams.energy0h, "gen" : 0 , "growth_step" : ExpParams.growth_step, "rsize" : ExpParams.rads[0], "vinit": ExpParams.energy0h}); //TODO different vamin and amin for haploids and diploids?
                 }
 
@@ -35 +35 @@
                 if (crossed == 1)
                 {
-                        crossover(parent, "vamin");
+                        crossover(parent, "minenergy");
                         crossed = 0;
                 }
@@ -61 +61 @@
 function onStep()
 {
-        var haploids = 0;
-        var diploids = 0;
-        var e_inc_h = 0.0;
-        var e_inc_d = 0.0;
-        var e_nut = 0.0;
-        var size_h = 0.0;
-        var size_d = 0.0;
-        var vmin_h = 0.0;
-        var vmin_d = 0.0;
-
-        for (var i = 0; i < Populations[0].size; i++)
-        {
-                var cr = Populations[0].get(i);
-                if (cr.user2["gen"] == 0)
-                {
-                        haploids += 1;
-                        e_inc_h += cr.energy;
-                        size_h += ExpParams.rads[0]; //TODO change of size dependent on gene
-                        vmin_h += cr.user1["vamin"];
-                }
-                else if (cr.user2["gen"] == 1)
-                {
-                        diploids += 1;
-                        e_inc_d += cr.energy;
-                        size_d += ExpParams.rads[1]; //TODO change of size dependent on gene
-                        vmin_d += cr.user1[0]["vamin"];
-                }
-        }
-
-        for (var i = 0; i < Populations[1].size; i++)
-        {
-                var cr = Populations[1].get(i);
-                e_nut += cr.energy;
-        }
-
-        if (haploids < 2 && diploids == 0)
-        {
-                Simulator.print("no more creatures, stopped.");
-                Simulator.stop();
-        }
-
-        var l1 = [haploids, diploids, Populations[1].size];
-        var l2 = [e_inc_h, e_inc_d, e_nut];
-        var dp = diploids;
-        var hp = haploids;
-        if (dp == 0) dp = 0.000001;
-        if (hp == 0) hp = 0.000001;
-        var l3 = [size_h / hp, vmin_h / hp, size_d / dp, vmin_d / dp];
 
         if (ExpParams.logging == 1)
         {
-                log(l1, "log.txt");
-                log(l2, "log2.txt");
-                log(l3, "log3.txt");
+                createStatistics();
         }
-        //food growth ---------------------------------------------
-        foodenergywaiting = foodenergywaiting + ExpParams.feedrate;
-        if (foodenergywaiting > ExpParams.feede0)
+
+        //nutrient growth ---------------------------------------------
+        nutrientenergywaiting = nutrientenergywaiting + ExpParams.feedrate;
+        if (nutrientenergywaiting > ExpParams.feede0)
         {
-                for (var i = 0; i < ExpParams.foodPop; i++)
+                for (var i = 0; i < ExpParams.nutrientPop; i++)
                 {
-                        addfood();
+                        addNutrient();
                 }
 
-                foodenergywaiting = 0.0;
+                nutrientenergywaiting = 0.0;
                 Simulator.checkpoint();
         }
@@ -183 +134 @@
 }
 
+function createStatistics()
+{
+        var haploids = 0;
+        var diploids = 0;
+        var e_inc_h = 0.0;
+        var e_inc_d = 0.0;
+        var e_nut = 0.0;
+        var size_h = 0.0;
+        var size_d = 0.0;
+        var vmin_h = 0.0;
+        var vmin_d = 0.0;
+
+        for (var i = 0; i < Populations[0].size; i++)
+        {
+                var cr = Populations[0].get(i);
+                if (cr.user2["gen"] == 0)
+                {
+                        haploids += 1;
+                        e_inc_h += cr.energy;
+                        size_h += ExpParams.rads[0]; //TODO change of size dependent on gene
+                        vmin_h += cr.user1["minenergy"];
+                }
+                else if (cr.user2["gen"] == 1)
+                {
+                        diploids += 1;
+                        e_inc_d += cr.energy;
+                        size_d += ExpParams.rads[1]; //TODO change of size dependent on gene
+                        vmin_d += cr.user1[0]["minenergy"];
+                }
+        }
+
+        for (var i = 0; i < Populations[1].size; i++)
+        {
+                var cr = Populations[1].get(i);
+                e_nut += cr.energy;
+        }
+
+        var l1 = [haploids, diploids, Populations[1].size];
+        var l2 = [e_inc_h, e_inc_d, e_nut];
+        var dp = diploids;
+        var hp = haploids;
+        if (dp == 0) dp = 0.000001;
+        if (hp == 0) hp = 0.000001;
+        var l3 = [size_h / hp, vmin_h / hp, size_d / dp, vmin_d / dp];
+
+        log(l1, "log_sizes.txt");
+        log(l2, "log_energy.txt");
+        log(l3, "log_sizes_v.txt");
+}
 // -------------------------------- step end --------------------------------