|
| 1 | +const PI = Math.PI; |
| 2 | +const SOLAR_MASS = 4 * PI * PI; |
| 3 | +const DAYS_PER_YEAR = 365.24; |
| 4 | + |
| 5 | +function Body(x, y, z, vx, vy, vz, mass) { |
| 6 | + this.x = x; |
| 7 | + this.y = y; |
| 8 | + this.z = z; |
| 9 | + this.vx = vx; |
| 10 | + this.vy = vy; |
| 11 | + this.vz = vz; |
| 12 | + this.mass = mass; |
| 13 | +} |
| 14 | + |
| 15 | +function Jupiter() { |
| 16 | + return new Body( |
| 17 | + 4.84143144246472090e+00, |
| 18 | + -1.16032004402742839e+00, |
| 19 | + -1.03622044471123109e-01, |
| 20 | + 1.66007664274403694e-03 * DAYS_PER_YEAR, |
| 21 | + 7.69901118419740425e-03 * DAYS_PER_YEAR, |
| 22 | + -6.90460016972063023e-05 * DAYS_PER_YEAR, |
| 23 | + 9.54791938424326609e-04 * SOLAR_MASS |
| 24 | + ); |
| 25 | +} |
| 26 | + |
| 27 | +function Saturn() { |
| 28 | + return new Body( |
| 29 | + 8.34336671824457987e+00, |
| 30 | + 4.12479856412430479e+00, |
| 31 | + -4.03523417114321381e-01, |
| 32 | + -2.76742510726862411e-03 * DAYS_PER_YEAR, |
| 33 | + 4.99852801234917238e-03 * DAYS_PER_YEAR, |
| 34 | + 2.30417297573763929e-05 * DAYS_PER_YEAR, |
| 35 | + 2.85885980666130812e-04 * SOLAR_MASS |
| 36 | + ); |
| 37 | +} |
| 38 | + |
| 39 | +function Uranus() { |
| 40 | + return new Body( |
| 41 | + 1.28943695621391310e+01, |
| 42 | + -1.51111514016986312e+01, |
| 43 | + -2.23307578892655734e-01, |
| 44 | + 2.96460137564761618e-03 * DAYS_PER_YEAR, |
| 45 | + 2.37847173959480950e-03 * DAYS_PER_YEAR, |
| 46 | + -2.96589568540237556e-05 * DAYS_PER_YEAR, |
| 47 | + 4.36624404335156298e-05 * SOLAR_MASS |
| 48 | + ); |
| 49 | +} |
| 50 | + |
| 51 | +function Neptune() { |
| 52 | + return new Body( |
| 53 | + 1.53796971148509165e+01, |
| 54 | + -2.59193146099879641e+01, |
| 55 | + 1.79258772950371181e-01, |
| 56 | + 2.68067772490389322e-03 * DAYS_PER_YEAR, |
| 57 | + 1.62824170038242295e-03 * DAYS_PER_YEAR, |
| 58 | + -9.51592254519715870e-05 * DAYS_PER_YEAR, |
| 59 | + 5.15138902046611451e-05 * SOLAR_MASS |
| 60 | + ); |
| 61 | +} |
| 62 | + |
| 63 | +function Sun() { |
| 64 | + return new Body(0.0, 0.0, 0.0, 0.0, 0.0, 0.0, SOLAR_MASS); |
| 65 | +} |
| 66 | + |
| 67 | +const bodies = Array(Sun(), Jupiter(), Saturn(), Uranus(), Neptune()); |
| 68 | + |
| 69 | +function offsetMomentum() { |
| 70 | + let px = 0; |
| 71 | + let py = 0; |
| 72 | + let pz = 0; |
| 73 | + const size = bodies.length; |
| 74 | + for (let i = 0; i < size; i++) { |
| 75 | + const body = bodies[i]; |
| 76 | + const mass = body.mass; |
| 77 | + px += body.vx * mass; |
| 78 | + py += body.vy * mass; |
| 79 | + pz += body.vz * mass; |
| 80 | + } |
| 81 | + |
| 82 | + const body = bodies[0]; |
| 83 | + body.vx = -px / SOLAR_MASS; |
| 84 | + body.vy = -py / SOLAR_MASS; |
| 85 | + body.vz = -pz / SOLAR_MASS; |
| 86 | +} |
| 87 | + |
| 88 | +function advance(dt) { |
| 89 | + const size = bodies.length; |
| 90 | + |
| 91 | + for (let i = 0; i < size; i++) { |
| 92 | + const bodyi = bodies[i]; |
| 93 | + let vxi = bodyi.vx; |
| 94 | + let vyi = bodyi.vy; |
| 95 | + let vzi = bodyi.vz; |
| 96 | + for (let j = i + 1; j < size; j++) { |
| 97 | + const bodyj = bodies[j]; |
| 98 | + const dx = bodyi.x - bodyj.x; |
| 99 | + const dy = bodyi.y - bodyj.y; |
| 100 | + const dz = bodyi.z - bodyj.z; |
| 101 | + |
| 102 | + const d2 = dx * dx + dy * dy + dz * dz; |
| 103 | + const mag = dt / (d2 * Math.sqrt(d2)); |
| 104 | + |
| 105 | + const massj = bodyj.mass; |
| 106 | + vxi -= dx * massj * mag; |
| 107 | + vyi -= dy * massj * mag; |
| 108 | + vzi -= dz * massj * mag; |
| 109 | + |
| 110 | + const massi = bodyi.mass; |
| 111 | + bodyj.vx += dx * massi * mag; |
| 112 | + bodyj.vy += dy * massi * mag; |
| 113 | + bodyj.vz += dz * massi * mag; |
| 114 | + } |
| 115 | + bodyi.vx = vxi; |
| 116 | + bodyi.vy = vyi; |
| 117 | + bodyi.vz = vzi; |
| 118 | + } |
| 119 | + |
| 120 | + for (let i = 0; i < size; i++) { |
| 121 | + const body = bodies[i]; |
| 122 | + body.x += dt * body.vx; |
| 123 | + body.y += dt * body.vy; |
| 124 | + body.z += dt * body.vz; |
| 125 | + } |
| 126 | +} |
| 127 | + |
| 128 | +function energy() { |
| 129 | + let e = 0; |
| 130 | + const size = bodies.length; |
| 131 | + |
| 132 | + for (let i = 0; i < size; i++) { |
| 133 | + const bodyi = bodies[i]; |
| 134 | + |
| 135 | + e += 0.5 * bodyi.mass * ( bodyi.vx * bodyi.vx + bodyi.vy * bodyi.vy + bodyi.vz * bodyi.vz ); |
| 136 | + |
| 137 | + for (let j = i + 1; j < size; j++) { |
| 138 | + const bodyj = bodies[j]; |
| 139 | + const dx = bodyi.x - bodyj.x; |
| 140 | + const dy = bodyi.y - bodyj.y; |
| 141 | + const dz = bodyi.z - bodyj.z; |
| 142 | + |
| 143 | + const distance = Math.sqrt(dx * dx + dy * dy + dz * dz); |
| 144 | + e -= (bodyi.mass * bodyj.mass) / distance; |
| 145 | + } |
| 146 | + } |
| 147 | + return e; |
| 148 | +} |
| 149 | + |
| 150 | +const n = +50000000; |
| 151 | + |
| 152 | +offsetMomentum(); |
| 153 | + |
| 154 | +console.log(energy().toFixed(9)); |
| 155 | +const start = Date.now(); |
| 156 | +for (let i = 0; i < n; i++) { |
| 157 | + advance(0.01); |
| 158 | +} |
| 159 | +const end = Date.now(); |
| 160 | +console.log(energy().toFixed(9)); |
| 161 | +console.log("elapsed:",end-start); |
0 commit comments