On 13/02/2013, at 12:57 PM, john skaller wrote:
> I have some more performance stuff coming.
>
> Now interestingly, I think I know why Felix slaughters C and Ocaml
> in nbody.
Bzzt .. wrong. Felix wasn't writing anything to memory.
Hmm.
> bodies.j.vel = bodies.j.vel + d * bimass * dt / (distance * d2);
This doesn't work. This
&bodies.j.val <- bodies.j.vel * ....
makes Felix give the same answer as Ocaml and C ..
and also take 5 min 30 secs (vs C/Ocaml 20s) to finish,
and Jupiter still doesn't move at all..
--
john skaller
[email protected]
http://felix-lang.org
------------------------------------------------------------------------------
Free Next-Gen Firewall Hardware Offer
Buy your Sophos next-gen firewall before the end March 2013
and get the hardware for free! Learn more.
http://p.sf.net/sfu/sophos-d2d-feb
_______________________________________________
Felix-language mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/felix-language
