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

Timestamp:

03/27/16 00:12:20 (8 years ago)

Author:

oriona

Message:

Method for definig new species added.

File:

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

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

@@ -28 +28 @@
 }
 
-function createForamGenotype(gen, species, chamber_num)
+function createForamMorphology(morphotype, gen, chamber_num)
 {
         var rad = getProperty(gen, "chamber_proculus"); 
-        var geno = "//0\np:" + chambers[species][0] + "sh=1,sx=" + rad + ",sy=" + rad + ",sz=" + rad + ", rz=3.14159265358979,vr=" + colors[gen];
-
-        chamber_num = Math.min(chamber_num, chambers[species].size - 1);
+        var geno = "//0\np:" + chambers[morphotype][0] + "sh=1,sx=" + rad + ",sy=" + rad + ",sz=" + rad + ", rz=3.14159265358979,vr=" + colors[gen];
+
+        chamber_num = Math.min(chamber_num, chambers[morphotype].size - 1);
 
         for (var i = 0; i < chamber_num; i++)
         {
                 rad = getProperty(gen, "chamber_proculus") + getProperty(gen, "chamber_difference") * (i + 1);
-                geno += "\n" + "p:" + chambers[species][i+1] + "sh=1,sx=" + rad + ",sy=" + rad + ",sz=" + rad + ",vr=" + colors[gen];
+                geno += "\n" + "p:" + chambers[morphotype][i+1] + "sh=1,sx=" + rad + ",sy=" + rad + ",sz=" + rad + ",vr=" + colors[gen];
         }
 
@@ -46 +46 @@
         }
 
-        if (species == 0) geno += "\nn:p=0,d=\"S\"";
+        if (morphotype == 0) geno += "\nn:p=0,d=\"S\"";
 
         return geno;
@@ -55 +55 @@
         if (mode->opt == 0) //initial
         {
-                mode->cr.data->genes = {"min_repro_energies" : [energyFromVolume(max_chamber_volume[0][getProperty(0, "min_repro_energ")],0), energyFromVolume(max_chamber_volume[1][getProperty(1, "min_repro_energ")],0)], "hibernation" : mode->species};
-                mode->cr.data->lifeparams = {"max_energy_level" : getProperty(0,"energies0"), "gen" : 0,  "hibernated" : 0, "species" : mode->species, "reproduce" : 0, "dir" : randomDir(), "chamber_growth" : -1, "division_time" : -1};
+                mode->cr.data->genes = String.deserialize(String.serialize(mode->genes)); 
+                mode->cr.data->lifeparams = {"max_energy_level" : mode->energy0, "gen" : 0,  "hibernated" : 0, "species" : mode->species, "reproduce" : 0, "dir" : randomDir(), "chamber_growth" : -1, "division_time" : -1};
         }
         else if (mode->opt  == 1) //child
         {
-                mode->cr.data->lifeparams = {"max_energy_level" : getProperty(1 - mode->parent_lifeparams->gen,"energies0"), "gen" : 1 - mode->parent_lifeparams->gen,  "hibernated" : 0, "species" : mode->parent_lifeparams->species, "reproduce" : 0, "dir" : randomDir(), "chamber_growth" : -1, "division_time" : -1};
-                mode->cr.data->genes = mode->parent_genes;
+                mode->cr.data->lifeparams = {"max_energy_level" : mode->energy0, "gen" : 1 - mode->parent_lifeparams->gen,  "hibernated" : 0, "species" : mode->parent_lifeparams->species, "reproduce" : 0, "dir" : randomDir(), "chamber_growth" : -1, "division_time" : -1};
+                mode->cr.data->genes = String.deserialize(String.serialize(mode->parent_genes));
         }
         else //grow
@@ -81 +81 @@
         //Simulator.print("parent: " + parent_energy + " result: " + result + " number " + number);
         return {"energy" : result, "number" : number};
+}
+
+function getEnergy0(radius)
+{
+        return energyFromVolume(micronsToFrams(radius),1);
 }
 
@@ -88 +93 @@
         if (clone == 1)
         {
-                var offspring = gametsDivision(parent.energy,getProperty(0,"energies0"));
+                var offspring = gametsDivision(parent.energy,getEnergy0(getGene(parent,"energies0",0)[0])); 
                 energy0 = offspring->energy;
                 number = offspring->number;
@@ -97 +102 @@
         else
         {
-                var offspring1 = gametsDivision(parent.energy,getProperty(1,"energies0"));
-                var offspring2 = gametsDivision(parent2.energy,getProperty(1,"energies0"));
+                var offspring1 = gametsDivision(parent.energy,getEnergy0(getGene(parent,"energies0", 0)[1])); 
+                var offspring2 = gametsDivision(parent2.energy,getEnergy0(getGene(parent2,"energies0", 0)[1])); 
                 energy0 = (offspring1->energy+offspring2->energy); 
                 number = ExpProperties.gametSuccessRate*(offspring1->number+offspring2->number)/2;
@@ -109 +114 @@
         for (var j = 0; j < number; j++)
         {
-                createOffspring(createForamGenotype(gen, parent.data->lifeparams->species, 0), energy0, new_genes, parent.data->lifeparams); 
+                createOffspring(createForamMorphology(getGene(parent, "morphotype", 0), gen, 0), energy0, new_genes, parent.data->lifeparams); 
         }
 }
@@ -115 +120 @@
 function reproduce_diploid(parent)
 {
-        var energy0 = getProperty(0,"energies0");
+        var energy0 =getEnergy0( getGene(parent,"energies0", 0)[0]); 
         var number = ((1 - (getProperty(parent.data->lifeparams->gen, "e_repro_cost"))) * parent.energy) / energy0;
 
@@ -132 +137 @@
                 for (var k = 0; k < 2; k++)
                 {
-                        createOffspring(createForamGenotype(1 - parent.data->lifeparams->gen, parent.data->lifeparams->species, 0), energy0, parent.data->genes[0], parent.data->lifeparams); 
+                        createOffspring(createForamMorphology(getGene(parent, "morphotype", 0), 1 - parent.data->lifeparams->gen, 0), energy0, parent.data->genes[0], parent.data->lifeparams); 
                 }
 
@@ -236 +241 @@
 }
 
-
-
 function foramReproduce(cr)
 {
-        var properEnergy = 0;
-        var reproduced = 0;
-
-        if (cr.data->lifeparams->gen == 0)
-        {
-                properEnergy = ( cr.energy >= cr.data->genes->min_repro_energies[cr.data->lifeparams->gen] );
-        }
-        else
-        {
-                properEnergy = ( cr.energy >= cr.data->genes[0]->min_repro_energies[cr.data->lifeparams->gen] ); //TODO gene selection
-        }
+        var properEnergy = cr.energy >= energyFromVolume(max_chamber_volume[0][getGene(cr, "min_repro_energies",0)[cr.data->lifeparams->gen]],0);
+        var reproduced = 0;     
 
         //if creature has proper energy
@@ -289 +283 @@
         cr.energy0 = energy;
         cr.energy = cr.energy0;
-        setGenotype({"cr" : cr, "parent_genes" : parent_genes, "parent_lifeparams" : parent_lifeparams, "opt" : 1});
+        setGenotype({"cr" : cr, "parent_genes" : parent_genes, "parent_lifeparams" : parent_lifeparams, "opt" : 1, "energy0" : cr.energy0});
         placeRandomlyNotColliding(cr);
 }