formulab 0.5.1 → 0.9.0

package/dist/electronics/index.js

@@ -1,3 +1,5 @@
+// Guards
+export { isOhmsLawInput } from './guards.js';
 // Functions
 export { ohmsLaw } from './ohmsLaw.js';
 export { reflowProfile, getPasteTypes } from './reflow.js';

package/dist/electronics/index.js.map

@@ -1 +1 @@
-{"version":3,"file":"index.js","sourceRoot":"","sources":["../../src/electronics/index.ts"],"names":[],"mappings":"AA6CA,YAAY;AACZ,OAAO,EAAE,OAAO,EAAE,MAAM,cAAc,CAAC;AACvC,OAAO,EAAE,aAAa,EAAE,aAAa,EAAE,MAAM,aAAa,CAAC;AAC3D,OAAO,EAAE,cAAc,EAAE,MAAM,eAAe,CAAC;AAC/C,OAAO,EAAE,OAAO,EAAE,MAAM,eAAe,CAAC;AACxC,OAAO,EAAE,WAAW,EAAE,MAAM,aAAa,CAAC;AAC1C,OAAO,EAAE,UAAU,EAAE,MAAM,YAAY,CAAC;AACxC,OAAO,EAAE,aAAa,EAAE,MAAM,UAAU,CAAC;AACzC,OAAO,EAAE,eAAe,EAAE,MAAM,gBAAgB,CAAC;AACjD,OAAO,EAAE,WAAW,EAAE,MAAM,UAAU,CAAC;AACvC,OAAO,EAAE,eAAe,EAAE,MAAM,cAAc,CAAC;AAC/C,OAAO,EAAE,UAAU,EAAE,MAAM,UAAU,CAAC"}
+{"version":3,"file":"index.js","sourceRoot":"","sources":["../../src/electronics/index.ts"],"names":[],"mappings":"AA6CA,SAAS;AACT,OAAO,EAAE,cAAc,EAAE,MAAM,aAAa,CAAC;AAE7C,YAAY;AACZ,OAAO,EAAE,OAAO,EAAE,MAAM,cAAc,CAAC;AACvC,OAAO,EAAE,aAAa,EAAE,aAAa,EAAE,MAAM,aAAa,CAAC;AAC3D,OAAO,EAAE,cAAc,EAAE,MAAM,eAAe,CAAC;AAC/C,OAAO,EAAE,OAAO,EAAE,MAAM,eAAe,CAAC;AACxC,OAAO,EAAE,WAAW,EAAE,MAAM,aAAa,CAAC;AAC1C,OAAO,EAAE,UAAU,EAAE,MAAM,YAAY,CAAC;AACxC,OAAO,EAAE,aAAa,EAAE,MAAM,UAAU,CAAC;AACzC,OAAO,EAAE,eAAe,EAAE,MAAM,gBAAgB,CAAC;AACjD,OAAO,EAAE,WAAW,EAAE,MAAM,UAAU,CAAC;AACvC,OAAO,EAAE,eAAe,EAAE,MAAM,cAAc,CAAC;AAC/C,OAAO,EAAE,UAAU,EAAE,MAAM,UAAU,CAAC"}

package/dist/energy/boilerEfficiency.d.ts

@@ -0,0 +1,13 @@
+import type { BoilerEfficiencyInput, BoilerEfficiencyResult } from './types.js';
+/**
+ * Calculate boiler thermal efficiency using the direct (input-output) method.
+ *
+ * η = (steamOutput × (steamEnthalpy − feedwaterEnthalpy)) / (fuelRate × fuelHeatValue) × 100
+ *
+ * Heat values are converted from kJ to kW by dividing by 3600 (kJ/h → kW).
+ *
+ * @param input - Boiler operating parameters
+ * @returns Boiler efficiency result with heat balance and optional annual metrics
+ */
+export declare function boilerEfficiency(input: BoilerEfficiencyInput): BoilerEfficiencyResult;
+//# sourceMappingURL=boilerEfficiency.d.ts.map

package/dist/energy/boilerEfficiency.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"boilerEfficiency.d.ts","sourceRoot":"","sources":["../../src/energy/boilerEfficiency.ts"],"names":[],"mappings":"AACA,OAAO,KAAK,EAAE,qBAAqB,EAAE,sBAAsB,EAAE,MAAM,YAAY,CAAC;AAEhF;;;;;;;;;GASG;AACH,wBAAgB,gBAAgB,CAAC,KAAK,EAAE,qBAAqB,GAAG,sBAAsB,CA8CrF"}

package/dist/energy/boilerEfficiency.js

@@ -0,0 +1,46 @@
+import { roundTo } from '../utils.js';
+/**
+ * Calculate boiler thermal efficiency using the direct (input-output) method.
+ *
+ * η = (steamOutput × (steamEnthalpy − feedwaterEnthalpy)) / (fuelRate × fuelHeatValue) × 100
+ *
+ * Heat values are converted from kJ to kW by dividing by 3600 (kJ/h → kW).
+ *
+ * @param input - Boiler operating parameters
+ * @returns Boiler efficiency result with heat balance and optional annual metrics
+ */
+export function boilerEfficiency(input) {
+    const { fuelRate, fuelHeatValue, steamOutput, steamEnthalpy, feedwaterEnthalpy, operatingHours, fuelCost, } = input;
+    if (fuelRate <= 0 || fuelHeatValue <= 0) {
+        return {
+            heatInput: 0,
+            heatOutput: 0,
+            heatLoss: 0,
+            efficiency: 0,
+            annualFuelCost: null,
+            annualHeatLoss: null,
+        };
+    }
+    // kJ/h → kW (÷ 3600)
+    const heatInput = (fuelRate * fuelHeatValue) / 3600;
+    const heatOutput = (steamOutput * (steamEnthalpy - feedwaterEnthalpy)) / 3600;
+    const heatLoss = heatInput - heatOutput;
+    const efficiency = (heatOutput / heatInput) * 100;
+    let annualFuelCost = null;
+    if (operatingHours !== undefined && fuelCost !== undefined) {
+        annualFuelCost = fuelRate * operatingHours * fuelCost;
+    }
+    let annualHeatLoss = null;
+    if (operatingHours !== undefined) {
+        annualHeatLoss = heatLoss * operatingHours;
+    }
+    return {
+        heatInput: roundTo(heatInput, 4),
+        heatOutput: roundTo(heatOutput, 4),
+        heatLoss: roundTo(heatLoss, 4),
+        efficiency: roundTo(efficiency, 4),
+        annualFuelCost: annualFuelCost !== null ? roundTo(annualFuelCost, 2) : null,
+        annualHeatLoss: annualHeatLoss !== null ? roundTo(annualHeatLoss, 4) : null,
+    };
+}
+//# sourceMappingURL=boilerEfficiency.js.map

