// Use this for initialization
void Start()
{
timeStream = TimeStream.Instance;
timer = 0;
maxTimer = 4;
credit = false;
}
// Use this for initialization
void Start()
{
timeStream = TimeStream.Instance;
date = timeStream.getTime();
audioData = GetComponent <AudioSource>();
}
// Use this for initialization
void Start()
{
timeStream = TimeStream.Instance;
startDate = timeStream.getTime();
age = initAge;
nbCoups = 0;
spriteRenderer = GetComponent <SpriteRenderer>();
if (spriteRenderer.sprite == null)
{
if (age == 0)
{
spriteRenderer.sprite = tree0;
}
else if (age == 1)
{
spriteRenderer.sprite = tree1;
}
else if (age == 2)
{
spriteRenderer.sprite = tree2;
}
else
{
spriteRenderer.sprite = tree3;
}
}
}
private void SafeRelease()
{
if (TimeStream != null)
{
TimeStream.Dispose();
}
TimeStream = null;
}
private void Start()
{
timeStream = TimeStream.Instance;
canGetWoodLog = true;
bonusState = false;
nbJumpsGift = 4;
stop = false;
}
// Use this for initialization
void Start()
{
timeStream = TimeStream.Instance;
isActivated = false;
spriteRenderer = GetComponent <SpriteRenderer>();
if (spriteRenderer.sprite == null)
{
spriteRenderer.sprite = axeSprite1;
}
}
// Use this for initialization
void Start()
{
timeStream = TimeStream.Instance;
spriteRenderer = GetComponent <SpriteRenderer>();
if (spriteRenderer.sprite == null)
{
spriteRenderer.sprite = riverSprite;
}
}
// Use this for initialization
void Start()
{
timeStream = TimeStream.Instance;
nbWood = 0;
spriteRenderer = GetComponent <SpriteRenderer>();
if (spriteRenderer.sprite == null)
{
spriteRenderer.sprite = atelier0;
}
}
public void Open(String cellFile)
{
if (_context.TimeStream != null)
{
Close();
}
_context.TimeStream = TimeStream.Open(cellFile);
_context.TimeStream.MoveTo(0);
_context.CurrentLayer = 0;
_context.CurrentTime = 0.0;
_context.Timer.Stop();
_context.Timer.Reset();
}
// Use this for initialization
void Start()
{
timeStream = TimeStream.Instance;
audioPause = false;
}
// Use this for initialization
void Start()
{
timeStream = TimeStream.Instance;
}
private void Awake()
{
instance = gameObject.GetComponent <TimeStream>();
}
private static bool movePoles = true; // move poles from top to sides
private static void Simulate(String testID, SimParams config, Random rnd)
{
// Prepare directory
if (Directory.Exists(testID))
{
Directory.Delete(testID, true);
}
Directory.CreateDirectory(testID);
Directory.SetCurrentDirectory(testID);
// Cache some values that we are goting to use for configuration
double R_Cell = config[SimParameter.Double.R_Cell];
double L_Poles = config[SimParameter.Double.L_Poles];
int N_MT_Total = config[SimParameter.Int.N_MT_Total];
double Dt = config[SimParameter.Double.Dt];
double Save_Freq_Macro = config[SimParameter.Double.Save_Freq_Macro];
double T_End = config[SimParameter.Double.T_End];
int layerCount = (int)(T_End / Save_Freq_Macro) + 1;
// Ready to start
int[] res_mts = new int[layerCount * 4 * cellsPerCase];
double[] res_angles = new double[layerCount * 2 * cellsPerCase];
Stopwatch sw = new Stopwatch();
int batchSize = Math.Max(1, Environment.ProcessorCount);
int simN = 0;
while (simN < cellsPerCase)
{
var parameters = new LaunchParameters();
var cellCount = Math.Min(batchSize, cellsPerCase - simN);
sw.Reset();
// Say greetings
if (cellCount == 1)
{
Console.WriteLine("Cell #" + (simN + 1).ToString());
}
else
{
Console.WriteLine(String.Format("Cells #{0}...#{1}", simN + 1, simN + cellCount));
}
Console.Write(" Simulating ... ");
// Update config in accordance with our flags
sw.Start();
parameters.Config = config;
config[SimParameter.Double.D_Trans] = hasDiff ? SimParams.GetDefaultValue(SimParameter.Double.D_Trans, false) : 0.0;
config[SimParameter.Double.D_Rot] = hasDiff ? SimParams.GetDefaultValue(SimParameter.Double.D_Rot, false) : 0.0;
// Set initial state
parameters.InitialStates = new InitialStates();
parameters.InitialStates.AddChromosomePair(0.0, 0.0, 1e-6, Math.PI / 2, 0.0, 0.0);
for (int i = 0; i < N_MT_Total * 2; i++)
{
double alpha = (rnd.NextDouble()) * Math.PI * 2;
double dx = rnd.NextDouble() * (i < N_MT_Total ? 1 : -1);
double dy = Math.Sqrt(1.0 - dx * dx) * Math.Cos(alpha);
double dz = Math.Sqrt(1.0 - dx * dx) * Math.Sin(alpha);
parameters.InitialStates.AddMT(i < N_MT_Total ? PoleType.Left : PoleType.Right,
dx, dy, dz,
0.0,
MTState.Polymerization);
}
// Set poles
if (movePoles)
{
parameters.PoleCoords = new PoleCoordinates();
double t = 180.0;
double pole_r = L_Poles / 2 * 1e-6;
for (int i = 0; i < t; i++)
{
double angle = i / (t - 1) * Math.PI / 2;
parameters.PoleCoords.AddRecord((double)i,
-pole_r * Math.Sin(angle), pole_r * Math.Cos(angle), 0.0,
pole_r * Math.Sin(angle), pole_r * Math.Cos(angle), 0.0);
}
}
// Simulate
parameters.Args.Mode = LaunchMode.New;
parameters.Args.Solver = new SimulatorConfig(SimulatorType.CPU);
parameters.Args.CellCount = cellCount;
parameters.Args.CellFile = "results.cell";
using (Launcher ml = new Launcher(".", "../Mitosis.exe", parameters))
{
var res = ml.StartAndWait();
if (res.ExitCode != 0 || !String.IsNullOrEmpty(res.Error))
{
throw new ApplicationException(String.Format("Simulation failed, error #{0}: {1}, {2}",
res.ExitCode, res.Error, res.Output));
}
}
sw.Stop();
Console.WriteLine(String.Format("{0} seconds", (int)sw.ElapsedMilliseconds / 1000));
// Analyze
Console.Write(" Counting ... ");
sw.Reset();
sw.Start();
for (int cellN = simN; cellN < simN + cellCount; cellN++)
{
var ts = TimeStream.Open(cellCount == 1
? parameters.Args.CellFile
: parameters.Args.CellFile.Replace("results.cell",
String.Format("results_{0}.cell",
cellN - simN)));
if (ts.LayerCount != layerCount)
{
throw new ApplicationException("Unexpected count of time layers");
}
ts.Reset();
int layerN = 0;
while (ts.MoveNext())
{
var cell = ts.Current.Cell;
var chr0 = cell.Chromosomes.ElementAt(0);
var chr1 = cell.Chromosomes.ElementAt(1);
// Count MTs
int c0l = 0, c0r = 0, c1l = 0, c1r = 0;
var c0 = chr0.BoundMTs;
foreach (var mt in c0)
{
if (mt.Pole.Type == PoleType.Left)
{
c0l += 1;
}
else
{
c0r += 1;
}
}
var c1 = chr1.BoundMTs;
foreach (var mt in c1)
{
if (mt.Pole.Type == PoleType.Left)
{
c1l += 1;
}
else
{
c1r += 1;
}
}
int offset = layerCount * 4 * cellN + layerN * 4;
res_mts[offset + 0] = c0l;
res_mts[offset + 1] = c0r;
res_mts[offset + 2] = c1l;
res_mts[offset + 3] = c1r;
// Count angles
if (cell.AreSpringsBroken)
{
throw new ApplicationException("Wrong setup, springs were broken!");
}
Vector3 dp = new Vector3(cell.GetPole(PoleType.Left).Position.X - cell.GetPole(PoleType.Right).Position.X,
cell.GetPole(PoleType.Left).Position.Y - cell.GetPole(PoleType.Right).Position.Y,
cell.GetPole(PoleType.Left).Position.Z - cell.GetPole(PoleType.Right).Position.Z);
if (dp.Length() == 0.0f)
{
dp = new Vector3(-1.0f, 0.0f, 0.0f);
}
else
{
dp = dp / dp.Length();
}
if (Math.Abs(dp.Y) + Math.Abs(dp.Z) > 1e-9)
{
throw new ApplicationException("Poles must be located symmetrically");
}
Vector3 dc_x = new Vector3((float)-chr0.Orientation[0, 0],
(float)-chr0.Orientation[1, 0],
(float)-chr0.Orientation[2, 0]);
if (dc_x.Length() == 0.0f)
{
dc_x = new Vector3(-1.0f, 0.0f, 0.0f);
}
else
{
dc_x = dc_x / dc_x.Length();
}
Vector3 dc_y = new Vector3((float)-chr0.Orientation[0, 1],
(float)-chr0.Orientation[1, 1],
(float)-chr0.Orientation[2, 1]);
if (dc_y.Length() == 0.0f)
{
dc_y = new Vector3(0.0f, 1.0f, 0.0f);
}
else
{
dc_y = dc_y / dc_y.Length();
}
offset = (layerCount * cellN + layerN) * 2;
res_angles[offset + 0] = Math.Acos(Math.Min(1.0, Math.Max(-1.0, Vector3.Dot(dp, dc_x))));
res_angles[offset + 1] = Math.Acos(Math.Min(1.0, Math.Max(-1.0, Vector3.Dot(dp, dc_y))));
layerN += 1;
}
ts.Dispose();
} // for cellN
sw.Stop();
Console.WriteLine(String.Format("{0} seconds", sw.ElapsedMilliseconds / 1000));
simN += cellCount;
} // while simN < cellsPerCase
// Save results
Console.WriteLine("Saving ...");
// A-results
using (var csv = File.CreateText("results_A.csv"))
{
csv.WriteLine("Connected MTs (pcs)");
csv.Write("Time (seconds),Averaged,Min,Max");
for (int i = 0; i < cellsPerCase; i++)
{
csv.Write(",Cell #" + (i + 1).ToString());
}
csv.WriteLine();
for (int layerN = 0; layerN < layerCount; layerN++)
{
StringBuilder str = new StringBuilder();
int min = Int32.MaxValue, max = Int32.MinValue, avg = 0;
for (int cellN = 0; cellN < cellsPerCase; cellN++)
{
int offset = layerN * 4 + layerCount * 4 * cellN;
int val = res_mts[offset] + res_mts[offset + 1] + res_mts[offset + 2] + res_mts[offset + 3];
str.Append(",");
str.Append(val.ToString());
min = Math.Min(min, val);
max = Math.Max(max, val);
avg += val;
}
csv.WriteLine(String.Format("{0},{1},{2},{3}{4}",
(Save_Freq_Macro * layerN).ToString(CultureInfo.InvariantCulture),
((double)avg / cellsPerCase).ToString(CultureInfo.InvariantCulture),
min, max, str.ToString()));
}
}
// B-results
using (var csv = File.CreateText("results_B.csv"))
{
csv.WriteLine("Kinetochore types (%)");
csv.WriteLine("Time (seconds),No MTs,Monotelic,Syntelic,Bioriented,Merotelic");
for (int layerN = 0; layerN < layerCount; layerN++)
{
int noMTs = 0, monotelic = 0, syntelic = 0, bioriented = 0, merotelic = 0;
for (int cellN = 0; cellN < cellsPerCase; cellN++)
{
int offset = layerN * 4 + layerCount * 4 * cellN;
int c0l = res_mts[offset + 0];
int c0r = res_mts[offset + 1];
int c1l = res_mts[offset + 2];
int c1r = res_mts[offset + 3];
int c0m = Math.Min(c0l, c0r);
int c0nm = Math.Max(c0l, c0r);
int c1m = Math.Min(c1l, c1r);
int c1nm = Math.Max(c1l, c1r);
int c0_sum = c0m + c0nm;
int c1_sum = c1m + c1nm;
int sum = c0_sum + c1_sum;
if (sum == 0)
{
noMTs += 1;
}
else if (c0m != 0 || c1m != 0)
{
merotelic += 1;
}
else if (c0_sum == 0 || c1_sum == 0)
{
monotelic += 1;
}
else if (Math.Min(c0l, c1l) != 0 || Math.Min(c0r, c1r) != 0)
{
syntelic += 1;
}
else if (c0m == 0 && c0nm != 0 && c1m == 0 && c1nm != 0)
{
bioriented += 1;
}
else
{
throw new ApplicationException("Internal error, unknown type of chromosome");
}
}
csv.WriteLine(String.Format("{0},{1},{2},{3},{4},{5}",
(Save_Freq_Macro * layerN).ToString(CultureInfo.InvariantCulture),
((double)noMTs / cellsPerCase).ToString(CultureInfo.InvariantCulture),
((double)monotelic / cellsPerCase).ToString(CultureInfo.InvariantCulture),
((double)syntelic / cellsPerCase).ToString(CultureInfo.InvariantCulture),
((double)bioriented / cellsPerCase).ToString(CultureInfo.InvariantCulture),
((double)merotelic / cellsPerCase).ToString(CultureInfo.InvariantCulture)));
}
}
// C-results
using (var csv = File.CreateText("results_C.csv"))
{
csv.WriteLine("Merotelic MTs (pcs)");
csv.Write("Time (seconds),Averaged,Min,Max");
for (int i = 0; i < cellsPerCase; i++)
{
csv.Write(",Cell #" + (i + 1).ToString());
}
csv.WriteLine();
for (int layerN = 0; layerN < layerCount; layerN++)
{
StringBuilder str = new StringBuilder();
int min = Int32.MaxValue, max = Int32.MinValue, avg = 0;
for (int cellN = 0; cellN < cellsPerCase; cellN++)
{
int offset = layerN * 4 + layerCount * 4 * cellN;
int val = Math.Min(res_mts[offset + 0], res_mts[offset + 1]) + Math.Min(res_mts[offset + 2], res_mts[offset + 3]);
str.Append(",");
str.Append(val.ToString());
min = Math.Min(min, val);
max = Math.Max(max, val);
avg += val;
}
csv.WriteLine(String.Format("{0},{1},{2},{3}{4}",
(Save_Freq_Macro * layerN).ToString(CultureInfo.InvariantCulture),
((double)avg / cellsPerCase).ToString(CultureInfo.InvariantCulture),
min, max, str.ToString()));
}
}
// D-results
using (var csv_x = File.CreateText("results_D_x.csv"))
using (var csv_y = File.CreateText("results_D_y.csv"))
{
var csvs = new StreamWriter[] { csv_x, csv_y };
for (int l = 0; l < csvs.Length; l++)
{
var csv = csvs[l];
csv.WriteLine("Angle, degrees");
csv.Write("Time (seconds),Averaged,Min,Max");
for (int i = 0; i < cellsPerCase; i++)
{
csv.Write(",Cell #" + (i + 1).ToString());
}
csv.WriteLine();
for (int layerN = 0; layerN < layerCount; layerN++)
{
StringBuilder str = new StringBuilder();
double min = Double.MaxValue, max = Double.MinValue, avg = 0.0;
for (int cellN = 0; cellN < cellsPerCase; cellN++)
{
double angle = res_angles[l + (layerN + layerCount * cellN) * 2];
if (angle > Math.PI)
{
angle = 2 * Math.PI - angle;
}
if (angle > Math.PI / 2)
{
angle = Math.PI - angle;
}
angle *= 180.0 / Math.PI;
str.Append(",");
str.Append(String.Format(CultureInfo.InvariantCulture, "{0:F1}", angle));
min = Math.Min(min, angle);
max = Math.Max(max, angle);
avg += angle;
}
csv.WriteLine(String.Format(CultureInfo.InvariantCulture,
"{0:F1},{1:F1},{2:F1},{3:F1}{4}",
Save_Freq_Macro * layerN,
(double)avg / cellsPerCase,
min, max, str.ToString()));
}
}
}
Directory.SetCurrentDirectory("..");
}
void Start()
{
axis = new Vector3(0, 0, 0);
timeStream = TimeStream.Instance;
}
private double GetEndTime()
{
return(TimeStream == null ? 0.0 : TimeStream.GetLayerTime(TimeStream.LayerCount - 1));
}
