@johnhenry/math 0.0.2 → 0.1.0

package/dist/index.d.ts

@@ -42,7 +42,7 @@ export { SpecialOperator } from "./SpecialOperator.ts";
 export * as Statistics from "./Statistics.ts";
 export { StringEvaluator } from "./StringEvaluator.ts";
 export { Structure, type StructureOptions } from "./Structure.ts";
-export { type Assumption, type BinaryFuncName, type CmpOp, DegenerateOdeError, type DifferentiationStep, type Expr, FUNCTION_NAMES, type FuncName, InfiniteSolutionsError, IntegrationSingularityError, NoClosedFormError, NonLinearSystemError, NotIntegrableError, NotSeparableError, type OdeClosedFormResult, ProductNotDifferentiableError, SeriesDivergesError, SingularSystemError, Symbolic, SystemDidNotConvergeError, UndeclaredVariableError, } from "./Symbolic.ts";
+export { type Assumption, type BinaryFuncName, type CmpOp, DegenerateOdeError, type DifferentiationStep, type Expr, FUNCTION_NAMES, type FuncName, InfiniteSolutionsError, IntegrationSingularityError, NoClosedFormError, NonLinearSystemError, NotIntegrableError, NotSeparableError, type OdeClosedFormResult, type ParseOptions, ProductNotDifferentiableError, SeriesDivergesError, SingularSystemError, Symbolic, SystemDidNotConvergeError, UndeclaredVariableError, } from "./Symbolic.ts";
 export { Type, TypeTag } from "./Type.ts";
 export { Utilities } from "./Utilities.ts";
 export { Vec4 } from "./Vec4.ts";

