iapws-if97 2.1.0 → 2.1.2

package/dist/saturation/common.js

@@ -1,11 +1,12 @@
-import { Region } from '../types.js';
-import { Pc, R2_T_MIN, RHOc } from '../constants.js';
+import { IF97Error, Region } from '../types.js';
+import { Pc, R2_T_MIN, RHOc, Tt } from '../constants.js';
 import { region1 } from '../regions/region1.js';
 import { region2 } from '../regions/region2.js';
 import { region3ByRhoT } from '../regions/region3.js';
 import { saturationPressure, saturationTemperature } from '../regions/region4.js';
 import { region3SatVolume } from '../regions/region3-subregions.js';
 import { newtonRaphson } from '../solvers/newton-raphson.js';
+import { normalizeRegion4Pressure, normalizeRegion4Temperature, } from './region4-boundaries.js';
 /** Tag a single-phase state as Region 4 with the given vapour quality. */
 function withRegion4Metadata(state, quality) {
     return { ...state, region: Region.Region4, quality };
@@ -19,16 +20,36 @@ function withRegion4Metadata(state, quality) {
  */
 function solveR3Density(p, T, x) {
     const v0 = region3SatVolume(p, T, x);
-    return newtonRaphson((rho) => region3ByRhoT(rho, T).pressure - p, 1 / v0);
+    return newtonRaphson((rho) => {
+        if (!Number.isFinite(rho) || rho <= 0) {
+            return Number.NaN;
+        }
+        return region3ByRhoT(rho, T).pressure - p;
+    }, 1 / v0);
 }
 /** Clamp vapour quality to [0, 1], snapping near-zero/one values to the boundary. */
 function clampQuality(x) {
+    if (!Number.isFinite(x))
+        return Number.NaN;
     if (x <= 1e-12)
         return 0;
     if (x >= 1 - 1e-12)
         return 1;
     return Math.max(0, Math.min(1, x));
 }
+function qualitySpanTolerance(liquidValue, vaporValue, mixtureValue) {
+    return 1e-12 * Math.max(1, Math.abs(liquidValue), Math.abs(vaporValue), Math.abs(mixtureValue));
+}
+export function rawQualityFromSaturationProperty(liquidValue, vaporValue, mixtureValue) {
+    if (!Number.isFinite(liquidValue) || !Number.isFinite(vaporValue) || !Number.isFinite(mixtureValue)) {
+        return Number.NaN;
+    }
+    const span = vaporValue - liquidValue;
+    if (Math.abs(span) <= qualitySpanTolerance(liquidValue, vaporValue, mixtureValue)) {
+        return Number.NaN;
+    }
+    return (mixtureValue - liquidValue) / span;
+}
 /**
  * Compute saturated liquid and vapour properties at a given pressure.
  *
@@ -38,24 +59,25 @@ function clampQuality(x) {
  * @param p - Saturation pressure [MPa]
  */
 export function saturationEndpointsAtPressure(p) {
-    const temperature = saturationTemperature(p);
-    if (Math.abs(p - Pc) < 1e-10) {
+    const pressure = normalizeRegion4Pressure(p);
+    const temperature = normalizeRegion4Temperature(saturationTemperature(pressure));
+    if (pressure === Pc) {
         const state = region3ByRhoT(RHOc, temperature);
-        return { pressure: p, temperature, liquid: state, vapor: state };
+        return { pressure, temperature, liquid: state, vapor: state };
     }
     if (temperature > R2_T_MIN) {
         return {
-            pressure: p,
+            pressure,
             temperature,
-            liquid: region3ByRhoT(solveR3Density(p, temperature, 0), temperature),
-            vapor: region3ByRhoT(solveR3Density(p, temperature, 1), temperature),
+            liquid: region3ByRhoT(solveR3Density(pressure, temperature, 0), temperature),
+            vapor: region3ByRhoT(solveR3Density(pressure, temperature, 1), temperature),
         };
     }
     return {
-        pressure: p,
+        pressure,
         temperature,
-        liquid: region1(p, temperature),
-        vapor: region2(p, temperature),
+        liquid: region1(pressure, temperature),
+        vapor: region2(pressure, temperature),
     };
 }
 /**
@@ -63,7 +85,17 @@ export function saturationEndpointsAtPressure(p) {
  * @param T - Saturation temperature [K]
  */
 export function saturationEndpointsAtTemperature(T) {
-    return saturationEndpointsAtPressure(saturationPressure(T));
+    const temperature = normalizeRegion4Temperature(T);
+    const pressure = normalizeRegion4Pressure(saturationPressure(temperature));
+    if (temperature < Tt) {
+        return {
+            pressure,
+            temperature,
+            liquid: region1(pressure, temperature),
+            vapor: region2(pressure, temperature),
+        };
+    }
+    return saturationEndpointsAtPressure(pressure);
 }
 /**
  * Compute vapour quality from any extensive saturation property.
@@ -75,7 +107,11 @@ export function saturationEndpointsAtTemperature(T) {
  * @param mixtureValue - Property value of the two-phase mixture
  */
 export function qualityFromSaturationProperty(liquidValue, vaporValue, mixtureValue) {
-    return clampQuality((mixtureValue - liquidValue) / (vaporValue - liquidValue));
+    const quality = rawQualityFromSaturationProperty(liquidValue, vaporValue, mixtureValue);
+    if (!Number.isFinite(quality)) {
+        throw new IF97Error('Cannot determine vapor quality from degenerate saturation endpoints');
+    }
+    return clampQuality(quality);
 }
 /**
  * Linearly mix saturated liquid and vapour properties at a given quality.
@@ -88,6 +124,9 @@ export function qualityFromSaturationProperty(liquidValue, vaporValue, mixtureVa
  */
 export function mixSaturationState(endpoints, qualityInput) {
     const quality = clampQuality(qualityInput);
+    if (!Number.isFinite(quality)) {
+        throw new IF97Error('Region 4 quality must be finite');
+    }
     const { liquid, vapor, pressure, temperature } = endpoints;
     if (quality === 0) {
         return withRegion4Metadata(liquid, 0);

package/dist/saturation/common.js.map

@@ -1 +1 @@
-{"version":3,"file":"common.js","sourceRoot":"","sources":["../../src/saturation/common.ts"],"names":[],"mappings":"AAQA,OAAO,EAAE,MAAM,EAAE,MAAM,aAAa,CAAC;AACrC,OAAO,EAAE,EAAE,EAAE,QAAQ,EAAE,IAAI,EAAE,MAAM,iBAAiB,CAAC;AACrD,OAAO,EAAE,OAAO,EAAE,MAAM,uBAAuB,CAAC;AAChD,OAAO,EAAE,OAAO,EAAE,MAAM,uBAAuB,CAAC;AAChD,OAAO,EAAE,aAAa,EAAE,MAAM,uBAAuB,CAAC;AACtD,OAAO,EAAE,kBAAkB,EAAE,qBAAqB,EAAE,MAAM,uBAAuB,CAAC;AAClF,OAAO,EAAE,gBAAgB,EAAE,MAAM,kCAAkC,CAAC;AACpE,OAAO,EAAE,aAAa,EAAE,MAAM,8BAA8B,CAAC;AAU7D,0EAA0E;AAC1E,SAAS,mBAAmB,CAAC,KAAsB,EAAE,OAAe;IAClE,OAAO,EAAE,GAAG,KAAK,EAAE,MAAM,EAAE,MAAM,CAAC,OAAO,EAAE,OAAO,EAAE,CAAC;AACvD,CAAC;AAED;;;;;;GAMG;AACH,SAAS,cAAc,CAAC,CAAS,EAAE,CAAS,EAAE,CAAQ;IACpD,MAAM,EAAE,GAAG,gBAAgB,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IACrC,OAAO,aAAa,CAClB,CAAC,GAAG,EAAE,EAAE,CAAC,aAAa,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,QAAQ,GAAG,CAAC,EAC3C,CAAC,GAAG,EAAE,CACP,CAAC;AACJ,CAAC;AAED,qFAAqF;AACrF,SAAS,YAAY,CAAC,CAAS;IAC7B,IAAI,CAAC,IAAI,KAAK;QAAE,OAAO,CAAC,CAAC;IACzB,IAAI,CAAC,IAAI,CAAC,GAAG,KAAK;QAAE,OAAO,CAAC,CAAC;IAC7B,OAAO,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;AACrC,CAAC;AAED;;;;;;;GAOG;AACH,MAAM,UAAU,6BAA6B,CAAC,CAAS;IACrD,MAAM,WAAW,GAAG,qBAAqB,CAAC,CAAC,CAAC,CAAC;IAE7C,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,KAAK,EAAE,CAAC;QAC7B,MAAM,KAAK,GAAG,aAAa,CAAC,IAAI,EAAE,WAAW,CAAC,CAAC;QAC/C,OAAO,EAAE,QAAQ,EAAE,CAAC,EAAE,WAAW,EAAE,MAAM,EAAE,KAAK,EAAE,KAAK,EAAE,KAAK,EAAE,CAAC;IACnE,CAAC;IAED,IAAI,WAAW,GAAG,QAAQ,EAAE,CAAC;QAC3B,OAAO;YACL,QAAQ,EAAE,CAAC;YACX,WAAW;YACX,MAAM,EAAE,aAAa,CAAC,cAAc,CAAC,CAAC,EAAE,WAAW,EAAE,CAAC,CAAC,EAAE,WAAW,CAAC;YACrE,KAAK,EAAE,aAAa,CAAC,cAAc,CAAC,CAAC,EAAE,WAAW,EAAE,CAAC,CAAC,EAAE,WAAW,CAAC;SACrE,CAAC;IACJ,CAAC;IAED,OAAO;QACL,QAAQ,EAAE,CAAC;QACX,WAAW;QACX,MAAM,EAAE,OAAO,CAAC,CAAC,EAAE,WAAW,CAAC;QAC/B,KAAK,EAAE,OAAO,CAAC,CAAC,EAAE,WAAW,CAAC;KAC/B,CAAC;AACJ,CAAC;AAED;;;GAGG;AACH,MAAM,UAAU,gCAAgC,CAAC,CAAS;IACxD,OAAO,6BAA6B,CAAC,kBAAkB,CAAC,CAAC,CAAC,CAAC,CAAC;AAC9D,CAAC;AAED;;;;;;;;GAQG;AACH,MAAM,UAAU,6BAA6B,CAC3C,WAAmB,EACnB,UAAkB,EAClB,YAAoB;IAEpB,OAAO,YAAY,CAAC,CAAC,YAAY,GAAG,WAAW,CAAC,GAAG,CAAC,UAAU,GAAG,WAAW,CAAC,CAAC,CAAC;AACjF,CAAC;AAED;;;;;;;;GAQG;AACH,MAAM,UAAU,kBAAkB,CAChC,SAA8B,EAC9B,YAAoB;IAEpB,MAAM,OAAO,GAAG,YAAY,CAAC,YAAY,CAAC,CAAC;IAC3C,MAAM,EAAE,MAAM,EAAE,KAAK,EAAE,QAAQ,EAAE,WAAW,EAAE,GAAG,SAAS,CAAC;IAE3D,IAAI,OAAO,KAAK,CAAC,EAAE,CAAC;QAClB,OAAO,mBAAmB,CAAC,MAAM,EAAE,CAAC,CAAC,CAAC;IACxC,CAAC;IAED,IAAI,OAAO,KAAK,CAAC,EAAE,CAAC;QAClB,OAAO,mBAAmB,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC;IACvC,CAAC;IAED,OAAO;QACL,MAAM,EAAE,MAAM,CAAC,OAAO;QACtB,QAAQ;QACR,WAAW;QACX,cAAc,EAAE,MAAM,CAAC,cAAc,GAAG,CAAC,CAAC,GAAG,OAAO,CAAC,GAAG,KAAK,CAAC,cAAc,GAAG,OAAO;QACtF,cAAc,EAAE,MAAM,CAAC,cAAc,GAAG,CAAC,CAAC,GAAG,OAAO,CAAC,GAAG,KAAK,CAAC,cAAc,GAAG,OAAO;QACtF,OAAO,EAAE,MAAM,CAAC,OAAO,GAAG,CAAC,CAAC,GAAG,OAAO,CAAC,GAAG,KAAK,CAAC,OAAO,GAAG,OAAO;QACjE,QAAQ,EAAE,MAAM,CAAC,QAAQ,GAAG,CAAC,CAAC,GAAG,OAAO,CAAC,GAAG,KAAK,CAAC,QAAQ,GAAG,OAAO;QACpE,EAAE,EAAE,IAAI;QACR,EAAE,EAAE,IAAI;QACR,YAAY,EAAE,IAAI;QAClB,OAAO;QACP,iBAAiB,EAAE,IAAI;QACvB,yBAAyB,EAAE,IAAI;KAChC,CAAC;AACJ,CAAC"}
+{"version":3,"file":"common.js","sourceRoot":"","sources":["../../src/saturation/common.ts"],"names":[],"mappings":"AAQA,OAAO,EAAE,SAAS,EAAE,MAAM,EAAE,MAAM,aAAa,CAAC;AAChD,OAAO,EAAE,EAAE,EAAE,QAAQ,EAAE,IAAI,EAAE,EAAE,EAAE,MAAM,iBAAiB,CAAC;AACzD,OAAO,EAAE,OAAO,EAAE,MAAM,uBAAuB,CAAC;AAChD,OAAO,EAAE,OAAO,EAAE,MAAM,uBAAuB,CAAC;AAChD,OAAO,EAAE,aAAa,EAAE,MAAM,uBAAuB,CAAC;AACtD,OAAO,EAAE,kBAAkB,EAAE,qBAAqB,EAAE,MAAM,uBAAuB,CAAC;AAClF,OAAO,EAAE,gBAAgB,EAAE,MAAM,kCAAkC,CAAC;AACpE,OAAO,EAAE,aAAa,EAAE,MAAM,8BAA8B,CAAC;AAC7D,OAAO,EACL,wBAAwB,EACxB,2BAA2B,GAC5B,MAAM,yBAAyB,CAAC;AAUjC,0EAA0E;AAC1E,SAAS,mBAAmB,CAAC,KAAsB,EAAE,OAAe;IAClE,OAAO,EAAE,GAAG,KAAK,EAAE,MAAM,EAAE,MAAM,CAAC,OAAO,EAAE,OAAO,EAAE,CAAC;AACvD,CAAC;AAED;;;;;;GAMG;AACH,SAAS,cAAc,CAAC,CAAS,EAAE,CAAS,EAAE,CAAQ;IACpD,MAAM,EAAE,GAAG,gBAAgB,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IACrC,OAAO,aAAa,CAClB,CAAC,GAAG,EAAE,EAAE;QACN,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC,EAAE,CAAC;YACtC,OAAO,MAAM,CAAC,GAAG,CAAC;QACpB,CAAC;QACD,OAAO,aAAa,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,QAAQ,GAAG,CAAC,CAAC;IAC5C,CAAC,EACD,CAAC,GAAG,EAAE,CACP,CAAC;AACJ,CAAC;AAED,qFAAqF;AACrF,SAAS,YAAY,CAAC,CAAS;IAC7B,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC;QAAE,OAAO,MAAM,CAAC,GAAG,CAAC;IAC3C,IAAI,CAAC,IAAI,KAAK;QAAE,OAAO,CAAC,CAAC;IACzB,IAAI,CAAC,IAAI,CAAC,GAAG,KAAK;QAAE,OAAO,CAAC,CAAC;IAC7B,OAAO,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;AACrC,CAAC;AAED,SAAS,oBAAoB,CAC3B,WAAmB,EACnB,UAAkB,EAClB,YAAoB;IAEpB,OAAO,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,WAAW,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,UAAU,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,YAAY,CAAC,CAAC,CAAC;AAClG,CAAC;AAED,MAAM,UAAU,gCAAgC,CAC9C,WAAmB,EACnB,UAAkB,EAClB,YAAoB;IAEpB,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,WAAW,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,UAAU,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,YAAY,CAAC,EAAE,CAAC;QACpG,OAAO,MAAM,CAAC,GAAG,CAAC;IACpB,CAAC;IAED,MAAM,IAAI,GAAG,UAAU,GAAG,WAAW,CAAC;IACtC,IAAI,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,IAAI,oBAAoB,CAAC,WAAW,EAAE,UAAU,EAAE,YAAY,CAAC,EAAE,CAAC;QAClF,OAAO,MAAM,CAAC,GAAG,CAAC;IACpB,CAAC;IAED,OAAO,CAAC,YAAY,GAAG,WAAW,CAAC,GAAG,IAAI,CAAC;AAC7C,CAAC;AAED;;;;;;;GAOG;AACH,MAAM,UAAU,6BAA6B,CAAC,CAAS;IACrD,MAAM,QAAQ,GAAG,wBAAwB,CAAC,CAAC,CAAC,CAAC;IAC7C,MAAM,WAAW,GAAG,2BAA2B,CAAC,qBAAqB,CAAC,QAAQ,CAAC,CAAC,CAAC;IAEjF,IAAI,QAAQ,KAAK,EAAE,EAAE,CAAC;QACpB,MAAM,KAAK,GAAG,aAAa,CAAC,IAAI,EAAE,WAAW,CAAC,CAAC;QAC/C,OAAO,EAAE,QAAQ,EAAE,WAAW,EAAE,MAAM,EAAE,KAAK,EAAE,KAAK,EAAE,KAAK,EAAE,CAAC;IAChE,CAAC;IAED,IAAI,WAAW,GAAG,QAAQ,EAAE,CAAC;QAC3B,OAAO;YACL,QAAQ;YACR,WAAW;YACX,MAAM,EAAE,aAAa,CAAC,cAAc,CAAC,QAAQ,EAAE,WAAW,EAAE,CAAC,CAAC,EAAE,WAAW,CAAC;YAC5E,KAAK,EAAE,aAAa,CAAC,cAAc,CAAC,QAAQ,EAAE,WAAW,EAAE,CAAC,CAAC,EAAE,WAAW,CAAC;SAC5E,CAAC;IACJ,CAAC;IAED,OAAO;QACL,QAAQ;QACR,WAAW;QACX,MAAM,EAAE,OAAO,CAAC,QAAQ,EAAE,WAAW,CAAC;QACtC,KAAK,EAAE,OAAO,CAAC,QAAQ,EAAE,WAAW,CAAC;KACtC,CAAC;AACJ,CAAC;AAED;;;GAGG;AACH,MAAM,UAAU,gCAAgC,CAAC,CAAS;IACxD,MAAM,WAAW,GAAG,2BAA2B,CAAC,CAAC,CAAC,CAAC;IACnD,MAAM,QAAQ,GAAG,wBAAwB,CAAC,kBAAkB,CAAC,WAAW,CAAC,CAAC,CAAC;IAE3E,IAAI,WAAW,GAAG,EAAE,EAAE,CAAC;QACrB,OAAO;YACL,QAAQ;YACR,WAAW;YACX,MAAM,EAAE,OAAO,CAAC,QAAQ,EAAE,WAAW,CAAC;YACtC,KAAK,EAAE,OAAO,CAAC,QAAQ,EAAE,WAAW,CAAC;SACtC,CAAC;IACJ,CAAC;IAED,OAAO,6BAA6B,CAAC,QAAQ,CAAC,CAAC;AACjD,CAAC;AAED;;;;;;;;GAQG;AACH,MAAM,UAAU,6BAA6B,CAC3C,WAAmB,EACnB,UAAkB,EAClB,YAAoB;IAEpB,MAAM,OAAO,GAAG,gCAAgC,CAAC,WAAW,EAAE,UAAU,EAAE,YAAY,CAAC,CAAC;IACxF,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,OAAO,CAAC,EAAE,CAAC;QAC9B,MAAM,IAAI,SAAS,CAAC,qEAAqE,CAAC,CAAC;IAC7F,CAAC;IACD,OAAO,YAAY,CAAC,OAAO,CAAC,CAAC;AAC/B,CAAC;AAED;;;;;;;;GAQG;AACH,MAAM,UAAU,kBAAkB,CAChC,SAA8B,EAC9B,YAAoB;IAEpB,MAAM,OAAO,GAAG,YAAY,CAAC,YAAY,CAAC,CAAC;IAC3C,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,OAAO,CAAC,EAAE,CAAC;QAC9B,MAAM,IAAI,SAAS,CAAC,iCAAiC,CAAC,CAAC;IACzD,CAAC;IACD,MAAM,EAAE,MAAM,EAAE,KAAK,EAAE,QAAQ,EAAE,WAAW,EAAE,GAAG,SAAS,CAAC;IAE3D,IAAI,OAAO,KAAK,CAAC,EAAE,CAAC;QAClB,OAAO,mBAAmB,CAAC,MAAM,EAAE,CAAC,CAAC,CAAC;IACxC,CAAC;IAED,IAAI,OAAO,KAAK,CAAC,EAAE,CAAC;QAClB,OAAO,mBAAmB,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC;IACvC,CAAC;IAED,OAAO;QACL,MAAM,EAAE,MAAM,CAAC,OAAO;QACtB,QAAQ;QACR,WAAW;QACX,cAAc,EAAE,MAAM,CAAC,cAAc,GAAG,CAAC,CAAC,GAAG,OAAO,CAAC,GAAG,KAAK,CAAC,cAAc,GAAG,OAAO;QACtF,cAAc,EAAE,MAAM,CAAC,cAAc,GAAG,CAAC,CAAC,GAAG,OAAO,CAAC,GAAG,KAAK,CAAC,cAAc,GAAG,OAAO;QACtF,OAAO,EAAE,MAAM,CAAC,OAAO,GAAG,CAAC,CAAC,GAAG,OAAO,CAAC,GAAG,KAAK,CAAC,OAAO,GAAG,OAAO;QACjE,QAAQ,EAAE,MAAM,CAAC,QAAQ,GAAG,CAAC,CAAC,GAAG,OAAO,CAAC,GAAG,KAAK,CAAC,QAAQ,GAAG,OAAO;QACpE,EAAE,EAAE,IAAI;QACR,EAAE,EAAE,IAAI;QACR,YAAY,EAAE,IAAI;QAClB,OAAO;QACP,iBAAiB,EAAE,IAAI;QACvB,yBAAyB,EAAE,IAAI;KAChC,CAAC;AACJ,CAAC"}

package/dist/saturation/region4-boundaries.d.ts

@@ -0,0 +1,14 @@
+export declare function normalizeRegion4Pressure(p: number): number;
+export declare function normalizeRegion4Temperature(T: number): number;
+export declare function clampRegion4TemperatureBelowCritical(T: number): number;
+export declare function isRegion4CriticalPressure(p: number): boolean;
+export declare function isRegion4CriticalTemperature(T: number): boolean;
+export declare function assertRegion4PressureStateInput(p: number, solverName: string): number;
+export declare function assertRegion4TemperatureStateInput(T: number, solverName: string): number;
+export declare function assertRegion4StateAllowed(pressure: number, temperature: number, solverName: string): void;
+export declare function isCriticalRegion4Enthalpy(h: number): boolean;
+export declare function isCriticalRegion4Entropy(s: number): boolean;
+export declare function assertCriticalRegion4PHInput(p: number, h: number, solverName: string): void;
+export declare function assertCriticalRegion4PSInput(p: number, s: number, solverName: string): void;
+export declare function assertCriticalRegion4HSInput(h: number, s: number, solverName: string): void;
+//# sourceMappingURL=region4-boundaries.d.ts.map

package/dist/saturation/region4-boundaries.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"region4-boundaries.d.ts","sourceRoot":"","sources":["../../src/saturation/region4-boundaries.ts"],"names":[],"mappings":"AAGA,wBAAgB,wBAAwB,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM,CAK1D;AAED,wBAAgB,2BAA2B,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM,CAK7D;AAED,wBAAgB,oCAAoC,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM,CAEtE;AAED,wBAAgB,yBAAyB,CAAC,CAAC,EAAE,MAAM,GAAG,OAAO,CAE5D;AAED,wBAAgB,4BAA4B,CAAC,CAAC,EAAE,MAAM,GAAG,OAAO,CAE/D;AAMD,wBAAgB,+BAA+B,CAAC,CAAC,EAAE,MAAM,EAAE,UAAU,EAAE,MAAM,GAAG,MAAM,CASrF;AAED,wBAAgB,kCAAkC,CAAC,CAAC,EAAE,MAAM,EAAE,UAAU,EAAE,MAAM,GAAG,MAAM,CASxF;AAED,wBAAgB,yBAAyB,CACvC,QAAQ,EAAE,MAAM,EAChB,WAAW,EAAE,MAAM,EACnB,UAAU,EAAE,MAAM,GACjB,IAAI,CAIN;AAED,wBAAgB,yBAAyB,CAAC,CAAC,EAAE,MAAM,GAAG,OAAO,CAE5D;AAED,wBAAgB,wBAAwB,CAAC,CAAC,EAAE,MAAM,GAAG,OAAO,CAE3D;AAED,wBAAgB,4BAA4B,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,UAAU,EAAE,MAAM,GAAG,IAAI,CAI3F;AAED,wBAAgB,4BAA4B,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,UAAU,EAAE,MAAM,GAAG,IAAI,CAI3F;AAED,wBAAgB,4BAA4B,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,UAAU,EAAE,MAAM,GAAG,IAAI,CAI3F"}

package/dist/saturation/region4-boundaries.js

@@ -0,0 +1,79 @@
+import * as C from '../constants.js';
+import { IF97Error, OutOfRangeError } from '../types.js';
+export function normalizeRegion4Pressure(p) {
+    if (Math.abs(p - C.P_MIN) <= C.REGION4_PRESSURE_TOLERANCE)
+        return C.P_MIN;
+    if (Math.abs(p - C.Pt) <= C.REGION4_PRESSURE_TOLERANCE)
+        return C.Pt;
+    if (Math.abs(p - C.Pc) <= C.REGION4_PRESSURE_TOLERANCE)
+        return C.Pc;
+    return p;
+}
+export function normalizeRegion4Temperature(T) {
+    if (Math.abs(T - C.T_MIN) <= C.REGION4_TEMPERATURE_TOLERANCE)
+        return C.T_MIN;
+    if (Math.abs(T - C.Tt) <= C.REGION4_TEMPERATURE_TOLERANCE)
+        return C.Tt;
+    if (Math.abs(T - C.Tc) <= C.REGION4_TEMPERATURE_TOLERANCE)
+        return C.Tc;
+    return T;
+}
+export function clampRegion4TemperatureBelowCritical(T) {
+    return Math.max(C.Tt, Math.min(T, C.Tc - C.REGION4_SUBCRITICAL_TEMPERATURE_MARGIN));
+}
+export function isRegion4CriticalPressure(p) {
+    return normalizeRegion4Pressure(p) === C.Pc;
+}
+export function isRegion4CriticalTemperature(T) {
+    return normalizeRegion4Temperature(T) === C.Tc;
+}
+function criticalPointMessage(solverName) {
+    return `${solverName} does not support the critical point P=${C.Pc} MPa, T=${C.Tc} K as a Region 4 saturation state`;
+}
+export function assertRegion4PressureStateInput(p, solverName) {
+    const pressure = normalizeRegion4Pressure(p);
+    if (pressure < C.Pt || pressure > C.Pc) {
+        throw new OutOfRangeError('Pressure', p, C.Pt, C.Pc);
+    }
+    if (pressure === C.Pc) {
+        throw new IF97Error(criticalPointMessage(solverName));
+    }
+    return pressure;
+}
+export function assertRegion4TemperatureStateInput(T, solverName) {
+    const temperature = normalizeRegion4Temperature(T);
+    if (temperature < C.Tt || temperature > C.Tc) {
+        throw new OutOfRangeError('Temperature', T, C.Tt, C.Tc);
+    }
+    if (temperature === C.Tc) {
+        throw new IF97Error(criticalPointMessage(solverName));
+    }
+    return temperature;
+}
+export function assertRegion4StateAllowed(pressure, temperature, solverName) {
+    if (isRegion4CriticalPressure(pressure) || isRegion4CriticalTemperature(temperature)) {
+        throw new IF97Error(criticalPointMessage(solverName));
+    }
+}
+export function isCriticalRegion4Enthalpy(h) {
+    return Math.abs(h - C.R3_H_CRT) <= C.REGION4_CRITICAL_HS_TOLERANCE;
+}
+export function isCriticalRegion4Entropy(s) {
+    return Math.abs(s - C.R3_S_CRT) <= C.REGION4_CRITICAL_HS_TOLERANCE;
+}
+export function assertCriticalRegion4PHInput(p, h, solverName) {
+    if (isRegion4CriticalPressure(p) && isCriticalRegion4Enthalpy(h)) {
+        throw new IF97Error(criticalPointMessage(solverName));
+    }
+}
+export function assertCriticalRegion4PSInput(p, s, solverName) {
+    if (isRegion4CriticalPressure(p) && isCriticalRegion4Entropy(s)) {
+        throw new IF97Error(criticalPointMessage(solverName));
+    }
+}
+export function assertCriticalRegion4HSInput(h, s, solverName) {
+    if (isCriticalRegion4Enthalpy(h) && isCriticalRegion4Entropy(s)) {
+        throw new IF97Error(criticalPointMessage(solverName));
+    }
+}
+//# sourceMappingURL=region4-boundaries.js.map

package/dist/saturation/region4-boundaries.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"region4-boundaries.js","sourceRoot":"","sources":["../../src/saturation/region4-boundaries.ts"],"names":[],"mappings":"AAAA,OAAO,KAAK,CAAC,MAAM,iBAAiB,CAAC;AACrC,OAAO,EAAE,SAAS,EAAE,eAAe,EAAE,MAAM,aAAa,CAAC;AAEzD,MAAM,UAAU,wBAAwB,CAAC,CAAS;IAChD,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC,0BAA0B;QAAE,OAAO,CAAC,CAAC,KAAK,CAAC;IAC1E,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,IAAI,CAAC,CAAC,0BAA0B;QAAE,OAAO,CAAC,CAAC,EAAE,CAAC;IACpE,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,IAAI,CAAC,CAAC,0BAA0B;QAAE,OAAO,CAAC,CAAC,EAAE,CAAC;IACpE,OAAO,CAAC,CAAC;AACX,CAAC;AAED,MAAM,UAAU,2BAA2B,CAAC,CAAS;IACnD,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC,6BAA6B;QAAE,OAAO,CAAC,CAAC,KAAK,CAAC;IAC7E,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,IAAI,CAAC,CAAC,6BAA6B;QAAE,OAAO,CAAC,CAAC,EAAE,CAAC;IACvE,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,IAAI,CAAC,CAAC,6BAA6B;QAAE,OAAO,CAAC,CAAC,EAAE,CAAC;IACvE,OAAO,CAAC,CAAC;AACX,CAAC;AAED,MAAM,UAAU,oCAAoC,CAAC,CAAS;IAC5D,OAAO,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,sCAAsC,CAAC,CAAC,CAAC;AACtF,CAAC;AAED,MAAM,UAAU,yBAAyB,CAAC,CAAS;IACjD,OAAO,wBAAwB,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC;AAC9C,CAAC;AAED,MAAM,UAAU,4BAA4B,CAAC,CAAS;IACpD,OAAO,2BAA2B,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC;AACjD,CAAC;AAED,SAAS,oBAAoB,CAAC,UAAkB;IAC9C,OAAO,GAAG,UAAU,0CAA0C,CAAC,CAAC,EAAE,WAAW,CAAC,CAAC,EAAE,mCAAmC,CAAC;AACvH,CAAC;AAED,MAAM,UAAU,+BAA+B,CAAC,CAAS,EAAE,UAAkB;IAC3E,MAAM,QAAQ,GAAG,wBAAwB,CAAC,CAAC,CAAC,CAAC;IAC7C,IAAI,QAAQ,GAAG,CAAC,CAAC,EAAE,IAAI,QAAQ,GAAG,CAAC,CAAC,EAAE,EAAE,CAAC;QACvC,MAAM,IAAI,eAAe,CAAC,UAAU,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC;IACvD,CAAC;IACD,IAAI,QAAQ,KAAK,CAAC,CAAC,EAAE,EAAE,CAAC;QACtB,MAAM,IAAI,SAAS,CAAC,oBAAoB,CAAC,UAAU,CAAC,CAAC,CAAC;IACxD,CAAC;IACD,OAAO,QAAQ,CAAC;AAClB,CAAC;AAED,MAAM,UAAU,kCAAkC,CAAC,CAAS,EAAE,UAAkB;IAC9E,MAAM,WAAW,GAAG,2BAA2B,CAAC,CAAC,CAAC,CAAC;IACnD,IAAI,WAAW,GAAG,CAAC,CAAC,EAAE,IAAI,WAAW,GAAG,CAAC,CAAC,EAAE,EAAE,CAAC;QAC7C,MAAM,IAAI,eAAe,CAAC,aAAa,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC;IAC1D,CAAC;IACD,IAAI,WAAW,KAAK,CAAC,CAAC,EAAE,EAAE,CAAC;QACzB,MAAM,IAAI,SAAS,CAAC,oBAAoB,CAAC,UAAU,CAAC,CAAC,CAAC;IACxD,CAAC;IACD,OAAO,WAAW,CAAC;AACrB,CAAC;AAED,MAAM,UAAU,yBAAyB,CACvC,QAAgB,EAChB,WAAmB,EACnB,UAAkB;IAElB,IAAI,yBAAyB,CAAC,QAAQ,CAAC,IAAI,4BAA4B,CAAC,WAAW,CAAC,EAAE,CAAC;QACrF,MAAM,IAAI,SAAS,CAAC,oBAAoB,CAAC,UAAU,CAAC,CAAC,CAAC;IACxD,CAAC;AACH,CAAC;AAED,MAAM,UAAU,yBAAyB,CAAC,CAAS;IACjD,OAAO,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC,6BAA6B,CAAC;AACrE,CAAC;AAED,MAAM,UAAU,wBAAwB,CAAC,CAAS;IAChD,OAAO,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC,6BAA6B,CAAC;AACrE,CAAC;AAED,MAAM,UAAU,4BAA4B,CAAC,CAAS,EAAE,CAAS,EAAE,UAAkB;IACnF,IAAI,yBAAyB,CAAC,CAAC,CAAC,IAAI,yBAAyB,CAAC,CAAC,CAAC,EAAE,CAAC;QACjE,MAAM,IAAI,SAAS,CAAC,oBAAoB,CAAC,UAAU,CAAC,CAAC,CAAC;IACxD,CAAC;AACH,CAAC;AAED,MAAM,UAAU,4BAA4B,CAAC,CAAS,EAAE,CAAS,EAAE,UAAkB;IACnF,IAAI,yBAAyB,CAAC,CAAC,CAAC,IAAI,wBAAwB,CAAC,CAAC,CAAC,EAAE,CAAC;QAChE,MAAM,IAAI,SAAS,CAAC,oBAAoB,CAAC,UAAU,CAAC,CAAC,CAAC;IACxD,CAAC;AACH,CAAC;AAED,MAAM,UAAU,4BAA4B,CAAC,CAAS,EAAE,CAAS,EAAE,UAAkB;IACnF,IAAI,yBAAyB,CAAC,CAAC,CAAC,IAAI,wBAAwB,CAAC,CAAC,CAAC,EAAE,CAAC;QAChE,MAAM,IAAI,SAAS,CAAC,oBAAoB,CAAC,UAAU,CAAC,CAAC,CAAC;IACxD,CAAC;AACH,CAAC"}

package/dist/saturation/region4-hs.d.ts

@@ -0,0 +1,3 @@
+import type { BasicProperties } from '../types.js';
+export declare function tryRegion4HSState(h: number, s: number): BasicProperties | null;
+//# sourceMappingURL=region4-hs.d.ts.map

package/dist/saturation/region4-hs.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"region4-hs.d.ts","sourceRoot":"","sources":["../../src/saturation/region4-hs.ts"],"names":[],"mappings":"AAEA,OAAO,KAAK,EAAE,eAAe,EAAE,MAAM,aAAa,CAAC;AAmMnD,wBAAgB,iBAAiB,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,eAAe,GAAG,IAAI,CAuD9E"}

package/dist/saturation/region4-hs.js

@@ -0,0 +1,165 @@
+import * as C from '../constants.js';
+import { bracketedNewton } from '../solvers/bracketed-newton.js';
+import { mixSaturationState, rawQualityFromSaturationProperty, saturationEndpointsAtTemperature, } from './common.js';
+import { assertRegion4StateAllowed, clampRegion4TemperatureBelowCritical, } from './region4-boundaries.js';
+const REGION4_QUALITY_TOLERANCE = 1e-6;
+function isClose(a, b) {
+    return Math.abs(a - b) <= 1e-9 * Math.max(1, Math.abs(a), Math.abs(b));
+}
+function normalizeEndpointQuality(quality) {
+    if (!Number.isFinite(quality)) {
+        return Number.NaN;
+    }
+    if (Math.abs(quality) <= REGION4_QUALITY_TOLERANCE) {
+        return 0;
+    }
+    if (Math.abs(quality - 1) <= REGION4_QUALITY_TOLERANCE) {
+        return 1;
+    }
+    return quality;
+}
+function isAdmissibleQuality(quality) {
+    return Number.isFinite(quality)
+        && quality >= -REGION4_QUALITY_TOLERANCE
+        && quality <= 1 + REGION4_QUALITY_TOLERANCE;
+}
+function entropyResidual(h, s, temperatureGuess) {
+    const endpoints = saturationEndpointsAtTemperature(clampRegion4TemperatureBelowCritical(temperatureGuess));
+    const quality = rawQualityFromSaturationProperty(endpoints.liquid.enthalpy, endpoints.vapor.enthalpy, h);
+    if (!Number.isFinite(quality)) {
+        return Number.NaN;
+    }
+    return endpoints.liquid.entropy * (1 - quality) + endpoints.vapor.entropy * quality - s;
+}
+function finalizeRegion4HSState(h, s, temperature) {
+    const endpoints = saturationEndpointsAtTemperature(temperature);
+    const qualityFromEnthalpy = rawQualityFromSaturationProperty(endpoints.liquid.enthalpy, endpoints.vapor.enthalpy, h);
+    const qualityFromEntropy = rawQualityFromSaturationProperty(endpoints.liquid.entropy, endpoints.vapor.entropy, s);
+    const normalizedQualityFromEnthalpy = normalizeEndpointQuality(qualityFromEnthalpy);
+    const normalizedQualityFromEntropy = normalizeEndpointQuality(qualityFromEntropy);
+    if (!isAdmissibleQuality(normalizedQualityFromEnthalpy) || !isAdmissibleQuality(normalizedQualityFromEntropy)) {
+        return null;
+    }
+    if (Math.abs(normalizedQualityFromEnthalpy - normalizedQualityFromEntropy) > REGION4_QUALITY_TOLERANCE) {
+        return null;
+    }
+    if (Math.abs(entropyResidual(h, s, temperature)) > C.REGION4_HS_RESIDUAL_TOLERANCE) {
+        return null;
+    }
+    try {
+        assertRegion4StateAllowed(endpoints.pressure, endpoints.temperature, 'Region4HS');
+    }
+    catch {
+        return null;
+    }
+    return mixSaturationState(endpoints, normalizedQualityFromEnthalpy);
+}
+function endpointStateForPhase(h, s, temperature, phase) {
+    let endpoints;
+    try {
+        endpoints = saturationEndpointsAtTemperature(temperature);
+    }
+    catch {
+        return null;
+    }
+    const endpoint = endpoints[phase];
+    if (!isClose(endpoint.enthalpy, h) || !isClose(endpoint.entropy, s)) {
+        return null;
+    }
+    try {
+        assertRegion4StateAllowed(endpoints.pressure, endpoints.temperature, 'Region4HS');
+    }
+    catch {
+        return null;
+    }
+    return mixSaturationState(endpoints, phase === 'liquid' ? 0 : 1);
+}
+function endpointEnthalpyResidual(h, temperatureGuess, phase) {
+    let endpoints;
+    try {
+        endpoints = saturationEndpointsAtTemperature(clampRegion4TemperatureBelowCritical(temperatureGuess));
+    }
+    catch {
+        return Number.NaN;
+    }
+    return endpoints[phase].enthalpy - h;
+}
+function tryRegion4HSEndpointState(h, s) {
+    const lower = C.Tt;
+    const upper = clampRegion4TemperatureBelowCritical(C.Tc);
+    for (const phase of ['liquid', 'vapor']) {
+        const lowerState = endpointStateForPhase(h, s, lower, phase);
+        if (lowerState !== null) {
+            return lowerState;
+        }
+        const upperState = endpointStateForPhase(h, s, upper, phase);
+        if (upperState !== null) {
+            return upperState;
+        }
+        let previousTemperature = lower;
+        let previousResidual = endpointEnthalpyResidual(h, lower, phase);
+        for (let i = 1; i <= C.REGION4_HS_BRACKET_SEGMENTS; i++) {
+            const temperature = lower + (upper - lower) * (i / C.REGION4_HS_BRACKET_SEGMENTS);
+            const residual = endpointEnthalpyResidual(h, temperature, phase);
+            if (!Number.isFinite(previousResidual) || !Number.isFinite(residual)) {
+                previousTemperature = temperature;
+                previousResidual = residual;
+                continue;
+            }
+            if (isClose(previousResidual, 0)) {
+                return endpointStateForPhase(h, s, previousTemperature, phase);
+            }
+            if (isClose(residual, 0)) {
+                return endpointStateForPhase(h, s, temperature, phase);
+            }
+            if (previousResidual * residual < 0) {
+                const root = bracketedNewton((candidateTemperature) => endpointEnthalpyResidual(h, candidateTemperature, phase), previousTemperature, temperature, (previousTemperature + temperature) / 2, { tolerance: C.REGION4_HS_RESIDUAL_TOLERANCE });
+                return endpointStateForPhase(h, s, root, phase);
+            }
+            previousTemperature = temperature;
+            previousResidual = residual;
+        }
+    }
+    return null;
+}
+export function tryRegion4HSState(h, s) {
+    const endpointState = tryRegion4HSEndpointState(h, s);
+    if (endpointState !== null) {
+        return endpointState;
+    }
+    const lower = C.Tt;
+    const upper = clampRegion4TemperatureBelowCritical(C.Tc);
+    const lowerState = finalizeRegion4HSState(h, s, lower);
+    if (lowerState !== null) {
+        return lowerState;
+    }
+    const upperState = finalizeRegion4HSState(h, s, upper);
+    if (upperState !== null) {
+        return upperState;
+    }
+    let previousTemperature = lower;
+    let previousResidual = entropyResidual(h, s, lower);
+    for (let i = 1; i <= C.REGION4_HS_BRACKET_SEGMENTS; i++) {
+        const temperature = lower + (upper - lower) * (i / C.REGION4_HS_BRACKET_SEGMENTS);
+        const residual = entropyResidual(h, s, temperature);
+        if (!Number.isFinite(previousResidual) || !Number.isFinite(residual)) {
+            previousTemperature = temperature;
+            previousResidual = residual;
+            continue;
+        }
+        if (isClose(previousResidual, 0)) {
+            return finalizeRegion4HSState(h, s, previousTemperature);
+        }
+        if (isClose(residual, 0)) {
+            return finalizeRegion4HSState(h, s, temperature);
+        }
+        if (previousResidual * residual < 0) {
+            const root = bracketedNewton((candidateTemperature) => entropyResidual(h, s, candidateTemperature), previousTemperature, temperature, (previousTemperature + temperature) / 2, { tolerance: C.REGION4_HS_RESIDUAL_TOLERANCE });
+            return finalizeRegion4HSState(h, s, root);
+        }
+        previousTemperature = temperature;
+        previousResidual = residual;
+    }
+    return null;
+}
+//# sourceMappingURL=region4-hs.js.map

package/dist/saturation/region4-hs.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"region4-hs.js","sourceRoot":"","sources":["../../src/saturation/region4-hs.ts"],"names":[],"mappings":"AAAA,OAAO,KAAK,CAAC,MAAM,iBAAiB,CAAC;AACrC,OAAO,EAAE,eAAe,EAAE,MAAM,gCAAgC,CAAC;AAEjE,OAAO,EACL,kBAAkB,EAClB,gCAAgC,EAChC,gCAAgC,GACjC,MAAM,aAAa,CAAC;AACrB,OAAO,EACL,yBAAyB,EACzB,oCAAoC,GACrC,MAAM,yBAAyB,CAAC;AAEjC,MAAM,yBAAyB,GAAG,IAAI,CAAC;AAGvC,SAAS,OAAO,CAAC,CAAS,EAAE,CAAS;IACnC,OAAO,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;AACzE,CAAC;AAED,SAAS,wBAAwB,CAAC,OAAe;IAC/C,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,OAAO,CAAC,EAAE,CAAC;QAC9B,OAAO,MAAM,CAAC,GAAG,CAAC;IACpB,CAAC;IACD,IAAI,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,IAAI,yBAAyB,EAAE,CAAC;QACnD,OAAO,CAAC,CAAC;IACX,CAAC;IACD,IAAI,IAAI,CAAC,GAAG,CAAC,OAAO,GAAG,CAAC,CAAC,IAAI,yBAAyB,EAAE,CAAC;QACvD,OAAO,CAAC,CAAC;IACX,CAAC;IACD,OAAO,OAAO,CAAC;AACjB,CAAC;AAED,SAAS,mBAAmB,CAAC,OAAe;IAC1C,OAAO,MAAM,CAAC,QAAQ,CAAC,OAAO,CAAC;WAC1B,OAAO,IAAI,CAAC,yBAAyB;WACrC,OAAO,IAAI,CAAC,GAAG,yBAAyB,CAAC;AAChD,CAAC;AAED,SAAS,eAAe,CAAC,CAAS,EAAE,CAAS,EAAE,gBAAwB;IACrE,MAAM,SAAS,GAAG,gCAAgC,CAChD,oCAAoC,CAAC,gBAAgB,CAAC,CACvD,CAAC;IACF,MAAM,OAAO,GAAG,gCAAgC,CAC9C,SAAS,CAAC,MAAM,CAAC,QAAQ,EACzB,SAAS,CAAC,KAAK,CAAC,QAAQ,EACxB,CAAC,CACF,CAAC;IAEF,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,OAAO,CAAC,EAAE,CAAC;QAC9B,OAAO,MAAM,CAAC,GAAG,CAAC;IACpB,CAAC;IAED,OAAO,SAAS,CAAC,MAAM,CAAC,OAAO,GAAG,CAAC,CAAC,GAAG,OAAO,CAAC,GAAG,SAAS,CAAC,KAAK,CAAC,OAAO,GAAG,OAAO,GAAG,CAAC,CAAC;AAC1F,CAAC;AAED,SAAS,sBAAsB,CAC7B,CAAS,EACT,CAAS,EACT,WAAmB;IAEnB,MAAM,SAAS,GAAG,gCAAgC,CAAC,WAAW,CAAC,CAAC;IAChE,MAAM,mBAAmB,GAAG,gCAAgC,CAC1D,SAAS,CAAC,MAAM,CAAC,QAAQ,EACzB,SAAS,CAAC,KAAK,CAAC,QAAQ,EACxB,CAAC,CACF,CAAC;IACF,MAAM,kBAAkB,GAAG,gCAAgC,CACzD,SAAS,CAAC,MAAM,CAAC,OAAO,EACxB,SAAS,CAAC,KAAK,CAAC,OAAO,EACvB,CAAC,CACF,CAAC;IACF,MAAM,6BAA6B,GAAG,wBAAwB,CAAC,mBAAmB,CAAC,CAAC;IACpF,MAAM,4BAA4B,GAAG,wBAAwB,CAAC,kBAAkB,CAAC,CAAC;IAElF,IAAI,CAAC,mBAAmB,CAAC,6BAA6B,CAAC,IAAI,CAAC,mBAAmB,CAAC,4BAA4B,CAAC,EAAE,CAAC;QAC9G,OAAO,IAAI,CAAC;IACd,CAAC;IAED,IAAI,IAAI,CAAC,GAAG,CAAC,6BAA6B,GAAG,4BAA4B,CAAC,GAAG,yBAAyB,EAAE,CAAC;QACvG,OAAO,IAAI,CAAC;IACd,CAAC;IAED,IAAI,IAAI,CAAC,GAAG,CAAC,eAAe,CAAC,CAAC,EAAE,CAAC,EAAE,WAAW,CAAC,CAAC,GAAG,CAAC,CAAC,6BAA6B,EAAE,CAAC;QACnF,OAAO,IAAI,CAAC;IACd,CAAC;IAED,IAAI,CAAC;QACH,yBAAyB,CAAC,SAAS,CAAC,QAAQ,EAAE,SAAS,CAAC,WAAW,EAAE,WAAW,CAAC,CAAC;IACpF,CAAC;IAAC,MAAM,CAAC;QACP,OAAO,IAAI,CAAC;IACd,CAAC;IAED,OAAO,kBAAkB,CAAC,SAAS,EAAE,6BAA6B,CAAC,CAAC;AACtE,CAAC;AAED,SAAS,qBAAqB,CAC5B,CAAS,EACT,CAAS,EACT,WAAmB,EACnB,KAAsB;IAEtB,IAAI,SAAS,CAAC;IACd,IAAI,CAAC;QACH,SAAS,GAAG,gCAAgC,CAAC,WAAW,CAAC,CAAC;IAC5D,CAAC;IAAC,MAAM,CAAC;QACP,OAAO,IAAI,CAAC;IACd,CAAC;IACD,MAAM,QAAQ,GAAG,SAAS,CAAC,KAAK,CAAC,CAAC;IAElC,IAAI,CAAC,OAAO,CAAC,QAAQ,CAAC,QAAQ,EAAE,CAAC,CAAC,IAAI,CAAC,OAAO,CAAC,QAAQ,CAAC,OAAO,EAAE,CAAC,CAAC,EAAE,CAAC;QACpE,OAAO,IAAI,CAAC;IACd,CAAC;IAED,IAAI,CAAC;QACH,yBAAyB,CAAC,SAAS,CAAC,QAAQ,EAAE,SAAS,CAAC,WAAW,EAAE,WAAW,CAAC,CAAC;IACpF,CAAC;IAAC,MAAM,CAAC;QACP,OAAO,IAAI,CAAC;IACd,CAAC;IAED,OAAO,kBAAkB,CAAC,SAAS,EAAE,KAAK,KAAK,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;AACnE,CAAC;AAED,SAAS,wBAAwB,CAC/B,CAAS,EACT,gBAAwB,EACxB,KAAsB;IAEtB,IAAI,SAAS,CAAC;IACd,IAAI,CAAC;QACH,SAAS,GAAG,gCAAgC,CAC1C,oCAAoC,CAAC,gBAAgB,CAAC,CACvD,CAAC;IACJ,CAAC;IAAC,MAAM,CAAC;QACP,OAAO,MAAM,CAAC,GAAG,CAAC;IACpB,CAAC;IACD,OAAO,SAAS,CAAC,KAAK,CAAC,CAAC,QAAQ,GAAG,CAAC,CAAC;AACvC,CAAC;AAED,SAAS,yBAAyB,CAChC,CAAS,EACT,CAAS;IAET,MAAM,KAAK,GAAG,CAAC,CAAC,EAAE,CAAC;IACnB,MAAM,KAAK,GAAG,oCAAoC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;IAEzD,KAAK,MAAM,KAAK,IAAI,CAAC,QAAQ,EAAE,OAAO,CAAU,EAAE,CAAC;QACjD,MAAM,UAAU,GAAG,qBAAqB,CAAC,CAAC,EAAE,CAAC,EAAE,KAAK,EAAE,KAAK,CAAC,CAAC;QAC7D,IAAI,UAAU,KAAK,IAAI,EAAE,CAAC;YACxB,OAAO,UAAU,CAAC;QACpB,CAAC;QAED,MAAM,UAAU,GAAG,qBAAqB,CAAC,CAAC,EAAE,CAAC,EAAE,KAAK,EAAE,KAAK,CAAC,CAAC;QAC7D,IAAI,UAAU,KAAK,IAAI,EAAE,CAAC;YACxB,OAAO,UAAU,CAAC;QACpB,CAAC;QAED,IAAI,mBAAmB,GAAG,KAAK,CAAC;QAChC,IAAI,gBAAgB,GAAG,wBAAwB,CAAC,CAAC,EAAE,KAAK,EAAE,KAAK,CAAC,CAAC;QAEjE,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,CAAC,2BAA2B,EAAE,CAAC,EAAE,EAAE,CAAC;YACxD,MAAM,WAAW,GAAG,KAAK,GAAG,CAAC,KAAK,GAAG,KAAK,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,2BAA2B,CAAC,CAAC;YAClF,MAAM,QAAQ,GAAG,wBAAwB,CAAC,CAAC,EAAE,WAAW,EAAE,KAAK,CAAC,CAAC;YAEjE,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,gBAAgB,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,QAAQ,CAAC,EAAE,CAAC;gBACrE,mBAAmB,GAAG,WAAW,CAAC;gBAClC,gBAAgB,GAAG,QAAQ,CAAC;gBAC5B,SAAS;YACX,CAAC;YAED,IAAI,OAAO,CAAC,gBAAgB,EAAE,CAAC,CAAC,EAAE,CAAC;gBACjC,OAAO,qBAAqB,CAAC,CAAC,EAAE,CAAC,EAAE,mBAAmB,EAAE,KAAK,CAAC,CAAC;YACjE,CAAC;YAED,IAAI,OAAO,CAAC,QAAQ,EAAE,CAAC,CAAC,EAAE,CAAC;gBACzB,OAAO,qBAAqB,CAAC,CAAC,EAAE,CAAC,EAAE,WAAW,EAAE,KAAK,CAAC,CAAC;YACzD,CAAC;YAED,IAAI,gBAAgB,GAAG,QAAQ,GAAG,CAAC,EAAE,CAAC;gBACpC,MAAM,IAAI,GAAG,eAAe,CAC1B,CAAC,oBAAoB,EAAE,EAAE,CAAC,wBAAwB,CAAC,CAAC,EAAE,oBAAoB,EAAE,KAAK,CAAC,EAClF,mBAAmB,EACnB,WAAW,EACX,CAAC,mBAAmB,GAAG,WAAW,CAAC,GAAG,CAAC,EACvC,EAAE,SAAS,EAAE,CAAC,CAAC,6BAA6B,EAAE,CAC/C,CAAC;gBACF,OAAO,qBAAqB,CAAC,CAAC,EAAE,CAAC,EAAE,IAAI,EAAE,KAAK,CAAC,CAAC;YAClD,CAAC;YAED,mBAAmB,GAAG,WAAW,CAAC;YAClC,gBAAgB,GAAG,QAAQ,CAAC;QAC9B,CAAC;IACH,CAAC;IAED,OAAO,IAAI,CAAC;AACd,CAAC;AAED,MAAM,UAAU,iBAAiB,CAAC,CAAS,EAAE,CAAS;IACpD,MAAM,aAAa,GAAG,yBAAyB,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IACtD,IAAI,aAAa,KAAK,IAAI,EAAE,CAAC;QAC3B,OAAO,aAAa,CAAC;IACvB,CAAC;IAED,MAAM,KAAK,GAAG,CAAC,CAAC,EAAE,CAAC;IACnB,MAAM,KAAK,GAAG,oCAAoC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;IACzD,MAAM,UAAU,GAAG,sBAAsB,CAAC,CAAC,EAAE,CAAC,EAAE,KAAK,CAAC,CAAC;IACvD,IAAI,UAAU,KAAK,IAAI,EAAE,CAAC;QACxB,OAAO,UAAU,CAAC;IACpB,CAAC;IAED,MAAM,UAAU,GAAG,sBAAsB,CAAC,CAAC,EAAE,CAAC,EAAE,KAAK,CAAC,CAAC;IACvD,IAAI,UAAU,KAAK,IAAI,EAAE,CAAC;QACxB,OAAO,UAAU,CAAC;IACpB,CAAC;IAED,IAAI,mBAAmB,GAAG,KAAK,CAAC;IAChC,IAAI,gBAAgB,GAAG,eAAe,CAAC,CAAC,EAAE,CAAC,EAAE,KAAK,CAAC,CAAC;IAEpD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,CAAC,2BAA2B,EAAE,CAAC,EAAE,EAAE,CAAC;QACxD,MAAM,WAAW,GAAG,KAAK,GAAG,CAAC,KAAK,GAAG,KAAK,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,2BAA2B,CAAC,CAAC;QAClF,MAAM,QAAQ,GAAG,eAAe,CAAC,CAAC,EAAE,CAAC,EAAE,WAAW,CAAC,CAAC;QAEpD,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,gBAAgB,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,QAAQ,CAAC,EAAE,CAAC;YACrE,mBAAmB,GAAG,WAAW,CAAC;YAClC,gBAAgB,GAAG,QAAQ,CAAC;YAC5B,SAAS;QACX,CAAC;QAED,IAAI,OAAO,CAAC,gBAAgB,EAAE,CAAC,CAAC,EAAE,CAAC;YACjC,OAAO,sBAAsB,CAAC,CAAC,EAAE,CAAC,EAAE,mBAAmB,CAAC,CAAC;QAC3D,CAAC;QAED,IAAI,OAAO,CAAC,QAAQ,EAAE,CAAC,CAAC,EAAE,CAAC;YACzB,OAAO,sBAAsB,CAAC,CAAC,EAAE,CAAC,EAAE,WAAW,CAAC,CAAC;QACnD,CAAC;QAED,IAAI,gBAAgB,GAAG,QAAQ,GAAG,CAAC,EAAE,CAAC;YACpC,MAAM,IAAI,GAAG,eAAe,CAC1B,CAAC,oBAAoB,EAAE,EAAE,CAAC,eAAe,CAAC,CAAC,EAAE,CAAC,EAAE,oBAAoB,CAAC,EACrE,mBAAmB,EACnB,WAAW,EACX,CAAC,mBAAmB,GAAG,WAAW,CAAC,GAAG,CAAC,EACvC,EAAE,SAAS,EAAE,CAAC,CAAC,6BAA6B,EAAE,CAC/C,CAAC;YACF,OAAO,sBAAsB,CAAC,CAAC,EAAE,CAAC,EAAE,IAAI,CAAC,CAAC;QAC5C,CAAC;QAED,mBAAmB,GAAG,WAAW,CAAC;QAClC,gBAAgB,GAAG,QAAQ,CAAC;IAC9B,CAAC;IAED,OAAO,IAAI,CAAC;AACd,CAAC"}

package/dist/saturation/two-phase.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"two-phase.d.ts","sourceRoot":"","sources":["../../src/saturation/two-phase.ts"],"names":[],"mappings":"AAAA;;;;;;;GAOG;AAEH,OAAO,KAAK,EAAE,eAAe,EAAE,MAAM,aAAa,CAAC;AAQnD;;;;GAIG;AACH,wBAAgB,OAAO,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,eAAe,CAS7D;AAID;;;;GAIG;AACH,wBAAgB,OAAO,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,eAAe,CAU7D"}
+{"version":3,"file":"two-phase.d.ts","sourceRoot":"","sources":["../../src/saturation/two-phase.ts"],"names":[],"mappings":"AAAA;;;;;;;GAOG;AAEH,OAAO,KAAK,EAAE,eAAe,EAAE,MAAM,aAAa,CAAC;AAanD;;;;GAIG;AACH,wBAAgB,OAAO,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,eAAe,CAY7D;AAID;;;;GAIG;AACH,wBAAgB,OAAO,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,eAAe,CAW7D"}