package/dist/energy/boilerEfficiency.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"boilerEfficiency.js","sourceRoot":"","sources":["../../src/energy/boilerEfficiency.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,OAAO,EAAE,MAAM,aAAa,CAAC;AAGtC;;;;;;;;;GASG;AACH,MAAM,UAAU,gBAAgB,CAAC,KAA4B;IAC3D,MAAM,EACJ,QAAQ,EACR,aAAa,EACb,WAAW,EACX,aAAa,EACb,iBAAiB,EACjB,cAAc,EACd,QAAQ,GACT,GAAG,KAAK,CAAC;IAEV,IAAI,QAAQ,IAAI,CAAC,IAAI,aAAa,IAAI,CAAC,EAAE,CAAC;QACxC,OAAO;YACL,SAAS,EAAE,CAAC;YACZ,UAAU,EAAE,CAAC;YACb,QAAQ,EAAE,CAAC;YACX,UAAU,EAAE,CAAC;YACb,cAAc,EAAE,IAAI;YACpB,cAAc,EAAE,IAAI;SACrB,CAAC;IACJ,CAAC;IAED,qBAAqB;IACrB,MAAM,SAAS,GAAG,CAAC,QAAQ,GAAG,aAAa,CAAC,GAAG,IAAI,CAAC;IACpD,MAAM,UAAU,GAAG,CAAC,WAAW,GAAG,CAAC,aAAa,GAAG,iBAAiB,CAAC,CAAC,GAAG,IAAI,CAAC;IAC9E,MAAM,QAAQ,GAAG,SAAS,GAAG,UAAU,CAAC;IACxC,MAAM,UAAU,GAAG,CAAC,UAAU,GAAG,SAAS,CAAC,GAAG,GAAG,CAAC;IAElD,IAAI,cAAc,GAAkB,IAAI,CAAC;IACzC,IAAI,cAAc,KAAK,SAAS,IAAI,QAAQ,KAAK,SAAS,EAAE,CAAC;QAC3D,cAAc,GAAG,QAAQ,GAAG,cAAc,GAAG,QAAQ,CAAC;IACxD,CAAC;IAED,IAAI,cAAc,GAAkB,IAAI,CAAC;IACzC,IAAI,cAAc,KAAK,SAAS,EAAE,CAAC;QACjC,cAAc,GAAG,QAAQ,GAAG,cAAc,CAAC;IAC7C,CAAC;IAED,OAAO;QACL,SAAS,EAAE,OAAO,CAAC,SAAS,EAAE,CAAC,CAAC;QAChC,UAAU,EAAE,OAAO,CAAC,UAAU,EAAE,CAAC,CAAC;QAClC,QAAQ,EAAE,OAAO,CAAC,QAAQ,EAAE,CAAC,CAAC;QAC9B,UAAU,EAAE,OAAO,CAAC,UAAU,EAAE,CAAC,CAAC;QAClC,cAAc,EAAE,cAAc,KAAK,IAAI,CAAC,CAAC,CAAC,OAAO,CAAC,cAAc,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI;QAC3E,cAAc,EAAE,cAAc,KAAK,IAAI,CAAC,CAAC,CAAC,OAAO,CAAC,cAAc,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI;KAC5E,CAAC;AACJ,CAAC"}

package/dist/energy/cusum.d.ts

@@ -0,0 +1,15 @@
+import type { CusumInput, CusumResult } from './types.js';
+/**
+ * CUSUM (Cumulative Sum) Control Chart — Page's Algorithm
+ *
+ * @formula
+ *   - C⁺ₙ = max(0, xₙ − (μ₀ + K) + C⁺ₙ₋₁)
+ *   - C⁻ₙ = max(0, (μ₀ − K) − xₙ + C⁻ₙ₋₁)
+ *   - Signal when C⁺ or C⁻ > H
+ *
+ * @reference Page, E.S. (1954). "Continuous inspection schemes"
+ * @reference ISO 50001:2018 — CUSUM for energy monitoring
+ * @reference Montgomery, D.C. "Introduction to Statistical Quality Control", Ch. 9
+ */
+export declare function cusum(input: CusumInput): CusumResult;
+//# sourceMappingURL=cusum.d.ts.map

package/dist/energy/cusum.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"cusum.d.ts","sourceRoot":"","sources":["../../src/energy/cusum.ts"],"names":[],"mappings":"AACA,OAAO,KAAK,EAAE,UAAU,EAAE,WAAW,EAAE,MAAM,YAAY,CAAC;AAE1D;;;;;;;;;;;GAWG;AACH,wBAAgB,KAAK,CAAC,KAAK,EAAE,UAAU,GAAG,WAAW,CA4DpD"}

package/dist/energy/cusum.js

@@ -0,0 +1,71 @@
+import { roundTo } from '../utils.js';
+/**
+ * CUSUM (Cumulative Sum) Control Chart — Page's Algorithm
+ *
+ * @formula
+ *   - C⁺ₙ = max(0, xₙ − (μ₀ + K) + C⁺ₙ₋₁)
+ *   - C⁻ₙ = max(0, (μ₀ − K) − xₙ + C⁻ₙ₋₁)
+ *   - Signal when C⁺ or C⁻ > H
+ *
+ * @reference Page, E.S. (1954). "Continuous inspection schemes"
+ * @reference ISO 50001:2018 — CUSUM for energy monitoring
+ * @reference Montgomery, D.C. "Introduction to Statistical Quality Control", Ch. 9
+ */
+export function cusum(input) {
+    const { values, target } = input;
+    const n = values.length;
+    if (n === 0) {
+        return {
+            cusumPositive: [],
+            cusumNegative: [],
+            signals: [],
+            isOutOfControl: false,
+            shiftDetected: 'none',
+        };
+    }
+    // Auto-calculate stdDev if not provided
+    let sigma = input.stdDev;
+    if (sigma == null) {
+        const mean = values.reduce((s, v) => s + v, 0) / n;
+        sigma = Math.sqrt(values.reduce((s, v) => s + (v - mean) ** 2, 0) / (n - 1));
+        if (!Number.isFinite(sigma) || sigma <= 0)
+            sigma = 1;
+    }
+    const K = input.allowance ?? sigma / 2;
+    const H = input.decisionInterval ?? 5 * sigma;
+    const cusumPositive = [];
+    const cusumNegative = [];
+    const signals = [];
+    let hasPositiveSignal = false;
+    let hasNegativeSignal = false;
+    let cPlus = 0;
+    let cMinus = 0;
+    for (let i = 0; i < n; i++) {
+        cPlus = Math.max(0, values[i] - (target + K) + cPlus);
+        cMinus = Math.max(0, (target - K) - values[i] + cMinus);
+        cusumPositive.push(roundTo(cPlus, 4));
+        cusumNegative.push(roundTo(cMinus, 4));
+        if (cPlus > H || cMinus > H) {
+            signals.push(i);
+            if (cPlus > H)
+                hasPositiveSignal = true;
+            if (cMinus > H)
+                hasNegativeSignal = true;
+        }
+    }
+    let shiftDetected = 'none';
+    if (hasPositiveSignal && hasNegativeSignal)
+        shiftDetected = 'both';
+    else if (hasPositiveSignal)
+        shiftDetected = 'positive';
+    else if (hasNegativeSignal)
+        shiftDetected = 'negative';
+    return {
+        cusumPositive,
+        cusumNegative,
+        signals,
+        isOutOfControl: signals.length > 0,
+        shiftDetected,
+    };
+}
+//# sourceMappingURL=cusum.js.map

package/dist/energy/cusum.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"cusum.js","sourceRoot":"","sources":["../../src/energy/cusum.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,OAAO,EAAE,MAAM,aAAa,CAAC;AAGtC;;;;;;;;;;;GAWG;AACH,MAAM,UAAU,KAAK,CAAC,KAAiB;IACrC,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,GAAG,KAAK,CAAC;IACjC,MAAM,CAAC,GAAG,MAAM,CAAC,MAAM,CAAC;IAExB,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC;QACZ,OAAO;YACL,aAAa,EAAE,EAAE;YACjB,aAAa,EAAE,EAAE;YACjB,OAAO,EAAE,EAAE;YACX,cAAc,EAAE,KAAK;YACrB,aAAa,EAAE,MAAM;SACtB,CAAC;IACJ,CAAC;IAED,wCAAwC;IACxC,IAAI,KAAK,GAAG,KAAK,CAAC,MAAM,CAAC;IACzB,IAAI,KAAK,IAAI,IAAI,EAAE,CAAC;QAClB,MAAM,IAAI,GAAG,MAAM,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC;QACnD,KAAK,GAAG,IAAI,CAAC,IAAI,CAAC,MAAM,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QAC7E,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC,IAAI,KAAK,IAAI,CAAC;YAAE,KAAK,GAAG,CAAC,CAAC;IACvD,CAAC;IAED,MAAM,CAAC,GAAG,KAAK,CAAC,SAAS,IAAI,KAAK,GAAG,CAAC,CAAC;IACvC,MAAM,CAAC,GAAG,KAAK,CAAC,gBAAgB,IAAI,CAAC,GAAG,KAAK,CAAC;IAE9C,MAAM,aAAa,GAAa,EAAE,CAAC;IACnC,MAAM,aAAa,GAAa,EAAE,CAAC;IACnC,MAAM,OAAO,GAAa,EAAE,CAAC;IAC7B,IAAI,iBAAiB,GAAG,KAAK,CAAC;IAC9B,IAAI,iBAAiB,GAAG,KAAK,CAAC;IAE9B,IAAI,KAAK,GAAG,CAAC,CAAC;IACd,IAAI,MAAM,GAAG,CAAC,CAAC;IAEf,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC3B,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,MAAM,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC;QACtD,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,MAAM,GAAG,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC;QAExD,aAAa,CAAC,IAAI,CAAC,OAAO,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC,CAAC;QACtC,aAAa,CAAC,IAAI,CAAC,OAAO,CAAC,MAAM,EAAE,CAAC,CAAC,CAAC,CAAC;QAEvC,IAAI,KAAK,GAAG,CAAC,IAAI,MAAM,GAAG,CAAC,EAAE,CAAC;YAC5B,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;YAChB,IAAI,KAAK,GAAG,CAAC;gBAAE,iBAAiB,GAAG,IAAI,CAAC;YACxC,IAAI,MAAM,GAAG,CAAC;gBAAE,iBAAiB,GAAG,IAAI,CAAC;QAC3C,CAAC;IACH,CAAC;IAED,IAAI,aAAa,GAAiC,MAAM,CAAC;IACzD,IAAI,iBAAiB,IAAI,iBAAiB;QAAE,aAAa,GAAG,MAAM,CAAC;SAC9D,IAAI,iBAAiB;QAAE,aAAa,GAAG,UAAU,CAAC;SAClD,IAAI,iBAAiB;QAAE,aAAa,GAAG,UAAU,CAAC;IAEvD,OAAO;QACL,aAAa;QACb,aAAa;QACb,OAAO;QACP,cAAc,EAAE,OAAO,CAAC,MAAM,GAAG,CAAC;QAClC,aAAa;KACd,CAAC;AACJ,CAAC"}