package/dist/index.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"index.d.ts","sourceRoot":"","sources":["../src/index.ts"],"names":[],"mappings":"AAAA;;;;GAIG;AAIH,OAAO,EAAE,SAAS,EAAE,MAAM,gBAAgB,CAAC;AAI3C,OAAO,EAAE,SAAS,EAAE,KAAK,QAAQ,EAAE,uBAAuB,EAAE,eAAe,EAAE,MAAM,gBAAgB,CAAC;AAEpG,OAAO,EAAE,aAAa,EAAE,MAAM,oBAAoB,CAAC;AACnD,OAAO,EACL,YAAY,EACZ,2BAA2B,EAC3B,YAAY,EACZ,OAAO,GACR,MAAM,8BAA8B,CAAC;AACtC,OAAO,KAAK,WAAW,MAAM,kBAAkB,CAAC;AAChD,OAAO,EAAE,aAAa,EAAE,MAAM,oBAAoB,CAAC;AACnD,OAAO,EAAE,KAAK,EAAE,MAAM,YAAY,CAAC;AACnC,OAAO,EAAE,OAAO,EAAE,MAAM,cAAc,CAAC;AACvC,OAAO,EACL,KAAK,sBAAsB,EAC3B,KAAK,oBAAoB,EACzB,aAAa,EACb,eAAe,EACf,KAAK,UAAU,GAChB,MAAM,oBAAoB,CAAC;AAC5B,OAAO,EAAE,UAAU,EAAE,MAAM,iBAAiB,CAAC;AAC7C,OAAO,EAAE,qBAAqB,EAAE,KAAK,kBAAkB,EAAE,MAAM,4BAA4B,CAAC;AAE5F,OAAO,EAAE,WAAW,EAAE,MAAM,kBAAkB,CAAC;AAC/C,OAAO,EAAE,UAAU,EAAE,MAAM,iBAAiB,CAAC;AAC7C,OAAO,EAAE,GAAG,EAAE,MAAM,UAAU,CAAC;AAC/B,OAAO,EACL,KAAK,wBAAwB,EAC7B,iBAAiB,EACjB,UAAU,IAAI,WAAW,EACzB,IAAI,IAAI,KAAK,EACb,yBAAyB,EACzB,uBAAuB,EACvB,eAAe,EACf,KAAK,sBAAsB,EAC3B,YAAY,IAAI,aAAa,GAC9B,MAAM,wBAAwB,CAAC;AAChC,OAAO,EAAE,QAAQ,EAAE,KAAK,KAAK,EAAE,WAAW,EAAE,MAAM,eAAe,CAAC;AAClE,OAAO,EAAE,KAAK,IAAI,EAAE,KAAK,EAAE,MAAM,YAAY,CAAC;AAC9C,OAAO,EACL,KAAK,IAAI,EACT,YAAY,EACZ,KAAK,QAAQ,EACb,KAAK,IAAI,EACT,KAAK,WAAW,EAChB,KAAK,IAAI,GACV,MAAM,mBAAmB,CAAC;AAE3B,OAAO,EACL,KAAK,cAAc,EACnB,KAAK,SAAS,EACd,UAAU,EACV,KAAK,MAAM,EACX,KAAK,WAAW,EAChB,KAAK,WAAW,EAChB,KAAK,WAAW,GACjB,MAAM,iBAAiB,CAAC;AACzB,OAAO,EAAE,WAAW,EAAE,MAAM,kBAAkB,CAAC;AAE/C,OAAO,EAAE,QAAQ,EAAE,MAAM,eAAe,CAAC;AACzC,OAAO,EAAE,QAAQ,EAAE,MAAM,eAAe,CAAC;AACzC,OAAO,EACL,oBAAoB,EACpB,KAAK,YAAY,EACjB,WAAW,EACX,YAAY,EACZ,KAAK,YAAY,EACjB,KAAK,cAAc,EACnB,WAAW,EACX,YAAY,EACZ,cAAc,GACf,MAAM,cAAc,CAAC;AACtB,OAAO,EAAE,KAAK,EAAE,MAAM,YAAY,CAAC;AACnC,OAAO,EACL,KAAK,WAAW,EAChB,KAAK,QAAQ,EACb,UAAU,EACV,KAAK,QAAQ,EACb,KAAK,SAAS,GACf,MAAM,iBAAiB,CAAC;AACzB,OAAO,EAAE,YAAY,EAAE,MAAM,mBAAmB,CAAC;AACjD,OAAO,EAAE,KAAK,wBAAwB,EAAE,SAAS,EAAE,KAAK,OAAO,EAAE,MAAM,gBAAgB,CAAC;AACxF,OAAO,EAAE,WAAW,EAAE,MAAM,kBAAkB,CAAC;AAC/C,OAAO,EAAE,OAAO,EAAE,MAAM,cAAc,CAAC;AACvC,OAAO,EAAE,cAAc,EAAE,eAAe,EAAE,kBAAkB,EAAE,MAAM,qBAAqB,CAAC;AAC1F,OAAO,EAAE,UAAU,EAAE,MAAM,iBAAiB,CAAC;AAC7C,OAAO,EAAE,QAAQ,EAAE,MAAM,eAAe,CAAC;AAEzC,OAAO,KAAK,QAAQ,MAAM,eAAe,CAAC;AAC1C,OAAO,EAAE,MAAM,EAAE,MAAM,aAAa,CAAC;AACrC,OAAO,EAAE,gBAAgB,EAAE,MAAM,uBAAuB,CAAC;AACzD,OAAO,EAAE,eAAe,EAAE,MAAM,sBAAsB,CAAC;AACvD,OAAO,KAAK,UAAU,MAAM,iBAAiB,CAAC;AAC9C,OAAO,EAAE,eAAe,EAAE,MAAM,sBAAsB,CAAC;AAEvD,OAAO,EAAE,SAAS,EAAE,KAAK,gBAAgB,EAAE,MAAM,gBAAgB,CAAC;AAClE,OAAO,EACL,KAAK,UAAU,EACf,KAAK,cAAc,EACnB,KAAK,KAAK,EACV,kBAAkB,EAClB,KAAK,mBAAmB,EACxB,KAAK,IAAI,EACT,cAAc,EACd,KAAK,QAAQ,EACb,sBAAsB,EACtB,2BAA2B,EAC3B,iBAAiB,EACjB,oBAAoB,EACpB,kBAAkB,EAClB,iBAAiB,EACjB,KAAK,mBAAmB,EACxB,6BAA6B,EAC7B,mBAAmB,EACnB,mBAAmB,EACnB,QAAQ,EACR,yBAAyB,EACzB,uBAAuB,GACxB,MAAM,eAAe,CAAC;AACvB,OAAO,EAAE,IAAI,EAAE,OAAO,EAAE,MAAM,WAAW,CAAC;AAE1C,OAAO,EAAE,SAAS,EAAE,MAAM,gBAAgB,CAAC;AAC3C,OAAO,EAAE,IAAI,EAAE,MAAM,WAAW,CAAC;AAEjC,OAAO,EAAE,MAAM,EAAE,MAAM,aAAa,CAAC;AACrC,OAAO,EAAE,cAAc,EAAE,MAAM,qBAAqB,CAAC;AAErD,OAAO,EAAE,KAAK,MAAM,EAAE,WAAW,EAAE,MAAM,kBAAkB,CAAC"}
+{"version":3,"file":"index.d.ts","sourceRoot":"","sources":["../src/index.ts"],"names":[],"mappings":"AAAA;;;;GAIG;AAIH,OAAO,EAAE,SAAS,EAAE,MAAM,gBAAgB,CAAC;AAI3C,OAAO,EAAE,SAAS,EAAE,KAAK,QAAQ,EAAE,uBAAuB,EAAE,eAAe,EAAE,MAAM,gBAAgB,CAAC;AAEpG,OAAO,EAAE,aAAa,EAAE,MAAM,oBAAoB,CAAC;AACnD,OAAO,EACL,YAAY,EACZ,2BAA2B,EAC3B,YAAY,EACZ,OAAO,GACR,MAAM,8BAA8B,CAAC;AACtC,OAAO,KAAK,WAAW,MAAM,kBAAkB,CAAC;AAChD,OAAO,EAAE,aAAa,EAAE,MAAM,oBAAoB,CAAC;AACnD,OAAO,EAAE,KAAK,EAAE,MAAM,YAAY,CAAC;AACnC,OAAO,EAAE,OAAO,EAAE,MAAM,cAAc,CAAC;AACvC,OAAO,EACL,KAAK,sBAAsB,EAC3B,KAAK,oBAAoB,EACzB,aAAa,EACb,eAAe,EACf,KAAK,UAAU,GAChB,MAAM,oBAAoB,CAAC;AAC5B,OAAO,EAAE,UAAU,EAAE,MAAM,iBAAiB,CAAC;AAC7C,OAAO,EAAE,qBAAqB,EAAE,KAAK,kBAAkB,EAAE,MAAM,4BAA4B,CAAC;AAE5F,OAAO,EAAE,WAAW,EAAE,MAAM,kBAAkB,CAAC;AAC/C,OAAO,EAAE,UAAU,EAAE,MAAM,iBAAiB,CAAC;AAC7C,OAAO,EAAE,GAAG,EAAE,MAAM,UAAU,CAAC;AAC/B,OAAO,EACL,KAAK,wBAAwB,EAC7B,iBAAiB,EACjB,UAAU,IAAI,WAAW,EACzB,IAAI,IAAI,KAAK,EACb,yBAAyB,EACzB,uBAAuB,EACvB,eAAe,EACf,KAAK,sBAAsB,EAC3B,YAAY,IAAI,aAAa,GAC9B,MAAM,wBAAwB,CAAC;AAChC,OAAO,EAAE,QAAQ,EAAE,KAAK,KAAK,EAAE,WAAW,EAAE,MAAM,eAAe,CAAC;AAClE,OAAO,EAAE,KAAK,IAAI,EAAE,KAAK,EAAE,MAAM,YAAY,CAAC;AAC9C,OAAO,EACL,KAAK,IAAI,EACT,YAAY,EACZ,KAAK,QAAQ,EACb,KAAK,IAAI,EACT,KAAK,WAAW,EAChB,KAAK,IAAI,GACV,MAAM,mBAAmB,CAAC;AAE3B,OAAO,EACL,KAAK,cAAc,EACnB,KAAK,SAAS,EACd,UAAU,EACV,KAAK,MAAM,EACX,KAAK,WAAW,EAChB,KAAK,WAAW,EAChB,KAAK,WAAW,GACjB,MAAM,iBAAiB,CAAC;AACzB,OAAO,EAAE,WAAW,EAAE,MAAM,kBAAkB,CAAC;AAE/C,OAAO,EAAE,QAAQ,EAAE,MAAM,eAAe,CAAC;AACzC,OAAO,EAAE,QAAQ,EAAE,MAAM,eAAe,CAAC;AACzC,OAAO,EACL,oBAAoB,EACpB,KAAK,YAAY,EACjB,WAAW,EACX,YAAY,EACZ,KAAK,YAAY,EACjB,KAAK,cAAc,EACnB,WAAW,EACX,YAAY,EACZ,cAAc,GACf,MAAM,cAAc,CAAC;AACtB,OAAO,EAAE,KAAK,EAAE,MAAM,YAAY,CAAC;AACnC,OAAO,EACL,KAAK,WAAW,EAChB,KAAK,QAAQ,EACb,UAAU,EACV,KAAK,QAAQ,EACb,KAAK,SAAS,GACf,MAAM,iBAAiB,CAAC;AACzB,OAAO,EAAE,YAAY,EAAE,MAAM,mBAAmB,CAAC;AACjD,OAAO,EAAE,KAAK,wBAAwB,EAAE,SAAS,EAAE,KAAK,OAAO,EAAE,MAAM,gBAAgB,CAAC;AACxF,OAAO,EAAE,WAAW,EAAE,MAAM,kBAAkB,CAAC;AAC/C,OAAO,EAAE,OAAO,EAAE,MAAM,cAAc,CAAC;AACvC,OAAO,EAAE,cAAc,EAAE,eAAe,EAAE,kBAAkB,EAAE,MAAM,qBAAqB,CAAC;AAC1F,OAAO,EAAE,UAAU,EAAE,MAAM,iBAAiB,CAAC;AAC7C,OAAO,EAAE,QAAQ,EAAE,MAAM,eAAe,CAAC;AAEzC,OAAO,KAAK,QAAQ,MAAM,eAAe,CAAC;AAC1C,OAAO,EAAE,MAAM,EAAE,MAAM,aAAa,CAAC;AACrC,OAAO,EAAE,gBAAgB,EAAE,MAAM,uBAAuB,CAAC;AACzD,OAAO,EAAE,eAAe,EAAE,MAAM,sBAAsB,CAAC;AACvD,OAAO,KAAK,UAAU,MAAM,iBAAiB,CAAC;AAC9C,OAAO,EAAE,eAAe,EAAE,MAAM,sBAAsB,CAAC;AAEvD,OAAO,EAAE,SAAS,EAAE,KAAK,gBAAgB,EAAE,MAAM,gBAAgB,CAAC;AAClE,OAAO,EACL,KAAK,UAAU,EACf,KAAK,cAAc,EACnB,KAAK,KAAK,EACV,kBAAkB,EAClB,KAAK,mBAAmB,EACxB,KAAK,IAAI,EACT,cAAc,EACd,KAAK,QAAQ,EACb,sBAAsB,EACtB,2BAA2B,EAC3B,iBAAiB,EACjB,oBAAoB,EACpB,kBAAkB,EAClB,iBAAiB,EACjB,KAAK,mBAAmB,EACxB,KAAK,YAAY,EACjB,6BAA6B,EAC7B,mBAAmB,EACnB,mBAAmB,EACnB,QAAQ,EACR,yBAAyB,EACzB,uBAAuB,GACxB,MAAM,eAAe,CAAC;AACvB,OAAO,EAAE,IAAI,EAAE,OAAO,EAAE,MAAM,WAAW,CAAC;AAE1C,OAAO,EAAE,SAAS,EAAE,MAAM,gBAAgB,CAAC;AAC3C,OAAO,EAAE,IAAI,EAAE,MAAM,WAAW,CAAC;AAEjC,OAAO,EAAE,MAAM,EAAE,MAAM,aAAa,CAAC;AACrC,OAAO,EAAE,cAAc,EAAE,MAAM,qBAAqB,CAAC;AAErD,OAAO,EAAE,KAAK,MAAM,EAAE,WAAW,EAAE,MAAM,kBAAkB,CAAC"}