package/dist/saturation/two-phase.js

@@ -7,9 +7,10 @@
  * Reference: IAPWS-IF97, Region 4 (Saturation Line)
  */
 import { OutOfRangeError } from '../types.js';
-import { Pc, Tc, Pt, T_MIN } from '../constants.js';
+import { assertFiniteNumber } from '../core/input-validation.js';
 import { saturationPressure } from '../regions/region4.js';
 import { mixSaturationState, saturationEndpointsAtPressure } from './common.js';
+import { assertRegion4PressureStateInput, assertRegion4TemperatureStateInput, normalizeRegion4Pressure, } from './region4-boundaries.js';
 // ─── solvePx ────────────────────────────────────────────────────────────────
 /**
  * Solve for thermodynamic state given pressure and vapour quality.
@@ -17,13 +18,12 @@ import { mixSaturationState, saturationEndpointsAtPressure } from './common.js';
  * @param x - Vapour quality [0, 1]
  */
 export function solvePx(p, x) {
+    assertFiniteNumber('Pressure', p);
+    assertFiniteNumber('Quality', x);
     if (x < 0 || x > 1) {
         throw new OutOfRangeError('Quality', x, 0, 1);
     }
-    if (p < Pt || p > Pc) {
-        throw new OutOfRangeError('Pressure', p, Pt, Pc);
-    }
-    return mixSaturationState(saturationEndpointsAtPressure(p), x);
+    return mixSaturationState(saturationEndpointsAtPressure(assertRegion4PressureStateInput(p, 'solvePx')), x);
 }
 // ─── solveTx ────────────────────────────────────────────────────────────────
 /**
@@ -32,13 +32,13 @@ export function solvePx(p, x) {
  * @param x - Vapour quality [0, 1]
  */
 export function solveTx(T, x) {
+    assertFiniteNumber('Temperature', T);
+    assertFiniteNumber('Quality', x);
     if (x < 0 || x > 1) {
         throw new OutOfRangeError('Quality', x, 0, 1);
     }
-    if (T < T_MIN || T > Tc) {
-        throw new OutOfRangeError('Temperature', T, T_MIN, Tc);
-    }
-    const p = saturationPressure(T);
+    const temperature = assertRegion4TemperatureStateInput(T, 'solveTx');
+    const p = normalizeRegion4Pressure(saturationPressure(temperature));
     return solvePx(p, x);
 }
 //# sourceMappingURL=two-phase.js.map