package/dist/energy/degreeDay.d.ts

@@ -0,0 +1,13 @@
+import type { DegreeDayInput, DegreeDayResult } from './types.js';
+/**
+ * Degree Day (HDD/CDD) Calculator — ISO 50001 EnB
+ *
+ * @formula
+ *   - HDD = Σ max(0, T_base_heating − T_daily)
+ *   - CDD = Σ max(0, T_daily − T_base_cooling)
+ *
+ * @reference ASHRAE Fundamentals Handbook — Degree Day Method
+ * @reference ISO 50001:2018 — Energy management systems
+ */
+export declare function degreeDay(input: DegreeDayInput): DegreeDayResult;
+//# sourceMappingURL=degreeDay.d.ts.map

package/dist/energy/degreeDay.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"degreeDay.d.ts","sourceRoot":"","sources":["../../src/energy/degreeDay.ts"],"names":[],"mappings":"AACA,OAAO,KAAK,EAAE,cAAc,EAAE,eAAe,EAAE,MAAM,YAAY,CAAC;AAElE;;;;;;;;;GASG;AACH,wBAAgB,SAAS,CAAC,KAAK,EAAE,cAAc,GAAG,eAAe,CAoChE"}

package/dist/energy/degreeDay.js

@@ -0,0 +1,46 @@
+import { roundTo } from '../utils.js';
+/**
+ * Degree Day (HDD/CDD) Calculator — ISO 50001 EnB
+ *
+ * @formula
+ *   - HDD = Σ max(0, T_base_heating − T_daily)
+ *   - CDD = Σ max(0, T_daily − T_base_cooling)
+ *
+ * @reference ASHRAE Fundamentals Handbook — Degree Day Method
+ * @reference ISO 50001:2018 — Energy management systems
+ */
+export function degreeDay(input) {
+    const { dailyTemps, baseHeating = 18, baseCooling = 24 } = input;
+    const totalDays = dailyTemps.length;
+    if (totalDays === 0) {
+        return { hdd: 0, cdd: 0, totalDays: 0, heatingDays: 0, coolingDays: 0, neutralDays: 0, avgTemp: 0 };
+    }
+    let hdd = 0;
+    let cdd = 0;
+    let heatingDays = 0;
+    let coolingDays = 0;
+    let sum = 0;
+    for (const temp of dailyTemps) {
+        sum += temp;
+        const hddContrib = Math.max(0, baseHeating - temp);
+        const cddContrib = Math.max(0, temp - baseCooling);
+        hdd += hddContrib;
+        cdd += cddContrib;
+        if (hddContrib > 0)
+            heatingDays++;
+        if (cddContrib > 0)
+            coolingDays++;
+    }
+    const neutralDays = totalDays - heatingDays - coolingDays;
+    const avgTemp = sum / totalDays;
+    return {
+        hdd: roundTo(hdd, 2),
+        cdd: roundTo(cdd, 2),
+        totalDays,
+        heatingDays,
+        coolingDays,
+        neutralDays,
+        avgTemp: roundTo(avgTemp, 2),
+    };
+}
+//# sourceMappingURL=degreeDay.js.map

package/dist/energy/degreeDay.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"degreeDay.js","sourceRoot":"","sources":["../../src/energy/degreeDay.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,OAAO,EAAE,MAAM,aAAa,CAAC;AAGtC;;;;;;;;;GASG;AACH,MAAM,UAAU,SAAS,CAAC,KAAqB;IAC7C,MAAM,EAAE,UAAU,EAAE,WAAW,GAAG,EAAE,EAAE,WAAW,GAAG,EAAE,EAAE,GAAG,KAAK,CAAC;IACjE,MAAM,SAAS,GAAG,UAAU,CAAC,MAAM,CAAC;IAEpC,IAAI,SAAS,KAAK,CAAC,EAAE,CAAC;QACpB,OAAO,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,SAAS,EAAE,CAAC,EAAE,WAAW,EAAE,CAAC,EAAE,WAAW,EAAE,CAAC,EAAE,WAAW,EAAE,CAAC,EAAE,OAAO,EAAE,CAAC,EAAE,CAAC;IACtG,CAAC;IAED,IAAI,GAAG,GAAG,CAAC,CAAC;IACZ,IAAI,GAAG,GAAG,CAAC,CAAC;IACZ,IAAI,WAAW,GAAG,CAAC,CAAC;IACpB,IAAI,WAAW,GAAG,CAAC,CAAC;IACpB,IAAI,GAAG,GAAG,CAAC,CAAC;IAEZ,KAAK,MAAM,IAAI,IAAI,UAAU,EAAE,CAAC;QAC9B,GAAG,IAAI,IAAI,CAAC;QACZ,MAAM,UAAU,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,WAAW,GAAG,IAAI,CAAC,CAAC;QACnD,MAAM,UAAU,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,GAAG,WAAW,CAAC,CAAC;QACnD,GAAG,IAAI,UAAU,CAAC;QAClB,GAAG,IAAI,UAAU,CAAC;QAClB,IAAI,UAAU,GAAG,CAAC;YAAE,WAAW,EAAE,CAAC;QAClC,IAAI,UAAU,GAAG,CAAC;YAAE,WAAW,EAAE,CAAC;IACpC,CAAC;IAED,MAAM,WAAW,GAAG,SAAS,GAAG,WAAW,GAAG,WAAW,CAAC;IAC1D,MAAM,OAAO,GAAG,GAAG,GAAG,SAAS,CAAC;IAEhC,OAAO;QACL,GAAG,EAAE,OAAO,CAAC,GAAG,EAAE,CAAC,CAAC;QACpB,GAAG,EAAE,OAAO,CAAC,GAAG,EAAE,CAAC,CAAC;QACpB,SAAS;QACT,WAAW;QACX,WAAW;QACX,WAAW;QACX,OAAO,EAAE,OAAO,CAAC,OAAO,EAAE,CAAC,CAAC;KAC7B,CAAC;AACJ,CAAC"}

package/dist/energy/heatPump.d.ts

@@ -0,0 +1,13 @@
+import type { HeatPumpInput, HeatPumpResult } from './types.js';
+/**
+ * Heat Pump COP Calculator
+ *
+ * @formula
+ *   - COP = heatingCapacity / (compressorPower + auxiliaryPower)
+ *   - COP_Carnot = T_sink / (T_sink − T_source)  (Kelvin)
+ *   - Efficiency = COP / COP_Carnot × 100
+ *
+ * @reference EN 14511 — Air conditioners, liquid chilling packages and heat pumps
+ */
+export declare function heatPump(input: HeatPumpInput): HeatPumpResult;
+//# sourceMappingURL=heatPump.d.ts.map