package/dist/index.js.map

@@ -1 +1 @@
-{"version":3,"file":"index.js","sourceRoot":"","sources":["../src/index.ts"],"names":[],"mappings":"AAAA;;;;GAIG;AAEH,sFAAsF;AACtF,4EAA4E;AAC5E,OAAO,EAAE,SAAS,EAAE,MAAM,gBAAgB,CAAC;AAC3C,2EAA2E;AAC3E,yEAAyE;AACzE,wEAAwE;AACxE,OAAO,EAAE,SAAS,EAAiB,uBAAuB,EAAE,eAAe,EAAE,MAAM,gBAAgB,CAAC;AACpG,uBAAuB;AACvB,OAAO,EAAE,aAAa,EAAE,MAAM,oBAAoB,CAAC;AACnD,OAAO,EACL,YAAY,EACZ,2BAA2B,EAC3B,YAAY,EACZ,OAAO,GACR,MAAM,8BAA8B,CAAC;AACtC,OAAO,KAAK,WAAW,MAAM,kBAAkB,CAAC;AAChD,OAAO,EAAE,aAAa,EAAE,MAAM,oBAAoB,CAAC;AACnD,OAAO,EAAE,KAAK,EAAE,MAAM,YAAY,CAAC;AACnC,OAAO,EAAE,OAAO,EAAE,MAAM,cAAc,CAAC;AACvC,OAAO,EAGL,aAAa,EACb,eAAe,GAEhB,MAAM,oBAAoB,CAAC;AAC5B,OAAO,EAAE,UAAU,EAAE,MAAM,iBAAiB,CAAC;AAC7C,OAAO,EAAE,qBAAqB,EAA2B,MAAM,4BAA4B,CAAC;AAC5F,wBAAwB;AACxB,OAAO,EAAE,WAAW,EAAE,MAAM,kBAAkB,CAAC;AAC/C,OAAO,EAAE,UAAU,EAAE,MAAM,iBAAiB,CAAC;AAC7C,OAAO,EAAE,GAAG,EAAE,MAAM,UAAU,CAAC;AAC/B,OAAO,EAEL,iBAAiB,EACjB,UAAU,IAAI,WAAW,EACzB,IAAI,IAAI,KAAK,EACb,yBAAyB,EACzB,uBAAuB,EACvB,eAAe,EAEf,YAAY,IAAI,aAAa,GAC9B,MAAM,wBAAwB,CAAC;AAChC,OAAO,EAAE,QAAQ,EAAc,WAAW,EAAE,MAAM,eAAe,CAAC;AAClE,OAAO,EAAa,KAAK,EAAE,MAAM,YAAY,CAAC;AAC9C,OAAO,EAEL,YAAY,GAKb,MAAM,mBAAmB,CAAC;AAC3B,wCAAwC;AACxC,OAAO,EAGL,UAAU,GAKX,MAAM,iBAAiB,CAAC;AACzB,OAAO,EAAE,WAAW,EAAE,MAAM,kBAAkB,CAAC;AAC/C,gCAAgC;AAChC,OAAO,EAAE,QAAQ,EAAE,MAAM,eAAe,CAAC;AACzC,OAAO,EAAE,QAAQ,EAAE,MAAM,eAAe,CAAC;AACzC,OAAO,EACL,oBAAoB,EAEpB,WAAW,EACX,YAAY,EAGZ,WAAW,EACX,YAAY,EACZ,cAAc,GACf,MAAM,cAAc,CAAC;AACtB,OAAO,EAAE,KAAK,EAAE,MAAM,YAAY,CAAC;AACnC,OAAO,EAGL,UAAU,GAGX,MAAM,iBAAiB,CAAC;AACzB,OAAO,EAAE,YAAY,EAAE,MAAM,mBAAmB,CAAC;AACjD,OAAO,EAAiC,SAAS,EAAgB,MAAM,gBAAgB,CAAC;AACxF,OAAO,EAAE,WAAW,EAAE,MAAM,kBAAkB,CAAC;AAC/C,OAAO,EAAE,OAAO,EAAE,MAAM,cAAc,CAAC;AACvC,OAAO,EAAE,cAAc,EAAE,eAAe,EAAE,kBAAkB,EAAE,MAAM,qBAAqB,CAAC;AAC1F,OAAO,EAAE,UAAU,EAAE,MAAM,iBAAiB,CAAC;AAC7C,OAAO,EAAE,QAAQ,EAAE,MAAM,eAAe,CAAC;AACzC,gBAAgB;AAChB,OAAO,KAAK,QAAQ,MAAM,eAAe,CAAC;AAC1C,OAAO,EAAE,MAAM,EAAE,MAAM,aAAa,CAAC;AACrC,OAAO,EAAE,gBAAgB,EAAE,MAAM,uBAAuB,CAAC;AACzD,OAAO,EAAE,eAAe,EAAE,MAAM,sBAAsB,CAAC;AACvD,OAAO,KAAK,UAAU,MAAM,iBAAiB,CAAC;AAC9C,OAAO,EAAE,eAAe,EAAE,MAAM,sBAAsB,CAAC;AACvD,kCAAkC;AAClC,OAAO,EAAE,SAAS,EAAyB,MAAM,gBAAgB,CAAC;AAClE,OAAO,EAIL,kBAAkB,EAGlB,cAAc,EAEd,sBAAsB,EACtB,2BAA2B,EAC3B,iBAAiB,EACjB,oBAAoB,EACpB,kBAAkB,EAClB,iBAAiB,EAEjB,6BAA6B,EAC7B,mBAAmB,EACnB,mBAAmB,EACnB,QAAQ,EACR,yBAAyB,EACzB,uBAAuB,GACxB,MAAM,eAAe,CAAC;AACvB,OAAO,EAAE,IAAI,EAAE,OAAO,EAAE,MAAM,WAAW,CAAC;AAC1C,qBAAqB;AACrB,OAAO,EAAE,SAAS,EAAE,MAAM,gBAAgB,CAAC;AAC3C,OAAO,EAAE,IAAI,EAAE,MAAM,WAAW,CAAC;AACjC,cAAc;AACd,OAAO,EAAE,MAAM,EAAE,MAAM,aAAa,CAAC;AACrC,OAAO,EAAE,cAAc,EAAE,MAAM,qBAAqB,CAAC;AACrD,iBAAiB;AACjB,OAAO,EAAe,WAAW,EAAE,MAAM,kBAAkB,CAAC"}
+{"version":3,"file":"index.js","sourceRoot":"","sources":["../src/index.ts"],"names":[],"mappings":"AAAA;;;;GAIG;AAEH,sFAAsF;AACtF,4EAA4E;AAC5E,OAAO,EAAE,SAAS,EAAE,MAAM,gBAAgB,CAAC;AAC3C,2EAA2E;AAC3E,yEAAyE;AACzE,wEAAwE;AACxE,OAAO,EAAE,SAAS,EAAiB,uBAAuB,EAAE,eAAe,EAAE,MAAM,gBAAgB,CAAC;AACpG,uBAAuB;AACvB,OAAO,EAAE,aAAa,EAAE,MAAM,oBAAoB,CAAC;AACnD,OAAO,EACL,YAAY,EACZ,2BAA2B,EAC3B,YAAY,EACZ,OAAO,GACR,MAAM,8BAA8B,CAAC;AACtC,OAAO,KAAK,WAAW,MAAM,kBAAkB,CAAC;AAChD,OAAO,EAAE,aAAa,EAAE,MAAM,oBAAoB,CAAC;AACnD,OAAO,EAAE,KAAK,EAAE,MAAM,YAAY,CAAC;AACnC,OAAO,EAAE,OAAO,EAAE,MAAM,cAAc,CAAC;AACvC,OAAO,EAGL,aAAa,EACb,eAAe,GAEhB,MAAM,oBAAoB,CAAC;AAC5B,OAAO,EAAE,UAAU,EAAE,MAAM,iBAAiB,CAAC;AAC7C,OAAO,EAAE,qBAAqB,EAA2B,MAAM,4BAA4B,CAAC;AAC5F,wBAAwB;AACxB,OAAO,EAAE,WAAW,EAAE,MAAM,kBAAkB,CAAC;AAC/C,OAAO,EAAE,UAAU,EAAE,MAAM,iBAAiB,CAAC;AAC7C,OAAO,EAAE,GAAG,EAAE,MAAM,UAAU,CAAC;AAC/B,OAAO,EAEL,iBAAiB,EACjB,UAAU,IAAI,WAAW,EACzB,IAAI,IAAI,KAAK,EACb,yBAAyB,EACzB,uBAAuB,EACvB,eAAe,EAEf,YAAY,IAAI,aAAa,GAC9B,MAAM,wBAAwB,CAAC;AAChC,OAAO,EAAE,QAAQ,EAAc,WAAW,EAAE,MAAM,eAAe,CAAC;AAClE,OAAO,EAAa,KAAK,EAAE,MAAM,YAAY,CAAC;AAC9C,OAAO,EAEL,YAAY,GAKb,MAAM,mBAAmB,CAAC;AAC3B,wCAAwC;AACxC,OAAO,EAGL,UAAU,GAKX,MAAM,iBAAiB,CAAC;AACzB,OAAO,EAAE,WAAW,EAAE,MAAM,kBAAkB,CAAC;AAC/C,gCAAgC;AAChC,OAAO,EAAE,QAAQ,EAAE,MAAM,eAAe,CAAC;AACzC,OAAO,EAAE,QAAQ,EAAE,MAAM,eAAe,CAAC;AACzC,OAAO,EACL,oBAAoB,EAEpB,WAAW,EACX,YAAY,EAGZ,WAAW,EACX,YAAY,EACZ,cAAc,GACf,MAAM,cAAc,CAAC;AACtB,OAAO,EAAE,KAAK,EAAE,MAAM,YAAY,CAAC;AACnC,OAAO,EAGL,UAAU,GAGX,MAAM,iBAAiB,CAAC;AACzB,OAAO,EAAE,YAAY,EAAE,MAAM,mBAAmB,CAAC;AACjD,OAAO,EAAiC,SAAS,EAAgB,MAAM,gBAAgB,CAAC;AACxF,OAAO,EAAE,WAAW,EAAE,MAAM,kBAAkB,CAAC;AAC/C,OAAO,EAAE,OAAO,EAAE,MAAM,cAAc,CAAC;AACvC,OAAO,EAAE,cAAc,EAAE,eAAe,EAAE,kBAAkB,EAAE,MAAM,qBAAqB,CAAC;AAC1F,OAAO,EAAE,UAAU,EAAE,MAAM,iBAAiB,CAAC;AAC7C,OAAO,EAAE,QAAQ,EAAE,MAAM,eAAe,CAAC;AACzC,gBAAgB;AAChB,OAAO,KAAK,QAAQ,MAAM,eAAe,CAAC;AAC1C,OAAO,EAAE,MAAM,EAAE,MAAM,aAAa,CAAC;AACrC,OAAO,EAAE,gBAAgB,EAAE,MAAM,uBAAuB,CAAC;AACzD,OAAO,EAAE,eAAe,EAAE,MAAM,sBAAsB,CAAC;AACvD,OAAO,KAAK,UAAU,MAAM,iBAAiB,CAAC;AAC9C,OAAO,EAAE,eAAe,EAAE,MAAM,sBAAsB,CAAC;AACvD,kCAAkC;AAClC,OAAO,EAAE,SAAS,EAAyB,MAAM,gBAAgB,CAAC;AAClE,OAAO,EAIL,kBAAkB,EAGlB,cAAc,EAEd,sBAAsB,EACtB,2BAA2B,EAC3B,iBAAiB,EACjB,oBAAoB,EACpB,kBAAkB,EAClB,iBAAiB,EAGjB,6BAA6B,EAC7B,mBAAmB,EACnB,mBAAmB,EACnB,QAAQ,EACR,yBAAyB,EACzB,uBAAuB,GACxB,MAAM,eAAe,CAAC;AACvB,OAAO,EAAE,IAAI,EAAE,OAAO,EAAE,MAAM,WAAW,CAAC;AAC1C,qBAAqB;AACrB,OAAO,EAAE,SAAS,EAAE,MAAM,gBAAgB,CAAC;AAC3C,OAAO,EAAE,IAAI,EAAE,MAAM,WAAW,CAAC;AACjC,cAAc;AACd,OAAO,EAAE,MAAM,EAAE,MAAM,aAAa,CAAC;AACrC,OAAO,EAAE,cAAc,EAAE,MAAM,qBAAqB,CAAC;AACrD,iBAAiB;AACjB,OAAO,EAAe,WAAW,EAAE,MAAM,kBAAkB,CAAC"}

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "@johnhenry/math",
-  "version": "0.0.2",
+  "version": "0.1.0",
   "publishConfig": {
     "access": "public",
     "provenance": true

package/src/DualNumber.ts

@@ -25,6 +25,16 @@ export class DualNumber {
     return new DualNumber(x, 1);
   }
 
+  /**
+   * Builds a dual number whose value may be outside the function's domain. A
+   * NaN value (e.g. `log` of a negative) makes the derivative NaN too, rather
+   * than leaving a finite-looking derivative (`x.deriv / x.value`) next to a
+   * meaningless value.
+   */
+  private static domain(value: number, deriv: number): DualNumber {
+    return new DualNumber(value, Number.isNaN(value) ? Number.NaN : deriv);
+  }
+
   private static lift(x: DualNumber | number): DualNumber {
     return x instanceof DualNumber ? x : new DualNumber(x, 0);
   }
@@ -62,7 +72,7 @@ export class DualNumber {
     // general formula entirely rather than let it hit that indeterminate
     // form (issue #74).
     if (n === 0) return new DualNumber(1, 0);
-    return new DualNumber(this.value ** n, n * this.value ** (n - 1) * this.deriv);
+    return DualNumber.domain(this.value ** n, n * this.value ** (n - 1) * this.deriv);
   }
 
   static sin(x: DualNumber): DualNumber {
@@ -80,11 +90,11 @@ export class DualNumber {
     return new DualNumber(e, e * x.deriv);
   }
   static log(x: DualNumber): DualNumber {
-    return new DualNumber(Math.log(x.value), x.deriv / x.value);
+    return DualNumber.domain(Math.log(x.value), x.deriv / x.value);
   }
   static sqrt(x: DualNumber): DualNumber {
     const s = Math.sqrt(x.value);
-    return new DualNumber(s, x.deriv / (2 * s));
+    return DualNumber.domain(s, x.deriv / (2 * s));
   }
 
   /** The derivative of `f` at `x`, evaluated exactly via dual numbers. */

package/src/Graph3DUtils.ts

@@ -59,6 +59,12 @@ export interface Placement3D {
 
 const v3 = (x: number, y: number, z = 0): Vec3 => ({ x, y, z });
 const asPoint = (p: Vector<number>): Vec3 => ({ x: p.x as number, y: p.y as number, z: p.z as number });
+const isFinitePoint = (p: Vec3): boolean => Number.isFinite(p.x) && Number.isFinite(p.y) && Number.isFinite(p.z);
+/** A triangle face, or `null` when any vertex is NaN / infinite (so the surface breaks there instead of drawing garbage). */
+const finiteFace = (a: Vector<number>, b: Vector<number>, c: Vector<number>): Face | null => {
+  const face: Face = [asPoint(a), asPoint(b), asPoint(c)];
+  return face.every(isFinitePoint) ? face : null;
+};
 
 export class Graph3DUtils {
   /** Sample `binaryOperation` over a 2D grid, returning a matrix of results. */
@@ -222,6 +228,8 @@ export class Graph3DUtils {
   /**
    * Build a surface mesh from a grid (matrix) of points, as two triangle
    * sweeps with distinct materials (bug fix: the first sweep now uses `alpha1`).
+   * Triangles touching a NaN or infinite vertex are omitted, so the surface
+   * breaks into pieces around them.
    */
   static pointMatrixToMesh3D(
     vector: Vector<Vector<Vector<number>>>,
@@ -239,11 +247,12 @@ export class Graph3DUtils {
     for (let i = 0; i < width - 1; i++) {
       for (let j = 0; j < height - 1; j++) {
         if (at(i, j) && vector[i + 1] && at(i, j + 1)) {
-          faces1.push([
-            asPoint(at(i, j) as Vector<number>),
-            asPoint(at(i + 1, j) as Vector<number>),
-            asPoint(at(i, j + 1) as Vector<number>),
-          ]);
+          const face = finiteFace(
+            at(i, j) as Vector<number>,
+            at(i + 1, j) as Vector<number>,
+            at(i, j + 1) as Vector<number>,
+          );
+          if (face) faces1.push(face);
         }
       }
     }
@@ -252,11 +261,12 @@ export class Graph3DUtils {
     for (let i = 1; i < width; i++) {
       for (let j = 1; j < height; j++) {
         if (at(i, j) && vector[i - 1] && at(i, j - 1)) {
-          faces2.push([
-            asPoint(at(i, j) as Vector<number>),
-            asPoint(at(i - 1, j) as Vector<number>),
-            asPoint(at(i, j - 1) as Vector<number>),
-          ]);
+          const face = finiteFace(
+            at(i, j) as Vector<number>,
+            at(i - 1, j) as Vector<number>,
+            at(i, j - 1) as Vector<number>,
+          );
+          if (face) faces2.push(face);
         }
       }
     }

package/src/GraphUtils.ts

@@ -98,7 +98,10 @@ export class GraphUtils {
     return [...points].map((p) => GraphUtils.barPlacement(p, barWidth));
   }
 
-  /** A poly-line path through the given points. */
+  /**
+   * A poly-line path through the given points. Points with a non-finite x or y
+   * (NaN, +/-Infinity) are skipped and split the line into separate segments.
+   */
   static vectorToCurve(
     points: Vector<Vector<number>>,
     thickness = 0,
@@ -112,13 +115,18 @@ export class GraphUtils {
   ): Path2D {
     const stroke: StrokeStyle = { thickness, color, alpha, pixelHinting, scaleMode, caps, joints, miterLimit };
     const commands: PathCommand[] = [];
-    if (points.length > 0) {
-      const first = point(points[0] as Vector<number>);
-      commands.push({ op: "moveTo", x: first.x, y: first.y });
-      for (let i = 1; i < points.length; i++) {
-        const pt = point(points[i] as Vector<number>);
-        commands.push({ op: "lineTo", x: pt.x, y: pt.y });
+    // A NaN or +/-Infinity coordinate (a pole, a domain error) breaks the path: the
+    // run before it ends, and the next finite point starts a new run with `moveTo`
+    // rather than a `lineTo` drawn through the gap.
+    let penDown = false;
+    for (let i = 0; i < points.length; i++) {
+      const pt = point(points[i] as Vector<number>);
+      if (!Number.isFinite(pt.x) || !Number.isFinite(pt.y)) {
+        penDown = false;
+        continue;
       }
+      commands.push({ op: penDown ? "lineTo" : "moveTo", x: pt.x, y: pt.y });
+      penDown = true;
     }
     return { stroke, commands };
   }

package/src/Structure.ts

@@ -272,9 +272,52 @@ export class Structure<T = unknown> {
     return new RangeError(`${op}: element ${String(x)} has no multiplicative inverse in this structure`);
   }
 
+  /**
+   * Determinant of a square matrix using only addition, subtraction and
+   * multiplication (Berkowitz's algorithm, O(n^4)), so it is valid over any
+   * commutative ring, including Z/nZ with composite `n` where a pivot may have
+   * no multiplicative inverse.
+   */
+  private determinantDivisionFree(alpha: SMat<T>): T {
+    const n = VectorUtils.width(alpha);
+    const a: T[][] = [...alpha].map((row) => [...(row as SVec<T>)]);
+    // p holds the characteristic-polynomial coefficients of the leading r x r block.
+    let p: T[] = [this.one];
+    for (let r = 1; r <= n; r++) {
+      const m = r - 1; // size of the leading block A_{r-1}
+      // Toeplitz column: [1, -a_rr, -R C, -R A C, -R A^2 C, ...]
+      const col: T[] = [this.one, this.negative(a[m][m] as T)];
+      let v: T[] = Array.from({ length: m }, (_, i) => a[i][m] as T); // C
+      for (let k = 2; k <= r; k++) {
+        let dot = this.zero;
+        for (let i = 0; i < m; i++) dot = this.add(dot, this.multiply(a[m][i] as T, v[i]));
+        col.push(this.negative(dot));
+        if (k < r) {
+          v = Array.from({ length: m }, (_, i) => {
+            let acc = this.zero;
+            for (let j = 0; j < m; j++) acc = this.add(acc, this.multiply(a[i][j] as T, v[j]));
+            return acc;
+          });
+        }
+      }
+      const next: T[] = [];
+      for (let i = 0; i <= r; i++) {
+        let acc = this.zero;
+        for (let j = 0; j < r && j <= i; j++) acc = this.add(acc, this.multiply(col[i - j], p[j]));
+        next.push(acc);
+      }
+      p = next;
+    }
+    // p = coefficients of det(xI - A); det(A) = (-1)^n * constant term.
+    const c = p[n];
+    return n % 2 === 0 ? c : this.negative(c);
+  }
+
   /**
    * Determinant via Gaussian elimination with row-swap pivoting, O(n^3)
-   * instead of the previous O(n!) cofactor (Laplace) expansion.
+   * instead of the previous O(n!) cofactor (Laplace) expansion. Over a ring where
+   * a needed pivot has no multiplicative inverse (e.g. Z/nZ, composite `n`) it
+   * falls back to a division-free method, so the result is always defined.
    */
   determinant(alpha: SMat<T>): T {
     if (VectorUtils.width(alpha) !== VectorUtils.height(alpha)) return this.zero;
@@ -297,7 +340,7 @@ export class Structure<T = unknown> {
         }
         if (swap === -1) {
           if (wasZero && mat.slice(col + 1).every((row) => this.equality(row[col] as T, this.zero))) return this.zero;
-          throw this.noInverse("determinant", mat[col][col] as T);
+          return this.determinantDivisionFree(alpha);
         }
         [mat[col], mat[swap]] = [mat[swap], mat[col]];
         sign = this.negative(sign);
@@ -331,8 +374,12 @@ export class Structure<T = unknown> {
    * operations, with a zero-pivot row swap (bug fix — the AS3 version had no
    * pivoting and divided by a zero pivot).
    *
-   * @throws RangeError if the matrix is not square, is singular, or (over a ring
-   *   such as Z/nZ with composite n) a pivot has no multiplicative inverse.
+   * Over a ring such as Z/nZ with composite `n`, when elimination hits a pivot
+   * with no inverse it falls back to `det⁻¹ · adj(A)`, which succeeds exactly
+   * when the determinant is a unit.
+   *
+   * @throws RangeError if the matrix is not square, is singular, or (over a ring)
+   *   its determinant has no multiplicative inverse.
    */
   invertMatrix(alpha: SMat<T>, checkDeterminant = true): SMat<T> {
     const n = VectorUtils.height(alpha);
@@ -341,7 +388,32 @@ export class Structure<T = unknown> {
     if (checkDeterminant && this.equality(this.determinant(alpha), this.zero)) {
       throw new RangeError("invertMatrix: matrix is singular");
     }
+    return this.invertByElimination(alpha) ?? this.invertByAdjugate(alpha);
+  }
+
+  /** `det⁻¹ · adj(A)`; needs no pivot inverses, only that det(A) is a unit. */
+  private invertByAdjugate(alpha: SMat<T>): SMat<T> {
+    const n = VectorUtils.height(alpha);
+    const rows: T[][] = [...alpha].map((row) => [...(row as SVec<T>)]);
+    const det = this.determinant(alpha);
+    if (this.equality(det, this.zero)) throw new RangeError("invertMatrix: matrix is singular");
+    if (!this.hasInverse(det)) throw this.noInverse("invertMatrix", det);
+    const detInv = this.reciprocal(det);
+    const out: T[][] = Array.from({ length: n }, () => new Array<T>(n));
+    for (let i = 0; i < n; i++) {
+      for (let j = 0; j < n; j++) {
+        const sub = rows.filter((_, r) => r !== i).map((row) => Vector.fromArray(row.filter((_, c) => c !== j)));
+        const minor = this.determinant(Vector.fromArray(sub));
+        const cofactor = (i + j) % 2 === 0 ? minor : this.negative(minor);
+        out[j][i] = this.multiply(detInv, cofactor); // adjugate is the transposed cofactor matrix
+      }
+    }
+    return Vector.fromArray(out.map((row) => Vector.fromArray(row)));
+  }
 
+  /** Gauss-Jordan; returns `null` if a pivot has no inverse (caller falls back to the adjugate). */
+  private invertByElimination(alpha: SMat<T>): SMat<T> | null {
+    const n = VectorUtils.height(alpha);
     const mat: T[][] = [...alpha].map((row) => [...(row as SVec<T>)]);
     const inv: T[][] = [...this.generateIdentity(n, n)].map((row) => [...(row as SVec<T>)]);
 
@@ -359,7 +431,7 @@ export class Structure<T = unknown> {
           if (wasZero && mat.slice(col + 1).every((row) => this.equality(row[col] as T, this.zero))) {
             throw new RangeError("invertMatrix: matrix is singular");
           }
-          throw this.noInverse("invertMatrix", mat[col][col] as T);
+          return null;
         }
         [mat[col], mat[swap]] = [mat[swap], mat[col]];
         [inv[col], inv[swap]] = [inv[swap], inv[col]];