package/dist/saturation/two-phase.js.map

@@ -1 +1 @@
-{"version":3,"file":"two-phase.js","sourceRoot":"","sources":["../../src/saturation/two-phase.ts"],"names":[],"mappings":"AAAA;;;;;;;GAOG;AAGH,OAAO,EAAE,eAAe,EAAE,MAAM,aAAa,CAAC;AAC9C,OAAO,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,KAAK,EAAE,MAAM,iBAAiB,CAAC;AACpD,OAAO,EAAE,kBAAkB,EAAE,MAAM,uBAAuB,CAAC;AAC3D,OAAO,EAAE,kBAAkB,EAAE,6BAA6B,EAAE,MAAM,aAAa,CAAC;AAEhF,+EAA+E;AAE/E;;;;GAIG;AACH,MAAM,UAAU,OAAO,CAAC,CAAS,EAAE,CAAS;IAC1C,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC;QACnB,MAAM,IAAI,eAAe,CAAC,SAAS,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IAChD,CAAC;IACD,IAAI,CAAC,GAAG,EAAE,IAAI,CAAC,GAAG,EAAE,EAAE,CAAC;QACrB,MAAM,IAAI,eAAe,CAAC,UAAU,EAAE,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;IACnD,CAAC;IAED,OAAO,kBAAkB,CAAC,6BAA6B,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;AACjE,CAAC;AAED,+EAA+E;AAE/E;;;;GAIG;AACH,MAAM,UAAU,OAAO,CAAC,CAAS,EAAE,CAAS;IAC1C,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC;QACnB,MAAM,IAAI,eAAe,CAAC,SAAS,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IAChD,CAAC;IACD,IAAI,CAAC,GAAG,KAAK,IAAI,CAAC,GAAG,EAAE,EAAE,CAAC;QACxB,MAAM,IAAI,eAAe,CAAC,aAAa,EAAE,CAAC,EAAE,KAAK,EAAE,EAAE,CAAC,CAAC;IACzD,CAAC;IAED,MAAM,CAAC,GAAG,kBAAkB,CAAC,CAAC,CAAC,CAAC;IAChC,OAAO,OAAO,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;AACvB,CAAC"}
+{"version":3,"file":"two-phase.js","sourceRoot":"","sources":["../../src/saturation/two-phase.ts"],"names":[],"mappings":"AAAA;;;;;;;GAOG;AAGH,OAAO,EAAE,eAAe,EAAE,MAAM,aAAa,CAAC;AAC9C,OAAO,EAAE,kBAAkB,EAAE,MAAM,6BAA6B,CAAC;AACjE,OAAO,EAAE,kBAAkB,EAAE,MAAM,uBAAuB,CAAC;AAC3D,OAAO,EAAE,kBAAkB,EAAE,6BAA6B,EAAE,MAAM,aAAa,CAAC;AAChF,OAAO,EACL,+BAA+B,EAC/B,kCAAkC,EAClC,wBAAwB,GACzB,MAAM,yBAAyB,CAAC;AAEjC,+EAA+E;AAE/E;;;;GAIG;AACH,MAAM,UAAU,OAAO,CAAC,CAAS,EAAE,CAAS;IAC1C,kBAAkB,CAAC,UAAU,EAAE,CAAC,CAAC,CAAC;IAClC,kBAAkB,CAAC,SAAS,EAAE,CAAC,CAAC,CAAC;IAEjC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC;QACnB,MAAM,IAAI,eAAe,CAAC,SAAS,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IAChD,CAAC;IAED,OAAO,kBAAkB,CACvB,6BAA6B,CAAC,+BAA+B,CAAC,CAAC,EAAE,SAAS,CAAC,CAAC,EAC5E,CAAC,CACF,CAAC;AACJ,CAAC;AAED,+EAA+E;AAE/E;;;;GAIG;AACH,MAAM,UAAU,OAAO,CAAC,CAAS,EAAE,CAAS;IAC1C,kBAAkB,CAAC,aAAa,EAAE,CAAC,CAAC,CAAC;IACrC,kBAAkB,CAAC,SAAS,EAAE,CAAC,CAAC,CAAC;IAEjC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC;QACnB,MAAM,IAAI,eAAe,CAAC,SAAS,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IAChD,CAAC;IAED,MAAM,WAAW,GAAG,kCAAkC,CAAC,CAAC,EAAE,SAAS,CAAC,CAAC;IACrE,MAAM,CAAC,GAAG,wBAAwB,CAAC,kBAAkB,CAAC,WAAW,CAAC,CAAC,CAAC;IACpE,OAAO,OAAO,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;AACvB,CAAC"}

