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 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 () { if (this.on) { this.on = false; Vec.obj.push (this); } } static neww (x = 0, y = 0, z = 0, temp = false) { let obj; if (Vec.obj.length > 0) { obj = Vec.obj.pop () .setup (x, y, z, temp); } else { obj = new Vec (x, y, z, temp); } return obj; } constructor (x = 0, y = 0, z = 0, temp = false) { this.setup (x, y, z, temp); } setup (x = 0, y = 0, z = 0, temp = false) { this.on = true; this.temp = temp; this.change (x, y, z); return this; } change (x = 0, y = 0, z = 0) { this.x = x; this.y = y; this.z = z; return this; } rad () { if (this.temp) {this.del ();} return Math.sqrt (this.x * this.x + this.y * this.y + this.z * this.z); } anx () { if (this.temp) {this.del ();} return Math.atan2 (this.y, this.x); } anz () { if (this.temp) {this.del ();} return Math.acos (this.z / this.rad ()); } unit () { let r = this.rad (); let res = Vec.neww (1); if (r != 0) { res .eql (this .div (r)); } if (this.temp) {this.del ();} return res; } eql (a) { this.change (a.x, a.y, a.z); if (a.temp) {a.del ();} return this; } op2 (a, x, y, z) { let res = Vec.neww (x, y, z, true); if (this.temp) {this.del ();} if (a.temp) {a.del ();} return res; } addSub (a, op) { return this.op2 (a, this.x + op * a.x, this.y + op * a.y, this.z + op * a.z); } add (a) { return this.addSub (a, 1); } sub (a) { return this.addSub (a, -1); } mulDiv (k, op) { if (op == -1) {k = 1 / k;} let res = Vec.neww (this.x * k, this.y * k, this.z * k, true); if (this.temp) {this.del ();} return res; } mul (k) { return this.mulDiv (k, 1); } div (k) { return this.mulDiv (k, -1); } addeql (a) { return this .eql (this .add (a)); } subeql (a) { return this .eql (this .sub (a)); } muleql (k) { return this .eql (this .mul (k)); } diveql (k) { return this .eql (this .div (k)); } sclMul (a) { if (this.temp) {this.del ();} if (a.temp) {a.del ();} return this.x * a.x + this.y * a.y + this.z * a.z; } vecMul (a) { return this.op2 (a, this.y * a.z - this.z * a.y, this.z * a.x - this.x * a.z, this.x * a.y - this.y * a.x); } ang (a) { return Math.acos (this.sclMul (a) / (this.rad () * a.rad ())); } rotate (axis, angle) { let az = Vec.neww () .eql (axis .vecMul (this)); if (az.rad () == 0) { az.del (); return this; } let axisu = Vec.neww (); axisu .eql (axis .unit ()); az .eql (axisu .mul (this .sclMul (axisu))); let ax = Vec.neww (); ax .eql (this .sub (az)); let r = ax.rad (); ax .diveql (r); let ay = Vec.neww (); ay .eql (axisu .vecMul (ax)); ax .muleql (r * Math.cos (angle)); ay .muleql (r * Math.sin (angle)); this .eql (ax .add (ay) .add (az)); ax.del (); ay.del (); az.del (); axisu.del (); return this; } } function solve2eq2var (a, b, c, d, e, f) { let det = a * e - d * b; let detx = c * e - f * b; let dety = a * f - d * c; let sol = Vec.neww (); sol.x = detx / det; sol.y = dety / det; return sol; } class TrfMat { static obj = []; del () { if (this.on) { this.on = false; TrfMat.obj.push (this); } } static neww (a = 1, b = 0, c = 0, d = 0, e = 1, f = 0, temp = false) { let obj; if (TrfMat.obj.length > 0) { obj = TrfMat.obj.pop () .setup (a, b, c, d, e, f, temp); } else { obj = new TrfMat (a, b, c, d, e, f, temp); } return obj; } constructor (a = 1, b = 0, c = 0, d = 0, e = 1, f = 0, temp = false) { this.mat = [0, 0, 0, 0, 0, 0]; this.setup (a, b, c, d, e, f, temp); } setup (a = 1, b = 0, c = 0, d = 0, e = 1, f = 0, temp = false) { this.on = true; this.temp = temp; this.change (a, b, c, d, e, f); return this; } change (a = 1, b = 0, c = 0, d = 0, e = 1, f = 0) { let m = this.mat; m[0] = a; m[1] = b; m[2] = c; m[3] = d; m[4] = e; m[5] = f; return this; } eql (a) { let m = a.mat; this.change (m[0], m[1], m[2], m[3], m[4], m[5]); if (a.temp) {a.del ();} return this; } mul (a) { let b = this.mat, c = a.mat; let res = TrfMat.neww ( b[0] * c[0] + b[1] * c[3], b[0] * c[1] + b[1] * c[4], b[0] * c[2] + b[1] * c[5] + b[2], b[3] * c[0] + b[4] * c[3], b[3] * c[1] + b[4] * c[4], b[3] * c[2] + b[4] * c[5] + b[5], true ); if (this.temp) {this.del ();} if (a.temp) {a.del ();} return res; } transformShape (a, b, c, d, e, f) { this .eql (TrfMat.neww (a, b, c, d, e, f, true) .mul (this)); } translateShape (x, y) { this.transformShape (1, 0, x, 0, 1, y); } scaleShape (kx, ky) { this.transformShape (kx, 0, 0, 0, ky, 0); } rotateShape (a) { let cos = Math.cos (a), sin = Math.sin (a); this.transformShape (cos, -sin, 0, sin, cos, 0); } translateUntwist (x, y) { this .eql (this .mul (TrfMat.neww (1, 0, -x, 0, 1, -y, true))); } scaleUntwist (kx, ky) { this .eql (this .mul (TrfMat.neww (1 / kx, 0, 0, 0, 1 / ky, 0, true))); } rotateUntwist (a) { let cos = Math.cos (-a), sin = Math.sin (-a); this .eql (this .mul (TrfMat.neww (cos, -sin, 0, sin, cos, 0, true))); } apply (c) { let m = this.mat; c.setTransform (m[0], m[3], m[1], m[4], m[2], m[5]); } save () { let res = TrfMat.neww (); res .eql (this); return res; } restore (mat, c) { this .eql (mat); this.apply (c); mat.del (); } poi (x, y) { let m = this.mat, p = Vec.neww (); p.x = m[0] * x + m[1] * y + m[2]; p.y = m[3] * x + m[4] * y + m[5]; return p; } revPoi (x, y) { let m = this.mat; let p = solve2eq2var (m[0], m[1], x - m[2], m[3], m[4], y - m[5]); return p; } det () { let m = this.mat; return m[0] * m[4] - m[3] * m[1]; } calFlip () { return Math.sign (this.det ()); } calScale () { return Math.sqrt (Math.abs (this.det ())); } deskew (x, y) { if (this.det () == 0) {return;} let p = this.revPoi (x, y); let k = this.calScale (); let m = this.mat; if (m[3] == 0) { m[3] += m[0]; m[4] += m[1];// skew } let s = (k - m[0]) / m[3]; m[0] += s * m[3]; m[1] += s * m[4];// skew s = -m[3] / m[0]; m[3] += s * m[0]; m[4] += s * m[1];// skew s = -m[1] / m[4]; m[0] += s * m[3]; m[1] += s * m[4];// skew let p1 = this.poi (p.x, p.y); m[2] += p1.x - x; m[5] += p1.y - y; p.del (); p1.del (); } } class Orbit { static obj = []; del () { if (this.on) { this.on = false; Orbit.obj.push (this); } } static neww () { let obj; if (Orbit.obj.length > 0) { obj = Orbit.obj.pop () .setup (); } else { obj = new Orbit (); } return obj; } constructor () { this.mat = TrfMat.neww (); this.r = Vec.neww (); this.v = Vec.neww (); this.L = Vec.neww (); this.setup (); } setup () { this.on = true; this.change (); return this; } change () { this.orbPln = null; this.cenPln = null; this.sizeA = 0; this.sizeB = 0; this.type = 0; this.p = 0; this.e = 0; this.a = 0; this.b = 0; this.c = 0; this.incline = 0; this.mat.change (); this.p2 = 0; this._2p = 0; this.e2_1 = 0; this.pe = 0; this.r.change (); this.v.change (); this.L.change (); this.rr = 0; return this; } update (plni, plnj) { if (plni.mass <= plnj.mass) { this.orbPln = plni; this.cenPln = plnj; } else { this.orbPln = plnj; this.cenPln = plni; } this.r .eql (this.orbPln.pos .sub (this.cenPln.pos)); this.r.y = -this.r.y; this.v .eql (this.orbPln.vel .sub (this.cenPln.vel)); this.v.y = -this.v.y; this.L .eql (this.r .vecMul (this.v)); let Lr = this.L.rad (); this.rr = this.r.rad (); if (Lr == 0) { this.type = 2; this.sizeA = -2 / this.rr; this.sizeB = 1; return; } let vr = this.v.rad (); let A = 6.67e-11 * (this.cenPln.mass + this.orbPln.mass); let E = vr * vr / 2 - A / this.rr; if (E >= 0) { this.type = 1; } else { this.type = 0; } this.a = -A / (2 * E); this.p = Lr * Lr / A; if (this.type == 0) { this.sizeA = -1 / this.a; this.sizeB = 1; } else if (this.type == 1) { this.sizeA = -this.p / this.a; this.sizeB = this.p; } this.update1 (this.orbPln.orbit, this.cenPln, 1); if (this.orbPln.mass == this.cenPln.mass) { this.update1 (this.cenPln.orbit, this.orbPln, -1); } } update1 (orb, pln, s) { if (orb.cenPln === null || this.sizeA < orb.sizeA || this.sizeA == orb.sizeA && this.sizeB < orb.sizeB) { orb.cenPln = pln; orb.sizeA = this.sizeA; orb.sizeB = this.sizeB; orb.type = this.type; orb.p = this.p; orb.a = this.a; orb.r .eql (this.r .mul (s)); orb.v .eql (this.v .mul (s)); orb.L .eql (this.L); orb.rr = this.rr; } } calAll () { this.e2_1 = -this.p / this.a; this.e = Math.sqrt (1 + this.e2_1); if (this.type == 0) { this.c = this.a * this.e; this.b = Math.sqrt (this.a * this.a - this.c * this.c); } else if (this.type == 1) { this.p2 = this.p * this.p; this._2p = 2 * this.p; this.pe = this.p * this.e; } let phi = Math.acos ((this.p / this.rr - 1) / this.e); if (this.r .sclMul (this.v) < 0) { phi = -phi; } this.r .rotate (this.L, -phi); let meter = 1 / st.scaleLength; if (this.type == 0) { this.a *= meter; this.b *= meter; this.c *= meter; } else if (this.type == 1) { this.p *= meter; this.p2 *= meter * meter; this._2p *= meter; this.pe *= meter; } let Oz = Vec.neww (0, 0, 1); let itsl = Vec.neww (); itsl .eql (this.L .vecMul (Oz)); this.mat.change (); this.incline = 1; if (itsl.rad () != 0) { let phi1 = this.L .ang (Oz); phi = Math.atan2 (-itsl.y, itsl.x); this.mat.rotateUntwist (-phi); this.mat.scaleUntwist (1, 1 / Math.cos (phi1)); this.mat.rotateUntwist (phi); this.incline = Math.abs (Math.cos (phi1)); this.L .rotate (itsl, phi1); this.r .rotate (itsl, phi1); } else if (this.L.z < 0) { phi = Math.atan2 (-this.r.y, this.r.x); this.mat.rotateUntwist (-phi); this.mat.scaleUntwist (1, -1); this.mat.rotateUntwist (phi); this.L .rotate (this.r, 180*deg); } phi = Math.atan2 (-this.r.y, this.r.x); this.mat.rotateUntwist (-phi); itsl.del (); Oz.del (); } calX (y) { if (this.e2_1 > 0) { return (this.pe - Math.sqrt (this.e2_1 * y * y + this.p2)) / this.e2_1; } else if (this.e2_1 == 0) { return (this.p2 - y * y) / this._2p; } } calPoi (x, y, dmin, dmax, s) { let res = Vec.neww (); let x0, y0, x1, y1, x2, y2, dx, dy, dr2, dr2min, dr2max; dr2min = dmin * dmin; dr2max = dmax * dmax; x0 = x; y0 = y; y2 = y + s * dmax * 2; x2 = this.calX (y2); while (Math.abs (x2 - x0) + Math.abs (y2 - y0) > 1e-3) { y1 = (y0 + y2) / 2; x1 = this.calX (y1); dx = x1 - x; dy = y1 - y; dr2 = dx * dx + dy * dy; if (dr2 < dr2min) { x0 = x1; y0 = y1; } else if (dr2 > dr2max) { x2 = x1; y2 = y1; } else { res.x = x1; res.y = y1; return res; } } } calMidPoi (x0, y0, x1, y1) { let k = this.e2_1 * x0 - this.pe; let res = solve2eq2var (y0, -y1, -x0 + x1, k, -this.e2_1 * x1 + this.pe, -y0 + y1); let poi = Vec.neww (y0 * res.x + x0, k * res.x + y0); res.del (); return poi; } 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); } static stdVelocity (a, mass0, mass1, radius) { return Math.sqrt (a * 6.67e-11 * (mass0 + mass1) / radius); } } class Planet { static obj = []; del () { if (this.on) { this.on = false; 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 obj; if (Planet.obj.length > 0) { obj = Planet.obj.pop () .setup (system, id, mass, type, radius, rx, ry, rz, vx, vy, vz, vRadius, red, green, blue); } else { obj = new Planet (system, id, mass, type, radius, rx, ry, rz, vx, vy, vz, vRadius, red, green, blue); } return obj; } 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.absorb = Vec.neww (); this.light = Vec.neww (); this.pos = Vec.neww (); this.vel = Vec.neww (); this.acc = Vec.neww (); this.orbit = Orbit.neww (); this.orbit.orbPln = this; this.vpos = Vec.neww (); this.color = Vec.neww (); this.colorRatio = Vec.neww (); this.setup (system, id, mass, type, radius, rx, ry, rz, vx, vy, vz, vRadius, red, green, blue); } // use for general case (without "made up" rules) setup (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.on = true; this.system = system; this.id = id; this.change (mass, type, radius, rx, ry, rz, vx, vy, vz, vRadius, red, green, blue); return this; } change (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.mass = mass; this.type = type; this.radius = radius; this.temp4 = 0; this.absorb.change (); this.light.change (); this.pos.change (rx, ry, rz); this.vel.change (vx, vy, vz); this.acc.change (); this.willDraw = true; this.orbit.change (); this.orbit.orbPln = this; this.vpos.change (); this.vRadius = vRadius; this.color.change (red, green, blue); this.colorRatio.change (); return this; } // use in this specific simulation setup1 (mass, rx, ry, rz, vx, vy, vz) { let s = this.system; this.id = s.totalPlanets; s.totalPlanets++; this.mass = mass; this.pos.change (rx, ry, rz); this.vel.change (vx, vy, vz); this.setupRule (); if (this.type == 1) { this.color.change (rdInt (0, 256), rdInt (0, 256), rdInt (0, 256)); } else { this.color.change (512, 256, 128); } let a; if (s.drawType == 2 && this.type == 1) { a = s.scaleColorRGBa; this.absorb.x = 1 - Math.exp ((this.color.x - 255) / a); this.absorb.y = 1 - Math.exp ((this.color.y - 255) / a); this.absorb.z = 1 - Math.exp ((this.color.z - 255) / a); } if (s.drawType == 1) { this.colorRatio .eql (this.color .div ((this.color.x + this.color.y + this.color.z) / 3)); this.absorb.x = rdFloat (0, 1); } return this; } merge (pln) { let s = this.system; let totMass = this.mass + pln.mass, a = this.mass / totMass, b = pln.mass / totMass; this.mass = totMass; this.setupRule (); if (s.drawType == 2 && this.type == 1) { this.absorb .muleql (a) .addeql (pln.absorb .mul (b)); } this.pos .muleql (a) .addeql (pln.pos .mul (b)); this.vel .muleql (a) .addeql (pln.vel .mul (b)); let sc; if (s.drawType == 2 && this.type == 1) { sc = s.scaleColorRGBa; this.color.x = Math.log (1 - this.absorb.x) * sc + 255; this.color.y = Math.log (1 - this.absorb.y) * sc + 255; this.color.z = Math.log (1 - this.absorb.z) * sc + 255; } else { this.color .muleql (a) .addeql (pln.color .mul (b)); } if (s.drawType == 1) { this.colorRatio .muleql (a) .addeql (pln.colorRatio .mul (b)); } let tr = s.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.light.change (s.calLightRed (temp), s.calLightGreen (temp), s.calLightBlue (temp)); } this.vRadius = s.calVradius (this.mass); } nextPosVel (time) { let dvel = Vec.neww (); dvel .eql (this.acc .mul (time)); this.pos .addeql (this.vel .add (dvel .div (2)) .mul (time)); this.vel .addeql (dvel); dvel.del (); } clearShine () { this.temp4 = 0; this.light.change (); } addShine (pln) { const dist = pln.pos .sub (this.pos) .rad (); const distRatio = pln.radius * pln.radius / (dist * dist); this.temp4 += pln.temp4 * distRatio / 4; switch (this.system.drawType) { case 1: this.light.x += 5.67e-8 * pln.temp4 * distRatio / 2 * (1 - this.absorb.x); break; case 2: this.light.x += pln.light.x * distRatio / 2 * (1 - this.absorb.x); this.light.y += pln.light.y * distRatio / 2 * (1 - this.absorb.y); this.light.z += pln.light.z * distRatio / 2 * (1 - this.absorb.z); } } calShine () { let s = this.system; let temp; switch (s.drawType) { case 1: this.light.x += 5.67e-8 * this.temp4 * this.absorb.x; break; case 2: temp = Math.sqrt (Math.sqrt (this.temp4)); this.light.x += s.calLightRed (temp) * this.absorb.x; this.light.y += s.calLightGreen (temp) * this.absorb.y; this.light.z += s.calLightBlue (temp) * this.absorb.z; } } // standard 2-body gravity orbit drawOrbit (c) { let orb = this.orbit; if (orb.cenPln === null || orb.orbPln === null) {return;} if (orb.type == 2) {return;} let pln = orb.cenPln; let m = this.system.cvMat; let rec = m.save (); c.save (); let color = filterToColor (this.color.x, this.color.y, this.color.z); c.strokeStyle = color; c.fillStyle = color; m.translateUntwist (-pln.vpos.x, -pln.vpos.y); orb.calAll (); c.lineWidth = 1 / Math.sqrt (orb.incline); m .eql (m .mul (orb.mat)); if (orb.type == 0) { m.translateUntwist (orb.c, 0); m.scaleUntwist (1, orb.a / orb.b); m.apply (c); c.beginPath (); c.arc (0, 0, orb.a, 0, 360*deg); c.stroke (); } else if (orb.type == 1) { m.apply (c); let poi, poi1, poi2; poi2 = Vec.neww (orb.calX (0), 0); for (let i = 0; i < 30; i++) { poi = poi2; poi2 = orb.calPoi (poi.x, poi.y, 9.5, 10.5, 1); poi1 = orb.calMidPoi (poi.x, poi.y, poi2.x, poi2.y); c.beginPath (); c.moveTo (poi.x, -poi.y); c.quadraticCurveTo (poi1.x, -poi1.y, poi2.x, -poi2.y); c.stroke (); c.beginPath (); c.moveTo (poi.x, poi.y); c.quadraticCurveTo (poi1.x, poi1.y, poi2.x, poi2.y); c.stroke (); poi.del (); poi1.del (); } poi2.del (); } c.restore (); m.restore (rec, c); } drawSimple (c) { let color = filterToColor (this.color.x, this.color.y, this.color.z); c.save (); c.lineWidth = 0; c.fillStyle = color; c.beginPath (); c.arc (this.vpos.x, this.vpos.y, this.vRadius, 0, 360*deg); c.fill (); c.restore (); } drawFusion (c) { let color, temp, s = this.system; switch (s.drawType) { case 1: color = s.calColorAll (5.67e-8 * this.temp4); color = filterToColor (color * this.colorRatio.x, color * this.colorRatio.y, color * this.colorRatio.z); break; case 2: temp = Math.sqrt (Math.sqrt (this.temp4)); color = s.calColorRGB (s.calLightRed (temp), s.calLightGreen (temp), s.calLightBlue (temp)); } c.save (); c.lineWidth = 0; c.fillStyle = color; c.beginPath (); c.arc (this.vpos.x, this.vpos.y, this.vRadius, 0, 360*deg); c.fill (); c.restore (); } drawRock (c) { let s = this.system, color; switch (s.drawType) { case 1: color = s.calColorAll (this.light.x); color = filterToColor (color * this.colorRatio.x, color * this.colorRatio.y, color * this.colorRatio.z); break; case 2: color = s.calColorRGB (this.light.x, this.light.y, this.light.z); } c.save (); c.lineWidth = 0; c.fillStyle = color; c.beginPath (); c.arc (this.vpos.x, this.vpos.y, this.vRadius, 0, 360*deg); 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 = canvas.getContext ("2d"); this.cvMat = TrfMat.neww (); this.canvas1 = canvas1; this.canvas1.width = width; this.canvas1.height = height; this.cvct1 = 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.drawAmount = 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)); } calVradius (mass) { return Math.pow (mass, this.scaleVradiusA) * this.scaleVradiusB; } scaleVradius (mass0, radius0, mass1, radius1) { let sol = solve2eq2var (Math.log (mass0), 1, Math.log (radius0), Math.log (mass1), 1, Math.log (radius1)); this.scaleVradiusA = sol.x; this.scaleVradiusB = Math.exp (sol.y); sol.del (); } calFusionTemp4 (mass) { return this.fusionTemp4 * Math.cbrt (mass); } // 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; } calColorAll (light) { return Math.log (light) * this.scaleColorAllA + this.scaleColorAllB; } scaleColorAll (light0, color0, light1, color1) { let sol = solve2eq2var (Math.log (light0), 1, color0, Math.log (light1), 1, color1); 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 = solve2eq2var (Math.log (minRed / this.maxLightRed), 1, color0, 0, 1, 255); this.scaleColorRGBa = sol.x; this.scaleColorRGBb = sol.y; sol.del (); } setup (cen2, plnAm, astAm, newAstAm, mergeType, mergeParam, moveRateType, moveRateParam, moveRateParam1, drawType, drawOrbit, drawAmount, bgStar, cenType) { this.mergeType = mergeType; this.mergeParam = mergeParam; this.moveRateType = moveRateType; this.moveRateParam = moveRateParam; this.moveRateParam1 = moveRateParam1; this.drawType = drawType; this.drawOrbit = drawOrbit; this.drawAmount = drawAmount; this.cenType = cenType; this.newAstAm = newAstAm; if (this.firstSetup) { alert ( "[Pause] / [Play] button to pause / play the simulation" + "\n[Zoom Unlock] / [Zoom Lock] button to turn on / off the zooming function" + "\n[Zoom Reset] button to reset the zooming" + "\nDouble tap a point to track 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, 15, 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 mass = 2e30; let pln = Planet.neww (this) .setup1 (mass, 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) .setup1 (mass, radius, 0, 0, 0, -Orbit.stdVelocity (1, this.cenPlanets[0].mass, mass, radius), 0); this.planets.push (pln); this.cenPlanets.push (pln); this.cenRatios.push (0); } this.inRange (); for (let i = 0; i < plnAm; i++) { this.stablePlanet (2e12, 1e22, 1e26, 0.2, 0.1); } for (let i = 0; i < astAm; i++) { this.stablePlanet (2e12, 1e20, 1e22, 0.3, 0.2); } 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); } } stablePlanet (radius, massa, massb, vrat, zRatio) { 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 (massa, massb); let cenMass = 0; for (let pln of this.cenPlanets) { cenMass += pln.mass; } let stbVelo = Orbit.stdVelocity (1, cenMass, mass, radius); let dv = stbVelo * vrat; this.planets.push (Planet.neww (this) .setup1 ( 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), 0 )); } outsideInPlanet (radius, massa, massb, a, b, vrat) { let x, y, z; do { x = rdFloat (-1, 1); y = rdFloat (-1, 1); z = rdFloat (-1, 1); } while (x * x + y * y + z * z > 1); let mass = rdFloat (massa, massb); let cenMass = 0; for (let pln of this.cenPlanets) { cenMass += pln.mass; } let stbVelo = Orbit.stdVelocity (rdFloat (a, b), cenMass, mass, radius); let dv = stbVelo * vrat; this.planets.push (Planet.neww (this) .setup1 ( mass, radius * x, radius * y, radius * z, stbVelo * -x + rdFloat (-dv, dv), stbVelo * -y + rdFloat (-dv, dv), stbVelo * -z + rdFloat (-dv, dv) )); } destroy (i) { let pln = this.planets [i]; if (this.tracking [1] === pln) { this.tracking [0] = false; } pln.del (); this.planets [i] = null; } 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.update1 (); if (newSec) { this.statStr = "
Draw: " + this.drawCount + " (1/s). Move: " + this.moveCount + " (1/s)."; this.moveCount = 0; this.drawCount = 0; } this.update2 (); this.proc = window.requestAnimationFrame ((t) => {st.update (t);}); } updateNow () { this.update1 (); this.update2 (); } update1 () { this.inRange (); this.draw (); } update2 () { doc ("stat") .innerHTML = "Year: " + (Math.floor (this.time / GravSys.secInYear) + 1970) + ". Planets: " + this.planets.length + " / " + this.totalPlanets + "." + this.statStr; let mat = this.cvMat, scdir, scale, ang, x, y; scale = mat.calScale (); ang = Math.atan2 (mat.mat[3], mat.mat[0]); x = mat.mat[2] - this.centerX; y = this.centerY - mat.mat[5]; if (scale >= 1) { scdir = "*"; } else { scdir = "/"; scale = 1 / scale; } doc ("stat") .innerHTML += "
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 (pln.pos.rad () / 1000), 3) + " (km). Vel.: " + Math.round (pln.vel.rad ()) + " (m/s). Temp.: " + Math.round (Math.sqrt (Math.sqrt (pln.temp4)) - 273.15) + " (^oC)."; } } move (timeDif) { let p = this.planets; for (let pln of p) { pln.nextPosVel (timeDif); pln.acc.change (); } 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.pos.x - plnj.pos.x; dy = plni.pos.y - plnj.pos.y; dz = plni.pos.z - plnj.pos.z; 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.pos.x - plnj.pos.x; dy = plni.pos.y - plnj.pos.y; dz = plni.pos.z - plnj.pos.z; 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.acc.x += gravj * -dx; plni.acc.y += gravj * -dy; plni.acc.z += gravj * -dz; plnj.acc.x += gravi * dx; plnj.acc.y += gravi * dy; plnj.acc.z += gravi * dz; let t; if (this.moveRateType == 2) { t = dist / (plni.vel .sub (plnj.vel) .rad () * timeDif); if (t < this.moveRateStat) { this.moveRateStat = t; }} }} let astAm = this.newAstAm * timeDif / this.scaleTime; let frac = astAm - Math.floor (astAm); if (frac > 0 && rdFloat (0, 1) < frac) { astAm = Math.ceil (astAm); } else { astAm = Math.floor (astAm); } for (let i = 0; i < astAm; i++) { this.outsideInPlanet (990 * this.scaleLength, 1e20, 1e22, 3, 10, 0.3); } 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 p = this.planets; let cpos = Vec.neww (), cvel = Vec.neww (), 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]; } cpos .addeql (cen[i].pos .mul (ratio)); cvel .addeql (cen[i].vel .mul (ratio)); cenRatioSum += ratio; } if (cenRatioSum != 0) { cpos .diveql (cenRatioSum); cvel .diveql (cenRatioSum); } for (let pln of p) { pln.pos .subeql (cpos); pln.vel .subeql (cvel); } cpos.del (); cvel.del (); for (let i = p.length - 1; i >= 0; i--) { let pln = p[i]; if (pln.pos.rad () > this.radius) { this.destroy (i); } } clrNull (p); } shineToDraw () { let p = this.planets; for (let pln of p) { if (!pln.willDraw || pln.type != 1) {continue;} pln.clearShine (); } for (let plni of p) { if (!plni.willDraw || plni.type !== 1) {continue;} for (let plnj of p) { if (plni === plnj || plnj.type !== 0) {continue;} plni.addShine (plnj); }} for (let pln of p) { if (!pln.willDraw || pln.type != 1) {continue;} pln.calShine (); } } draw () { let c = this.cvct; let p = this.planets; for (let pln of p) { pln.willDraw = false; pln.vpos .eql (pln.pos .div (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.vpos.z > plnj.vpos.z || plni.vpos.z == plnj.vpos.z && plni.mass < plnj.mass) { p[i] = plnj; p[j] = plni; plni = plnj; }}} for (let i = p.length - 1, j = this.drawAmount - 1; i >= 0 && j >= 0; i--, j--) { let pln = p[i]; let mat, x, y, r; mat = this.cvMat; let poi = mat.poi (pln.vpos.x, pln.vpos.y); x = poi.x; y = poi.y; r = mat.calScale () * pln.vRadius; poi.del (); 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) { for (let pln of p) { pln.orbit.cenPln = null; } let orbit = Orbit.neww (); 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]; orbit.update (plni, plnj); }} orbit.del (); for (let pln of p) { if (!pln.willDraw) {continue;} pln.drawOrbit (c); }} if (this.drawType == 1 || this.drawType == 2) { this.shineToDraw (); } for (let pln of p) { 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++; } zoomReset () { this.cvMat.change (1, 0, this.centerX, 0, 1, this.centerY); this.cvMat.apply (this.cvct); } translate (x, y) { this.cvMat.translateShape (x, y); this.cvMat.apply (this.cvct); } scale (k) { this.cvMat.scaleShape (k, k); this.cvMat.apply (this.cvct); } rotate (a) { this.cvMat.rotateShape (a); this.cvMat.apply (this.cvct); } track () { let tr = this.tracking, pln = tr[1]; let poi = this.cvMat.poi (pln.vpos.x, pln.vpos.y); this.translate (tr[2] - poi.x, tr[3] - poi.y); poi.del (); } 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, s2 = side / 2, s6 = side / 6; c.beginPath (); c.moveTo (s2, s6); c.lineTo (s2, s2); c.lineTo (s6, s2); c.moveTo (-s6, s2); c.lineTo (-s2, s2); c.lineTo (-s2, s6); c.moveTo (-s2, -s6); c.lineTo (-s2, -s2); c.lineTo (-s6, -s2); c.moveTo (s6, -s2); c.lineTo (s2, -s2); c.lineTo (s2, -s6); c.stroke (); c.restore (); } trackClosest (x, y) { let m = this.cvMat; let poi = m.revPoi (x, y); x = poi.x; y = poi.y; poi.del (); 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.vpos.x - x, pln.vpos.y - y); if (t < mind2) {mind2 = t; index = i;} } let tr = this.tracking; tr[0] = true; tr[1] = this.planets [index]; let pln = tr[1]; poi = m.poi (pln.vpos.x, pln.vpos.y); tr[2] = poi.x; tr[3] = poi.y; poi.del (); } trackRandom () { let tr = this.tracking, p = this.planets; tr[0] = true; tr[1] = p [rdInt (0, p.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) { let s = this.gravSys; if (!this.zoomLocked && this.zoomState == 1) { this.zoom (ev); } if (s.tracking [0]) { this.updateTrack (); } if (s.pausing) {s.updateNow ();} } 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; } 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.rotate (a); 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.scale (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; } track (ev) { let s = this.gravSys; let trTime = ev.timeStamp / 1000 - this.trackTime; if (ev.targetTouches.length == 1) { if (trTime > 0 && trTime <= 0.2) { s.trackClosest ((this.trackX + ev.targetTouches [0].pageX) / 2, (this.trackY + ev.targetTouches [0].pageY) / 2); this.trackTime = -Infinity; if (s.pausing) {s.updateNow ();} } else { this.trackTime = ev.timeStamp / 1000; this.trackX = ev.targetTouches [0].pageX; this.trackY = ev.targetTouches [0].pageY; }} } updateTrack () { let s = this.gravSys, m = s.cvMat, tr = s.tracking, pln = tr[1]; let poi = m.poi (pln.vpos.x, pln.vpos.y); tr[2] = poi.x; tr[3] = poi.y; poi.del (); } } 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, 5, 1, 1, 0, 300, 0, 2, false, 30, true, 0); }); docWhen ("mode2", "click", () => { doc ("menu").style.display = "none"; st.test = [false]; st.setup (false, 10, 30, 5, 1, 1, 2, 500, 10, 1, true, 20, false, 1); }); docWhen ("mode3", "click", () => { doc ("menu").style.display = "none"; st.test = [false]; st.setup (false, 0, 0, 2, 3, 100, 0, 1000, 0, 0, true, 20, 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, 0, 0, 3, 1, 1, 2, 500, 10, 1, true, 20, false, 1); // test orbit /*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.updateNow ();} }); docWhen ("trackStop", "click", () => { st.tracking [0] = false; }); docWhen ("trackRandom", "click", () => { st.trackRandom (); }); };

Planets System (3D)