package/dist/energy/heatPump.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"heatPump.d.ts","sourceRoot":"","sources":["../../src/energy/heatPump.ts"],"names":[],"mappings":"AACA,OAAO,KAAK,EAAE,aAAa,EAAE,cAAc,EAAE,MAAM,YAAY,CAAC;AAEhE;;;;;;;;;GASG;AACH,wBAAgB,QAAQ,CAAC,KAAK,EAAE,aAAa,GAAG,cAAc,CA+C7D"}

package/dist/energy/heatPump.js

@@ -0,0 +1,48 @@
+import { roundTo } from '../utils.js';
+/**
+ * Heat Pump COP Calculator
+ *
+ * @formula
+ *   - COP = heatingCapacity / (compressorPower + auxiliaryPower)
+ *   - COP_Carnot = T_sink / (T_sink − T_source)  (Kelvin)
+ *   - Efficiency = COP / COP_Carnot × 100
+ *
+ * @reference EN 14511 — Air conditioners, liquid chilling packages and heat pumps
+ */
+export function heatPump(input) {
+    const { sourceTemp, sinkTemp, heatingCapacity, compressorPower, auxiliaryPower = 0, operatingHours, electricityRate, boilerEfficiency, fuelCost, } = input;
+    const totalPower = compressorPower + auxiliaryPower;
+    const cop = totalPower > 0 ? heatingCapacity / totalPower : 0;
+    // Carnot COP (theoretical maximum)
+    const tSinkK = sinkTemp + 273.15;
+    const tSourceK = sourceTemp + 273.15;
+    const deltaT = tSinkK - tSourceK;
+    const copCarnot = deltaT > 0 ? tSinkK / deltaT : 0;
+    const efficiency = copCarnot > 0 ? roundTo((cop / copCarnot) * 100, 2) : 0;
+    // Annual calculations
+    let annualElectricity = null;
+    let annualElecCost = null;
+    let annualFuelCost = null;
+    let annualSavings = null;
+    if (operatingHours != null) {
+        annualElectricity = roundTo(totalPower * operatingHours, 2);
+        if (electricityRate != null) {
+            annualElecCost = roundTo(annualElectricity * electricityRate, 2);
+            if (boilerEfficiency != null && fuelCost != null && boilerEfficiency > 0) {
+                const heatDelivered = heatingCapacity * operatingHours;
+                annualFuelCost = roundTo((heatDelivered / boilerEfficiency) * fuelCost, 2);
+                annualSavings = roundTo(annualFuelCost - annualElecCost, 2);
+            }
+        }
+    }
+    return {
+        cop: roundTo(cop, 3),
+        copCarnot: roundTo(copCarnot, 3),
+        efficiency,
+        annualElectricity,
+        annualElecCost,
+        annualFuelCost,
+        annualSavings,
+    };
+}
+//# sourceMappingURL=heatPump.js.map

package/dist/energy/heatPump.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"heatPump.js","sourceRoot":"","sources":["../../src/energy/heatPump.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,OAAO,EAAE,MAAM,aAAa,CAAC;AAGtC;;;;;;;;;GASG;AACH,MAAM,UAAU,QAAQ,CAAC,KAAoB;IAC3C,MAAM,EACJ,UAAU,EAAE,QAAQ,EAAE,eAAe,EAAE,eAAe,EACtD,cAAc,GAAG,CAAC,EAAE,cAAc,EAAE,eAAe,EACnD,gBAAgB,EAAE,QAAQ,GAC3B,GAAG,KAAK,CAAC;IAEV,MAAM,UAAU,GAAG,eAAe,GAAG,cAAc,CAAC;IACpD,MAAM,GAAG,GAAG,UAAU,GAAG,CAAC,CAAC,CAAC,CAAC,eAAe,GAAG,UAAU,CAAC,CAAC,CAAC,CAAC,CAAC;IAE9D,mCAAmC;IACnC,MAAM,MAAM,GAAG,QAAQ,GAAG,MAAM,CAAC;IACjC,MAAM,QAAQ,GAAG,UAAU,GAAG,MAAM,CAAC;IACrC,MAAM,MAAM,GAAG,MAAM,GAAG,QAAQ,CAAC;IACjC,MAAM,SAAS,GAAG,MAAM,GAAG,CAAC,CAAC,CAAC,CAAC,MAAM,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC;IAEnD,MAAM,UAAU,GAAG,SAAS,GAAG,CAAC,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC,GAAG,GAAG,SAAS,CAAC,GAAG,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAE3E,sBAAsB;IACtB,IAAI,iBAAiB,GAAkB,IAAI,CAAC;IAC5C,IAAI,cAAc,GAAkB,IAAI,CAAC;IACzC,IAAI,cAAc,GAAkB,IAAI,CAAC;IACzC,IAAI,aAAa,GAAkB,IAAI,CAAC;IAExC,IAAI,cAAc,IAAI,IAAI,EAAE,CAAC;QAC3B,iBAAiB,GAAG,OAAO,CAAC,UAAU,GAAG,cAAc,EAAE,CAAC,CAAC,CAAC;QAE5D,IAAI,eAAe,IAAI,IAAI,EAAE,CAAC;YAC5B,cAAc,GAAG,OAAO,CAAC,iBAAiB,GAAG,eAAe,EAAE,CAAC,CAAC,CAAC;YAEjE,IAAI,gBAAgB,IAAI,IAAI,IAAI,QAAQ,IAAI,IAAI,IAAI,gBAAgB,GAAG,CAAC,EAAE,CAAC;gBACzE,MAAM,aAAa,GAAG,eAAe,GAAG,cAAc,CAAC;gBACvD,cAAc,GAAG,OAAO,CAAC,CAAC,aAAa,GAAG,gBAAgB,CAAC,GAAG,QAAQ,EAAE,CAAC,CAAC,CAAC;gBAC3E,aAAa,GAAG,OAAO,CAAC,cAAc,GAAG,cAAc,EAAE,CAAC,CAAC,CAAC;YAC9D,CAAC;QACH,CAAC;IACH,CAAC;IAED,OAAO;QACL,GAAG,EAAE,OAAO,CAAC,GAAG,EAAE,CAAC,CAAC;QACpB,SAAS,EAAE,OAAO,CAAC,SAAS,EAAE,CAAC,CAAC;QAChC,UAAU;QACV,iBAAiB;QACjB,cAAc;QACd,cAAc;QACd,aAAa;KACd,CAAC;AACJ,CAAC"}

package/dist/energy/index.d.ts

@@ -5,5 +5,13 @@ export { motorEfficiency } from './motorEfficiency.js';
 export { pfCorrection } from './pfCorrection.js';
 export { powerCost } from './powerCost.js';
 export { vfdSavings } from './vfdSavings.js';
-export type { SolarOutputInput, SolarOutputResult, CarbonFootprintInput, CarbonFootprintResult, CompressedAirCostInput, CompressedAirCostResult, MotorEfficiencyInput, MotorEfficiencyResult, PfCorrectionInput, PfCorrectionResult, PowerCostInput, PowerCostResult, VfdSavingsInput, VfdSavingsResult, } from './types.js';
+export { boilerEfficiency } from './boilerEfficiency.js';
+export { transformerLoss } from './transformerLoss.js';
+export { insulationRoi } from './insulationRoi.js';
+export { ledRoi } from './ledRoi.js';
+export { heatPump } from './heatPump.js';
+export { degreeDay } from './degreeDay.js';
+export { windOutput } from './windOutput.js';
+export { cusum } from './cusum.js';
+export type { SolarOutputInput, SolarOutputResult, CarbonFootprintInput, CarbonFootprintResult, CompressedAirCostInput, CompressedAirCostResult, MotorEfficiencyInput, MotorEfficiencyResult, PfCorrectionInput, PfCorrectionResult, PowerCostInput, PowerCostResult, VfdSavingsInput, VfdSavingsResult, BoilerEfficiencyInput, BoilerEfficiencyResult, TransformerLossInput, TransformerLossResult, InsulationRoiInput, InsulationRoiResult, LedRoiInput, LedRoiResult, HeatPumpInput, HeatPumpResult, DegreeDayInput, DegreeDayResult, WindOutputInput, WindOutputResult, CusumInput, CusumResult, } from './types.js';
 //# sourceMappingURL=index.d.ts.map