package/dist/saturation.d.ts

@@ -1,5 +1,10 @@
 /**
- * Public saturation helper exports.
+ * Public low-level saturation helper exports.
+ *
+ * These functions expose the mathematical Region 4 boundary relations.
+ * They normalize values near Pt/Pc/Tt/Tc, but they are not state constructors:
+ * high-level solvers such as solvePx/solveTx apply stricter triple-point and
+ * critical-point rules.
  */
 export { saturationPressure, saturationTemperature } from './regions/region4.js';
 //# sourceMappingURL=saturation.d.ts.map

package/dist/saturation.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"saturation.d.ts","sourceRoot":"","sources":["../src/saturation.ts"],"names":[],"mappings":"AAAA;;GAEG;AAEH,OAAO,EAAE,kBAAkB,EAAE,qBAAqB,EAAE,MAAM,sBAAsB,CAAC"}
+{"version":3,"file":"saturation.d.ts","sourceRoot":"","sources":["../src/saturation.ts"],"names":[],"mappings":"AAAA;;;;;;;GAOG;AAEH,OAAO,EAAE,kBAAkB,EAAE,qBAAqB,EAAE,MAAM,sBAAsB,CAAC"}

package/dist/saturation.js

@@ -1,5 +1,10 @@
 /**
- * Public saturation helper exports.
+ * Public low-level saturation helper exports.
+ *
+ * These functions expose the mathematical Region 4 boundary relations.
+ * They normalize values near Pt/Pc/Tt/Tc, but they are not state constructors:
+ * high-level solvers such as solvePx/solveTx apply stricter triple-point and
+ * critical-point rules.
  */
 export { saturationPressure, saturationTemperature } from './regions/region4.js';
 //# sourceMappingURL=saturation.js.map

