source code
// The Computer Language Benchmarks Game
// https://salsa.debian.org/benchmarksgame-team/benchmarksgame/
//
// ported from C# version with local variable changes by Anthony Lloyd
#nowarn "9"
[<Literal>]
let N = 5
[<Literal>]
let DT = 0.01
[<Literal>]
let DaysPeryear = 365.24
[<Literal>]
let Solarmass = 39.4784176043574
[<Literal>]
let Jx = 4.84143144246472090e+00
[<Literal>]
let Jy = -1.16032004402742839e+00
[<Literal>]
let Jz = -1.03622044471123109e-01
[<Literal>]
let Jvx = 1.66007664274403694e-03
[<Literal>]
let Jvy = 7.69901118419740425e-03
[<Literal>]
let Jvz = -6.90460016972063023e-05
[<Literal>]
let Jmass = 9.54791938424326609e-04
[<Literal>]
let Sx = 8.34336671824457987e+00
[<Literal>]
let Sy = 4.12479856412430479e+00
[<Literal>]
let Sz = -4.03523417114321381e-01
[<Literal>]
let Svx = -2.76742510726862411e-03
[<Literal>]
let Svy = 4.99852801234917238e-03
[<Literal>]
let Svz = 2.30417297573763929e-05
[<Literal>]
let Smass = 2.85885980666130812e-04
[<Literal>]
let Ux = 1.28943695621391310e+01
[<Literal>]
let Uy = -1.51111514016986312e+01
[<Literal>]
let Uz = -2.23307578892655734e-01
[<Literal>]
let Uvx = 2.96460137564761618e-03
[<Literal>]
let Uvy = 2.37847173959480950e-03
[<Literal>]
let Uvz = -2.96589568540237556e-05
[<Literal>]
let Umass = 4.36624404335156298e-05
[<Literal>]
let Nx = 1.53796971148509165e+01
[<Literal>]
let Ny = -2.59193146099879641e+01
[<Literal>]
let Nz = 1.79258772950371181e-01
[<Literal>]
let Nvx = 2.68067772490389322e-03
[<Literal>]
let Nvy = 1.62824170038242295e-03
[<Literal>]
let Nvz = -9.51592254519715870e-05
[<Literal>]
let Nmass = 5.15138902046611451e-05
open Microsoft.FSharp.NativeInterop
[<Struct>]
type NBody =
val mutable X: float
val mutable Y: float
val mutable Z: float
val mutable VX: float
val mutable VY: float
val mutable VZ: float
val mutable Mass : float
new(x,y,z,vx,vy,vz,mass) = {X=x;Y=y;Z=z;VX=vx;VY=vy;VZ=vz;Mass=mass}
[<EntryPoint>]
let main args =
let ptrBody = NativePtr.stackalloc<NBody> N
// sun
NBody (
0.0,
0.0,
0.0,
(Jvx * Jmass + Svx * Smass + Uvx * Umass + Nvx * Nmass) * -DaysPeryear,
(Jvy * Jmass + Svy * Smass + Uvy * Umass + Nvy * Nmass) * -DaysPeryear,
(Jvz * Jmass + Svz * Smass + Uvz * Umass + Nvz * Nmass) * -DaysPeryear,
Solarmass
) |> NativePtr.set ptrBody 0
NBody ( // jupiter
Jx,
Jy,
Jz,
Jvx * DaysPeryear,
Jvy * DaysPeryear,
Jvz * DaysPeryear,
Jmass * Solarmass
) |> NativePtr.set ptrBody 1
NBody ( // saturn
Sx,
Sy,
Sz,
Svx * DaysPeryear,
Svy * DaysPeryear,
Svz * DaysPeryear,
Smass * Solarmass
) |> NativePtr.set ptrBody 2
NBody ( // uranus
Ux,
Uy,
Uz,
Uvx * DaysPeryear,
Uvy * DaysPeryear,
Uvz * DaysPeryear,
Umass * Solarmass
) |> NativePtr.set ptrBody 3
NBody ( // neptune
Nx,
Ny,
Nz,
Nvx * DaysPeryear,
Nvy * DaysPeryear,
Nvz * DaysPeryear,
Nmass * Solarmass
) |> NativePtr.set ptrBody 4
let inline sqr x = x * x
let inline energy() =
let mutable e = 0.0
for i = 0 to N-1 do
let bi = NativePtr.get ptrBody i
e <- e + 0.5 * bi.Mass * (sqr bi.VX + sqr bi.VY + sqr bi.VZ)
for j = i+1 to N-1 do
let bj = NativePtr.get ptrBody j
e <- e - bi.Mass * bj.Mass /
sqrt(sqr(bi.X-bj.X) + sqr(bi.Y-bj.Y) + sqr(bi.Z-bj.Z))
e
energy().ToString("F9") |> stdout.WriteLine
let inline getD2 dx dy dz =
let d2 = sqr dx + sqr dy + sqr dz
d2 * sqrt d2
let inline moveNative i j =
NativePtr.toNativeInt ptrBody + (7n * nativeint i + j) * 8n
|> NativePtr.ofNativeInt<float>
let inline getNative i j = moveNative i j |> NativePtr.read
let inline setNative i j v = NativePtr.write (moveNative i j) v
let inline addNative i j v =
let p = moveNative i j
NativePtr.read p + v |> NativePtr.write p
let inline subNative i j v =
let p = moveNative i j
NativePtr.read p - v |> NativePtr.write p
let mutable advancements = if args.Length=0 then 1000 else int args.[0]
while (advancements <- advancements - 1; advancements>=0) do
for i = 0 to N-1 do
let biX = getNative i 0n
let biY = getNative i 1n
let biZ = getNative i 2n
let mutable biVX = getNative i 3n
let mutable biVY = getNative i 4n
let mutable biVZ = getNative i 5n
let biMass = getNative i 6n
for j = i+1 to N-1 do
let dx = getNative j 0n - biX
let dy = getNative j 1n - biY
let dz = getNative j 2n - biZ
let mag = DT / getD2 dx dy dz
dx * biMass * mag |> subNative j 3n
dy * biMass * mag |> subNative j 4n
dz * biMass * mag |> subNative j 5n
let bjMass = getNative j 6n * mag
biVX <- biVX + dx * bjMass
biVY <- biVY + dy * bjMass
biVZ <- biVZ + dz * bjMass
biVX * DT |> addNative i 0n
biVY * DT |> addNative i 1n
biVZ * DT |> addNative i 2n
setNative i 3n biVX
setNative i 4n biVY
setNative i 5n biVZ
energy().ToString("F9") |> stdout.WriteLine
exit 0
notes, command-line, and program output
NOTES:
64-bit Ubuntu quad core
.NET SDK 8.0.200
Host Version: 8.0.2
Commit: 1381d5ebd2
<ServerGarbageCollection>true
F# 8.0
Wed, 06 Mar 2024 20:43:56 GMT
MAKE:
cp nbody.fsharpcore-3.fsharpcore Program.fs
cp Include/fsharpcore/tmp.fsproj .
mkdir obj
cp Include/fsharpcore/project.assets.json ./obj
/usr/bin/dotnet build -c Release --no-restore --no-self-contained -r linux-x64
MSBuild version 17.9.4+90725d08d for .NET
tmp -> /home/dunham/all-benchmarksgame/benchmarksgame_i53330/nbody/tmp/bin/Release/net8.0/linux-x64/tmp.dll
Build succeeded.
0 Warning(s)
0 Error(s)
Time Elapsed 00:00:07.26
9.55s to complete and log all make actions
COMMAND LINE:
./bin/Release/net8.0/linux-x64/tmp 50000000
PROGRAM OUTPUT:
-0.169075164
-0.169059907