<div dir="ltr"><div><div>I'd be interested in helping also and I might be of some assistance on this one. I'm currently working on a thorn that handles particles and uses boost's serialization library to keep them with their local grids. It currently 'works' but is in need of more testing and bug fixes before it'll be ready for science runs. It's also worth nothing that I abandoned using grid array's some time ago (for the reasons Erik mentioned) in favor of my own vector's of particles. <br>
<br></div>If there's a discussion on google hangouts, let me know as I'd be very interested in talking about this.<br><br></div>-Matt<br></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Wed, Jan 15, 2014 at 10:45 AM, Erik Schnetter <span dir="ltr"><<a href="mailto:schnetter@cct.lsu.edu" target="_blank">schnetter@cct.lsu.edu</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="im">On Jan 15, 2014, at 9:05 , <a href="mailto:sruehe@astrophysik.uni-kiel.de">sruehe@astrophysik.uni-kiel.de</a> wrote:<br>
<br>
> Hello EinsteinToolkit-Community,<br>
><br>
> first of all i want thank you very much for the fast and helpful answer to<br>
> my last Email.<br>
><br>
> I mentioned i'm at the beginning of my phd. So before i'm wasting too much<br>
> time on working at something an other already knows, i want to ask some<br>
> perhaps a little bit naively sounding questions.<br>
><br>
> I have studied smoothed particle hydrodynamics on a fixed spacetime, so i<br>
> get the idea devolping a sph based thorn. I found no publications, so is<br>
> perhaps anybody already working on this problem, or has one shown that<br>
> this isn't possible?<br>
> Because this is an community based development, it is quite interesting<br>
> whether it could be an useful thorn.<br>
><br>
> I'm not sure whether my questions are answerable, but perhaps you can help<br>
> me a second time.<br>
<br>
</div>Stefan<br>
<br>
Some years ago I developed an SPH thorn for Newtonian hydrodynamics, and coupled this to a grid-based elliptic solver for self-gravity. This code was able to evolve e.g. single stars. This was a proof of concept; I did not push this into the direction new physics.<br>
<br>
A GR SPH implementation would be easier since it does not need a Newtonian solver. Instead, one would need to define how to set a grid-base T_ab (for spacetime evolution, if desired), and to either evaluate or interpolate the metric at the particle locations. For both, Cactus offers supporting infrastructure.<br>
<br>
The missing parts in Cactus have to do with handling particles. Cactus offers parallelized 1D arrays where particles can be stored. It does not yet offer methods to re-distribute particles between different processes. Of course, these can be added, and I would be very interesting in helping out.<br>
<br>
Would you have time for a Google Hangout in the coming days to chat? I'm located in the EST time zone, six hours behind you, so your late afternoon would be a good time for me.<br>
<span class="HOEnZb"><font color="#888888"><br>
-erik<br>
<br>
--<br>
Erik Schnetter <<a href="mailto:schnetter@cct.lsu.edu">schnetter@cct.lsu.edu</a>><br>
<a href="http://www.perimeterinstitute.ca/personal/eschnetter/" target="_blank">http://www.perimeterinstitute.ca/personal/eschnetter/</a><br>
<br>
My email is as private as my paper mail. I therefore support encrypting<br>
and signing email messages. Get my PGP key from <a href="http://pgp.mit.edu/" target="_blank">http://pgp.mit.edu/</a>.<br>
<br>
</font></span><br>_______________________________________________<br>
Users mailing list<br>
<a href="mailto:Users@einsteintoolkit.org">Users@einsteintoolkit.org</a><br>
<a href="http://lists.einsteintoolkit.org/mailman/listinfo/users" target="_blank">http://lists.einsteintoolkit.org/mailman/listinfo/users</a><br>
<br></blockquote></div><br></div>