package/dist/energy/index.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"index.d.ts","sourceRoot":"","sources":["../../src/energy/index.ts"],"names":[],"mappings":"AACA,OAAO,EAAE,WAAW,EAAE,MAAM,kBAAkB,CAAC;AAC/C,OAAO,EAAE,eAAe,EAAE,MAAM,sBAAsB,CAAC;AACvD,OAAO,EAAE,iBAAiB,EAAE,MAAM,wBAAwB,CAAC;AAC3D,OAAO,EAAE,eAAe,EAAE,MAAM,sBAAsB,CAAC;AACvD,OAAO,EAAE,YAAY,EAAE,MAAM,mBAAmB,CAAC;AACjD,OAAO,EAAE,SAAS,EAAE,MAAM,gBAAgB,CAAC;AAC3C,OAAO,EAAE,UAAU,EAAE,MAAM,iBAAiB,CAAC;AAG7C,YAAY,EAEV,gBAAgB,EAChB,iBAAiB,EAEjB,oBAAoB,EACpB,qBAAqB,EAErB,sBAAsB,EACtB,uBAAuB,EAEvB,oBAAoB,EACpB,qBAAqB,EAErB,iBAAiB,EACjB,kBAAkB,EAElB,cAAc,EACd,eAAe,EAEf,eAAe,EACf,gBAAgB,GACjB,MAAM,YAAY,CAAC"}
+{"version":3,"file":"index.d.ts","sourceRoot":"","sources":["../../src/energy/index.ts"],"names":[],"mappings":"AACA,OAAO,EAAE,WAAW,EAAE,MAAM,kBAAkB,CAAC;AAC/C,OAAO,EAAE,eAAe,EAAE,MAAM,sBAAsB,CAAC;AACvD,OAAO,EAAE,iBAAiB,EAAE,MAAM,wBAAwB,CAAC;AAC3D,OAAO,EAAE,eAAe,EAAE,MAAM,sBAAsB,CAAC;AACvD,OAAO,EAAE,YAAY,EAAE,MAAM,mBAAmB,CAAC;AACjD,OAAO,EAAE,SAAS,EAAE,MAAM,gBAAgB,CAAC;AAC3C,OAAO,EAAE,UAAU,EAAE,MAAM,iBAAiB,CAAC;AAC7C,OAAO,EAAE,gBAAgB,EAAE,MAAM,uBAAuB,CAAC;AACzD,OAAO,EAAE,eAAe,EAAE,MAAM,sBAAsB,CAAC;AACvD,OAAO,EAAE,aAAa,EAAE,MAAM,oBAAoB,CAAC;AACnD,OAAO,EAAE,MAAM,EAAE,MAAM,aAAa,CAAC;AACrC,OAAO,EAAE,QAAQ,EAAE,MAAM,eAAe,CAAC;AACzC,OAAO,EAAE,SAAS,EAAE,MAAM,gBAAgB,CAAC;AAC3C,OAAO,EAAE,UAAU,EAAE,MAAM,iBAAiB,CAAC;AAC7C,OAAO,EAAE,KAAK,EAAE,MAAM,YAAY,CAAC;AAGnC,YAAY,EAEV,gBAAgB,EAChB,iBAAiB,EAEjB,oBAAoB,EACpB,qBAAqB,EAErB,sBAAsB,EACtB,uBAAuB,EAEvB,oBAAoB,EACpB,qBAAqB,EAErB,iBAAiB,EACjB,kBAAkB,EAElB,cAAc,EACd,eAAe,EAEf,eAAe,EACf,gBAAgB,EAEhB,qBAAqB,EACrB,sBAAsB,EAEtB,oBAAoB,EACpB,qBAAqB,EAErB,kBAAkB,EAClB,mBAAmB,EAEnB,WAAW,EACX,YAAY,EAEZ,aAAa,EACb,cAAc,EAEd,cAAc,EACd,eAAe,EAEf,eAAe,EACf,gBAAgB,EAEhB,UAAU,EACV,WAAW,GACZ,MAAM,YAAY,CAAC"}

package/dist/energy/index.js

@@ -6,4 +6,12 @@ export { motorEfficiency } from './motorEfficiency.js';
 export { pfCorrection } from './pfCorrection.js';
 export { powerCost } from './powerCost.js';
 export { vfdSavings } from './vfdSavings.js';
+export { boilerEfficiency } from './boilerEfficiency.js';
+export { transformerLoss } from './transformerLoss.js';
+export { insulationRoi } from './insulationRoi.js';
+export { ledRoi } from './ledRoi.js';
+export { heatPump } from './heatPump.js';
+export { degreeDay } from './degreeDay.js';
+export { windOutput } from './windOutput.js';
+export { cusum } from './cusum.js';
 //# sourceMappingURL=index.js.map

package/dist/energy/index.js.map

@@ -1 +1 @@
-{"version":3,"file":"index.js","sourceRoot":"","sources":["../../src/energy/index.ts"],"names":[],"mappings":"AAAA,yBAAyB;AACzB,OAAO,EAAE,WAAW,EAAE,MAAM,kBAAkB,CAAC;AAC/C,OAAO,EAAE,eAAe,EAAE,MAAM,sBAAsB,CAAC;AACvD,OAAO,EAAE,iBAAiB,EAAE,MAAM,wBAAwB,CAAC;AAC3D,OAAO,EAAE,eAAe,EAAE,MAAM,sBAAsB,CAAC;AACvD,OAAO,EAAE,YAAY,EAAE,MAAM,mBAAmB,CAAC;AACjD,OAAO,EAAE,SAAS,EAAE,MAAM,gBAAgB,CAAC;AAC3C,OAAO,EAAE,UAAU,EAAE,MAAM,iBAAiB,CAAC"}
+{"version":3,"file":"index.js","sourceRoot":"","sources":["../../src/energy/index.ts"],"names":[],"mappings":"AAAA,yBAAyB;AACzB,OAAO,EAAE,WAAW,EAAE,MAAM,kBAAkB,CAAC;AAC/C,OAAO,EAAE,eAAe,EAAE,MAAM,sBAAsB,CAAC;AACvD,OAAO,EAAE,iBAAiB,EAAE,MAAM,wBAAwB,CAAC;AAC3D,OAAO,EAAE,eAAe,EAAE,MAAM,sBAAsB,CAAC;AACvD,OAAO,EAAE,YAAY,EAAE,MAAM,mBAAmB,CAAC;AACjD,OAAO,EAAE,SAAS,EAAE,MAAM,gBAAgB,CAAC;AAC3C,OAAO,EAAE,UAAU,EAAE,MAAM,iBAAiB,CAAC;AAC7C,OAAO,EAAE,gBAAgB,EAAE,MAAM,uBAAuB,CAAC;AACzD,OAAO,EAAE,eAAe,EAAE,MAAM,sBAAsB,CAAC;AACvD,OAAO,EAAE,aAAa,EAAE,MAAM,oBAAoB,CAAC;AACnD,OAAO,EAAE,MAAM,EAAE,MAAM,aAAa,CAAC;AACrC,OAAO,EAAE,QAAQ,EAAE,MAAM,eAAe,CAAC;AACzC,OAAO,EAAE,SAAS,EAAE,MAAM,gBAAgB,CAAC;AAC3C,OAAO,EAAE,UAAU,EAAE,MAAM,iBAAiB,CAAC;AAC7C,OAAO,EAAE,KAAK,EAAE,MAAM,YAAY,CAAC"}

package/dist/energy/insulationRoi.d.ts

