[](https://gleam.run/) [](https://www.erlang.org/) [](https://gleam.run/news/multi-target-compilation/) [](https://hex.pm/packages/viva_math) [](https://en.wikipedia.org/wiki/PAD_emotional_state_model) [](./test) [](./CHANGELOG.md) [](./LICENSE) --- *"Numbers feel. Catastrophes are smooth. The BEAM thinks in PAD."*

[!IMPORTANT] viva_math IS NOT JUST A MATH LIB. It is the kernel of VIVA's affective computing stack — PAD emotional space, Cusp catastrophe for mood transitions, Friston's Free Energy Principle, attractor dynamics, and information theory — all expressed as small, type-safe Gleam functions.

Built on gleam_community_maths, targeting Erlang + JavaScript.

🎯 Overview

Core mathematical foundations for VIVA — a sentient digital life research project. The library models emotional dynamics as a dynamical system: PAD state vectors, Cusp catastrophe for sudden mood shifts, Free Energy for interoception, and Shannon-family entropy for affective complexity.

The math is small, deliberate, and grounded in real papers (Mehrabian 1996, Thom 1972, Friston 2010, Shannon 1948).

⚡ Quick Start

gleam add viva_math

import viva_math/attractor
import viva_math/cusp
import viva_math/free_energy
import viva_math/vector

pub fn main() {
  // PAD emotional state: Pleasure / Arousal / Dominance
  let state = vector.pad(-0.3, 0.7, -0.2)

  // Nearest discrete emotion
  let emotion = attractor.classify_emotion(state)
  // -> "fear"

  // Bistability check (sudden mood shift possible?)
  let params = cusp.from_arousal_dominance(0.7, -0.2)
  let volatile = cusp.is_bistable(params)
  // -> True

  // Prediction error vs expected state
  let expected = vector.pad(0.0, 0.0, 0.0)
  let fe = free_energy.compute_state(expected, state, expected, 0.1)
  // fe.feeling -> Surprised | Alarmed
}

🏗️ Architecture

   ┌──────────────────────────────────────────────────────────┐
   │                Gleam application code                    │
   │              viva_math/{common,vector,...}               │
   └────────┬─────────────────────────────────────────────────┘
            │
   ┌────────▼─────────────────────────────────────────────────┐
   │                   viva_math modules                      │
   │                                                          │
   │   ┌────────┐  ┌────────┐  ┌────────┐  ┌──────────────┐   │
   │   │ common │  │ vector │  │  cusp  │  │ free_energy  │   │
   │   │ clamp  │  │ Vec3   │  │ Thom   │  │  Friston     │   │
   │   │ sigmoid│  │ PAD ℝ³ │  │ 1972   │  │  2010        │   │
   │   │ noise  │  │ ops    │  │ stoch. │  │  homeostasis │   │
   │   └────────┘  └────────┘  └────────┘  └──────────────┘   │
   │                                                          │
   │   ┌────────────┐                  ┌──────────────────┐   │
   │   │ attractor  │                  │     entropy      │   │
   │   │ Mehrabian  │                  │  Shannon · KL    │   │
   │   │   1996     │                  │  JS · Rényi      │   │
   │   │ 8 emotions │                  │  hybrid affect.  │   │
   │   └────────────┘                  └──────────────────┘   │
   └────────┬─────────────────────────────────────────────────┘
            │
   ┌────────▼─────────────────────────────────────────────────┐
   │            gleam_community_maths (trig, stats, …)        │
   └──────────────────────────────────────────────────────────┘

🧬 Theoretical Background

PAD Model — Mehrabian (1996)

Emotions live as points in 3D vector space:

• Pleasure[-1, 1] — sadness ↔ joy
• Arousal[-1, 1] — calm ↔ excitement
• Dominance[-1, 1] — submission ↔ control

Cusp Catastrophe — Thom (1972)

Sudden mood transitions modeled by the potential:

V(x) = x⁴/4 + αx²/2 + βx

When arousal pushes α < 0 and the discriminant Δ > 0, the system goes bistable — tiny perturbations trigger discrete mood jumps. The stochastic variant adds Wiener noise (dV/dx + σξ(t)) integrated via Euler-Maruyama.

Free Energy Principle — Friston (2010)

Agents minimize "surprise" via prediction:

F ≈ Prediction_Error² + Complexity

Low free energy → predictions match reality (homeostasis). High free energy → significant mismatch (alarm).

Attractor Dynamics

Eight basic emotions form attractors in PAD space:

Information Theory — Shannon (1948) + extensions

Shannon entropy H(p), KL divergence D(p‖q), Jensen-Shannon JS(p,q), Rényi H_α, and a hybrid affective entropyH_hybrid = αH₁ + (1-α)H₂ for mixed emotional states.

🎨 Design Principles

📚 Public API Highlights

Vec3 / PAD space

import viva_math/vector

let v = vector.pad(0.5, -0.2, 0.7)
let mag = vector.magnitude(v)
let unit = vector.normalize(v)
let dist = vector.euclidean(v, vector.pad(0.0, 0.0, 0.0))

Cusp catastrophe — deterministic + stochastic

import viva_math/cusp

let params = cusp.from_arousal_dominance(0.7, -0.2)
let v = cusp.potential(params, 0.3)         // V(x)
let g = cusp.gradient(params, 0.3)          // dV/dx

// Stochastic mood walk
let trajectory =
cusp.simulate_stochastic(
cusp.StochasticCuspParams(params, sigma: 0.05, seed: 42),
start: 0.0,
steps: 1000,
dt: 0.01,
)

Free Energy

import viva_math/free_energy
import viva_math/vector

let expected = vector.pad(0.0, 0.0, 0.0)
let actual   = vector.pad(-0.3, 0.7, -0.2)
let state    = free_energy.compute_state(expected, actual, expected, 0.1)
// state.feeling -> Calm | Surprised | Alarmed
// state.free_energy -> Float

Entropy and divergence

import viva_math/entropy

let h = entropy.shannon([0.2, 0.3, 0.5])
let d = entropy.kl_divergence([0.2, 0.8], [0.5, 0.5])
let js = entropy.jensen_shannon([0.2, 0.8], [0.5, 0.5])
let r = entropy.renyi([0.2, 0.3, 0.5], alpha: 2.0)

🗺️ Roadmap

🤝 Contributing

git checkout -b feature/your-feature
gleam test                  # 58 tests
gleam format --check src test
gleam docs build

See CONTRIBUTING.md for guidelines.

📖 References

• Mehrabian (1996) — Pleasure-arousal-dominance: A general framework
• Thom (1972) — Structural Stability and Morphogenesis
• Friston (2010) — The free-energy principle: a unified brain theory?
• Grasman et al. (2009) — Fitting the Cusp Catastrophe in R
• Oravecz et al. (2009) — Ornstein-Uhlenbeck Process in Affective Dynamics
• Shannon (1948) — A Mathematical Theory of Communication

🌌 VIVA Ecosystem

**Star if math should feel like something ⭐** [](https://github.com/gabrielmaialva33/viva_math) *Created by Gabriel Maia · MIT License*