package/dist/saturation.js.map

@@ -1 +1 @@
-{"version":3,"file":"saturation.js","sourceRoot":"","sources":["../src/saturation.ts"],"names":[],"mappings":"AAAA;;GAEG;AAEH,OAAO,EAAE,kBAAkB,EAAE,qBAAqB,EAAE,MAAM,sBAAsB,CAAC"}
+{"version":3,"file":"saturation.js","sourceRoot":"","sources":["../src/saturation.ts"],"names":[],"mappings":"AAAA;;;;;;;GAOG;AAEH,OAAO,EAAE,kBAAkB,EAAE,qBAAqB,EAAE,MAAM,sBAAsB,CAAC"}

package/dist/solvers/nelder-mead.d.ts

@@ -1,7 +1,3 @@
-/**
- * Nelder-Mead (downhill simplex) minimizer.
- * Minimizes f(x) where x is an array of numbers.
- */
 export interface NelderMeadOptions {
     maxIterations?: number;
     nonZeroDelta?: number;

package/dist/solvers/nelder-mead.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"nelder-mead.d.ts","sourceRoot":"","sources":["../../src/solvers/nelder-mead.ts"],"names":[],"mappings":"AAAA;;;GAGG;AAEH,MAAM,WAAW,iBAAiB;IAChC,aAAa,CAAC,EAAE,MAAM,CAAC;IACvB,YAAY,CAAC,EAAE,MAAM,CAAC;IACtB,SAAS,CAAC,EAAE,MAAM,CAAC;IACnB,aAAa,CAAC,EAAE,MAAM,CAAC;IACvB,YAAY,CAAC,EAAE,MAAM,CAAC;IACtB,GAAG,CAAC,EAAE,MAAM,CAAC;IACb,GAAG,CAAC,EAAE,MAAM,CAAC;IACb,GAAG,CAAC,EAAE,MAAM,CAAC;IACb,KAAK,CAAC,EAAE,MAAM,CAAC;CAChB;AAED,MAAM,WAAW,gBAAgB;IAC/B,CAAC,EAAE,MAAM,EAAE,CAAC;IACZ,EAAE,EAAE,MAAM,CAAC;CACZ;AAED,wBAAgB,UAAU,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,MAAM,EAAE,KAAK,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,IAAI,CAAC,EAAE,iBAAiB,GAAG,gBAAgB,CAiG/G"}
+{"version":3,"file":"nelder-mead.d.ts","sourceRoot":"","sources":["../../src/solvers/nelder-mead.ts"],"names":[],"mappings":"AAMA,MAAM,WAAW,iBAAiB;IAChC,aAAa,CAAC,EAAE,MAAM,CAAC;IACvB,YAAY,CAAC,EAAE,MAAM,CAAC;IACtB,SAAS,CAAC,EAAE,MAAM,CAAC;IACnB,aAAa,CAAC,EAAE,MAAM,CAAC;IACvB,YAAY,CAAC,EAAE,MAAM,CAAC;IACtB,GAAG,CAAC,EAAE,MAAM,CAAC;IACb,GAAG,CAAC,EAAE,MAAM,CAAC;IACb,GAAG,CAAC,EAAE,MAAM,CAAC;IACb,KAAK,CAAC,EAAE,MAAM,CAAC;CAChB;AAED,MAAM,WAAW,gBAAgB;IAC/B,CAAC,EAAE,MAAM,EAAE,CAAC;IACZ,EAAE,EAAE,MAAM,CAAC;CACZ;AAED,wBAAgB,UAAU,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,MAAM,EAAE,KAAK,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,IAAI,CAAC,EAAE,iBAAiB,GAAG,gBAAgB,CAqH/G"}