@@ -0,0 +1,14 @@
+import type { InsulationRoiInput, InsulationRoiResult } from './types.js';
+/**
+ * Calculate insulation ROI for pipes, vessels, or surfaces.
+ *
+ * bareHeatLoss = surfaceCoefficient × area × ΔT
+ * insulatedHeatLoss = area × ΔT / (1/surfaceCoefficient + thickness/k)
+ * annualEnergySaved = (bare − insulated) × hours / 1000 / boilerEfficiency
+ * payback = installationCost / annualCostSaved
+ *
+ * @param input - Insulation parameters
+ * @returns Insulation ROI result with heat savings and payback
+ */
+export declare function insulationRoi(input: InsulationRoiInput): InsulationRoiResult;
+//# sourceMappingURL=insulationRoi.d.ts.map

package/dist/energy/insulationRoi.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"insulationRoi.d.ts","sourceRoot":"","sources":["../../src/energy/insulationRoi.ts"],"names":[],"mappings":"AACA,OAAO,KAAK,EAAE,kBAAkB,EAAE,mBAAmB,EAAE,MAAM,YAAY,CAAC;AAE1E;;;;;;;;;;GAUG;AACH,wBAAgB,aAAa,CAAC,KAAK,EAAE,kBAAkB,GAAG,mBAAmB,CAuD5E"}

package/dist/energy/insulationRoi.js

@@ -0,0 +1,52 @@
+import { roundTo } from '../utils.js';
+/**
+ * Calculate insulation ROI for pipes, vessels, or surfaces.
+ *
+ * bareHeatLoss = surfaceCoefficient × area × ΔT
+ * insulatedHeatLoss = area × ΔT / (1/surfaceCoefficient + thickness/k)
+ * annualEnergySaved = (bare − insulated) × hours / 1000 / boilerEfficiency
+ * payback = installationCost / annualCostSaved
+ *
+ * @param input - Insulation parameters
+ * @returns Insulation ROI result with heat savings and payback
+ */
+export function insulationRoi(input) {
+    const { surfaceArea, tempDifference, insulationK, insulationThickness, surfaceCoefficient = 10, operatingHours, energyCost, boilerEfficiency = 0.8, installationCost, } = input;
+    if (surfaceArea <= 0 || tempDifference <= 0 || insulationK <= 0 || insulationThickness <= 0) {
+        return {
+            bareHeatLoss: 0,
+            insulatedHeatLoss: 0,
+            heatSaved: 0,
+            heatLossReduction: 0,
+            annualEnergySaved: 0,
+            annualCostSaved: 0,
+            paybackPeriod: null,
+        };
+    }
+    // Bare surface heat loss: Q = h × A × ΔT (W)
+    const bareHeatLoss = surfaceCoefficient * surfaceArea * tempDifference;
+    // Insulated heat loss: Q = A × ΔT / (1/h + thickness_m / k)
+    const thicknessM = insulationThickness / 1000; // mm → m
+    const thermalResistance = (1 / surfaceCoefficient) + (thicknessM / insulationK);
+    const insulatedHeatLoss = (surfaceArea * tempDifference) / thermalResistance;
+    const heatSaved = bareHeatLoss - insulatedHeatLoss;
+    const heatLossReduction = bareHeatLoss > 0 ? (heatSaved / bareHeatLoss) * 100 : 0;
+    // Annual energy saved (kWh) = heat saved (W) × hours / 1000 / boiler efficiency
+    // Divide by boiler efficiency because saved heat = saved fuel input
+    const annualEnergySaved = (heatSaved * operatingHours) / 1000 / boilerEfficiency;
+    const annualCostSaved = annualEnergySaved * energyCost;
+    let paybackPeriod = null;
+    if (installationCost !== undefined && installationCost > 0 && annualCostSaved > 0) {
+        paybackPeriod = installationCost / annualCostSaved;
+    }
+    return {
+        bareHeatLoss: roundTo(bareHeatLoss, 4),
+        insulatedHeatLoss: roundTo(insulatedHeatLoss, 4),
+        heatSaved: roundTo(heatSaved, 4),
+        heatLossReduction: roundTo(heatLossReduction, 4),
+        annualEnergySaved: roundTo(annualEnergySaved, 4),
+        annualCostSaved: roundTo(annualCostSaved, 2),
+        paybackPeriod: paybackPeriod !== null ? roundTo(paybackPeriod, 4) : null,
+    };
+}
+//# sourceMappingURL=insulationRoi.js.map

package/dist/energy/insulationRoi.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"insulationRoi.js","sourceRoot":"","sources":["../../src/energy/insulationRoi.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,OAAO,EAAE,MAAM,aAAa,CAAC;AAGtC;;;;;;;;;;GAUG;AACH,MAAM,UAAU,aAAa,CAAC,KAAyB;IACrD,MAAM,EACJ,WAAW,EACX,cAAc,EACd,WAAW,EACX,mBAAmB,EACnB,kBAAkB,GAAG,EAAE,EACvB,cAAc,EACd,UAAU,EACV,gBAAgB,GAAG,GAAG,EACtB,gBAAgB,GACjB,GAAG,KAAK,CAAC;IAEV,IAAI,WAAW,IAAI,CAAC,IAAI,cAAc,IAAI,CAAC,IAAI,WAAW,IAAI,CAAC,IAAI,mBAAmB,IAAI,CAAC,EAAE,CAAC;QAC5F,OAAO;YACL,YAAY,EAAE,CAAC;YACf,iBAAiB,EAAE,CAAC;YACpB,SAAS,EAAE,CAAC;YACZ,iBAAiB,EAAE,CAAC;YACpB,iBAAiB,EAAE,CAAC;YACpB,eAAe,EAAE,CAAC;YAClB,aAAa,EAAE,IAAI;SACpB,CAAC;IACJ,CAAC;IAED,6CAA6C;IAC7C,MAAM,YAAY,GAAG,kBAAkB,GAAG,WAAW,GAAG,cAAc,CAAC;IAEvE,4DAA4D;IAC5D,MAAM,UAAU,GAAG,mBAAmB,GAAG,IAAI,CAAC,CAAC,SAAS;IACxD,MAAM,iBAAiB,GAAG,CAAC,CAAC,GAAG,kBAAkB,CAAC,GAAG,CAAC,UAAU,GAAG,WAAW,CAAC,CAAC;IAChF,MAAM,iBAAiB,GAAG,CAAC,WAAW,GAAG,cAAc,CAAC,GAAG,iBAAiB,CAAC;IAE7E,MAAM,SAAS,GAAG,YAAY,GAAG,iBAAiB,CAAC;IACnD,MAAM,iBAAiB,GAAG,YAAY,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,SAAS,GAAG,YAAY,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;IAElF,gFAAgF;IAChF,oEAAoE;IACpE,MAAM,iBAAiB,GAAG,CAAC,SAAS,GAAG,cAAc,CAAC,GAAG,IAAI,GAAG,gBAAgB,CAAC;IACjF,MAAM,eAAe,GAAG,iBAAiB,GAAG,UAAU,CAAC;IAEvD,IAAI,aAAa,GAAkB,IAAI,CAAC;IACxC,IAAI,gBAAgB,KAAK,SAAS,IAAI,gBAAgB,GAAG,CAAC,IAAI,eAAe,GAAG,CAAC,EAAE,CAAC;QAClF,aAAa,GAAG,gBAAgB,GAAG,eAAe,CAAC;IACrD,CAAC;IAED,OAAO;QACL,YAAY,EAAE,OAAO,CAAC,YAAY,EAAE,CAAC,CAAC;QACtC,iBAAiB,EAAE,OAAO,CAAC,iBAAiB,EAAE,CAAC,CAAC;QAChD,SAAS,EAAE,OAAO,CAAC,SAAS,EAAE,CAAC,CAAC;QAChC,iBAAiB,EAAE,OAAO,CAAC,iBAAiB,EAAE,CAAC,CAAC;QAChD,iBAAiB,EAAE,OAAO,CAAC,iBAAiB,EAAE,CAAC,CAAC;QAChD,eAAe,EAAE,OAAO,CAAC,eAAe,EAAE,CAAC,CAAC;QAC5C,aAAa,EAAE,aAAa,KAAK,IAAI,CAAC,CAAC,CAAC,OAAO,CAAC,aAAa,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI;KACzE,CAAC;AACJ,CAAC"}

