function doc (id) { return document.getElementById (id); } function nodeWhen (node, ev, func) { node.addEventListener (ev, func); } function docWhen (id, ev, func) { nodeWhen (doc (id), ev, func); } function clr () { doc ("output") .innerHTML = "Output:
"; } function wr (a = "") { doc ("output") .innerHTML += a; } function wrln (a = "") {wr (a + "
");} function lnwr (a = "") {wr ("
" + a);} function rdFloat (a, b) { return Math.random () * (b - a) + a; } function rdInt (a, b) { return Math.floor (rdFloat (a, b)); } function clrNull (arr) { let i = 0, j = 0, n = arr.length; for (; j < n; j++) { if (arr[j] !== null) { if (i < j) { arr[i] = arr[j]; } i++; } } arr.length = i; } function square (a, b = 0, c = 0) { return a * a + b * b + c * c; } function abs (a, b = 0, c = 0) { return Math.sqrt (square (a, b, c)); } function floatStr (a, b) { let c, d, e, f, g, h; if (a == 0) {return "0";} if (a > 0) {c = "";} else {c = "-"; a *= -1;} g = Math.floor (Math.log10 (a)); d = Math.round (a * Math.pow (10, b - 1 - g)); f = d % Math.pow (10, b - 1); d = Math.floor (d / Math.pow (10, b - 1)); e = "."; h = Math.pow (10, b - 1); while (h > 1) { h /= 10; if (f < h) {e += "0";} else {e += f; break;} } return c + d + e + "*10^" + g; } let deg = Math.PI / 180, rad = 180 / Math.PI; function stdRad (a) { let _2pi = 2 * Math.PI; a = a % _2pi; if (a < 0) {a += _2pi;} return a; } function filterToColor (r, g, b) { r = Math.round (r); g = Math.round (g); b = Math.round (b); if (r < 0) {r = 0;} else if (r > 255) {r = 255;} if (g < 0) {g = 0;} else if (g > 255) {g = 255;} if (b < 0) {b = 0;} else if (b > 255) {b = 255;} return toColor (r, g, b); } function toColor (r, g, b) { let color = "#"; color += toColor1 (r); color += toColor1 (g); color += toColor1 (b); return color; } function toColor1 (ce) { let m, n, c = ""; n = ce % 16; m = (ce - n) / 16; c += toColor2 (m); c += toColor2 (n); return c; } function toColor2 (a) { if (a < 10) { return "" + a; } else { switch (a) { case 10: return "a"; case 11: return "b"; case 12: return "c"; case 13: return "d"; case 14: return "e"; case 15: return "f"; } } } class Vec { static obj = []; del () { Vec.obj.push (this); } static neww (x = 0, y = 0, z = 0) { let vec; if (Vec.obj.length > 0) { vec = Vec.obj.pop () .setup (x, y, z); } else { vec = new Vec (x, y, z); } return vec; } constructor (x = 0, y = 0, z = 0) { this.setup (x, y, z); } setup (x = 0, y = 0, z = 0) { this.temp = false; this.x = x; this.y = y; this.z = z; return this; } rad () { return abs (this.x, this.y, this.z); } eql (a) { this.x = a.x; this.y = a.y; this.z = a.z; if (a.temp) { a.del (); } return this; } add (a) { let r = Vec.neww (this.x + a.x, this.y + a.y, this.z + a.z); r.temp = true; if (this.temp) {this.del ();} if (a.temp) {a.del ();} return r; } sub (a) { let r = Vec.neww (this.x - a.x, this.y - a.y, this.z - a.z); r.temp = true; if (this.temp) {this.del ();} if (a.temp) {a.del ();} return r; } mul (k) { let r = Vec.neww (this.x * k, this.y * k, this.z * k); r.temp = true; if (this.temp) {this.del ();} return r; } div (k) { let r = Vec.neww (this.x / k, this.y / k, this.z / k); r.temp = true; if (this.temp) {this.del ();} return r; } sclMul (a) { return this.x * a.x + this.y * a.y + this.z * a.z; } vecMul (a) { let r = Vec.neww (this.y * a.z - this.z * a.y, this.z * a.x - this.x * a.z, this.x * a.y - this.y * a.x); r.temp = true; if (this.temp) {this.del ();} if (b.temp) {b.del ();} return r; } ang (a) { return Math.acos (this.sclMul (a) / (this.rad () * a.rad ())); } rotate (axis, alpha) { let bz = Vec.neww (); bz .eql (axis .div (axis.rad ())); bz .eql (bz .mul (a .sclMul (bz))); let bx = Vec.neww (); bx .eql (a .sub (bz)); let r = bx.rad (); bx .eql (bx .div (r)); let by = Vec.neww (); by .eql (bz .vecMul (bx) .div (by.rad ())); bx .eql (bx .mul (r * Math.cos (alpha))); by .eql (by .mul (r * Math.sin (alpha))); a .eql (bx .add (by) .add (bz)); bx.del (); by.del (); bz.del (); } } function solve2eq2var (a, b, c, d, e, f, sol) { let det = a * e - d * b; let detx = c * e - f * b; let dety = a * f - d * c; sol.x = detx / det; sol.y = dety / det; } class Planet { static obj = []; del () { Planet.obj.push (this); } static neww (system = null, id = 0, mass = 0, type = 0, radius = 0, rx = 0, ry = 0, rz = 0, vx = 0, vy = 0, vz = 0, vRadius = 0, red = 0, green = 0, blue = 0) { let pln; if (Planet.obj.length > 0) { pln = Planet.obj.pop () .constr (system, id, mass, type, radius, rx, ry, rz, vx, vy, vz, vRadius, red, green, blue); } else { pln = new Planet (system, id, mass, type, radius, rx, ry, rz, vx, vy, vz, vRadius, red, green, blue); } return pln; } constructor (system = null, id = 0, mass = 0, type = 0, radius = 0, rx = 0, ry = 0, rz = 0, vx = 0, vy = 0, vz = 0, vRadius = 0, red = 0, green = 0, blue = 0) { this.constr (system, id, mass, type, radius, rx, ry, rz, vx, vy, vz, vRadius, red, green, blue); } // use for general case (without "made up" rules) constr (system = null, id = 0, mass = 0, type = 0, radius = 0, rx = 0, ry = 0, rz = 0, vx = 0, vy = 0, vz = 0, vRadius = 0, red = 0, green = 0, blue = 0) { this.system = system; this.id = id; this.mass = mass; this.type = type; this.radius = radius; this.temp4 = 0; this.absorb = 0; this.absorbRed = 0; this.absorbGreen = 0; this.absorbBlue = 0; this.lightRed = 0; this.lightGreen = 0; this.lightBlue = 0; if (this.lightSecBound === undefined) { this.lightSecBound = []; this.lightSecRed = []; this.lightSecGreen = []; this.lightSecBlue = []; } else { this.lightSecBound.length = 0; this.lightSecRed.length = 0; this.lightSecGreen.length = 0; this.lightSecBlue.length = 0; } this.rx = rx; this.ry = ry; this.rz = rz; this.vx = vx; this.vy = vy; this.vz = vz; this.ax = 0; this.ay = 0; this.az = 0; this.willDraw = true; this.orbitPln = [null, 0, 0]; this.vrx = 0; this.vry = 0; this.vrz = 0; this.vRadius = vRadius; this.red = red; this.green = green; this.blue = blue; this.redRatio = 0; this.greenRatio = 0; this.blueRatio = 0; return this; } // use in this specific simulation setup (mass, rx, ry, rz, vx, vy, vz) { this.id = this.system.totalPlanets; this.system.totalPlanets++; this.mass = mass; this.rx = rx; this.ry = ry; this.rz = rz; this.vx = vx; this.vy = vy; this.vz = vz; this.setupRule (); if (this.type == 1) { this.red = rdInt (0, 256); this.green = rdInt (0, 256); this.blue = rdInt (0, 256); } else { this.red = 512; this.green = 256; this.blue = 128; } let rgb, a; if (this.system.drawType == 2 && this.type == 1) { a = this.system.scaleColorRGBA; this.absorbRed = 1 - Math.exp ((this.red - 255) / a); this.absorbGreen = 1 - Math.exp ((this.green - 255) / a); this.absorbBlue = 1 - Math.exp ((this.blue - 255) / a); this.absorb = (this.absorbRed + this.absorbGreen + this.absorbBlue) / 3; } if (this.system.drawType == 1) { rgb = (this.red + this.green + this.blue) / 3; this.redRatio = this.red / rgb; this.greenRatio = this.green / rgb; this.blueRatio = this.blue / rgb; this.absorb = rdFloat (0, 1); } return this; } merge (pln) { let totMass = this.mass + pln.mass, a = this.mass / totMass, b = pln.mass / totMass; this.mass = totMass; this.setupRule (); if (this.system.drawType == 2 && this.type == 1) { this.absorb = this.absorb * a + pln.absorb * b; this.absorbRed = this.absorbRed * a + pln.absorbRed * b; this.absorbGreen = this.absorbGreen * a + pln.absorbGreen * b; this.absorbBlue = this.absorbBlue * a + pln.absorbBlue * b; } this.rx = this.rx * a + pln.rx * b; this.ry = this.ry * a + pln.ry * b; this.rz = this.rz * a + pln.rz * b; this.vx = this.vx * a + pln.vx * b; this.vy = this.vy * a + pln.vy * b; this.vz = this.vz * a + pln.vz * b; let sc; if (this.system.drawType == 2 && this.type == 1) { sc = this.system.scaleColorRGBA; this.red = Math.log (1 - this.absorbRed) * sc + 255; this.green = Math.log (1 - this.absorbGreen) * sc + 255; this.blue = Math.log (1 - this.absorbBlue) * sc + 255; } else { this.red = this.red * a + pln.red * b; this.green = this.green * a + pln.green * b; this.blue = this.blue * a + pln.blue * b; } if (this.system.drawType == 1) { this.redRatio = this.redRatio * a + pln.redRatio * b; this.greenRatio = this.greenRatio * a + pln.greenRatio * b; this.blueRatio = this.blueRatio * a + pln.blueRatio * b; } let tr = this.system.tracking; if (tr [0] && tr [1] === pln) { tr [1] = this; } return this; } // some common rule setupRule () { let s = this.system; this.type = this.mass >= 1.6e29 ? 0 : 1; let denst = this.type == 0 ? s.fusionDensity : s.rockDensity; this.radius = GravSys.calRadius (denst, this.mass); let temp; if (this.type == 0) { this.temp4 = s.calFusionTemp4 (this.mass); temp = Math.sqrt (Math.sqrt (this.temp4)); this.lightRed = s.calLightRed (temp); this.lightGreen = s.calLightGreen (temp); this.lightBlue = s.calLightBlue (temp); } this.vRadius = s.calVRadius (this.mass); } nextPosVel (time) { let dvx = this.ax * time, dvy = this.ay * time, dvz = this.az * time, drx = (this.vx + dvx / 2) * time, dry = (this.vy + dvy / 2) * time, drz = (this.vz + dvz / 2) * time; this.rx += drx; this.ry += dry; this.rz += drz; this.vx += dvx; this.vy += dvy; this.vz += dvz; } clearShine () { const mode = this.system.drawType; this.temp4 = 0; this.lightSecBound.length = 0; switch (mode) { case 2: this.lightSecRed.length = 0; this.lightSecBlue.length = 0; case 1: this.lightSecGreen.length = 0; } } addShine (pln) { const mode = this.system.drawType; const dx = pln.rx - this.rx; const dy = pln.ry - this.ry; const dz = pln.rz - this.rz; const angle = Math.atan2 (dy, dx); let sta = stdRad (angle - 90*deg), end = stdRad (angle + 90*deg), sm, lg; if (sta < end) {sm = sta; lg = end;} else {sm = end; lg = sta;} const distRatio = pln.radius * pln.radius / square (dx, dy, dz); this.temp4 += pln.temp4 * distRatio / 4; const light = 5.67e-8 * pln.temp4 * distRatio / 2 * (1 - this.absorb); const lRed = pln.lightRed * distRatio / 2 * (1 - this.absorbRed); const lGreen = pln.lightGreen * distRatio / 2 * (1 - this.absorbGreen); const lBlue = pln.lightBlue * distRatio / 2 * (1 - this.absorbBlue); const bnd = this.lightSecBound, red = this.lightSecRed, green = this.lightSecGreen, blue = this.lightSecBlue; let id, len; if (bnd.length == 0) { bnd.splice (0, 0, sm, lg); if (sm == sta) { switch (mode) { case 1: green.splice (0, 0, light, 0); break; case 2: red.splice (0, 0, lRed, 0); green.splice (0, 0, lGreen, 0); blue.splice (0, 0, lBlue, 0); } } else { switch (mode) { case 1: green.splice (0, 0, 0, light); break; case 2: red.splice (0, 0, 0, lRed); green.splice (0, 0, 0, lGreen); blue.splice (0, 0, 0, lBlue); } } } else { len = bnd.length; for (let i = 0; i < len; i++) { if (bnd [i] >= sm) { id = i; break; } } bnd.splice (id + len / 2, 0, lg); bnd.splice (id, 0, sm); switch (mode) { case 1: green.splice (id + len / 2, 0, green [(id + len / 2 - 1) % len]); green.splice (id, 0, green [(id + len) % (len + 1)]); break; case 2: red.splice (id + len / 2, 0, red [(id + len / 2 - 1) % len]); red.splice (id, 0, red [(id + len) % (len + 1)]); green.splice (id + len / 2, 0, green [(id + len / 2 - 1) % len]); green.splice (id, 0, green [(id + len) % (len + 1)]); blue.splice (id + len / 2, 0, blue [(id + len / 2 - 1) % len]); blue.splice (id, 0, blue [(id + len) % (len + 1)]); } len = bnd.length; if (bnd [id] == end) { id += len / 2; } switch (mode) { case 1: for (let i = 0; i < len / 2; i++) { green [(id + i) % len] += light; } break; case 2: for (let i = 0; i < len / 2; i++) { red [(id + i) % len] += lRed; green [(id + i) % len] += lGreen; blue [(id + i) % len] += lBlue; }} } } calShine () { let st = this.system; let mode = st.drawType; let lRed, lGreen, lBlue, temp, light; switch (mode) { case 1: light = 5.67e-8 * this.temp4 * this.absorb; break; case 2: temp = Math.sqrt (Math.sqrt (this.temp4)); lRed = st.calLightRed (temp) * this.absorbRed; lGreen = st.calLightGreen (temp) * this.absorbGreen; lBlue = st.calLightBlue (temp) * this.absorbBlue; } const red = this.lightSecRed, green = this.lightSecGreen, blue = this.lightSecBlue, len = green.length; switch (mode) { case 1: for (let i = 0; i < len; i++) { green [i] += light; } break; case 2: for (let i = 0; i < len; i++) { red [i] += lRed; green [i] += lGreen; blue [i] += lBlue; }} } // standard 2-body gravity orbit drawOrbit (c) { c.save (); c.lineWidth = 1; let color = filterToColor (this.red, this.green, this.blue); c.strokeStyle = color; let pln = this.orbitPln [0]; c.translate (pln.vrx, pln.vry); let r = Vec.neww (this.rx - pln.rx, -(this.ry - pln.ry), this.rz - pln.rz); let v = Vec.neww (this.vx - pln.vx, -(this.vy - pln.vy), this.vz - pln.vz); let L = Vec.neww (); L .eql (r .vecMul (v)); let L2 = L.rad (); let mode = 0; if (L2 == 0) { return; mode = 2; } L2 *= L2; let rr = r.rad (); let A = 6.67e-11 * (pln.mass + this.mass); let E = v.rad (); E = E * E / 2 - A / rr; if (mode == 0 && E >= 0) { mode=1; } let p = L2 / A; let a = -A / (2 * E); let ELA = -p / a; let ecc = Math.sqrt (1 + ELA); let cenDif = a * ecc; let b = Math.sqrt (a * a - cenDif * cenDif); let phi = Math.acos ((p / rr - 1) / ecc); let r_v = r .sclMul (v); if (r_v >= 0) { phi = -phi; } let rmin = Vec.neww () .eql (r); rmin .rotate (L, phi); let Oz = Vec.neww (0, 0, 1); let itsl = Vec.neww () .eql (L .vecMul (Oz)); let itslr = itsl.rad (); if (itslr !== 0) { itsl .eql (itsl .div (itslr)); phi = Math.atan2 (itsl.y, itsl.x); c.rotate (phi); let phi1 = L .ang (Oz); L .rotate (itsl, phi1); rmin .rotate (itsl, phi1); c.scale (1, Math.cos (phi1)); c.rotate (-phi); } else if (L.z < 0) { phi = Math.atan2 (rmin.y, rmin.x); c.rotate (phi); L .rotate (rmin, 180*deg); c.scale (1, -1); c.rotate (-phi); } phi = Math.atan2 (rmin.y, rmin.x); c.rotate (phi); let meter = 1 / this.system.scaleLength; a *= meter; b *= meter; cenDif *= meter; p *= meter; if (mode == 0) { c.translate (-cenDif, 0); c.scale (1, b / a); c.beginPath (); c.arc (0, 0, a, 0, 360*deg); c.stroke (); } else if (mode == 1) { let x0, y0, x1, y1, crX, crY, eq00, eq01, eq10, eq11, eq02, eq12, er, r2; er = ecc * p; r2 = rr * p; let sqrt, yl, yu, xm, ym, l2, lx, ly; y0 = 0; sqrt = Math.sqrt (ELA * y0 * y0 + r2); x0 = (er - sqrt) / ELA; yl = 0; yu = 10; do { ym = (yl + yu) / 2; sqrt = Math.sqrt (ELA * ym * ym + r2); xm = (er - sqrt) / ELA; lx = xm - x0; ly = ym - y0; l2 = lx * lx + ly * ly; if (l2 < 16) {yl = ym;} if (l2 > 36) {yu = ym;} } while (l2 < 16 || l2 > 36); x0 = xm; y0 = ym; for (let i = 0; i < 15; i++) { eq00 = -y0; eq10 = sqrt; yl = y0; yu = y0 + 20; do { ym = (yl + yu) / 2; sqrt = Math.sqrt (ELA * ym * ym + r2); xm = (er - sqrt) / ELA; lx = xm - x0; ly = ym - y0; l2 = lx * lx + ly * ly; if (l2 < 81) {yl = ym;} if (l2 > 121) {yu = ym;} } while (l2 < 81 || l2 > 121); x1 = xm; y1 = ym; eq01 = y1; eq11 = -sqrt; eq02 = x1 - x0; eq12 = y1 - y0; let det0 = eq02 * eq11 - eq12 * eq01; let det1 = eq00*eq12 - eq10*eq02; let det = eq00*eq11-eq10*eq01; let t0 = det0/det; let t1 = det1/det; crX = eq00*t0+x0; crY = eq10*t0+y0; c.beginPath (); c.moveTo (x0, y0); c.quadraticCurveTo (crX, crY, x1, y1); c.moveTo (x0, -y0); c.quadraticCurveTo (crX, -crY, x1, -y1); c.stroke (); if (L > 0) { c.beginPath (); c.arc (x1, y1, 1.5, 0, 360*deg); c.fill (); c.beginPath (); c.arc (x0, -y0, 1.5, 0, 360*deg); c.fill (); } else { c.beginPath (); c.arc (x0, y0, 1.5, 0, 360*deg); c.fill (); c.beginPath (); c.arc (x1, -y1, 1.5, 0, 360*deg); c.fill (); } yl = y1; yu = y1 + 20; do { ym = (yl + yu) / 2; sqrt = Math.sqrt (ELA * ym * ym + r2); xm = (er - sqrt) / ELA; lx = xm - x1; ly = ym - y1; l2 = lx * lx + ly * ly; if (l2 < 81) {yl = ym;} if (l2 > 121) {yu = ym;} } while (l2 < 81 || l2 > 121); x0 = xm; y0 = ym; }} r.del (); rmin.del (); v.del (); L.del (); itsl.del (); Oz.del (); c.restore (); } drawSimple (c) { c.save (); c.lineWidth = 0; c.fillStyle = filterToColor (this.red, this.green, this.blue); c.beginPath (); c.arc (this.vrx, this.vry, this.vRadius, 0, 360*deg); c.fill (); c.restore (); } drawFusion (c) { let color, temp, s = this.system, mode = s.drawType; switch (mode) { case 1: color = s.calColorAll (5.67e-8 * this.temp4); color = filterToColor (color * this.redRatio, color * this.greenRatio, color * this.blueRatio); break; case 2: temp = Math.sqrt (Math.sqrt (this.temp4)); color = s.calColorRGB (s.calLightRed (temp), s.calLightGreen (temp), st.calLightBlue (temp)); } c.save (); c.lineWidth = 0; c.fillStyle = color; c.beginPath (); c.arc (this.vrx, this.vry, this.vRadius, 0, 360*deg); c.fill (); c.restore (); } drawRock (c) { let n = this.lightSecBound.length, s = this.system, mode = s.drawType; c.save (); c.lineWidth = 0; for (let i = 0; i < n; i++) { let color; switch (mode) { case 1: color = s.calColorAll (this.lightSecGreen [i]); color = filterToColor (color * this.redRatio, color * this.greenRatio, color * this.blueRatio); break; case 2: color = s.calColorRGB (this.lightSecRed [i], this.lightSecGreen [i], this.lightSecBlue [i]); } c.fillStyle = color; c.beginPath (); c.arc (this.vrx, this.vry, this.vRadius, this.lightSecBound [i] - 0.05, this.lightSecBound [(i + 1) % n] + 0.05); c.lineTo (this.vrx, this.vry); c.closePath (); c.fill (); } c.restore (); } } class GravSys { static secInYear = 60 * 60 * 24 * 365.25; constructor (canvas, canvas1, width, height) { this.canvas = canvas; this.canvas.width = width; this.canvas.height = height; this.cvct = this.canvas.getContext ("2d"); this.cvMat = [1, 0, 0, 0, 1, 0]; this.canvas1 = canvas1; this.canvas1.width = width; this.canvas1.height = height; this.cvct1 = this.canvas1.getContext ("2d"); this.proc = null; this.pausing = false; this.firstSetup = true; this.mergeType = 0; this.mergeParam = 0; this.moveRateType = 0; this.moveRateParam = 0; this.moveRateParam1 = 0; this.moveRateStat = 0; this.drawType = 0; this.drawOrbit = false; this.drawParam1 = 0; this.cenType = 0; this.planets = []; this.cenPlanets = []; this.cenRatios = []; this.newAstAm = 0; this.radius = 0; this.time = 0; this.timeDif = 0; this.moveTime = 0; this.rockDensity = 0; this.fusionDensity = 0; this.fusionTemp4 = 0; this.maxLightRed = 0; this.maxLightGreen = 0; this.maxLightBlue = 0; this.lightFreq0 = 0; this.lightFreq1 = 0; this.lightFreq2 = 0; this.lightFreq3 = 0; this.centerX = this.canvas.width / 2; this.centerY = this.canvas.height / 2; this.vtime = 0; this.vtimeDif = 0; this.tracking = [false, null]; this.scaleLength = 0; this.scaleTime = 0; this.scaleVRadiusA = 0; this.scaleVRadiusB = 0; this.scaleColorAllA = 0; this.scaleColorAllB = 0; this.scaleColorRGBA = 0; this.scaleColorRGBB = 0; this.drawCount = 0; this.moveCount = 0; this.test = null; this.stat = null; this.statStr = ""; this.totalPlanets = 0; } static calDensity (mass, radius) { return 3 * mass / (4 * Math.PI * radius * radius * radius); } static calRadius (denst, mass) { return Math.cbrt (3 * mass / (4 * Math.PI * denst)); } static stdVelocity (a, mass0, mass1, radius) { return Math.sqrt (a * 6.67e-11 * (mass0 + mass1) / radius); } // integral from freq0 to freq1 of 2pi * f^2/c^2 * hf / (e^(hf/kT) - 1) * (1 - |f - f0| / deltaf); static _2pih_c2 = 2 * Math.PI * 6.625e-34 / (3e8 * 3e8); static h_k = 6.625e-34 / 1.38e-23; static lightSpectrum (temp, freq0, freq1) { let integ = 0; const n = 10; const freqCen = (freq0 + freq1) / 2; const freqDif = (freq1 - freq0) / n; const freqRange = (freq1 - freq0) / 2; const h_kT = GravSys.h_k / temp; let sum = 0; for (let i = 2; i < n; i += 2) { const freq = freq0 + freqDif * i; sum += freq * freq * freq / (Math.exp (h_kT * freq) - 1) * (1 - Math.abs (freq - freqCen) / freqRange); } integ = sum * 2; sum = 0; for (let i = 1; i < n; i += 2) { const freq = freq0 + freqDif * i; sum += freq * freq * freq / (Math.exp (h_kT * freq) - 1) * (1 - Math.abs (freq - freqCen) / freqRange); } integ += sum * 4; integ *= GravSys._2pih_c2 * freqDif / 3; return integ; } calVRadius (mass) { return Math.pow (mass, this.scaleVRadiusA) * this.scaleVRadiusB; } scaleVRadius (mass0, radius0, mass1, radius1) { let sol = Vec.neww (); solve2eq2var (Math.log (mass0), 1, Math.log (radius0), Math.log (mass1), 1, Math.log (radius1), sol); this.scaleVRadiusA = sol.x; this.scaleVRadiusB = Math.exp (sol.y); sol.del (); } calFusionTemp4 (mass) { return this.fusionTemp4 * Math.cbrt (mass); } calColorAll (light) { return Math.log (light) * this.scaleColorAllA + this.scaleColorAllB; } scaleColorAll (light0, color0, light1, color1) { let sol = Vec.neww (); solve2eq2var (Math.log (light0), 1, color0, Math.log (light1), 1, color1, sol); this.scaleColorAllA = sol.x; this.scaleColorAllB = sol.y; sol.del (); } calLightRed (temp) { return GravSys.lightSpectrum (temp, this.lightFreq0, this.lightFreq1); } calLightGreen (temp) { return GravSys.lightSpectrum (temp, this.lightFreq1, this.lightFreq2); } calLightBlue (temp) { return GravSys.lightSpectrum (temp, this.lightFreq2, this.lightFreq3); } calColorRGB (red, green, blue) { red /= this.maxLightRed; green /= this.maxLightGreen; blue /= this.maxLightBlue; let mx = Math.max (red, green, blue); if (mx > 1) { red /= mx; green /= mx; blue /= mx; } red = Math.log (red) * this.scaleColorRGBA + this.scaleColorRGBB; green = Math.log (green) * this.scaleColorRGBA + this.scaleColorRGBB; blue = Math.log (blue) * this.scaleColorRGBA + this.scaleColorRGBB; return filterToColor (red, green, blue); } scaleColorRGB (temp0, color0) { const temp1 = Math.sqrt (Math.sqrt (this.calFusionTemp4 (2e30))); let distRatio = 7e8 / 1.5e11; distRatio *= distRatio; this.maxLightRed = this.calLightRed (temp1) * distRatio; this.maxLightGreen = this.calLightGreen (temp1) * distRatio; this.maxLightBlue = this.calLightBlue (temp1) * distRatio; let minRed = this.calLightRed (temp0); let sol = Vec.neww (); solve2eq2var (Math.log (minRed / this.maxLightRed), 1, color0, 0, 1, 255, sol); this.scaleColorRGBA = sol.x; this.scaleColorRGBB = sol.y; sol.del (); } setup (cen2, plnAm, astAm, newAstAm, mergeType, mergeParam, moveRateType, moveRateParam, moveRateParam1, drawType, drawOrbit, drawParam1, bgStar, cenType) { this.mergeType = mergeType; this.mergeParam = mergeParam; this.moveRateType = moveRateType; this.moveRateParam = moveRateParam; this.moveRateParam1 = moveRateParam1; this.drawType = drawType; this.drawOrbit = drawOrbit; this.drawParam1 = drawParam1; this.cenType = cenType; this.newAstAm = newAstAm; if (this.firstSetup) { alert ( "\"Pause\" / \"Play\" button to pause / play the simulation" + "\n\"Space Unlock\" / \"Space Lock\" button to turn on / off the zooming function" + "\n\"Space Reset\" button to reset the zooming" + "\nDouble tap a point to tracking the closest planet" + "\n\"Track Stop\" button to stop the tracking" + "\n\"Track Random\" button to track a random planet" ); } this.scaleLength = 5e9; // meter/pixel this.scaleTime = 3e7; // model sec / sec this.radius = 1000 * this.scaleLength; this.rockDensity = GravSys.calDensity (6e24, 6.4e6); // earth this.fusionDensity = GravSys.calDensity (2e30, 7e8); // sun this.lightFreq0 = 3e8 / 760e-9; this.lightFreq1 = this.lightFreq0 *4/3; this.lightFreq2 = this.lightFreq0 *5/3; this.lightFreq3 = this.lightFreq0 *2; this.scaleVRadius (2e30, 20, 6e24, 5); // sun and earth this.fusionTemp4 = 1.4e3 * 1.5e11 * 1.5e11 / (7e8 * 7e8 * 5.67e-8 * Math.cbrt (2e30)); let temp4; switch (this.drawType) { case 1: temp4 = 30; temp4 *= temp4; temp4 *= temp4; this.scaleColorAll (5.67e-8 * this.calFusionTemp4 (2e30), 255, 5.67e-8 * temp4, 0); break; case 2: this.scaleColorRGB (700, 0); } this.time = Math.round (new Date () . getTime () / 1000); for (let i = this.planets.length - 1; i >= 0; i--) { this.destroy (i); } this.planets.length = 0; this.cenPlanets.length = 0; this.cenRatios.length = 0; this.totalPlanets = 0; this.tracking [0] = false; this.tracking [1] = null; this.statStr = ""; if (bgStar) { this.backgroundStar (); } else { this.cvct1.clearRect (0, 0, this.canvas1.width, this.canvas1.height); } let pln = Planet.neww (this) .setup (2e30, 0, 0, 0, 0, 0, 0); this.planets.push (pln); this.cenPlanets.push (pln); this.cenRatios.push (100); if (cen2) { let radius = 4e11; let mass = 2e29; pln = Planet.neww (this) .setup (mass, radius, 0, 0, 0, -GravSys.stdVelocity (1, this.cenPlanets[0].mass, mass, radius), 0); this.planets.push (pln); this.cenPlanets.push (pln); this.cenRatios.push (0); } for (let i = 0; i < plnAm; i++) { let radius = 2e12; let x, y; do { x = rdFloat (-radius, radius); y = rdFloat (-radius, radius); } while (x * x + y * y > radius * radius); radius = abs (x, y); let alpha = Math.atan2 (y, x); let mass = rdFloat (1e22, 1e26); let stbVelo = GravSys.stdVelocity (1, this.cenPlanets[0].mass, mass, radius); let dv = stbVelo / 5; let zRatio = 0.1*0; this.planets.push (Planet.neww (this) .setup ( mass, x, y, radius * rdFloat (-zRatio, zRatio), stbVelo * Math.cos (alpha - 90*deg) + rdFloat (-dv, dv), stbVelo * Math.sin (alpha - 90*deg) + rdFloat (-dv, dv), stbVelo * rdFloat (-zRatio, zRatio) )); } for (let i = 0; i < astAm; i++) { let radius = 2e12; let x, y; do { x = rdFloat (-radius, radius); y = rdFloat (-radius, radius); } while (x * x + y * y > radius * radius); radius = abs (x, y); let alpha = Math.atan2 (y, x); let mass = rdFloat (1e20, 1e22); let stbVelo = GravSys.stdVelocity (1, this.cenPlanets[0].mass, mass, radius); let dv = stbVelo / 3; let zRatio = 0.2*0; this.planets.push (Planet.neww (this) .setup ( mass, x, y, radius * rdFloat (-zRatio, zRatio), stbVelo * Math.cos (alpha - 90*deg) + rdFloat (-dv, dv), stbVelo * Math.sin (alpha - 90*deg) + rdFloat (-dv, dv), stbVelo * rdFloat (-zRatio, zRatio) )); } switch (this.moveRateType) { case 0: case 2: this.moveTime = this.scaleTime / this.moveRateParam; break; case 1: this.moveTime = this.scaleTime / 60; } this.zoomReset (); doc ("back") .style.display = ""; doc ("pause") .innerHTML = "Pause"; if (this.firstSetup) { this.firstSetup = false; } this.proc = window.requestAnimationFrame ((t) => {st.update (t);}); } backgroundStar () { let c = this.cvct1; c.fillStyle = "#fff"; let wid = this.canvas1.width, hei = this.canvas1.height, dist = 10; for (let i = (wid / dist - 1) * (hei / dist - 1); i > 0; i--) { c.fillRect (rdInt (0, wid), rdInt (0, hei), 1, 1); } dist = 20; for (let i = (wid / dist - 1) * (hei / dist - 1); i > 0; i--) { c.fillRect (rdInt (0, wid - 1), rdInt (0, hei - 1), 2, 2); } dist = 40; for (let i = (wid / dist - 1) * (hei / dist - 1); i > 0; i--) { c.fillRect (rdInt (0, wid - 2), rdInt (0, hei - 2), 3, 3); } } newAsteroid () { let radius = 990 * this.scaleLength; let alpha = rdFloat (0, 360*deg); let cos = Math.cos (alpha); let sin = Math.sin (alpha); let mass = rdFloat (1e20, 1e22); let stbVelo = GravSys.stdVelocity (3, this.cenPlanets[0].mass, mass, radius); let dv = stbVelo / 3; let zRatio = 0.5*0; this.planets.push (Planet.neww (this) .setup ( mass, radius * cos, radius * sin, radius * rdFloat (-zRatio, zRatio), stbVelo * -cos + rdFloat (-dv, dv), stbVelo * -sin + rdFloat (-dv, dv), stbVelo * rdFloat (-zRatio, zRatio) )); } destroy (i) { let pln = this.planets [i]; if (this.tracking [1] === pln) { this.tracking [0] = false; } pln.del (); this.planets [i] = null; } setOrbit (cenPln, pln, plane, dir, ascend, incline, periArg, peri, eccen, anomaly, epoch, time) { let a = Vec.neww () .eql (dir); a .rotate (plane, ascend); let b = Vec.neww () .eql (a); let c = Vec.neww () .eql (plane .vecMul (a)); c .rotate (a, incline); let d = Vec.neww () .eql (plane); d .rotate (a, incline); b .rotate (d, periArg); } pause () { if (this.pausing) {return;} window.cancelAnimationFrame (this.proc); this.pausing = true; doc ("pause") .innerHTML = "Play"; } play () { if (!this.pausing) {return;} doc ("pause") .innerHTML = "Pause"; this.proc = window.requestAnimationFrame ((t) => {st.update (t);}); } togglePause () { if (this.pausing) { this.play (); } else { this.pause (); } } update (curTime) { let newSec; if (Math.floor (curTime / 1000) > Math.floor (this.vtime)) { newSec = true; } else { newSec = false; } this.vtimeDif = curTime / 1000 - this.vtime; this.vtime += this.vtimeDif; if (this.pausing) { this.pausing = false; this.timeDif = this.scaleTime / 60; } else if (this.vtimeDif <= 0 || this.vtimeDif > 1 / 15) { this.timeDif = this.scaleTime / 60; } else { this.timeDif = this.vtimeDif * this.scaleTime; } this.time += this.timeDif; if (this.moveRateType == 1) { if (1 / this.vtimeDif < 50) { this.moveTime *= 3; } else { this.moveTime *= 0.8; } } let tmDif = this.timeDif; while (tmDif >= this.moveTime) { tmDif -= this.moveTime; this.move (this.moveTime); } if (tmDif > 0) { this.time -= tmDif; } this.inRange (); this.draw (); if (newSec) { this.statStr = "Year: " + (Math.floor (this.time / GravSys.secInYear) + 1970) + ". Planets: " + this.planets.length + " / " + this.totalPlanets + ".
Draw: " + this.drawCount + " (1/s). Move: " + this.moveCount + " (1/s)."; this.moveCount = 0; this.drawCount = 0; } let scdir, scale, ang, x, y; scale = abs (this.cvMat[0], this.cvMat[1]); ang = Math.atan2 (this.cvMat[3], this.cvMat[0]); x = this.cvMat[2] - this.centerX; y = this.centerY - this.cvMat[5]; if (scale >= 1) {scdir = "*";} else {scdir = "/"; scale = 1 / scale;} doc ("stat") .innerHTML = this.statStr + "
Scale: " + scdir + Math.round (scale) + ". Angle: " + Math.round (ang*rad) + ". X: " + Math.round (x) + ". Y: " + Math.round (y) + "."; if (this.tracking [0]) { let pln = this.tracking [1]; doc ("stat") .innerHTML += "
ID: " + pln.id + ". Dist.: " + floatStr (Math.round (abs (pln.rx, pln.ry, pln.rz) / 1000), 3) + " (km). Vel.: " + Math.round (abs (pln.vx, pln.vy, pln.vz)) + " (m/s). Temp.: " + Math.round (Math.sqrt (Math.sqrt (pln.temp4)) - 273.15) + " (^oC)."; } this.proc = window.requestAnimationFrame ((t) => {st.update (t);}); } move (timeDif) { let p = this.planets; for (let i = p.length - 1; i >= 0; i--) { let pln = p [i]; pln.nextPosVel (timeDif); pln.ax = 0; pln.ay = 0; pln.az = 0; } if (this.mergeType != 0) { for (let i = p.length - 1; i >= 1; i--) { let plni = p [i]; for (let j = i - 1; j >= 0; j--) { let plnj = p [j]; let dx, dy, dz, minDist; dx = plni.rx - plnj.rx; dy = plni.ry - plnj.ry; dz = plni.rz - plnj.rz; switch (this.mergeType) { case 1: minDist = (plni.vRadius + plnj.vRadius) * this.scaleLength / this.mergeParam; break; case 2: minDist = this.mergeParam; break; case 3: minDist = this.mergeParam * Math.max (plni.radius, plnj.radius); } if (dx * dx + dy * dy + dz * dz < minDist * minDist) { if (plni.mass > plnj.mass || plni.mass == plnj.mass && plni.id > plnj.id) { plni.merge (plnj); this.destroy (j); p [i] = null; p [j] = plni; } else { plnj.merge (plni); this.destroy (i); } break; } }} clrNull (p); } this.moveRateStat = Infinity; for (let i = p.length - 1; i >= 1; i--) { let plni = p [i]; for (let j = i - 1; j >= 0; j--) { let plnj = p [j]; let dx, dy, dz, dist2, dist, grav, gravi, gravj; dx = plni.rx - plnj.rx; dy = plni.ry - plnj.ry; dz = plni.rz - plnj.rz; dist2 = dx * dx + dy * dy + dz * dz; dist = Math.sqrt (dist2); grav = 6.67e-11 / (dist2 * dist); gravi = grav * plni.mass; gravj = grav * plnj.mass; plni.ax += gravj * -dx; plni.ay += gravj * -dy; plni.az += gravj * -dz; plnj.ax += gravi * dx; plnj.ay += gravi * dy; plnj.az += gravi * dz; let t; if (this.moveRateType == 2) { t = dist / (abs (plni.vx - plnj.vx, plni.vy - plnj.vy, plni.vz - plnj.vz) * timeDif); if (t < this.moveRateStat) { this.moveRateStat = t; }} }} let astAm = this.newAstAm * timeDif / this.scaleTime; while (astAm >= 1) { this.newAsteroid (); astAm--; } if (astAm > 0 && rdFloat (0, 1) < astAm) { this.newAsteroid (); } if (this.moveRateType == 2) { if (this.moveRateStat < this.moveRateParam1) { this.moveTime = this.scaleTime * this.moveRateStat / (this.moveRateParam * this.moveRateParam1); } else { this.moveTime = this.scaleTime / this.moveRateParam; } } this.moveCount++; } inRange () { let cenrx = 0, cenry = 0, cenrz = 0, cenvx = 0, cenvy = 0, cenvz = 0, cenRatioSum = 0, cen = this.cenPlanets; for (let i = cen.length - 1; i >= 0; i--) { let ratio; switch (this.cenType) { case 0: ratio = cen [i].mass; break; case 1: ratio = this.cenRatios[i]; } cenrx += cen [i].rx * ratio; cenry += cen [i].ry * ratio; cenrz += cen [i].rz * ratio; cenvx += cen [i].vx * ratio; cenvy += cen [i].vy * ratio; cenvz += cen [i].vz * ratio; cenRatioSum += ratio; } if (cenRatioSum != 0) { cenrx /= cenRatioSum; cenry /= cenRatioSum; cenrz /= cenRatioSum; cenvx /= cenRatioSum; cenvy /= cenRatioSum; cenvz /= cenRatioSum; } for (let i = this.planets.length - 1; i >= 0; i--) { let pln = this.planets [i]; pln.rx -= cenrx; pln.ry -= cenry; pln.rz -= cenrz; pln.vx -= cenvx; pln.vy -= cenvy; pln.vz -= cenvz; } for (let i = this.planets.length - 1; i >= 0; i--) { let pln = this.planets [i]; if (abs (pln.rx, pln.ry, pln.rz) > this.radius) { this.destroy (i); } } clrNull (this.planets); } shineToDraw () { let plns = this.planets; for (let i = plns.length - 1; i >= 0; i--) { let pln = plns [i]; if (!pln.willDraw || pln.type != 1) {continue;} pln.clearShine (); } for (let i = plns.length - 1; i >= 0; i--) { let plni = plns [i]; if (!plni.willDraw || plni.type !== 1) {continue;} for (let j = plns.length - 1; j >= 0; j--) { let plnj = plns [j]; if (i === j || plnj.type !== 0) {continue;} plni.addShine (plnj); }} for (let i = plns.length - 1; i >= 0; i--) { let pln = plns [i]; if (!pln.willDraw || pln.type != 1) {continue;} pln.calShine (); } } draw () { let c = this.cvct; let p = this.planets; for (let i = p.length - 1; i >= 0; i--) { let pln = p [i]; pln.willDraw = false; pln.vrx = pln.rx / this.scaleLength; pln.vry = pln.ry / this.scaleLength; pln.vrz = pln.rz / this.scaleLength; } if (this.tracking [0]) { this.track (); } for (let i = p.length - 1; i > 0; i--) { let plni = p [i]; for (let j = i - 1; j >= 0; j--) { let plnj = p [j]; if (plni.vrz > plnj.vrz || plni.vrz == plnj.vrz && plni.mass < plnj.mass) { p [i] = plnj; p [j] = plni; plni = plnj; }}} for (let i = p.length - 1, j = this.drawParam1 - 1; i >= 0 && j >= 0; i--, j--) { let pln = p [i]; let mat, x, y, r; mat = this.cvMat; x = mat[0]*pln.vrx + mat[1]*pln.vry + mat[2]; y = mat[3]*pln.vrx + mat[4]*pln.vry + mat[5]; r = abs (mat[0], mat[1]) * pln.vRadius; if (x + r > 0 && x - r < this.canvas.width && y + r > 0 && y - r < this.canvas.height) { pln.willDraw = true; } else { pln.willDraw = false; j++; } } c.save (); c.setTransform (1, 0, 0, 1, 0, 0); c.clearRect (0, 0, this.canvas.width, this.canvas.height); c.restore (); c.save (); c.lineWidth = 100; c.strokeStyle = "#fff"; c.beginPath (); c.arc (0, 0, 1050, 0, 360*deg); c.stroke (); let cos = (6 + Math.sqrt (255)) / 22; let sin = Math.sqrt (1 - cos * cos); c.lineWidth = 1; c.fillStyle = "#fff"; c.beginPath (); c.moveTo (0, -1200); c.lineTo (1100 * sin, -1100 * cos); c.lineTo (-1100 * sin, -1100 * cos); c.closePath (); c.fill (); c.restore (); if (this.drawOrbit) { let cenPln = this.cenPlanets [0]; for (let i = p.length - 1; i >= 0; i--) { let pln = p [i]; if (!pln.willDraw) {continue;} pln.orbitPln [0] = null; } for (let i = p.length - 1; i >= 1; i--) { let plni = p [i]; if (!plni.willDraw) {continue;} for (let j = i - 1; j >= 0; j--) { let plnj = p [j]; if (!plnj.willDraw) {continue;} let r = Vec.neww (plni.rx - plnj.rx, -(plni.ry - plnj.ry), plni.rz - plnj.rz); let v = Vec.neww (plni.vx - plnj.vx, -(plni.vy - plnj.vy), plni.vz - plnj.vz); let L = Vec.neww (); L .eql (r .vecMul (v)); let L2 = L.rad (); if (L2 == 0) { continue; } L2 *= L2; let rr = r.rad (); let A = 6.67e-11 * (plni.mass + plnj.mass); let E = v.rad (); E = E * E / 2 - A / rr; let p = L2 / A; let a = -A / (2 * E); let ELA = -p / a; let sa, sb; if (ELA < 0) { sa = ELA / p; sb = 1; } else { sa = ELA; sb = p; } if (plni.orbitPln [0] === null || sa < plni.orbitPln [1] || sa == plni.orbitPln [1] && sb < plni.orbitPln [2]) { plni.orbitPln [0] = plnj; plni.orbitPln [1] = sa; plni.orbitPln [2] = sb; } if (plnj.orbitPln [0] === null || sa < plnj.orbitPln [1] || sa == plnj.orbitPln [1] && sb < plnj.orbitPln [2]) { plnj.orbitPln [0] = plni; plnj.orbitPln [1] = sa; plnj.orbitPln [2] = sb; } }} for (let i = p.length - 1; i >= 0; i--) { let pln = p [i]; if (!pln.willDraw) {continue;} pln.drawOrbit (c); }} if (this.drawType == 1 || this.drawType == 2) { this.shineToDraw (); } for (let i = p.length - 1; i >= 0; i--) { let pln = p [i]; if (!pln.willDraw) {continue;} if (this.drawType == 0) { pln.drawSimple (c); } else { switch (pln.type) { case 0: pln.drawFusion (c); break; case 1: pln.drawRock (c); }} } if (this.tracking [0]) { this.drawTrack (); } this.drawCount++; } static trfmata = [0, 0, 0, 0, 0, 0]; static trfmatb = [0, 0, 0, 0, 0, 0]; static setTrfMat (mat, a, b, c, d, e, f) { mat [0] = a; mat [1] = b; mat [2] = c; mat [3] = d; mat [4] = e; mat [5] = f; } static setTrfMat1 (mat, mat1) { GravSys.setTrfMat (mat, mat1 [0], mat1 [1], mat1 [2], mat1 [3], mat1 [4], mat1 [5]); } static mulTrfMat (mat, mat1, mat2) { mat [0] = mat1 [0] * mat2 [0] + mat1 [1] * mat2 [3]; mat [1] = mat1 [0] * mat2 [1] + mat1 [1] * mat2 [4]; mat [2] = mat1 [0] * mat2 [2] + mat1 [1] * mat2 [5] + mat1 [2]; mat [3] = mat1 [3] * mat2 [0] + mat1 [4] * mat2 [3]; mat [4] = mat1 [3] * mat2 [1] + mat1 [4] * mat2 [4]; mat [5] = mat1 [3] * mat2 [2] + mat1 [4] * mat2 [5] + mat1 [5]; } sealTrfMat () { let m = this.cvMat; this.cvct.setTransform (m [0], m [3], m [1], m [4], m [2], m [5]); } setTransform (a, b, c, d, e, f) { GravSys.setTrfMat (this.cvMat, a, b, c, d, e, f) ; this.sealTrfMat (); } transform (a, b, c, d, e, f) { GravSys.setTrfMat (GravSys.trfmata, a, b, c, d, e, f); GravSys.mulTrfMat (GravSys.trfmatb, GravSys.trfmata, this.cvMat); GravSys.setTrfMat1 (this.cvMat, GravSys.trfmatb); this.sealTrfMat (); } zoomReset () { this.setTransform (1, 0, this.centerX, 0, 1, this.centerY); } translate (x, y) { this.transform (1, 0, x, 0, 1, y); } rotate0 (a) { let alpha = a*deg, cos = Math.cos (alpha), sin = Math.sin (alpha); this.transform (cos, -sin, 0, sin, cos, 0); } scale0 (k) { this.transform (k, 0, 0, 0, k, 0); } rotateCenter (a) { this.translate (-this.centerX, -this.centerY); this.rotate0 (a); this.translate (this.centerX, this.centerY); } scaleCenter (k) { this.translate (-this.centerX, -this.centerY); this.scale0 (k); this.translate (this.centerX, this.centerY); } recTransform (a, b, c, d, e, f) { GravSys.setTrfMat (GravSys.trfmata, a, b, c, d, e, f); GravSys.mulTrfMat (GravSys.trfmatb, this.cvMat, GravSys.trfmata); GravSys.setTrfMat1 (this.cvMat, GravSys.trfmatb); } recTranslate (x, y) { this.recTransform (1, 0, x, 0, 1, y); } recRotate (a) { let alpha = a*deg, cos = Math.cos (alpha), sin = Math.sin (alpha); this.recTransform (cos, -sin, 0, sin, cos, 0); } recScale (k) { this.recTransform (k, 0, 0, 0, k, 0); } track () { let mat = this.cvMat, tr = this.tracking, pln, plnx, plny; pln = tr [1]; plnx = mat [0] * pln.vrx + mat [1] * pln.vry + mat [2]; plny = mat [3] * pln.vrx + mat [4] * pln.vry + mat [5]; this.translate (tr [2] - plnx, tr [3] - plny); } drawTrack () { let c = this.cvct; c.save (); c.setTransform (1, 0, 0, 1, this.tracking [2], this.tracking [3]); c.lineWidth = 2; c.strokeStyle = "#fff"; let side = 50; c.beginPath (); c.moveTo (side / 2, side / 6); c.lineTo (side / 2, side / 2); c.lineTo (side / 6, side / 2); c.moveTo (-side / 6, side / 2); c.lineTo (-side / 2, side / 2); c.lineTo (-side / 2, side / 6); c.moveTo (-side / 2, -side / 6); c.lineTo (-side / 2, -side / 2); c.lineTo (-side / 6, -side / 2); c.moveTo (side / 6, -side / 2); c.lineTo (side / 2, -side / 2); c.lineTo (side / 2, -side / 6); c.stroke (); c.restore (); } trackClosest (x, y) { let mat = this.cvMat, a = mat [0], b = mat [1], c = x - mat [2], d = mat [3], e = mat [4], f = y - mat [5]; let det = a*e - d*b; let detA = c*e - f*b; let detB = a*f - d*c; x = detA / det; y = detB / det; let mind2 = Infinity, index = 0; for (let i = this.planets.length - 1; i >= 0; i--) { let pln = this.planets [i]; if (!pln.willDraw) {continue;} let t = square (pln.vrx - x, pln.vry - y); if (t < mind2) {mind2 = t; index = i;} } let tr = this.tracking; tr [0] = true; tr [1] = this.planets [index]; let pln = tr [1]; tr [2] = a * pln.vrx + b * pln.vry + mat [2]; tr [3] = d * pln.vrx + e * pln.vry + mat [5]; } trackRandom () { let tr = this.tracking, plns = this.planets; tr [0] = true; tr [1] = plns [rdInt (0, plns.length)]; tr [2] = this.centerX; tr [3] = this.centerY; } } class Control { constructor (gravSys) { this.gravSys = gravSys; let cv = this.gravSys.canvas; nodeWhen (cv, "touchstart", (ev) => { ct.touchStart (ev); }); nodeWhen (cv, "touchmove", (ev) => { ct.touchMove (ev); }); nodeWhen (cv, "touchend", (ev) => { ct.touchEnd (ev); }); nodeWhen (cv, "touchcancel", (ev) => { ct.touchCancel (ev); }); this.zoomLocked = false; this.zoomState = 0; this.zoomPoint0id = 0; this.zoomPoint0x = 0; this.zoomPoint0y = 0; this.zoomPoint1id = 0; this.zoomPoint1x = 0; this.zoomPoint1y = 0; this.trackTime = -Infinity; this.trackX = 0; this.trackY = 0; } touchStart (ev) { if (!this.zoomLocked) { this.zoomPointUpdate (ev); } this.track (ev); } touchMove (ev) { if (!this.zoomLocked && this.zoomState === 1) { this.zoom (ev); } if (this.gravSys.tracking [0]) { this.updateTrack (); } } touchEnd (ev) { if (!this.zoomLocked) { this.zoomPointUpdate (ev); } } touchCancel (ev) { if (!this.zoomLocked) { this.zoomPointUpdate (ev); } } zoomLock () { this.zoomLocked = !this.zoomLocked; this.zoomState = 0; if (this.zoomLocked) { doc ("zoomLock") .innerHTML = "Zoom Unlock"; } else { doc ("zoomLock") .innerHTML = "Zoom Lock"; } } zoomPointUpdate (ev) { let tch = ev.targetTouches; switch (tch.length) { case 2: this.zoomPoint1id = tch [1].identifier; this.zoomPoint1x = tch [1].pageX; this.zoomPoint1y = tch [1].pageY; case 1: this.zoomPoint0id = tch [0].identifier; this.zoomPoint0x = tch [0].pageX; this.zoomPoint0y = tch [0].pageY; } this.zoomState = 1; } zoom (ev) { let tch = ev.targetTouches; switch (tch.length) { case 1: this.zoom1point (tch [0].pageX, tch [0].pageY); break; case 2: if (tch [0].identifier == this.zoomPoint0id) { this.zoom2point (tch [0].pageX, tch [0].pageY, tch [1].pageX, tch [1].pageY); } else { this.zoom2point (tch [1].pageX, tch [1].pageY, tch [0].pageX, tch [0].pageY); }} } zoom1point (x, y) { let s = this.gravSys; if (x == this.zoomPoint0x && y == this.zoomPoint0y) {return;} s.translate (x - this.zoomPoint0x, y - this.zoomPoint0y); this.zoomPoint0x = x; this.zoomPoint0y = y; if (s.pausing) {s.draw ();} } zoom2point (x0, y0, x1, y1) { let s = this.gravSys; if (x0 == this.zoomPoint0x && y0 == this.zoomPoint0y && x1 == this.zoomPoint1x && y1 == this.zoomPoint1y) {return;} let x, y, a, k; a = Math.atan2 (y1 - y0, x1 - x0) - Math.atan2 (this.zoomPoint1y - this.zoomPoint0y, this.zoomPoint1x - this.zoomPoint0x); s.rotate0 (a * rad); let cos = Math.cos (a), sin = Math.sin (a); x = this.zoomPoint0x * cos - this.zoomPoint0y * sin; y = this.zoomPoint0x * sin + this.zoomPoint0y * cos; k = abs (x1 - x0, y1 - y0) / abs (this.zoomPoint1x - this.zoomPoint0x, this.zoomPoint1y - this.zoomPoint0y); s.scale0 (k); x *= k; y *= k; x = x0 - x; y = y0 - y; s.translate (x, y); this.zoomPoint0x = x0; this.zoomPoint0y = y0; this.zoomPoint1x = x1; this.zoomPoint1y = y1; if (s.pausing) {s.draw ();} } track (ev) { let trTime = ev.timeStamp / 1000 - this.trackTime; if (ev.targetTouches.length === 1) { if (trTime > 0 && trTime <= 0.2) { this.gravSys.trackClosest ((this.trackX + ev.targetTouches [0].pageX) / 2, (this.trackY + ev.targetTouches [0].pageY) / 2); this.trackTime = -Infinity; } else { this.trackTime = ev.timeStamp / 1000; this.trackX = ev.targetTouches [0].pageX; this.trackY = ev.targetTouches [0].pageY; }} } updateTrack () { let s = this.gravSys, mat = s.cvMat, tr = s.tracking, pln = tr [1]; tr [2] = mat [0] * pln.vrx + mat [1] * pln.vry + mat [2]; tr [3] = mat [3] * pln.vrx + mat [4] * pln.vry + mat [5]; } } let st, ct, stPausingWhenMenuOn = false; window.onload = () => { st = new GravSys (doc ("maincv"), doc ("bgcv"), 360, 600); ct = new Control (st); docWhen ("mode1", "click", () => { doc ("menu").style.display = "none"; st.test = [false]; st.setup (true, 30, 100, 0.5, 1, 5, 0, 300, 0, 2, false, 30, true, 0); }); docWhen ("mode2", "click", () => { doc ("menu").style.display = "none"; st.test = [false]; st.setup (false, 10, 30, 1, 1, 1, 2, 500, 10, 1, true, 20, false, 1); }); docWhen ("mode3", "click", () => { doc ("menu").style.display = "none"; st.test = [false]; st.setup (true, 10, 50, 3, 3, 100, 0, 100, 0, 1, true, 30, false, 0); }); docWhen ("mode4", "click", () => { doc ("menu").style.display = "none"; st.test = [false]; st.setup (false, 10, 30, 0, 3, 10, 2, 1000, 10000, 0, true, 30, false, 0); st.scaleTime = 3e6; }); docWhen ("test", "click", () => { doc ("menu").style.display = "none"; st.test = [true]; st.setup (false, 2, 6, 0, 1, 1, 2, 500, 10, 1, true, 20, false, 1); /*st.setUp ( true, // has second star 30, // planets amount 30, // asteroids amount 3, // new asteroids per second 1, // mode of merging, 0 to not merge, 1 to merge when images collide, 2 to merge when distance is less than a specified distance, 3 to merge when distance is less than a specified times of radiuses 0, // parameter for merging 0, // mode of raising move rate, 0 to not raise, 1 to raise until dropping of draw rate, 2 to raise when distance is less than a specified times of moving distances 100, // default move rate 100, // parameter for raising move rate 2, // mode of drawing, 0 to draw planets as a simple circle, 1 and 2 to draw with light shine from stars true, // draw or not the orbit 30, // number of planets to draw true, // draw or not the background star 0, // mode of choosing the center point );*/ }); docWhen ("back", "click", () => { doc ("menu").style.display = "none"; if (!stPausingWhenMenuOn) { st.play (); } }); doc ("back").style.display = "none"; docWhen ("menuOn", "click", () => { stPausingWhenMenuOn = st.pausing; if (!stPausingWhenMenuOn) { st.pause (); } doc ("menu").style.display = ""; }); docWhen ("pause", "click", () => { st.togglePause (); }); docWhen ("zoomLock", "click", () => { ct.zoomLock (); }); docWhen ("zoomReset", "click", () => { st.zoomReset (); if (st.pausing) {st.draw ();} }); docWhen ("trackStop", "click", () => { st.tracking [0] = false; }); docWhen ("trackRandom", "click", () => { st.trackRandom (); }); };

Planets System