driftlet 0.5.0 → 0.6.0

package/types/solver.d.ts

@@ -44,6 +44,7 @@ export class Solver {
     jR: Float64Array<any>;
     time: number;
     lastDt: number;
+    atSteady: boolean;
     contactFlux: {
         left: Float64Array<any>;
         right: Float64Array<any>;
@@ -155,8 +156,11 @@ export class Solver {
     _segmentMixing(s: any, bL: any, bR: any, Dm: any, h: any, mat: any): void;
     mixWork: Float64Array<any>[];
     _dc(g: any, k: any, d: any): void;
+    bulkRate(rx: any, g: any): number;
+    bulkOneWay(rx: any, g: any): number;
     _srhRate(law: any, cn: any, cp: any, a: any): {
         rate: number;
+        oneWay: number;
         Cn: number;
         Cp: number;
         Ca: number;
@@ -243,12 +247,15 @@ export class Solver {
         push: (values: any) => any;
     };
     _errorNorm(ref: any): number;
+    _warmFrom(src: any): void;
+    solvedV: any;
     _snapshot(): {
         u: Float64Array<ArrayBuffer>;
         uLo: Float64Array<ArrayBuffer>;
         cOld: Float64Array<ArrayBuffer>;
         time: number;
         lastDt: number;
+        atSteady: boolean;
         contactDStart: {
             left: number;
             right: number;
@@ -309,7 +316,7 @@ export class Solver {
      * The clock is not advanced, and open-system conservation bookkeeping restarts here.
      */
     solveSteady(opts?: {}): any;
-    solvedV: any;
+    _steadyFromHere(opts: any): any;
     _generationContinuation(opts: any, r: any): {
         converged: boolean;
         steps: any;

package/types/stack.d.ts

@@ -19,6 +19,7 @@ export function combine(...parts: {
  * @param {object} [more]
  */
 export function layer(material: string, length: number, more?: object): any;
+export function drive(d: any, path: any): any;
 /**
  * An ohmic contact: each carrier in equilibrium with a metal at the terminal voltage, the
  * first being the terminal species, and every other sharing its voltage (offset 0, as holes and

package/types/types.d.ts

@@ -510,12 +510,23 @@ export type Solution = {
         current: number;
     }[];
     interfaces: {
+        left: string;
+        right: string;
         dipole: number;
         sheetCharge: number;
         D: number;
         N: Record<string, number>;
         rates: number[];
     }[];
+    /**
+     * each bulk reaction's
+     * forward rate per node (mol/(m³·s), NaN where it doesn't run), and per region and in all (mol/(m²·s))
+     */
+    bulkReactions: {
+        rate: Float64Array;
+        regions: number[];
+        total: number;
+    }[];
     /**
      * total charge in the device, C/m²
      */
@@ -538,6 +549,10 @@ export type Solution = {
      * s
      */
     time: number;
+    /**
+     * whether this is a converged steady state (solve()), not a transient's
+     */
+    steady: boolean;
     /**
      * advance(): whether tEnd was reached (false if the frame budget ran out)
      */