Thanks for using Compiler Explorer
Sponsors
Jakt
C++
Ada
Analysis
Android Java
Android Kotlin
Assembly
C
C3
Carbon
C++ (Circle)
CIRCT
Clean
CMake
CMakeScript
COBOL
C++ for OpenCL
MLIR
Cppx
Cppx-Blue
Cppx-Gold
Cpp2-cppfront
Crystal
C#
CUDA C++
D
Dart
Elixir
Erlang
Fortran
F#
Go
Haskell
HLSL
Hook
Hylo
ispc
Java
Julia
Kotlin
LLVM IR
LLVM MIR
Modula-2
Nim
Objective-C
Objective-C++
OCaml
OpenCL C
Pascal
Pony
Python
Racket
Ruby
Rust
Snowball
Scala
Solidity
Spice
Swift
LLVM TableGen
Toit
TypeScript Native
V
Vala
Visual Basic
Zig
Javascript
GIMPLE
swift source #1
Output
Compile to binary object
Link to binary
Execute the code
Intel asm syntax
Demangle identifiers
Verbose demangling
Filters
Unused labels
Library functions
Directives
Comments
Horizontal whitespace
Debug intrinsics
Compiler
x86-64 swiftc 3.1.1
x86-64 swiftc 4.0.2
x86-64 swiftc 4.0.3
x86-64 swiftc 4.1
x86-64 swiftc 4.1.1
x86-64 swiftc 4.1.2
x86-64 swiftc 4.2
x86-64 swiftc 5.0
x86-64 swiftc 5.1
x86-64 swiftc 5.10
x86-64 swiftc 5.2
x86-64 swiftc 5.3
x86-64 swiftc 5.4
x86-64 swiftc 5.5
x86-64 swiftc 5.6
x86-64 swiftc 5.7
x86-64 swiftc 5.8
x86-64 swiftc 5.9
x86-64 swiftc nightly
Options
Source code
/// The Computer Language Benchmarks Game /// https://salsa.debian.org/benchmarksgame-team/benchmarksgame/ /// /// contributed by Stuart Carnie /// /// derived from numerous versions, including Swift and Go /// /// This code requires hardware supporting /// CPU feature AVX2, using SIMD4 types. import Foundation let SOLAR_MASS = 4 * Double.pi * Double.pi let DAYS_PER_YEAR = 365.24 public typealias Vec4 = SIMD4<Double> func dot<V : SIMD>(_ a: V, _ b: V) -> V.Scalar where V.Scalar : FloatingPoint { var total = V.Scalar(0) for i in a.indices { total = total.addingProduct(a[i], b[i]) } return total } struct Body { public var pos: Vec4 public var vel: Vec4 public let mass: Double init() { let zero = Vec4(0, 0, 0, 0) self.init(pos: zero, vel: zero, mass: 0) } init(pos: Vec4, vel: Vec4, mass: Double) { self.pos = pos self.vel = vel self.mass = mass } } // ensures loops are static so they can be unrolled explicitly at compile-time let planetsCount = 5 struct System { var bodies: [Body] init(_ bodies: [Body]) { self.bodies = bodies start() } func energy() -> Double { var energy = 0.0 for i in 0..<planetsCount { let bodyi = bodies[i] let posi = bodyi.pos let veli = bodyi.vel energy += 0.5 * bodyi.mass * dot(veli, veli) for j in i+1..<planetsCount { let bodyj = bodies[j] let posj = bodyj.pos let d = posi - posj let distance = sqrt(dot(d, d)) energy -= (bodyi.mass * bodyj.mass) / distance } } return energy } mutating func step(dt: Double) { for i in 0..<planetsCount { let bodyi = bodies[i] for j in i+1..<planetsCount { let bodyj = bodies[j] let d = bodyi.pos - bodyj.pos let dSquared = dot(d, d) let distance = sqrt(dSquared) let mag = dt / (dSquared * distance) let nveli = bodies[i].vel - d * SIMD4(repeating: bodyj.mass) * mag bodies[i].vel = SIMD4(nveli.x, nveli.y, nveli.z, nveli.w) let nvelj = bodies[j].vel + d * SIMD4(repeating: bodyi.mass) * mag bodies[j].vel = SIMD4(nvelj.x, nvelj.y, nvelj.z, nvelj.w) } } for i in 0..<planetsCount { let body = bodies[i] var pos = body.pos pos += dt * body.vel bodies[i].pos = SIMD4(pos.x, pos.y, pos.z, pos.w) } } // Adjusts momentum of the sun private mutating func start() { var p = SIMD4<Double>() for i in 0..<planetsCount { let body = bodies[i] let mass = SIMD4(repeating: body.mass) let vel = body.vel p += vel * mass } let solarMass = SIMD4(repeating: SOLAR_MASS) let vel = -p / solarMass bodies[0].vel = SIMD4(vel.x, vel.y, vel.z, vel.w) } } var sun = Body( pos: SIMD4(x: 0.0, y: 0.0, z: 0.0, w: 0), vel: SIMD4(x: 0.0, y: 0.0, z: 0.0, w: 0), mass: SOLAR_MASS) var jupiter = Body( pos: SIMD4( x: 4.8414314424647209, y: -1.16032004402742839, z: -0.103622044471123109, w: 0), vel: SIMD4( x: 1.66007664274403694e-03 * DAYS_PER_YEAR, y: 7.69901118419740425e-03 * DAYS_PER_YEAR, z: -6.90460016972063023e-05 * DAYS_PER_YEAR, w: 0), mass: 9.54791938424326609e-04 * SOLAR_MASS) var saturn = Body( pos: SIMD4( x: 8.34336671824457987, y: 4.12479856412430479, z: -4.03523417114321381e-01, w: 0), vel: SIMD4( x: -2.76742510726862411e-03 * DAYS_PER_YEAR, y: 4.99852801234917238e-03 * DAYS_PER_YEAR, z: 2.30417297573763929e-05 * DAYS_PER_YEAR, w: 0), mass: 2.85885980666130812e-04 * SOLAR_MASS ) var uranus = Body( pos: SIMD4( x: 1.28943695621391310e+01, y: -1.51111514016986312e+01, z: -2.23307578892655734e-01, w: 0), vel: SIMD4( x: 2.96460137564761618e-03 * DAYS_PER_YEAR, y: 2.37847173959480950e-03 * DAYS_PER_YEAR, z: -2.96589568540237556e-05 * DAYS_PER_YEAR, w: 0), mass: 4.36624404335156298e-05 * SOLAR_MASS ) var neptune = Body( pos: SIMD4( x: 1.53796971148509165e+01, y: -2.59193146099879641e+01, z: 1.79258772950371181e-01, w: 0), vel: SIMD4( x: 2.68067772490389322e-03 * DAYS_PER_YEAR, y: 1.62824170038242295e-03 * DAYS_PER_YEAR, z: -9.51592254519715870e-05 * DAYS_PER_YEAR, w: 0), mass: 5.15138902046611451e-05 * SOLAR_MASS ) var system = System([sun, jupiter, saturn, uranus, neptune]) let n: Int if CommandLine.argc > 1 { n = Int(CommandLine.arguments[1]) ?? 1000 } else { n = 1000 } print(String(format: "%.9f", system.energy())) for _ in 0..<n { system.step(dt: 0.01) } print(String(format: "%.9f", system.energy()))
Become a Patron
Sponsor on GitHub
Donate via PayPal
Source on GitHub
Mailing list
Installed libraries
Wiki
Report an issue
How it works
Contact the author
CE on Mastodon
About the author
Statistics
Changelog
Version tree