package/dist/energy/ledRoi.d.ts

@@ -0,0 +1,14 @@
+import type { LedRoiInput, LedRoiResult } from './types.js';
+/**
+ * Calculate LED lighting retrofit ROI.
+ *
+ * energySaved = fixtureCount × (oldWatts − newWatts) × hours / 1000
+ * costSaved = energySaved × rate
+ * payback = totalInvestment / costSaved
+ * co2Saved = energySaved × co2Factor
+ *
+ * @param input - LED retrofit parameters
+ * @returns LED ROI result with energy savings, cost savings, and payback
+ */
+export declare function ledRoi(input: LedRoiInput): LedRoiResult;
+//# sourceMappingURL=ledRoi.d.ts.map

package/dist/energy/ledRoi.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"ledRoi.d.ts","sourceRoot":"","sources":["../../src/energy/ledRoi.ts"],"names":[],"mappings":"AACA,OAAO,KAAK,EAAE,WAAW,EAAE,YAAY,EAAE,MAAM,YAAY,CAAC;AAE5D;;;;;;;;;;GAUG;AACH,wBAAgB,MAAM,CAAC,KAAK,EAAE,WAAW,GAAG,YAAY,CAkDvD"}

package/dist/energy/ledRoi.js

@@ -0,0 +1,49 @@
+import { roundTo } from '../utils.js';
+/**
+ * Calculate LED lighting retrofit ROI.
+ *
+ * energySaved = fixtureCount × (oldWatts − newWatts) × hours / 1000
+ * costSaved = energySaved × rate
+ * payback = totalInvestment / costSaved
+ * co2Saved = energySaved × co2Factor
+ *
+ * @param input - LED retrofit parameters
+ * @returns LED ROI result with energy savings, cost savings, and payback
+ */
+export function ledRoi(input) {
+    const { fixtureCount, oldWatts, newWatts, operatingHours, electricityRate, fixtureCost = 0, installationCost = 0, co2Factor = 0.5, } = input;
+    if (fixtureCount <= 0 || operatingHours <= 0) {
+        return {
+            oldAnnualEnergy: 0,
+            newAnnualEnergy: 0,
+            annualEnergySaved: 0,
+            energyReduction: 0,
+            annualCostSaved: 0,
+            totalInvestment: 0,
+            paybackPeriod: null,
+            co2Saved: 0,
+        };
+    }
+    const oldAnnualEnergy = (fixtureCount * oldWatts * operatingHours) / 1000;
+    const newAnnualEnergy = (fixtureCount * newWatts * operatingHours) / 1000;
+    const annualEnergySaved = oldAnnualEnergy - newAnnualEnergy;
+    const energyReduction = oldAnnualEnergy > 0 ? (annualEnergySaved / oldAnnualEnergy) * 100 : 0;
+    const annualCostSaved = annualEnergySaved * electricityRate;
+    const totalInvestment = (fixtureCost * fixtureCount) + installationCost;
+    let paybackPeriod = null;
+    if (totalInvestment > 0 && annualCostSaved > 0) {
+        paybackPeriod = totalInvestment / annualCostSaved;
+    }
+    const co2Saved = annualEnergySaved * co2Factor;
+    return {
+        oldAnnualEnergy: roundTo(oldAnnualEnergy, 4),
+        newAnnualEnergy: roundTo(newAnnualEnergy, 4),
+        annualEnergySaved: roundTo(annualEnergySaved, 4),
+        energyReduction: roundTo(energyReduction, 4),
+        annualCostSaved: roundTo(annualCostSaved, 2),
+        totalInvestment: roundTo(totalInvestment, 2),
+        paybackPeriod: paybackPeriod !== null ? roundTo(paybackPeriod, 4) : null,
+        co2Saved: roundTo(co2Saved, 4),
+    };
+}
+//# sourceMappingURL=ledRoi.js.map

package/dist/energy/ledRoi.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"ledRoi.js","sourceRoot":"","sources":["../../src/energy/ledRoi.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,OAAO,EAAE,MAAM,aAAa,CAAC;AAGtC;;;;;;;;;;GAUG;AACH,MAAM,UAAU,MAAM,CAAC,KAAkB;IACvC,MAAM,EACJ,YAAY,EACZ,QAAQ,EACR,QAAQ,EACR,cAAc,EACd,eAAe,EACf,WAAW,GAAG,CAAC,EACf,gBAAgB,GAAG,CAAC,EACpB,SAAS,GAAG,GAAG,GAChB,GAAG,KAAK,CAAC;IAEV,IAAI,YAAY,IAAI,CAAC,IAAI,cAAc,IAAI,CAAC,EAAE,CAAC;QAC7C,OAAO;YACL,eAAe,EAAE,CAAC;YAClB,eAAe,EAAE,CAAC;YAClB,iBAAiB,EAAE,CAAC;YACpB,eAAe,EAAE,CAAC;YAClB,eAAe,EAAE,CAAC;YAClB,eAAe,EAAE,CAAC;YAClB,aAAa,EAAE,IAAI;YACnB,QAAQ,EAAE,CAAC;SACZ,CAAC;IACJ,CAAC;IAED,MAAM,eAAe,GAAG,CAAC,YAAY,GAAG,QAAQ,GAAG,cAAc,CAAC,GAAG,IAAI,CAAC;IAC1E,MAAM,eAAe,GAAG,CAAC,YAAY,GAAG,QAAQ,GAAG,cAAc,CAAC,GAAG,IAAI,CAAC;IAC1E,MAAM,iBAAiB,GAAG,eAAe,GAAG,eAAe,CAAC;IAC5D,MAAM,eAAe,GAAG,eAAe,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,iBAAiB,GAAG,eAAe,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;IAC9F,MAAM,eAAe,GAAG,iBAAiB,GAAG,eAAe,CAAC;IAE5D,MAAM,eAAe,GAAG,CAAC,WAAW,GAAG,YAAY,CAAC,GAAG,gBAAgB,CAAC;IAExE,IAAI,aAAa,GAAkB,IAAI,CAAC;IACxC,IAAI,eAAe,GAAG,CAAC,IAAI,eAAe,GAAG,CAAC,EAAE,CAAC;QAC/C,aAAa,GAAG,eAAe,GAAG,eAAe,CAAC;IACpD,CAAC;IAED,MAAM,QAAQ,GAAG,iBAAiB,GAAG,SAAS,CAAC;IAE/C,OAAO;QACL,eAAe,EAAE,OAAO,CAAC,eAAe,EAAE,CAAC,CAAC;QAC5C,eAAe,EAAE,OAAO,CAAC,eAAe,EAAE,CAAC,CAAC;QAC5C,iBAAiB,EAAE,OAAO,CAAC,iBAAiB,EAAE,CAAC,CAAC;QAChD,eAAe,EAAE,OAAO,CAAC,eAAe,EAAE,CAAC,CAAC;QAC5C,eAAe,EAAE,OAAO,CAAC,eAAe,EAAE,CAAC,CAAC;QAC5C,eAAe,EAAE,OAAO,CAAC,eAAe,EAAE,CAAC,CAAC;QAC5C,aAAa,EAAE,aAAa,KAAK,IAAI,CAAC,CAAC,CAAC,OAAO,CAAC,aAAa,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI;QACxE,QAAQ,EAAE,OAAO,CAAC,QAAQ,EAAE,CAAC,CAAC;KAC/B,CAAC;AACJ,CAAC"}

package/dist/energy/transformerLoss.d.ts

@@ -0,0 +1,13 @@
+import type { TransformerLossInput, TransformerLossResult } from './types.js';
+/**
+ * Calculate transformer losses and efficiency.
+ *
+ * totalLoss = coreLoss + copperLoss × loadFraction²
+ * efficiency = (outputPower / (outputPower + totalLoss)) × 100
+ * optimalLoad = √(coreLoss / copperLoss)
+ *
+ * @param input - Transformer parameters
+ * @returns Transformer loss result with efficiency and optimal loading
+ */
+export declare function transformerLoss(input: TransformerLossInput): TransformerLossResult;
+//# sourceMappingURL=transformerLoss.d.ts.map

