<html dir="ltr">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style id="owaParaStyle" type="text/css"></style>
</head>
<body ocsi="0" fpstyle="1">
<div style="direction: ltr;font-family: Tahoma;color: #000000;font-size: 10pt;">Hi Matt,<br>
Many thanks for your reply. I had thought the intracellular walker placement may be an issue, but even setting up simulations with very few walkers doesn't appear to help - none of them progress. I will start by taking a look at the checks you suggest!<br>
<br>
Best wishes,<br>
Damien<br>
<div style="font-family: Times New Roman; color: #000000; font-size: 16px">
<hr tabindex="-1">
<div style="direction: ltr;" id="divRpF963500"><font color="#000000" face="Tahoma" size="2"><b>From:</b> camino-users-bounces@www.nitrc.org [camino-users-bounces@www.nitrc.org] on behalf of Matt Budde [matt.budde@yahoo.com]<br>
<b>Sent:</b> 16 July 2012 19:56<br>
<b>To:</b> Damien McHugh; camino-users@www.nitrc.org<br>
<b>Subject:</b> Re: [Camino-users] ply files in Monte Carlo simulations<br>
</font><br>
</div>
<div></div>
<div>
<div style="color:#000; background-color:#fff; font-family:times new roman,new york,times,serif; font-size:12pt">
<div><span>Hi Damien,</span></div>
<div><span><br>
</span></div>
<div><span>I spent considerable time a couple of years ago running some simulations with ply files. I can give you my insights, but recognize that things could have changed since then.</span></div>
<div><span><br>
</span></div>
<div><span>Regarding intra/extra cellular spaces with mesh substrates: Whether a point within the substrate is intracellular or extracellular is determined by whether a vector passing through the point crosses the mesh an even (intracellular) or odd/zero (extracellular)
number of times. The complication is that if there are odd edges, overlaps, or gaps in your mesh, things can get very complicated. With the '-initial intra' flag, an initial position is selected at random within the entire substrate, and if it happens to
be intracellular according to the vector crossing, it is kept, if it is extracellular, it is discarded and the process continues until all your initial points are determined. If only a small portion of your substrate is intracellular, it can take considerably
longer to randomly select points within this space. Hence, that might be the situation you are seeing. If you start with a much smaller number of walkers, does the intra condition eventually get going? </span></div>
<div><span><br>
</span></div>
<div><span>I also did several checks along the way to ensure that things were working appropriately. The first was to plot the initial positions of the walkers compared to the mesh to ensure that these were in the right locations. (I wrote some matlab routines
to read and plot my ply files and overlay the walker data stored in the trajectory files from camino). Early on, I had issues with units and visualizing the data helped me get things scaled appropriately, including resolving an issue with intra/extra points.
The second check was to plot the full </span><span style="font-size:12pt">random walk </span><span style="font-size:12pt">history for at least a small number of walkers to ensure that they stayed within their compartment for the duration of the simulation
(with permeability set to zero). I did have the problem initially that sometimes a walker would cross over the mesh instead of being reflected by it. The exact solution escapes me now, but I think it had to do with either a rounding or floating point precision
issue. In any case, I had modified and recompiled the code to resolve that issue. For what it's worth, I was using camino version 841.</span></div>
<div><span style="font-size:12pt"><br>
</span></div>
<div><span style="font-size:12pt">Feel free to email me directly if you have any followup questions,</span></div>
<div><span style="font-size:12pt"><br>
</span></div>
<div><span style="font-size:12pt">-matt</span></div>
<div><span style="font-size:12pt"><br>
</span></div>
<div><span style="font-size:12pt"><br>
</span></div>
<div><br>
</div>
<div style="font-family:'times new roman','new york',times,serif; font-size:12pt">
<div style="font-family:'times new roman','new york',times,serif; font-size:12pt">
<div dir="ltr"><font face="Arial" size="2">
<hr size="1">
<b><span style="font-weight:bold">From:</span></b> Damien McHugh <damien.mchugh-2@postgrad.manchester.ac.uk><br>
<b><span style="font-weight:bold">To:</span></b> "camino-users@www.nitrc.org" <camino-users@www.nitrc.org>
<br>
<b><span style="font-weight:bold">Sent:</span></b> Monday, July 16, 2012 9:14 AM<br>
<b><span style="font-weight:bold">Subject:</span></b> [Camino-users] ply files in Monte Carlo simulations<br>
</font></div>
<br>
<div id="yiv90054068"><style type="text/css">
<!--
#yiv90054068 p
{margin-top:0;
margin-bottom:0}
-->
BODY {direction: ltr;font-family: Tahoma;color: #000000;font-size: 10pt;}P {margin-top:0;margin-bottom:0;}BODY {scrollbar-base-color:undefined;scrollbar-highlight-color:undefined;scrollbar-darkshadow-color:undefined;scrollbar-track-color:undefined;scrollbar-arrow-color:undefined}BODY {scrollbar-base-color:undefined;scrollbar-highlight-color:undefined;scrollbar-darkshadow-color:undefined;scrollbar-track-color:undefined;scrollbar-arrow-color:undefined}BODY {scrollbar-base-color:undefined;scrollbar-highlight-color:undefined;scrollbar-darkshadow-color:undefined;scrollbar-track-color:undefined;scrollbar-arrow-color:undefined}</style>
<div>
<div style="direction:ltr; font-family:Tahoma; color:rgb(0,0,0); font-size:10pt">
Hello all,<br>
I am having trouble running simulations of intra-cellular diffusion when using .ply files as substrates. If I set '-initial extra' the simulation will run, with the terminal output showing, for example:<br>
<br>
Jul 10, 2012 2:53:34 PM simulation.DiffusionSimulation <init><br>
INFO: running simulation: 100000 walkers, 1000 timesteps, p= 0.0<br>
Jul 10, 2012 2:53:34 PM simulation.DiffusionSimulation <init><br>
INFO: dynamics duration 0.13<br>
Jul 10, 2012 2:53:34 PM simulation.DiffusionSimulation initialiseWalkers<br>
INFO: writing initial walker positions<br>
Jul 10, 2012 2:53:34 PM simulation.DiffusionSimulation initialiseWalkers<br>
INFO: walkr positions written<br>
0.0% <br>
<br>
However, if I use '-initial intra' I get the following:<br>
<br>
INFO: running simulation: 100000 walkers, 1000 timesteps, p= 0.0<br>
Jul 10, 2012 3:08:52 PM simulation.DiffusionSimulation <init><br>
INFO: dynamics duration 0.13<br>
<br>
The simulation does not progress from here and 0.0% doesn't print to the terminal. Any advice regarding the use of ply files is much appreciated.
<br>
<br>
Many thanks,<br>
Damien<br>
<br>
</div>
</div>
</div>
<br>
_______________________________________________<br>
Camino-users mailing list<br>
<a href="mailto:Camino-users@www.nitrc.org" target="_blank">Camino-users@www.nitrc.org</a><br>
<a href="http://www.nitrc.org/mailman/listinfo/camino-users" target="_blank">http://www.nitrc.org/mailman/listinfo/camino-users</a><br>
<br>
<br>
</div>
</div>
</div>
</div>
</div>
</div>
</body>
</html>