@jjlmoya/utils-science 1.25.0 → 1.27.0

package/src/tool/natural-selection-drift/natural-selection-drift.css

@@ -0,0 +1,429 @@
+.ns-app {
+  --ns-bg: #f7fafc;
+  --ns-text: #0f172a;
+  --ns-muted: rgba(15, 23, 42, 0.62);
+  --ns-panel: rgba(255, 255, 255, 0.34);
+  --ns-border: rgba(255, 255, 255, 0.18);
+  --ns-track: rgba(15, 23, 42, 0.12);
+  --ns-thumb: #fff;
+  --ns-accent: #16a34a;
+
+  padding: 1.25rem;
+  color: var(--ns-text);
+}
+
+.theme-light .ns-app {
+  --ns-bg: #f8fbff;
+  --ns-text: #0f172a;
+  --ns-muted: rgba(15, 23, 42, 0.58);
+  --ns-panel: rgba(255, 255, 255, 0.72);
+  --ns-border: rgba(148, 163, 184, 0.22);
+  --ns-track: rgba(15, 23, 42, 0.14);
+  --ns-thumb: #fff;
+  --ns-accent: #0f766e;
+}
+
+.theme-dark .ns-app {
+  --ns-bg: #020617;
+  --ns-text: #f8fafc;
+  --ns-muted: rgba(255, 255, 255, 0.7);
+  --ns-panel: rgba(255, 255, 255, 0.03);
+  --ns-border: rgba(255, 255, 255, 0.08);
+  --ns-track: rgba(255, 255, 255, 0.1);
+  --ns-thumb: #fff;
+  --ns-accent: #86efac;
+}
+
+.ns-card {
+  position: relative;
+  min-height: min(84vh, 920px);
+  overflow: hidden;
+  padding: clamp(16px, 4vw, 24px);
+  border-radius: 32px;
+  background:
+    radial-gradient(circle at 50% 50%, rgba(34, 197, 94, 0.08), transparent 34%),
+    radial-gradient(circle at 15% 10%, rgba(56, 189, 248, 0.09), transparent 24%),
+    var(--ns-bg);
+  border: 1px solid var(--ns-border);
+  box-shadow: 0 30px 90px rgba(2, 6, 23, 0.18);
+}
+
+.theme-light .ns-card {
+  background:
+    radial-gradient(circle at 50% 50%, rgba(45, 212, 191, 0.16), transparent 32%),
+    radial-gradient(circle at 15% 10%, rgba(59, 130, 246, 0.12), transparent 22%),
+    linear-gradient(180deg, rgba(255, 255, 255, 0.96), rgba(244, 248, 255, 0.94)),
+    var(--ns-bg);
+  box-shadow: 0 28px 80px rgba(148, 163, 184, 0.22);
+}
+
+.ns-canvas {
+  position: absolute;
+  inset: 0;
+  width: 100%;
+  height: 100%;
+  display: block;
+  z-index: 1;
+}
+
+.ns-canvas-badge {
+  position: absolute;
+  top: clamp(16px, 3vw, 24px);
+  right: clamp(16px, 3vw, 24px);
+  z-index: 2;
+  padding: 0.4rem 0.7rem;
+  border-radius: 999px;
+  background: rgba(2, 6, 23, 0.75);
+  border: 1px solid rgba(255, 255, 255, 0.12);
+  color: rgba(255, 255, 255, 0.92);
+  font-size: 11px;
+  font-weight: 700;
+  text-shadow: 0 1px 2px rgba(0, 0, 0, 0.35);
+  pointer-events: none;
+}
+
+.theme-dark .ns-canvas-badge {
+  background: rgba(255, 255, 255, 0.08);
+}
+
+.theme-light .ns-canvas-badge {
+  background: rgba(255, 255, 255, 0.88);
+  border-color: rgba(148, 163, 184, 0.2);
+  color: #0f172a;
+  text-shadow: none;
+  box-shadow: 0 10px 30px rgba(148, 163, 184, 0.16);
+}
+
+.ns-top-hud {
+  position: absolute;
+  left: max(360px, clamp(16px, 3vw, 24px));
+  right: clamp(120px, 18vw, 220px);
+  top: clamp(16px, 3vw, 24px);
+  z-index: 2;
+  display: grid;
+  grid-template-columns: repeat(2, minmax(0, 1fr));
+  gap: 1rem;
+  pointer-events: none;
+}
+
+.ns-top-hud-item {
+  background: transparent;
+  border: 0;
+  padding: 0;
+}
+
+.ns-top-hud-item span {
+  display: block;
+  color: var(--ns-muted);
+  font-size: 11px;
+  letter-spacing: 0.08em;
+  text-transform: uppercase;
+}
+
+.ns-top-hud-item strong {
+  display: block;
+  margin-top: 8px;
+  font-size: clamp(1.1rem, 1.6vw, 1.6rem);
+  line-height: 1;
+  color: var(--ns-text);
+}
+
+.ns-console {
+  position: absolute;
+  left: clamp(16px, 4vw, 24px);
+  top: clamp(16px, 4vw, 24px);
+  bottom: clamp(16px, 4vw, 24px);
+  z-index: 2;
+  width: min(340px, calc(100% - 2rem));
+  display: flex;
+  flex-direction: column;
+  gap: 1rem;
+  padding: clamp(16px, 2.5vw, 24px);
+  border-radius: 24px;
+  background: rgba(0, 0, 0, 0.02);
+  border: 1px solid rgba(255, 255, 255, 0.08);
+  backdrop-filter: blur(16px);
+}
+
+.theme-dark .ns-console {
+  background: rgba(255, 255, 255, 0.03);
+}
+
+.theme-light .ns-console {
+  background: rgba(255, 255, 255, 0.84);
+  border-color: rgba(148, 163, 184, 0.18);
+  box-shadow: 0 18px 50px rgba(148, 163, 184, 0.14);
+}
+
+.ns-console-header h3 {
+  margin: 0.2rem 0 0;
+  font-size: 1.15rem;
+  line-height: 1.1;
+}
+
+.ns-console-kicker {
+  display: inline-block;
+  font-size: 0.72rem;
+  letter-spacing: 0.18em;
+  text-transform: uppercase;
+  color: var(--ns-muted);
+}
+
+.ns-slider-row {
+  position: relative;
+  display: grid;
+  gap: 0.45rem;
+}
+
+.ns-slider-row label {
+  display: flex;
+  align-items: center;
+  justify-content: space-between;
+  gap: 1rem;
+  font-size: 0.9rem;
+  font-weight: 600;
+  color: var(--ns-muted);
+}
+
+.ns-slider-value {
+  position: static;
+  padding: 0;
+  border-radius: 0;
+  background: transparent;
+  color: var(--ns-accent);
+  font-size: 0.78rem;
+  line-height: 1;
+  font-variant-numeric: tabular-nums;
+  opacity: 1;
+  transition: none;
+  pointer-events: none;
+}
+
+.theme-dark .ns-slider-value {
+  background: transparent;
+  color: var(--ns-accent);
+}
+
+.ns-slider-value.is-visible {
+  opacity: 1;
+}
+
+.ns-console input[type="range"] {
+  width: 100%;
+  appearance: none;
+  -webkit-appearance: none;
+  height: 18px;
+  background: transparent;
+  margin: 0;
+}
+
+.ns-console input[type="range"]::-webkit-slider-runnable-track {
+  height: 1px;
+  background: var(--ns-track);
+  border-radius: 999px;
+}
+
+.ns-console input[type="range"]::-moz-range-track {
+  height: 1px;
+  background: var(--ns-track);
+  border-radius: 999px;
+}
+
+.ns-console input[type="range"]::-webkit-slider-thumb {
+  -webkit-appearance: none;
+  appearance: none;
+  width: 16px;
+  height: 16px;
+  margin-top: -8px;
+  border-radius: 999px;
+  border: 0;
+  background: var(--ns-thumb);
+  box-shadow: 0 2px 4px rgba(0, 0, 0, 0.2);
+}
+
+.ns-console input[type="range"]::-moz-range-thumb {
+  width: 16px;
+  height: 16px;
+  border-radius: 999px;
+  border: 0;
+  background: var(--ns-thumb);
+  box-shadow: 0 2px 4px rgba(0, 0, 0, 0.2);
+}
+
+.ns-run-btn {
+  margin-top: auto;
+  height: 40px;
+  display: inline-flex;
+  align-items: center;
+  justify-content: center;
+  gap: 0.5rem;
+  border: 0;
+  border-radius: 999px;
+  padding: 0.9rem 1rem;
+  background: linear-gradient(135deg, rgba(22, 163, 74, 0.95), rgba(56, 189, 248, 0.92));
+  color: #001018;
+  font-weight: 800;
+  box-shadow: 0 16px 30px rgba(34, 197, 94, 0.16);
+  transition: transform 0.18s ease, box-shadow 0.18s ease, filter 0.18s ease;
+}
+
+.ns-run-btn:hover {
+  transform: translateY(-1px);
+  filter: saturate(1.06);
+}
+
+.ns-run-btn:active {
+  transform: translateY(1px) scale(0.99);
+}
+
+.ns-hud {
+  position: absolute;
+  left: max(360px, clamp(16px, 3vw, 24px));
+  right: clamp(16px, 3vw, 24px);
+  bottom: clamp(16px, 3vw, 24px);
+  z-index: 2;
+  display: grid;
+  grid-template-columns: repeat(2, minmax(0, 1fr));
+  gap: 1rem;
+  pointer-events: none;
+}
+
+.ns-hud-item {
+  padding: 0.5rem 0;
+  background: transparent;
+  border: 0;
+}
+
+.ns-hud-item span {
+  display: block;
+  color: var(--ns-muted);
+  font-size: 11px;
+  letter-spacing: 0.08em;
+  text-transform: uppercase;
+}
+
+.ns-hud-item strong {
+  display: block;
+  margin-top: 8px;
+  font-size: clamp(1.4rem, 2vw, 2.5rem);
+  line-height: 1;
+  color: var(--ns-text);
+}
+
+.ns-alleles {
+  position: absolute;
+  right: clamp(16px, 3vw, 24px);
+  top: clamp(72px, 8vw, 108px);
+  z-index: 2;
+  display: grid;
+  grid-template-columns: minmax(180px, 226px);
+  gap: 0.3rem;
+  pointer-events: none;
+}
+
+.ns-allele-chip {
+  --chip-accent: rgba(255, 255, 255, 0.24);
+  --rank-scale: 0.5;
+
+  display: grid;
+  grid-template-columns: 1fr auto;
+  align-items: center;
+  gap: 0.35rem;
+  padding: 0.36rem 0.5rem 0.36rem 0.78rem;
+  border-radius: 14px;
+  background: rgba(2, 6, 23, 0.42);
+  border: 1px solid rgba(255, 255, 255, 0.08);
+  backdrop-filter: blur(14px);
+  position: relative;
+  overflow: hidden;
+  transition: transform 0.18s ease, opacity 0.18s ease;
+}
+
+.ns-allele-chip::before {
+  content: "";
+  position: absolute;
+  inset: 0;
+  width: 4px;
+  background: var(--chip-accent);
+}
+
+.ns-allele-chip.allele-1 { --chip-accent: linear-gradient(180deg, #7dd3fc, #0284c7); }
+.ns-allele-chip.allele-2 { --chip-accent: linear-gradient(180deg, #86efac, #16a34a); }
+.ns-allele-chip.allele-3 { --chip-accent: linear-gradient(180deg, #f9a8d4, #db2777); }
+.ns-allele-chip.allele-4 { --chip-accent: linear-gradient(180deg, #fde68a, #f59e0b); }
+.ns-allele-chip.allele-5 { --chip-accent: linear-gradient(180deg, #c4b5fd, #7c3aed); }
+.ns-allele-chip.allele-6 { --chip-accent: linear-gradient(180deg, #fdba74, #ea580c); }
+.ns-allele-chip.allele-7 { --chip-accent: linear-gradient(180deg, #6ee7b7, #059669); }
+.ns-allele-chip.allele-8 { --chip-accent: linear-gradient(180deg, #fca5a5, #ef4444); }
+
+.ns-allele-chip.is-updated {
+  animation: chipPulse 180ms ease-out;
+}
+
+.theme-dark .ns-allele-chip {
+  background: rgba(255, 255, 255, 0.03);
+}
+
+.theme-light .ns-allele-chip {
+  background: rgba(255, 255, 255, 0.78);
+  border-color: rgba(148, 163, 184, 0.18);
+  box-shadow: 0 12px 28px rgba(148, 163, 184, 0.12);
+}
+
+.ns-allele-chip span {
+  font-size: 8px;
+  letter-spacing: 0.08em;
+  text-transform: uppercase;
+  color: var(--ns-muted);
+}
+
+.ns-allele-chip strong {
+  font-size: 0.92rem;
+  line-height: 1;
+  color: var(--ns-text);
+  text-shadow: 0 0 10px rgba(255,255,255,0.04);
+  justify-self: end;
+}
+
+.ns-allele-bar {
+  grid-column: 1 / -1;
+  height: 2px;
+  border-radius: 999px;
+  background: rgba(255,255,255,0.06);
+  overflow: hidden;
+}
+
+.ns-allele-bar i {
+  display: block;
+  width: calc(var(--rank-scale) * 100%);
+  height: 100%;
+  border-radius: inherit;
+  background: var(--chip-accent);
+  box-shadow: 0 0 18px rgba(255,255,255,0.18);
+  transition: width 180ms ease;
+}
+
+@keyframes chipPulse {
+  0% { transform: translateY(0) scale(1); }
+  50% { transform: translateY(-1px) scale(1.02); }
+  100% { transform: translateY(0) scale(1); }
+}
+
+@media (min-width: 900px) {
+  .ns-console {
+    width: 320px;
+  }
+}
+
+@media (max-width: 840px) {
+  .ns-top-hud,
+  .ns-hud {
+    left: clamp(16px, 3vw, 24px);
+    right: clamp(16px, 3vw, 24px);
+  }
+
+  .ns-alleles {
+    top: 124px;
+    grid-template-columns: minmax(0, 1fr);
+  }
+}

package/src/tool/natural-selection-drift/render.ts

@@ -0,0 +1,219 @@
+import type { Particle } from './simulation';
+
+type SimConfig = {
+  population: number;
+  mutationRate: number;
+  selectionPressure: number;
+  driftIntensity: number;
+  alleleCount: number;
+  innovationRate: number;
+};
+
+type Outcome = { diversity: number };
+
+const pickFactor = (allele: number) => {
+  if (allele === 0) return 1.1;
+  if (allele === 1) return 0.95;
+  return 0.8;
+};
+
+const pickFertility = (allele: number, base: number) => {
+  if (allele === 0) return base * 1.08;
+  if (allele === 1) return base;
+  return base * 0.96;
+};
+const pickChildAllele = (parentAllele: number, mutation: number, innovationRate: number, alleleCount: number) => {
+  const roll = Math.random();
+  if (roll < mutation * 0.6) return (parentAllele + 1 + Math.floor(Math.random() * 2)) % alleleCount;
+  if (roll < mutation * 0.6 + innovationRate) return alleleCount + Math.floor(Math.random() * 3);
+  return parentAllele;
+};
+
+function clampPosition(value: number) {
+  return Math.max(0, Math.min(1, value));
+}
+
+function updateNeighborCohesion(particle: Particle, neighbors: { x: number; y: number }[]) {
+  let cohesionX = 0;
+  let cohesionY = 0;
+  neighbors.forEach((other) => {
+    cohesionX += other.x - particle.x;
+    cohesionY += other.y - particle.y;
+  });
+  return { cohesionX, cohesionY };
+}
+
+function moveParticle(args: {
+  particle: Particle;
+  cohesionX: number;
+  cohesionY: number;
+  config: SimConfig;
+  index: number;
+  squeeze: number;
+}) {
+  const { particle, cohesionX, cohesionY, config, index, squeeze } = args;
+  const centerPull = 0.00018 + config.selectionPressure * 0.0012;
+  const jitter = Math.sin(index + performance.now() * 0.0007) * 0.00025 + (Math.random() - 0.5) * 0.00045;
+  particle.vx += (0.5 - particle.x) * centerPull + jitter + cohesionX * (1 - Math.max(0.08, config.driftIntensity)) * 0.0011 + (Math.random() - 0.5) * config.driftIntensity * 0.0008;
+  particle.vy += (0.5 - particle.y) * centerPull + jitter + cohesionY * (1 - Math.max(0.08, config.driftIntensity)) * 0.0011 + (Math.random() - 0.5) * config.driftIntensity * 0.0008;
+  particle.vx *= 0.99;
+  particle.vy *= 0.99;
+  particle.x += particle.vx;
+  particle.y += particle.vy;
+
+  const dx = particle.x - 0.5;
+  const dy = particle.y - 0.5;
+  const distance = Math.sqrt(dx * dx + dy * dy);
+  if (distance > 0.52 * squeeze) {
+    const scale = (0.52 * squeeze) / distance;
+    particle.x = 0.5 + dx * scale;
+    particle.y = 0.5 + dy * scale;
+    particle.vx *= -0.25;
+    particle.vy *= -0.25;
+  }
+}
+
+function mutateAndCull(particle: Particle, config: SimConfig) {
+  if (Math.random() < config.mutationRate * 0.02) {
+    particle.allele = (particle.allele + 1 + Math.floor(Math.random() * 2)) % Math.max(2, Number(config.alleleCount || 3));
+    particle.vx += (Math.random() - 0.5) * 0.006;
+    particle.vy += (Math.random() - 0.5) * 0.006;
+  }
+  const agePressure = Math.min(1.35, 0.55 + particle.age * 0.01);
+  const killChance = config.selectionPressure * 0.00006 * agePressure * pickFactor(particle.allele);
+  if (Math.random() < killChance) {
+    particle.alive = false;
+  }
+}
+
+function gatherAlleleCentroids(particles: Particle[]) {
+  const groups = new Map<number, { x: number; y: number; count: number }>();
+  particles.forEach((particle) => {
+    if (!particle.alive) return;
+    const group = groups.get(particle.allele) || { x: 0, y: 0, count: 0 };
+    group.x += particle.x;
+    group.y += particle.y;
+    group.count += 1;
+    groups.set(particle.allele, group);
+  });
+  return groups;
+}
+
+function createBirthFromParent(parent: Particle, config: SimConfig, mutationPush: number) {
+  const childAllele = pickChildAllele(parent.allele, config.mutationRate, config.innovationRate, Math.max(2, Number(config.alleleCount || 3)));
+  return {
+    x: clampPosition(parent.x + (Math.random() - 0.5) * 0.02),
+    y: clampPosition(parent.y + (Math.random() - 0.5) * 0.02),
+    vx: parent.vx * 0.72 + (Math.random() - 0.5) * mutationPush * 0.008,
+    vy: parent.vy * 0.72 + (Math.random() - 0.5) * mutationPush * 0.008,
+    size: Math.max(1.1, parent.size * (0.9 + Math.random() * 0.18)),
+    allele: childAllele,
+    alive: true,
+    hue: parent.hue,
+    energy: Math.max(0.25, Math.min(0.95, parent.energy + (Math.random() - 0.5) * 0.2)),
+    bias: parent.bias,
+    age: 0,
+  } satisfies Particle;
+}
+
+function createFallbackBirth(parent: Particle) {
+  return {
+    x: clampPosition(parent.x + (Math.random() - 0.5) * 0.02),
+    y: clampPosition(parent.y + (Math.random() - 0.5) * 0.02),
+    vx: parent.vx * 0.7,
+    vy: parent.vy * 0.7,
+    size: Math.max(1.1, parent.size * 0.95),
+    allele: parent.allele,
+    alive: true,
+    hue: parent.hue,
+    energy: parent.energy,
+    bias: parent.bias,
+    age: 0,
+  } satisfies Particle;
+}
+
+function reproduce(particles: Particle[], config: SimConfig, diversityFactor: number) {
+  const survivors = particles.filter((particle) => particle.alive);
+  const mutationPush = Math.max(0.1, config.mutationRate * 1.6);
+  const reproductionRate = 0.035 + config.selectionPressure * 0.02 + diversityFactor * 0.012;
+  const minBirths = Math.max(1, Math.round(survivors.length * 0.06));
+  const targetBirths = Math.min(Math.max(2, minBirths), Math.max(6, Math.round(config.population * 0.08)));
+  const births: Particle[] = [];
+
+  for (let i = 0; i < survivors.length && births.length < targetBirths; i += 1) {
+    const parent = survivors[i];
+    if (Math.random() > pickFertility(parent.allele, reproductionRate)) continue;
+    births.push(createBirthFromParent(parent, config, mutationPush));
+  }
+
+  while (births.length < minBirths && survivors.length > 0) {
+    const parent = survivors[births.length % survivors.length];
+    births.push(createFallbackBirth(parent));
+  }
+
+  return births;
+}
+
+function trimPopulation(particles: Particle[], population: number) {
+  if (particles.length <= population * 1.25) return particles;
+  return particles.sort((a, b) => {
+    if (a.alive !== b.alive) return a.alive ? -1 : 1;
+    return a.age - b.age;
+  }).slice(-Math.round(population * 1.25));
+}
+
+export function evolveParticles(particles: Particle[], config: SimConfig, outcome: Outcome) {
+  const diversityFactor = Math.max(0.08, outcome.diversity / 100);
+  const squeeze = 1 - config.selectionPressure * 0.3;
+  const groups = gatherAlleleCentroids(particles);
+
+  particles.forEach((particle, index) => {
+    particle.age += 1;
+    if (!particle.alive) return;
+    const centroid = groups.get(particle.allele);
+    const neighbors = centroid && centroid.count > 1 ? [{ x: centroid.x / centroid.count, y: centroid.y / centroid.count }] : [];
+    const { cohesionX, cohesionY } = updateNeighborCohesion(particle, neighbors);
+    moveParticle({ particle, cohesionX, cohesionY, config, index, squeeze });
+    mutateAndCull(particle, config);
+  });
+
+  const births = reproduce(particles, config, diversityFactor);
+  if (births.length > 0) {
+    const limit = Math.min(births.length, Math.max(4, Math.round(config.population * 0.08)));
+    for (let i = 0; i < limit; i += 1) {
+      particles.push(births[i]);
+    }
+  }
+  return trimPopulation(particles, config.population);
+}
+
+export function drawBackground(ctx: CanvasRenderingContext2D, rect: DOMRect, pressure: number) {
+  const background = ctx.createRadialGradient(rect.width / 2, rect.height / 2, 40, rect.width / 2, rect.height / 2, rect.width * 0.65);
+  background.addColorStop(0, `rgba(56, 189, 248, ${0.08 + pressure * 0.12})`);
+  background.addColorStop(0.45, 'rgba(34, 197, 94, 0.08)');
+  background.addColorStop(1, 'rgba(2, 6, 23, 0.14)');
+  ctx.fillStyle = background;
+  ctx.fillRect(0, 0, rect.width, rect.height);
+  const centerGlow = ctx.createRadialGradient(rect.width / 2, rect.height / 2, 20, rect.width / 2, rect.height / 2, rect.width * 0.34);
+  centerGlow.addColorStop(0, `rgba(255,255,255,${0.02 + pressure * 0.04})`);
+  centerGlow.addColorStop(1, 'rgba(255,255,255,0)');
+  ctx.fillStyle = centerGlow;
+  ctx.fillRect(0, 0, rect.width, rect.height);
+}
+
+export function drawParticles(ctx: CanvasRenderingContext2D, rect: DOMRect, particles: Particle[], atmosphere: { pressure: number; drift: number }) {
+  const { pressure, drift } = atmosphere;
+  particles.forEach((particle) => {
+    if (!particle.alive) return;
+    const x = particle.x * rect.width;
+    const y = particle.y * rect.height;
+    const radius = particle.size * (0.75 + pressure * 0.35);
+    ctx.beginPath();
+    ctx.fillStyle = `hsla(${particle.hue}, 86%, ${62 + particle.energy * 14}%, ${0.55 + particle.energy * 0.35})`;
+    ctx.shadowBlur = 10 + drift * 24;
+    ctx.shadowColor = `hsla(${particle.hue}, 86%, 66%, 0.22)`;
+    ctx.arc(x, y, radius, 0, Math.PI * 2);
+    ctx.fill();
+  });
+  ctx.shadowBlur = 0;
+}