package/dist/energy/transformerLoss.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"transformerLoss.d.ts","sourceRoot":"","sources":["../../src/energy/transformerLoss.ts"],"names":[],"mappings":"AACA,OAAO,KAAK,EAAE,oBAAoB,EAAE,qBAAqB,EAAE,MAAM,YAAY,CAAC;AAE9E;;;;;;;;;GASG;AACH,wBAAgB,eAAe,CAAC,KAAK,EAAE,oBAAoB,GAAG,qBAAqB,CAmElF"}

package/dist/energy/transformerLoss.js

@@ -0,0 +1,63 @@
+import { roundTo } from '../utils.js';
+/**
+ * Calculate transformer losses and efficiency.
+ *
+ * totalLoss = coreLoss + copperLoss × loadFraction²
+ * efficiency = (outputPower / (outputPower + totalLoss)) × 100
+ * optimalLoad = √(coreLoss / copperLoss)
+ *
+ * @param input - Transformer parameters
+ * @returns Transformer loss result with efficiency and optimal loading
+ */
+export function transformerLoss(input) {
+    const { ratedCapacity, coreLoss, copperLoss, loadFraction, powerFactor = 0.85, operatingHours, energyCost, } = input;
+    if (ratedCapacity <= 0 || coreLoss < 0 || copperLoss < 0) {
+        return {
+            outputPower: 0,
+            totalLoss: 0,
+            coreLossAtLoad: 0,
+            copperLossAtLoad: 0,
+            efficiency: 0,
+            optimalLoadFraction: 0,
+            annualLossEnergy: null,
+            annualLossCost: null,
+        };
+    }
+    // Output power in W = kVA × loadFraction × powerFactor × 1000
+    const outputPower = ratedCapacity * loadFraction * powerFactor * 1000;
+    // Core loss is constant regardless of load
+    const coreLossAtLoad = coreLoss;
+    // Copper loss scales with square of load fraction
+    const copperLossAtLoad = copperLoss * loadFraction * loadFraction;
+    const totalLoss = coreLossAtLoad + copperLossAtLoad;
+    // Efficiency = output / (output + losses) × 100
+    let efficiency = 0;
+    if (outputPower + totalLoss > 0) {
+        efficiency = (outputPower / (outputPower + totalLoss)) * 100;
+    }
+    // Optimal load fraction where core loss = copper loss at load
+    // coreLoss = copperLoss × k² → k = √(coreLoss / copperLoss)
+    let optimalLoadFraction = 0;
+    if (copperLoss > 0) {
+        optimalLoadFraction = Math.sqrt(coreLoss / copperLoss);
+    }
+    let annualLossEnergy = null;
+    let annualLossCost = null;
+    if (operatingHours !== undefined) {
+        annualLossEnergy = (totalLoss * operatingHours) / 1000; // W·h → kWh
+        if (energyCost !== undefined) {
+            annualLossCost = annualLossEnergy * energyCost;
+        }
+    }
+    return {
+        outputPower: roundTo(outputPower, 4),
+        totalLoss: roundTo(totalLoss, 4),
+        coreLossAtLoad: roundTo(coreLossAtLoad, 4),
+        copperLossAtLoad: roundTo(copperLossAtLoad, 4),
+        efficiency: roundTo(efficiency, 4),
+        optimalLoadFraction: roundTo(optimalLoadFraction, 4),
+        annualLossEnergy: annualLossEnergy !== null ? roundTo(annualLossEnergy, 4) : null,
+        annualLossCost: annualLossCost !== null ? roundTo(annualLossCost, 2) : null,
+    };
+}
+//# sourceMappingURL=transformerLoss.js.map

package/dist/energy/transformerLoss.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"transformerLoss.js","sourceRoot":"","sources":["../../src/energy/transformerLoss.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,OAAO,EAAE,MAAM,aAAa,CAAC;AAGtC;;;;;;;;;GASG;AACH,MAAM,UAAU,eAAe,CAAC,KAA2B;IACzD,MAAM,EACJ,aAAa,EACb,QAAQ,EACR,UAAU,EACV,YAAY,EACZ,WAAW,GAAG,IAAI,EAClB,cAAc,EACd,UAAU,GACX,GAAG,KAAK,CAAC;IAEV,IAAI,aAAa,IAAI,CAAC,IAAI,QAAQ,GAAG,CAAC,IAAI,UAAU,GAAG,CAAC,EAAE,CAAC;QACzD,OAAO;YACL,WAAW,EAAE,CAAC;YACd,SAAS,EAAE,CAAC;YACZ,cAAc,EAAE,CAAC;YACjB,gBAAgB,EAAE,CAAC;YACnB,UAAU,EAAE,CAAC;YACb,mBAAmB,EAAE,CAAC;YACtB,gBAAgB,EAAE,IAAI;YACtB,cAAc,EAAE,IAAI;SACrB,CAAC;IACJ,CAAC;IAED,8DAA8D;IAC9D,MAAM,WAAW,GAAG,aAAa,GAAG,YAAY,GAAG,WAAW,GAAG,IAAI,CAAC;IAEtE,2CAA2C;IAC3C,MAAM,cAAc,GAAG,QAAQ,CAAC;IAEhC,kDAAkD;IAClD,MAAM,gBAAgB,GAAG,UAAU,GAAG,YAAY,GAAG,YAAY,CAAC;IAElE,MAAM,SAAS,GAAG,cAAc,GAAG,gBAAgB,CAAC;IAEpD,gDAAgD;IAChD,IAAI,UAAU,GAAG,CAAC,CAAC;IACnB,IAAI,WAAW,GAAG,SAAS,GAAG,CAAC,EAAE,CAAC;QAChC,UAAU,GAAG,CAAC,WAAW,GAAG,CAAC,WAAW,GAAG,SAAS,CAAC,CAAC,GAAG,GAAG,CAAC;IAC/D,CAAC;IAED,8DAA8D;IAC9D,4DAA4D;IAC5D,IAAI,mBAAmB,GAAG,CAAC,CAAC;IAC5B,IAAI,UAAU,GAAG,CAAC,EAAE,CAAC;QACnB,mBAAmB,GAAG,IAAI,CAAC,IAAI,CAAC,QAAQ,GAAG,UAAU,CAAC,CAAC;IACzD,CAAC;IAED,IAAI,gBAAgB,GAAkB,IAAI,CAAC;IAC3C,IAAI,cAAc,GAAkB,IAAI,CAAC;IACzC,IAAI,cAAc,KAAK,SAAS,EAAE,CAAC;QACjC,gBAAgB,GAAG,CAAC,SAAS,GAAG,cAAc,CAAC,GAAG,IAAI,CAAC,CAAC,YAAY;QACpE,IAAI,UAAU,KAAK,SAAS,EAAE,CAAC;YAC7B,cAAc,GAAG,gBAAgB,GAAG,UAAU,CAAC;QACjD,CAAC;IACH,CAAC;IAED,OAAO;QACL,WAAW,EAAE,OAAO,CAAC,WAAW,EAAE,CAAC,CAAC;QACpC,SAAS,EAAE,OAAO,CAAC,SAAS,EAAE,CAAC,CAAC;QAChC,cAAc,EAAE,OAAO,CAAC,cAAc,EAAE,CAAC,CAAC;QAC1C,gBAAgB,EAAE,OAAO,CAAC,gBAAgB,EAAE,CAAC,CAAC;QAC9C,UAAU,EAAE,OAAO,CAAC,UAAU,EAAE,CAAC,CAAC;QAClC,mBAAmB,EAAE,OAAO,CAAC,mBAAmB,EAAE,CAAC,CAAC;QACpD,gBAAgB,EAAE,gBAAgB,KAAK,IAAI,CAAC,CAAC,CAAC,OAAO,CAAC,gBAAgB,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI;QACjF,cAAc,EAAE,cAAc,KAAK,IAAI,CAAC,CAAC,CAAC,OAAO,CAAC,cAAc,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI;KAC5E,CAAC;AACJ,CAAC"}