Improvements to dropsim stuff.

Author: Nishimura-Katsuo

analyze-simulations.js

@@ -0,0 +1,192 @@
+const fs = require('fs');
+
+const ensureDir = (dir) => {
+  if (!fs.existsSync(dir)) {
+    fs.mkdirSync(dir, { recursive: true });
+  }
+};
+
+const safeName = (value) => value.toString().replace(/[^A-Za-z0-9._ \-\[\]\(\)]+/g, '_');
+
+const mean = (values) => values.reduce((sum, v) => sum + v, 0) / (values.length || 1);
+
+const stddev = (values, avg) => {
+  if (values.length < 2) {
+    return 0;
+  }
+  const variance = values.reduce((sum, v) => sum + Math.pow(v - avg, 2), 0) / (values.length - 1);
+  return Math.sqrt(variance);
+};
+
+const normalSample = (avg, deviation) => {
+  if (deviation <= 0) {
+    return avg;
+  }
+  let u = 0;
+  let v = 0;
+  while (u === 0) {
+    u = Math.random();
+  }
+  while (v === 0) {
+    v = Math.random();
+  }
+  const z = Math.sqrt(-2.0 * Math.log(u)) * Math.cos(2.0 * Math.PI * v);
+  return avg + deviation * z;
+};
+
+const buildDistribution = (values, samplesCount, binCount) => {
+  const avg = mean(values);
+  const deviation = stddev(values, avg);
+  
+  const spread = deviation > 0 ? deviation * 4 : Math.max(1, Math.abs(avg) * 0.1);
+  const minX = avg - spread;
+  const maxX = avg + spread;
+  const bins = new Array(binCount).fill(0);
+  const step = (maxX - minX) / binCount;
+
+  for (let i = 0; i < samplesCount; i++) {
+    const sample = normalSample(avg, deviation);
+    const clamped = Math.max(minX, Math.min(maxX, sample));
+    const idx = Math.min(binCount - 1, Math.max(0, Math.floor((clamped - minX) / step)));
+    bins[idx] += 1;
+  }
+
+  const distribution = bins.map((count, index) => {
+    const x = minX + step * (index + 0.5);
+    const y = count / samplesCount;
+    return { x, y };
+  });
+
+  const maxY = distribution.reduce((maxValue, point) => Math.max(maxValue, point.y), 0);
+
+  return {
+    avg,
+    deviation,
+    minX,
+    maxX,
+    minY: 0,
+    maxY,
+    distribution
+  };
+};
+
+const formatNumber = (value) => {
+  if (!Number.isFinite(value)) {
+    return '0';
+  }
+
+  let digits = Math.log10(Math.abs(value));
+
+  if (digits < -2) {
+    return value.toExponential(2);
+  }
+  return value.toFixed(2);
+};
+
+const datapoints = {};
+const maxseq = {};
+const runs = {};
+const totalsdir = __dirname + '/simulation-totals/';
+const graphsdir = __dirname + '/simulation-graphs/';
+
+ensureDir(totalsdir);
+ensureDir(graphsdir);
+
+for (let dir of [totalsdir, graphsdir]) {
+  for (let file of fs.readdirSync(dir)) {
+    if (file.endsWith('.json') || file.endsWith('.svg')) {
+        fs.unlinkSync(dir + file);
+    }
+  }
+}
+
+for (let file of fs.readdirSync(__dirname + '/simulations/')) {
+  if (file.endsWith('.json')) {
+    const data = JSON.parse(fs.readFileSync(__dirname + '/simulations/' + file));
+
+    if (data.runs > 0) {
+      datapoints[data.tc] = datapoints[data.tc] || {};
+      datapoints[data.tc][data.playermod] = datapoints[data.tc][data.playermod] || {};
+      maxseq[data.tc] = maxseq[data.tc] || {};
+      maxseq[data.tc][data.playermod] = maxseq[data.tc][data.playermod] || 0;
+      maxseq[data.tc][data.playermod] = Math.max(maxseq[data.tc][data.playermod], data.seq);
+
+      runs[data.tc] = runs[data.tc] || {};
+      runs[data.tc][data.playermod] = runs[data.tc][data.playermod] || 0;
+      runs[data.tc][data.playermod] += data.runs;
+  
+      for (let item of data.drops) {
+        let [itc, count, magic, rare, set, unique] = item;
+        datapoints[data.tc][data.playermod][itc] = datapoints[data.tc][data.playermod][itc] || [];
+        datapoints[data.tc][data.playermod][itc][data.seq] = count / data.runs;
+      }
+    }
+  }
+}
+
+for (let tc in maxseq) {
+  for (let playermod in maxseq[tc]) {
+    for (let itc in datapoints[tc][playermod]) {
+      let data = datapoints[tc][playermod][itc].map(v => v || 0), runcount = runs[tc][playermod];
+
+      const samplesCount = Math.max(5000, data.length * 200);
+      const binCount = 60;
+      const dist = buildDistribution(data, samplesCount, binCount);
+      const safeBase = safeName(`${tc} [${playermod}][${itc}]`);
+      const totalsPath = totalsdir + safeBase + '.json';
+      const graphPath = graphsdir + safeBase + '.svg';
+
+      const totalsPayload = {
+        tc,
+        playermod: Number(playermod),
+        itc,
+        runs: runcount,
+        samples: data.length,
+        mean: dist.avg,
+        stddev: dist.deviation,
+        minX: dist.minX,
+        maxX: dist.maxX,
+        minY: dist.minY,
+        maxY: dist.maxY,
+        distribution: dist.distribution
+      };
+
+      fs.writeFileSync(totalsPath, JSON.stringify(totalsPayload, null, 2));
+
+      const plotWidth = 700;
+      const plotHeight = 350;
+      const gutterTop = 40;
+      const gutterBottom = 30;
+      const gutterX = 80;
+      const width = gutterX + plotWidth + gutterX;
+      const height = gutterTop + plotHeight + gutterBottom;
+      const boxX = -gutterX;
+      const boxY = -gutterTop;
+      const rangeX = dist.maxX - dist.minX || 1;
+      const rangeY = dist.maxY - dist.minY || 1;
+
+      const points = dist.distribution.map((point, index) => {
+        const x = ((point.x - dist.minX) / rangeX) * plotWidth;
+        const y = plotHeight - ((point.y - dist.minY) / rangeY) * plotHeight;
+        return `${index === 0 ? 'M' : 'L'} ${x.toFixed(2)} ${y.toFixed(2)}`;
+      }).join(' ');
+
+      const svg = [
+        `<svg xmlns="http://www.w3.org/2000/svg" width="${width}" height="${height}" viewBox="${boxX} ${boxY} ${width} ${height}">`,
+        `<rect x="${boxX}" y="${boxY}" width="${width}" height="${height}" fill="#11161c"/>`,
+        `<rect x="0" y="0" width="${plotWidth}" height="${plotHeight}" fill="#0b0f14" stroke="#2a323d" stroke-width="2"/>`,
+        `<line x1="0" y1="${plotHeight}" x2="${plotWidth}" y2="${plotHeight}" stroke="#2a323d" stroke-width="2"/>`,
+        `<line x1="0" y1="0" x2="0" y2="${plotHeight}" stroke="#2a323d" stroke-width="2"/>`,
+        `<path d="${points}" fill="none" stroke="#8ad1ff" stroke-width="2"/>`,
+        `<text x="${(width) / 2 - gutterX}" y="${-gutterTop / 2 + 5}" fill="#e6f0ff" font-size="18" text-anchor="middle">${tc} | players ${playermod * 2 - 1} | ${itc} | ${formatNumber(dist.avg)} ± ${formatNumber(dist.deviation)}</text>`,
+        `<text x="0" y="${plotHeight + gutterBottom / 2 + 5}" fill="#d7e3f2" font-size="14" text-anchor="start">${formatNumber(dist.minX)}</text>`,
+        `<text x="${plotWidth}" y="${plotHeight + gutterBottom / 2 + 5}" fill="#d7e3f2" font-size="14" text-anchor="end">${formatNumber(dist.maxX)}</text>`,
+        `<text x="-10" y="${plotHeight + 4}" fill="#d7e3f2" font-size="14" text-anchor="end">${formatNumber(dist.minY)}</text>`,
+        `<text x="-10" y="4" fill="#d7e3f2" font-size="14" text-anchor="end">${formatNumber(dist.maxY)}</text>`,
+        `</svg>`
+      ].join('');
+
+      fs.writeFileSync(graphPath, svg);
+    }
+  }
+}

build.sh

@@ -1,2 +1,3 @@
 #!/bin/bash
-g++ dropsim.cpp -o dropsim -std=c++17 -O3
+g++ dropsim.cpp -o dropsim -std=c++17 -O3
+g++ dropsim.cpp -o dropsimbase -D"BASETC=1" -std=c++17 -O3