@mojir/dvala 0.0.1 → 0.0.2

package/dist/modules/math.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"math.js","sources":["../../src/utils/debug/getCodeMarker.ts","../../src/errors.ts","../../src/constants/constants.ts","../../src/typeGuards/number.ts","../../src/typeGuards/annotatedArrays.ts","../../src/utils/arity.ts","../../src/builtin/modules/math/index.ts"],"sourcesContent":["import type { SourceCodeInfo } from '../../tokenizer/token'\n\nexport function getCodeMarker(sourceCodeInfo: SourceCodeInfo): string {\n  if (!sourceCodeInfo.position || !sourceCodeInfo.code)\n    return ''\n\n  const leftPadding = sourceCodeInfo.position.column - 1\n  const rightPadding = sourceCodeInfo.code.length - leftPadding - 1\n  return `${' '.repeat(Math.max(leftPadding, 0))}^${' '.repeat(Math.max(rightPadding, 0))}`\n}\n","import { getCodeMarker } from '../src/utils/debug/getCodeMarker'\nimport type { Arr } from './interface'\nimport type { SourceCodeInfo } from './tokenizer/token'\n\nfunction getDvalaErrorMessage(message: string, sourceCodeInfo?: SourceCodeInfo) {\n  if (!sourceCodeInfo) {\n    return message\n  }\n  const location = `${sourceCodeInfo.position.line}:${sourceCodeInfo.position.column}`\n  const filePathLine = sourceCodeInfo.filePath\n    ? `\\n${sourceCodeInfo.filePath}:${location}`\n    : `\\nLocation ${location}`\n  const codeLine = `\\n${sourceCodeInfo.code}`\n  const codeMarker = `\\n${getCodeMarker(sourceCodeInfo)}`\n  return `${message}${filePathLine}${codeLine}${codeMarker}`\n}\n\nexport class RecurSignal extends Error {\n  public params: Arr\n  constructor(params: Arr) {\n    super(`recur, params: ${params}`)\n    Object.setPrototypeOf(this, RecurSignal.prototype)\n    this.name = 'RecurSignal'\n    this.params = params\n  }\n}\n\nexport class DvalaError extends Error {\n  public readonly sourceCodeInfo?: SourceCodeInfo\n  public readonly shortMessage: string\n  constructor(err: unknown, sourceCodeInfo: SourceCodeInfo | undefined) {\n    const message = err instanceof Error\n      ? err.message\n      : `${err}`\n\n    super(getDvalaErrorMessage(message, sourceCodeInfo))\n    this.shortMessage = message\n    this.sourceCodeInfo = sourceCodeInfo\n    Object.setPrototypeOf(this, DvalaError.prototype)\n    this.name = 'DvalaError'\n  }\n\n  public getCodeMarker(): string | undefined {\n    return this.sourceCodeInfo && getCodeMarker(this.sourceCodeInfo)\n  }\n}\n\nexport class UserDefinedError extends DvalaError {\n  public userMessage: string\n  constructor(message: string, sourceCodeInfo?: SourceCodeInfo) {\n    super(message, sourceCodeInfo)\n    this.userMessage = message\n    Object.setPrototypeOf(this, UserDefinedError.prototype)\n    this.name = 'UserDefinedError'\n  }\n}\n\nexport class AssertionError extends DvalaError {\n  constructor(message: string | Error, sourceCodeInfo?: SourceCodeInfo) {\n    super(message, sourceCodeInfo)\n    Object.setPrototypeOf(this, AssertionError.prototype)\n    this.name = 'AssertionError'\n  }\n}\n\nexport class UndefinedSymbolError extends DvalaError {\n  public symbol: string\n  constructor(symbolName: string, sourceCodeInfo?: SourceCodeInfo) {\n    const message = `Undefined symbol '${symbolName}'.`\n    super(message, sourceCodeInfo)\n    this.symbol = symbolName\n    Object.setPrototypeOf(this, UndefinedSymbolError.prototype)\n    this.name = 'UndefinedSymbolError'\n  }\n}\n\nexport function isDvalaError(error: unknown): error is DvalaError {\n  return error instanceof DvalaError\n}\n","export const NodeTypes = {\n  Number: 1,\n  String: 2,\n  NormalExpression: 3,\n  SpecialExpression: 4,\n  UserDefinedSymbol: 5,\n  NormalBuiltinSymbol: 6,\n  SpecialBuiltinSymbol: 7,\n  ReservedSymbol: 8,\n  Binding: 9,\n  Spread: 10,\n} as const\n\nconst NodeTypesSet = new Set(Object.values(NodeTypes))\n\nexport type NodeType = typeof NodeTypes[keyof typeof NodeTypes]\n\nexport function getNodeTypeName(type: NodeType): keyof typeof NodeTypes {\n  return Object.keys(NodeTypes).find(key => NodeTypes[key as keyof typeof NodeTypes] === type) as keyof typeof NodeTypes\n}\n\n// TODO, is this needed?\nexport function isNodeType(type: unknown): type is NodeType {\n  return typeof type === 'number' && NodeTypesSet.has(type as NodeType)\n}\n\nconst functionTypes = [\n  'UserDefined',\n  'Partial',\n  'Comp',\n  'Constantly',\n  'Juxt',\n  'Complement',\n  'EveryPred',\n  'SomePred',\n  'Fnull',\n  'Builtin',\n  'SpecialBuiltin',\n  'NativeJsFunction',\n  'Module',\n] as const\n\nconst functionTypeSet = new Set(functionTypes)\n\nexport type FunctionType = typeof functionTypes[number]\n\nexport function isFunctionType(type: unknown): type is FunctionType {\n  return typeof type === 'string' && functionTypeSet.has(type as FunctionType)\n}\n","import { DvalaError } from '../errors'\nimport type { SourceCodeInfo } from '../tokenizer/token'\nimport { valueToString } from '../utils/debug/debugTools'\nimport { getSourceCodeInfo } from '../utils/debug/getSourceCodeInfo'\n\ntype SignOptions =\n  | {\n    positive?: true\n    negative?: never\n    nonPositive?: never\n    nonNegative?: never\n    zero?: never\n    nonZero?: never\n  }\n  | {\n    positive?: never\n    negative?: true\n    nonPositive?: never\n    nonNegative?: never\n    zero?: never\n    nonZero?: never\n  }\n  | {\n    positive?: never\n    negative?: never\n    nonPositive?: true\n    nonNegative?: never\n    zero?: never\n    nonZero?: never\n  }\n  | {\n    positive?: never\n    negative?: never\n    nonPositive?: never\n    nonNegative?: true\n    zero?: never\n    nonZero?: never\n  }\n  | {\n    positive?: never\n    negative?: never\n    nonPositive?: never\n    nonNegative?: never\n    zero?: true\n    nonZero?: never\n  }\n  | {\n    positive?: never\n    negative?: never\n    nonPositive?: never\n    nonNegative?: never\n    zero?: never\n    nonZero?: true\n  }\n\ntype GtOptions =\n  | {\n    gt?: number\n    gte?: never\n  }\n  | {\n    gt?: never\n    gte?: number\n  }\n\ntype LtOptions =\n  | {\n    lt?: number\n    lte?: never\n  }\n  | {\n    lt?: never\n    lte?: number\n  }\n\ntype NumberOptions = {\n  integer?: true\n  finite?: true\n} & SignOptions &\nGtOptions &\nLtOptions\n\nfunction getRangeString(options: NumberOptions): string {\n  const hasUpperAndLowerBound\n    = (typeof options.gt === 'number' || typeof options.gte === 'number')\n      && (typeof options.lt === 'number' || typeof options.lte === 'number')\n  if (hasUpperAndLowerBound) {\n    return `${typeof options.gt === 'number' ? `${options.gt} < n ` : `${options.gte} <= n `}${\n      typeof options.lt === 'number' ? `< ${options.lt}` : `<= ${options.lte}`\n    }`\n  }\n  else if (typeof options.gt === 'number' || typeof options.gte === 'number') {\n    return `${typeof options.gt === 'number' ? `n > ${options.gt}` : `n >= ${options.gte}`}`\n  }\n  else if (typeof options.lt === 'number' || typeof options.lte === 'number') {\n    return `${typeof options.lt === 'number' ? `n < ${options.lt}` : `n <= ${options.lte}`}`\n  }\n  else { return '' }\n}\n\nfunction getSignString(options: NumberOptions): string {\n  return options.positive\n    ? 'positive'\n    : options.negative\n      ? 'negative'\n      : options.nonNegative\n        ? 'non negative'\n        : options.nonPositive\n          ? 'non positive'\n          : options.nonZero\n            ? 'non zero'\n            : ''\n}\n\nfunction getNumberTypeName(options: NumberOptions): string {\n  if (options.zero)\n    return 'zero'\n\n  const sign = getSignString(options)\n  const numberType = options.integer ? 'integer' : 'number'\n  const finite = options.finite ? 'finite' : ''\n  const range = getRangeString(options)\n\n  return [sign, finite, numberType, range].filter(x => !!x).join(' ')\n}\n\nexport function isNumber(value: unknown, options: NumberOptions = {}): value is number {\n  if (typeof value !== 'number')\n    return false\n\n  if (Number.isNaN(value))\n    return false\n\n  if (options.integer && !Number.isInteger(value))\n    return false\n\n  if (options.finite && !Number.isFinite(value))\n    return false\n\n  if (options.zero && value !== 0)\n    return false\n\n  if (options.nonZero && value === 0)\n    return false\n\n  if (options.positive && value <= 0)\n    return false\n\n  if (options.negative && value >= 0)\n    return false\n\n  if (options.nonPositive && value > 0)\n    return false\n\n  if (options.nonNegative && value < 0)\n    return false\n\n  if (typeof options.gt === 'number' && value <= options.gt)\n    return false\n\n  if (typeof options.gte === 'number' && value < options.gte)\n    return false\n\n  if (typeof options.lt === 'number' && value >= options.lt)\n    return false\n\n  if (typeof options.lte === 'number' && value > options.lte)\n    return false\n\n  return true\n}\n\nexport function assertNumber(\n  value: unknown,\n  sourceCodeInfo?: SourceCodeInfo,\n  options: NumberOptions = {},\n): asserts value is number {\n  if (!isNumber(value, options)) {\n    throw new DvalaError(\n      `Expected ${getNumberTypeName(options)}, got ${valueToString(value)}.`,\n      getSourceCodeInfo(value, sourceCodeInfo),\n    )\n  }\n}\n\nexport function asNumber(\n  value: unknown,\n  sourceCodeInfo: SourceCodeInfo | undefined,\n  options: NumberOptions = {},\n): number {\n  assertNumber(value, sourceCodeInfo, options)\n  return value\n}\n","import { DvalaError } from '../errors'\nimport type { Any } from '../interface'\nimport type { SourceCodeInfo } from '../tokenizer/token'\nimport { isNumber } from './number'\n\nconst annotatedArrays = new WeakSet<unknown[]>()\nconst vectors = new WeakSet<unknown[]>()\nconst notVectors = new WeakSet<unknown[]>()\nconst matrices = new WeakSet<unknown[]>()\nconst notMatrices = new WeakSet<unknown[]>()\nconst grids = new WeakSet<unknown[]>()\nconst notGrids = new WeakSet<unknown[]>()\n\nexport function annotate<T>(value: T): T {\n  if (!Array.isArray(value)) {\n    return value\n  }\n  if (annotatedArrays.has(value)) {\n    return value\n  }\n  isVector(value)\n  if (!isMatrix(value)) {\n    isGrid(value)\n  }\n\n  return value\n}\nexport function isVector(vector: unknown): vector is number[] {\n  if (!Array.isArray(vector)) {\n    return false\n  }\n\n  if (vectors.has(vector)) {\n    return true\n  }\n  if (notVectors.has(vector)) {\n    return false\n  }\n\n  if (vector.every(elem => isNumber(elem))) {\n    annotatedArrays.add(vector)\n    vectors.add(vector)\n    return true\n  }\n  notVectors.add(vector)\n  return false\n}\n\nexport function assertVector(vector: unknown, sourceCodeInfo: SourceCodeInfo | undefined): asserts vector is number[] {\n  if (!isVector(vector)) {\n    throw new DvalaError(`Expected a vector, but got ${vector}`, sourceCodeInfo)\n  }\n}\n\nexport function is2dVector(vector: unknown): vector is [number, number] {\n  if (!isVector(vector)) {\n    return false\n  }\n  return vector.length === 2\n}\nexport function assert2dVector(vector: unknown, sourceCodeInfo: SourceCodeInfo | undefined): asserts vector is [number, number] {\n  if (!is2dVector(vector)) {\n    throw new DvalaError(`Expected a 2d vector, but got ${vector}`, sourceCodeInfo)\n  }\n}\n\nexport function is3dVector(vector: unknown): vector is [number, number, number] {\n  if (!isVector(vector)) {\n    return false\n  }\n  return vector.length === 3\n}\nexport function assert3dVector(vector: unknown, sourceCodeInfo: SourceCodeInfo | undefined): asserts vector is [number, number, number] {\n  if (!is3dVector(vector)) {\n    throw new DvalaError(`Expected a 3d vector, but got ${vector}`, sourceCodeInfo)\n  }\n}\n\nexport function assertNonEmptyVector(vector: unknown, sourceCodeInfo: SourceCodeInfo | undefined): asserts vector is number[] {\n  assertVector(vector, sourceCodeInfo)\n  if (vector.length === 0) {\n    throw new DvalaError(`Expected a non empty vector, but got ${vector}`, sourceCodeInfo)\n  }\n}\n\nexport function isGrid(grid: unknown): grid is unknown[][] {\n  if (!Array.isArray(grid)) {\n    return false\n  }\n  if (grids.has(grid)) {\n    return true\n  }\n  if (notGrids.has(grid)) {\n    return false\n  }\n  if (grid.length === 0) {\n    notGrids.add(grid)\n    return false\n  }\n  if (!Array.isArray(grid[0])) {\n    notGrids.add(grid)\n    return false\n  }\n  const nbrOfCols = grid[0].length\n  for (const row of grid.slice(1)) {\n    if (!Array.isArray(row)) {\n      notGrids.add(grid)\n      return false\n    }\n    if (row.length !== nbrOfCols) {\n      notGrids.add(grid)\n      return false\n    }\n  }\n  annotatedArrays.add(grid)\n  grids.add(grid)\n  return true\n}\n\nexport function assertGrid(grid: unknown, sourceCodeInfo: SourceCodeInfo | undefined): asserts grid is Any[][] {\n  if (!isGrid(grid)) {\n    throw new DvalaError(`Expected a grid, but got ${grid}`, sourceCodeInfo)\n  }\n}\n\nexport function isMatrix(matrix: unknown): matrix is number[][] {\n  if (!Array.isArray(matrix)) {\n    return false\n  }\n  if (matrices.has(matrix)) {\n    return true\n  }\n  if (notMatrices.has(matrix)) {\n    return false\n  }\n  if (matrix.length === 0) {\n    notMatrices.add(matrix)\n    return false\n  }\n  if (!Array.isArray(matrix[0]) || matrix[0].length === 0) {\n    notMatrices.add(matrix)\n    return false\n  }\n  const nbrOfCols = matrix[0].length\n  for (const row of matrix) {\n    if (!Array.isArray(row) || row.length !== nbrOfCols || row.some(cell => !isNumber(cell))) {\n      notMatrices.add(matrix)\n      return false\n    }\n  }\n  annotatedArrays.add(matrix)\n  grids.add(matrix)\n  matrices.add(matrix)\n  return true\n}\n\nexport function assertMatrix(matrix: unknown, sourceCodeInfo: SourceCodeInfo | undefined): asserts matrix is number[][] {\n  if (!isMatrix(matrix)) {\n    throw new DvalaError(`Expected a matrix, but got ${matrix}`, sourceCodeInfo)\n  }\n}\n\nexport function assertSquareMatrix(matrix: unknown, sourceCodeInfo: SourceCodeInfo | undefined): asserts matrix is number[][] {\n  if (!isMatrix(matrix)) {\n    throw new DvalaError(`Expected a matrix, but got ${matrix}`, sourceCodeInfo)\n  }\n  if (matrix.length !== matrix[0]!.length) {\n    throw new DvalaError(`Expected square matrix, but got ${matrix.length} and ${matrix[0]!.length}`, sourceCodeInfo)\n  }\n}\n\nexport function isSquareMatrix(matrix: unknown): matrix is number[][] {\n  if (!isMatrix(matrix)) {\n    return false\n  }\n  if (matrix.length !== matrix[0]!.length) {\n    return false\n  }\n  return true\n}\n","import type { Arity } from '../builtin/interface'\nimport { DvalaError } from '../errors'\nimport type { FunctionLike } from '../parser/types'\nimport type { SourceCodeInfo } from '../tokenizer/token'\nimport { isColl } from '../typeGuards/dvala'\nimport { valueToString } from './debug/debugTools'\n\nexport function arityAccepts(arity: Arity, nbrOfParams: number): boolean {\n  const { min, max } = arity\n  if (typeof min === 'number' && nbrOfParams < min) {\n    return false\n  }\n  if (typeof max === 'number' && nbrOfParams > max) {\n    return false\n  }\n  return true\n}\n\nexport function arityAcceptsMin(arity: Arity, nbrOfParams: number): boolean {\n  const { min } = arity\n  if (typeof min === 'number' && nbrOfParams < min) {\n    return false\n  }\n  return true\n}\n\nexport function getCommonArityFromFunctions(params: FunctionLike[]): Arity | null {\n  return params.reduce((acc: Arity | null, param): Arity | null => {\n    if (acc === null) {\n      return null\n    }\n    const arity: Arity = (typeof param === 'number' || isColl(param)) ? toFixedArity(1) : param.arity\n    const { min: aMin, max: aMax } = arity\n    const { min: bMin, max: bMax } = acc\n    const min = typeof aMin === 'number' && typeof bMin === 'number'\n      ? Math.max(aMin, bMin)\n      : typeof aMin === 'number' ? aMin : typeof bMin === 'number' ? bMin : undefined\n    const max = typeof aMax === 'number' && typeof bMax === 'number'\n      ? Math.min(aMax, bMax)\n      : typeof aMax === 'number' ? aMax : typeof bMax === 'number' ? bMax : undefined\n\n    if (typeof min === 'number' && typeof max === 'number' && min > max) {\n      return null\n    }\n\n    return { min, max }\n  }, {})\n}\n\nexport function getArityFromFunction(param: FunctionLike): Arity {\n  return (typeof param === 'number' || isColl(param)) ? toFixedArity(1) : param.arity\n}\n\nexport function assertNumberOfParams(arity: Arity, length: number, sourceCodeInfo: SourceCodeInfo | undefined): void {\n  const { min, max } = arity\n  if (typeof min === 'number' && length < min) {\n    throw new DvalaError(\n      `Wrong number of arguments, expected at least ${min}, got ${valueToString(length)}.`,\n      sourceCodeInfo,\n    )\n  }\n\n  if (typeof max === 'number' && length > max) {\n    throw new DvalaError(\n      `Wrong number of arguments, expected at most ${max}, got ${valueToString(length)}.`,\n      sourceCodeInfo,\n    )\n  }\n}\n\nexport function canBeOperator(count: Arity): boolean {\n  if (typeof count.max === 'number' && count.max < 2) {\n    return false\n  }\n\n  if (typeof count.min === 'number' && count.min > 2) {\n    return false\n  }\n\n  return true\n}\n\nexport function toFixedArity(arity: number): Arity {\n  return { min: arity, max: arity }\n}\n","import type { SourceCodeInfo } from '../../../tokenizer/token'\nimport { isMatrix, isVector } from '../../../typeGuards/annotatedArrays'\nimport { isNumber } from '../../../typeGuards/number'\nimport { toFixedArity } from '../../../utils/arity'\nimport { DvalaError } from '../../../errors'\nimport type { BuiltinNormalExpressions } from '../../interface'\nimport type { DvalaModule } from '../interface'\n\ntype NumberVectorOrMatrix = number | number[] | number[][]\n\nfunction getNumberVectorOrMatrixOperation(\n  param: unknown,\n  sourceCodeInfo: SourceCodeInfo | undefined,\n):\n  | ['number', number]\n  | ['vector', number[]]\n  | ['matrix', number[][]] {\n  if (isVector(param)) {\n    return ['vector', param]\n  }\n  if (isMatrix(param)) {\n    return ['matrix', param]\n  }\n  if (!isNumber(param)) {\n    throw new DvalaError(`Invalid parameter type: ${typeof param}`, sourceCodeInfo)\n  }\n  return ['number', param]\n}\n\nfunction unaryMathOp(\n  fn: (val: number) => number,\n): (params: unknown[], sourceCodeInfo: SourceCodeInfo | undefined) => NumberVectorOrMatrix {\n  return ([param], sourceCodeInfo) => {\n    const [operation, operand] = getNumberVectorOrMatrixOperation(param, sourceCodeInfo)\n    if (operation === 'number') {\n      return fn(operand)\n    }\n    else if (operation === 'vector') {\n      return operand.map(val => fn(val))\n    }\n    else {\n      return operand.map(row => row.map(val => fn(val)))\n    }\n  }\n}\n\nconst mathUtilsFunctions: BuiltinNormalExpressions = {\n  'sin': {\n    evaluate: unaryMathOp(val => Math.sin(val)),\n    arity: toFixedArity(1),\n    docs: {\n      category: 'math',\n      returns: { type: ['number', 'vector', 'matrix'] },\n      args: {\n        x: { type: ['number', 'vector', 'matrix'] },\n      },\n      variants: [{ argumentNames: ['x'] }],\n      description: 'The `sin` function computes the sine of an angle (in radians), working on `numbers` and element-wise on `vectors` and `matrices`. When applied to collections, it returns the sine of each element while preserving the original structure.',\n      seeAlso: ['math.asin', 'math.sinh', 'math.cos', 'math.tan', 'math.to-rad'],\n      examples: [\n        `let { sin } = import(math);\nsin(0)`,\n        `let { sin } = import(math);\nsin(1)`,\n        `let { sin } = import(math);\nsin(PI)`,\n        `let { sin } = import(math);\nsin(-0.5)`,\n        `let { sin } = import(math);\nsin([1, 2, 3])`,\n        `let { sin } = import(math);\nsin([[1, 2], [3, 4]])`,\n      ],\n    },\n  },\n  'asin': {\n    evaluate: unaryMathOp(val => Math.asin(val)),\n    arity: toFixedArity(1),\n    docs: {\n      category: 'math',\n      returns: { type: ['number', 'vector', 'matrix'] },\n      args: {\n        x: { type: ['number', 'vector', 'matrix'] },\n      },\n      variants: [{ argumentNames: ['x'] }],\n      description: 'The `asin` function computes the arcsine (inverse sine) of a `number` in radians, working on `numbers` and element-wise on `vectors` and `matrices`. When applied to collections, it returns the arcsine of each element while preserving the original structure.',\n      seeAlso: ['math.sin', 'math.asinh', 'math.acos', 'math.atan'],\n      examples: [\n        `let { asin } = import(math);\nasin(0)`,\n        `let { asin } = import(math);\nasin(1)`,\n        `let { asin } = import(math);\nasin(-0.5)`,\n        `let { asin } = import(math);\nasin([1, 2, 3])`,\n        `let { asin } = import(math);\nasin([[1, 2], [3, 4]])`,\n      ],\n    },\n  },\n  'sinh': {\n    evaluate: unaryMathOp(val => Math.sinh(val)),\n    arity: toFixedArity(1),\n    docs: {\n      category: 'math',\n      returns: { type: ['number', 'vector', 'matrix'] },\n      args: {\n        x: { type: ['number', 'vector', 'matrix'] },\n      },\n      variants: [{ argumentNames: ['x'] }],\n      description: 'The `sinh` function computes the hyperbolic sine of a `number`, working on `numbers` and element-wise on `vectors` and `matrices`. When applied to collections, it returns the hyperbolic sine of each element while preserving the original structure.',\n      seeAlso: ['math.asinh', 'math.sin', 'math.cosh', 'math.tanh'],\n      examples: [\n        `let { sinh } = import(math);\nsinh(0)`,\n        `let { sinh } = import(math);\nsinh(1)`,\n        `let { sinh } = import(math);\nsinh(-0.5)`,\n        `let { sinh } = import(math);\nsinh([0.1, 0.2, 0.3])`,\n        `let { sinh } = import(math);\nsinh([[0.1, 0.2], [0.3, 0.4]])`,\n      ],\n    },\n  },\n  'asinh': {\n    evaluate: unaryMathOp(val => Math.asinh(val)),\n    arity: toFixedArity(1),\n    docs: {\n      category: 'math',\n      returns: { type: ['number', 'vector', 'matrix'] },\n      args: {\n        x: { type: ['number', 'vector', 'matrix'] },\n      },\n      variants: [{ argumentNames: ['x'] }],\n      description: 'The `asinh` function computes the inverse hyperbolic sine of a `number`, working on `numbers` and element-wise on `vectors` and `matrices`. When applied to collections, it returns the inverse hyperbolic sine of each element while preserving the original structure.',\n      seeAlso: ['math.sinh', 'math.asin', 'math.acosh', 'math.atanh'],\n      examples: [\n        `let { asinh } = import(math);\nasinh(10)`,\n        `let { asinh } = import(math);\nasinh(90)`,\n        `let { asinh } = import(math);\nasinh(50)`,\n        `let { asinh } = import(math);\nasinh([10, 20, 30])`,\n        `let { asinh } = import(math);\nasinh([[10, 20], [30, 40]])`,\n      ],\n    },\n  },\n  'cos': {\n    evaluate: unaryMathOp(val => Math.cos(val)),\n    arity: toFixedArity(1),\n    docs: {\n      category: 'math',\n      returns: { type: ['number', 'vector', 'matrix'] },\n      args: {\n        x: { type: ['number', 'vector', 'matrix'] },\n      },\n      variants: [{ argumentNames: ['x'] }],\n      description: 'The `cos` function computes the cosine of an angle (in radians), working on `numbers` and element-wise on `vectors` and `matrices`. When applied to collections, it returns the cosine of each element while preserving the original structure.',\n      seeAlso: ['math.acos', 'math.cosh', 'math.sin', 'math.tan', 'math.to-rad'],\n      examples: [\n        `let { cos } = import(math);\ncos(0)`,\n        `let { cos } = import(math);\ncos(1)`,\n        `let { cos } = import(math);\ncos(PI)`,\n        `let { cos } = import(math);\ncos(-0.5)`,\n        `let { cos } = import(math);\ncos([1, 2, 3])`,\n        `let { cos } = import(math);\ncos([[1, 2], [3, 4]])`,\n      ],\n    },\n  },\n  'acos': {\n    evaluate: unaryMathOp(val => Math.acos(val)),\n    arity: toFixedArity(1),\n    docs: {\n      category: 'math',\n      returns: { type: ['number', 'vector', 'matrix'] },\n      args: {\n        x: { type: ['number', 'vector', 'matrix'] },\n      },\n      variants: [{ argumentNames: ['x'] }],\n      description: 'The `acos` function computes the arccosine (inverse cosine) of a `number` in radians, working on `numbers` and element-wise on `vectors` and `matrices`. When applied to collections, it returns the arccosine of each element while preserving the original structure.',\n      seeAlso: ['math.cos', 'math.acosh', 'math.asin', 'math.atan'],\n      examples: [\n        `let { acos } = import(math);\nacos(0)`,\n        `let { acos } = import(math);\nacos(1)`,\n        `let { acos } = import(math);\nacos(-0.5)`,\n        `let { acos } = import(math);\nacos([0.1, 0.2, 0.3])`,\n        `let { acos } = import(math);\nacos([[0.1, 0.2], [0.3, 0.4]])`,\n      ],\n    },\n  },\n  'cosh': {\n    evaluate: unaryMathOp(val => Math.cosh(val)),\n    arity: toFixedArity(1),\n    docs: {\n      category: 'math',\n      returns: { type: ['number', 'vector', 'matrix'] },\n      args: {\n        x: { type: ['number', 'vector', 'matrix'] },\n      },\n      variants: [{ argumentNames: ['x'] }],\n      description: 'The `cosh` function computes the hyperbolic cosine of a `number`, working on `numbers` and element-wise on `vectors` and `matrices`. When applied to collections, it returns the hyperbolic cosine of each element while preserving the original structure.',\n      seeAlso: ['math.acosh', 'math.cos', 'math.sinh', 'math.tanh'],\n      examples: [\n        `let { cosh } = import(math);\ncosh(0)`,\n        `let { cosh } = import(math);\ncosh(1)`,\n        `let { cosh } = import(math);\ncosh(-0.5)`,\n        `let { cosh } = import(math);\ncosh([0.1, 0.2, 0.3])`,\n        `let { cosh } = import(math);\ncosh([[0.1, 0.2], [0.3, 0.4]])`,\n      ],\n    },\n  },\n  'acosh': {\n    evaluate: unaryMathOp(val => Math.acosh(val)),\n    arity: toFixedArity(1),\n    docs: {\n      category: 'math',\n      returns: { type: ['number', 'vector', 'matrix'] },\n      args: {\n        x: { type: ['number', 'vector', 'matrix'] },\n      },\n      variants: [{ argumentNames: ['x'] }],\n      description: 'The `acosh` function computes the inverse hyperbolic cosine of a `number`, working on `numbers` and element-wise on `vectors` and `matrices`. When applied to collections, it returns the inverse hyperbolic cosine of each element while preserving the original structure.',\n      seeAlso: ['math.cosh', 'math.acos', 'math.asinh', 'math.atanh'],\n      examples: [\n        `let { acosh } = import(math);\nacosh(1)`,\n        `let { acosh } = import(math);\nacosh(2)`,\n        `let { acosh } = import(math);\nacosh(100)`,\n        `let { acosh } = import(math);\nacosh(50)`,\n        `let { acosh } = import(math);\nacosh([1, 2, 3])`,\n        `let { acosh } = import(math);\nacosh([[1, 2], [3, 4]])`,\n      ],\n    },\n  },\n  'tan': {\n    evaluate: unaryMathOp(val => Math.tan(val)),\n    arity: toFixedArity(1),\n    docs: {\n      category: 'math',\n      returns: { type: ['number', 'vector', 'matrix'] },\n      args: {\n        x: { type: ['number', 'vector', 'matrix'] },\n      },\n      variants: [{ argumentNames: ['x'] }],\n      description: 'The `tan` function computes the tangent of an angle (in radians), working on `numbers` and element-wise on `vectors` and `matrices`. When applied to collections, it returns the tangent of each element while preserving the original structure.',\n      seeAlso: ['math.atan', 'math.tanh', 'math.sin', 'math.cos', 'math.to-rad'],\n      examples: [\n        `let { tan } = import(math);\ntan(0)`,\n        `let { tan } = import(math);\ntan(1)`,\n        `let { tan } = import(math);\ntan(PI)`,\n        `let { tan } = import(math);\ntan(-0.5)`,\n        `let { tan } = import(math);\ntan([1, 2, 3])`,\n        `let { tan } = import(math);\ntan([[1, 2], [3, 4]])`,\n      ],\n    },\n  },\n  'atan': {\n    evaluate: unaryMathOp(val => Math.atan(val)),\n    arity: toFixedArity(1),\n    docs: {\n      category: 'math',\n      returns: { type: ['number', 'vector', 'matrix'] },\n      args: {\n        x: { type: ['number', 'vector', 'matrix'] },\n      },\n      variants: [{ argumentNames: ['x'] }],\n      description: 'The `atan` function computes the arctangent (inverse tangent) of a `number` in radians, working on `numbers` and element-wise on `vectors` and `matrices`. When applied to collections, it returns the arctangent of each element while preserving the original structure.',\n      seeAlso: ['math.tan', 'math.atanh', 'math.asin', 'math.acos'],\n      examples: [\n        `let { atan } = import(math);\natan(0)`,\n        `let { atan } = import(math);\natan(1)`,\n        `let { atan } = import(math);\natan(-0.5)`,\n        `let { atan } = import(math);\natan([0.1, 0.2, 0.3])`,\n        `let { atan } = import(math);\natan([[0.1, 0.2], [0.3, 0.4]])`,\n      ],\n    },\n  },\n  'tanh': {\n    evaluate: unaryMathOp(val => Math.tanh(val)),\n    arity: toFixedArity(1),\n    docs: {\n      category: 'math',\n      returns: { type: ['number', 'vector', 'matrix'] },\n      args: {\n        x: { type: ['number', 'vector', 'matrix'] },\n      },\n      variants: [{ argumentNames: ['x'] }],\n      description: 'The `tanh` function computes the hyperbolic tangent of a `number`, working on `numbers` and element-wise on `vectors` and `matrices`. When applied to collections, it returns the hyperbolic tangent of each element while preserving the original structure.',\n      seeAlso: ['math.atanh', 'math.tan', 'math.sinh', 'math.cosh'],\n      examples: [\n        `let { tanh } = import(math);\ntanh(0)`,\n        `let { tanh } = import(math);\ntanh(1)`,\n        `let { tanh } = import(math);\ntanh(-0.5)`,\n        `let { tanh } = import(math);\ntanh(50)`,\n      ],\n    },\n  },\n  'atanh': {\n    evaluate: unaryMathOp(val => Math.atanh(val)),\n    arity: toFixedArity(1),\n    docs: {\n      category: 'math',\n      returns: { type: ['number', 'vector', 'matrix'] },\n      args: {\n        x: { type: ['number', 'vector', 'matrix'] },\n      },\n      variants: [{ argumentNames: ['x'] }],\n      description: 'The `atanh` function computes the inverse hyperbolic tangent of a `number`, working on `numbers` and element-wise on `vectors` and `matrices`. When applied to collections, it returns the inverse hyperbolic tangent of each element while preserving the original structure.',\n      seeAlso: ['math.tanh', 'math.atan', 'math.asinh', 'math.acosh'],\n      examples: [\n        `let { atanh } = import(math);\natanh(0)`,\n        `let { atanh } = import(math);\natanh(0.9)`,\n        `let { atanh } = import(math);\natanh(-0.5)`,\n        `let { atanh } = import(math);\natanh([0.1, 0.2, 0.3])`,\n        `let { atanh } = import(math);\natanh([[0.1, 0.2], [0.3, 0.4]])`,\n      ],\n    },\n  },\n  'ln': {\n    evaluate: unaryMathOp(val => Math.log(val)),\n    arity: toFixedArity(1),\n    docs: {\n      category: 'math',\n      returns: { type: ['number', 'vector', 'matrix'] },\n      args: {\n        x: { type: ['number', 'vector', 'matrix'] },\n      },\n      variants: [{ argumentNames: ['x'] }],\n      description: 'The `ln` function computes the natural logarithm (base `e`) of a `number`, working on `numbers` and element-wise on `vectors` and `matrices`. When applied to collections, it returns the natural logarithm of each element while preserving the original structure.',\n      seeAlso: ['math.log2', 'math.log10', '^'],\n      examples: [\n        `let { ln } = import(math);\nln(0.01)`,\n        `let { ln } = import(math);\nln(2.5)`,\n        `let { ln } = import(math);\nln(E)`,\n        `let { ln } = import(math);\nln([1, 2, 3])`,\n        `let { ln } = import(math);\nln([[1, 2], [3, 4]])`,\n      ],\n    },\n  },\n  'log2': {\n    evaluate: unaryMathOp(val => Math.log2(val)),\n    arity: toFixedArity(1),\n    docs: {\n      category: 'math',\n      returns: { type: ['number', 'vector', 'matrix'] },\n      args: {\n        x: { type: ['number', 'vector', 'matrix'] },\n      },\n      variants: [{ argumentNames: ['x'] }],\n      description: 'The `log2` function computes the base `2` logarithm of a `number`, working on `numbers` and element-wise on `vectors` and `matrices`. When applied to collections, it returns the base-2 logarithm of each element while preserving the original structure.',\n      seeAlso: ['math.ln', 'math.log10'],\n      examples: [\n        `let { log2 } = import(math);\nlog2(0.01)`,\n        `let { log2 } = import(math);\nlog2(2 ^ 12)`,\n        `let { log2 } = import(math);\nlog2(2.5)`,\n        `let { log2 } = import(math);\nlog2([1, 2, 3])`,\n        `let { log2 } = import(math);\nlog2([[1, 2], [3, 4]])`,\n      ],\n    },\n  },\n  'log10': {\n    evaluate: unaryMathOp(val => Math.log10(val)),\n    arity: toFixedArity(1),\n    docs: {\n      category: 'math',\n      returns: { type: ['number', 'vector', 'matrix'] },\n      args: {\n        x: { type: ['number', 'vector', 'matrix'] },\n      },\n      variants: [{ argumentNames: ['x'] }],\n      description: 'The `log10` function computes the base `10` logarithm of a `number`, working on `numbers` and element-wise on `vectors` and `matrices`. When applied to collections, it returns the base-10 logarithm of each element while preserving the original structure.',\n      seeAlso: ['math.ln', 'math.log2'],\n      examples: [\n        `let { log10 } = import(math);\nlog10(0.01)`,\n        `let { log10 } = import(math);\nlog10(10 ^ 12)`,\n        `let { log10 } = import(math);\nlog10(2.5)`,\n        `let { log10 } = import(math);\nlog10([1, 2, 3])`,\n        `let { log10 } = import(math);\nlog10([[1, 2], [3, 4]])`,\n      ],\n    },\n  },\n  'to-rad': {\n    evaluate: unaryMathOp(val => (val * Math.PI) / 180),\n    arity: toFixedArity(1),\n    docs: {\n      category: 'math',\n      returns: { type: ['number', 'vector', 'matrix'] },\n      args: {\n        x: { type: ['number', 'vector', 'matrix'] },\n      },\n      variants: [{ argumentNames: ['x'] }],\n      description: 'The `to-rad` function converts an angle from degrees to radians, working on `numbers` and element-wise on `vectors` and `matrices`. When applied to collections, it converts each element while preserving the original structure.',\n      seeAlso: ['math.to-deg', 'math.sin', 'math.cos', 'math.tan'],\n      examples: [\n        `let { to-rad } = import(math);\nto-rad(0)`,\n        `let { to-rad } = import(math);\nto-rad(90)`,\n        `let { to-rad } = import(math);\nto-rad(180)`,\n        `let { to-rad } = import(math);\nto-rad(360)`,\n        `let { to-rad } = import(math);\nto-rad([0, 90, 180])`,\n        `let { to-rad } = import(math);\nto-rad([[0, 90], [180, 360]])`,\n      ],\n    },\n  },\n  'to-deg': {\n    evaluate: unaryMathOp(val => (val * 180) / Math.PI),\n    arity: toFixedArity(1),\n    docs: {\n      category: 'math',\n      returns: { type: ['number', 'vector', 'matrix'] },\n      args: {\n        x: { type: ['number', 'vector', 'matrix'] },\n      },\n      variants: [{ argumentNames: ['x'] }],\n      description: 'The `to-deg` function converts an angle from radians to degrees, working on `numbers` and element-wise on `vectors` and `matrices`. When applied to collections, it converts each element while preserving the original structure.',\n      seeAlso: ['math.to-rad'],\n      examples: [\n        `let { to-deg } = import(math);\nto-deg(0)`,\n        `let { to-deg } = import(math);\nto-deg(PI)`,\n        `let { to-deg } = import(math);\nto-deg(PI / 2)`,\n        `let { to-deg } = import(math);\nto-deg(3 * PI / 2)`,\n        `let { to-deg } = import(math);\nto-deg([0, PI, PI / 2])`,\n        `let { to-deg } = import(math);\nto-deg([[0, PI], [PI / 2, 3 * PI / 2]])`,\n      ],\n    },\n  },\n}\n\nexport const mathUtilsModule: DvalaModule = {\n  name: 'math',\n  functions: mathUtilsFunctions,\n}\n"],"names":["getCodeMarker","sourceCodeInfo","position","code","leftPadding","column","rightPadding","length","repeat","Math","max","DvalaError","Error","shortMessage","constructor","err","message","super","location","line","filePath","getDvalaErrorMessage","this","Object","setPrototypeOf","prototype","name","Set","values","Number","String","NormalExpression","SpecialExpression","UserDefinedSymbol","NormalBuiltinSymbol","SpecialBuiltinSymbol","ReservedSymbol","Binding","Spread","isNumber","value","options","isNaN","integer","isInteger","finite","isFinite","zero","nonZero","positive","negative","nonPositive","nonNegative","gt","gte","lt","lte","annotatedArrays","WeakSet","vectors","notVectors","matrices","notMatrices","grids","toFixedArity","arity","min","getNumberVectorOrMatrixOperation","param","vector","Array","isArray","has","every","elem","add","matrix","nbrOfCols","row","some","cell","isMatrix","unaryMathOp","fn","operation","operand","map","val","mathUtilsModule","functions","sin","evaluate","docs","category","returns","type","args","x","variants","argumentNames","description","seeAlso","examples","asin","sinh","asinh","cos","acos","cosh","acosh","tan","atan","tanh","atanh","ln","log","log2","log10","PI"],"mappings":"aAEM,SAAUA,EAAcC,GAC5B,IAAKA,EAAeC,WAAaD,EAAeE,KAC9C,MAAO,GAET,MAAMC,EAAcH,EAAeC,SAASG,OAAS,EAC/CC,EAAeL,EAAeE,KAAKI,OAASH,EAAc,EAChE,MAAO,GAAG,IAAII,OAAOC,KAAKC,IAAIN,EAAa,OAAO,IAAII,OAAOC,KAAKC,IAAIJ,EAAc,KACtF,CCkBM,MAAOK,UAAmBC,MACdX,eACAY,aAChB,WAAAC,CAAYC,EAAcd,GACxB,MAAMe,EAAUD,aAAeH,MAC3BG,EAAIC,QACJ,GAAGD,IAEPE,MA/BJ,SAA8BD,EAAiBf,GAC7C,IAAKA,EACH,OAAOe,EAET,MAAME,EAAW,GAAGjB,EAAeC,SAASiB,QAAQlB,EAAeC,SAASG,SAM5E,MAAO,GAAGW,IALWf,EAAemB,SAChC,KAAKnB,EAAemB,YAAYF,IAChC,cAAcA,QACIjB,EAAeE,SACbH,EAAcC,IAExC,CAoBUoB,CAAqBL,EAASf,IACpCqB,KAAKT,aAAeG,EACpBM,KAAKrB,eAAiBA,EACtBsB,OAAOC,eAAeF,KAAMX,EAAWc,WACvCH,KAAKI,KAAO,YACb,CAEM,aAAA1B,GACL,OAAOsB,KAAKrB,gBAAkBD,EAAcsB,KAAKrB,eAClD,EC/BkB,IAAI0B,IAAIJ,OAAOK,OAbX,CACvBC,OAAQ,EACRC,OAAQ,EACRC,iBAAkB,EAClBC,kBAAmB,EACnBC,kBAAmB,EACnBC,oBAAqB,EACrBC,qBAAsB,EACtBC,eAAgB,EAChBC,QAAS,EACTC,OAAQ,eCoHMC,EAASC,EAAgBC,EAAyB,IAChE,MAAqB,iBAAVD,KAGPX,OAAOa,MAAMF,OAGbC,EAAQE,UAAYd,OAAOe,UAAUJ,QAGrCC,EAAQI,SAAWhB,OAAOiB,SAASN,QAGnCC,EAAQM,MAAkB,IAAVP,OAGhBC,EAAQO,SAAqB,IAAVR,OAGnBC,EAAQQ,UAAYT,GAAS,OAG7BC,EAAQS,UAAYV,GAAS,OAG7BC,EAAQU,aAAeX,EAAQ,OAG/BC,EAAQW,aAAeZ,EAAQ,OAGT,iBAAfC,EAAQY,IAAmBb,GAASC,EAAQY,QAG5B,iBAAhBZ,EAAQa,KAAoBd,EAAQC,EAAQa,SAG7B,iBAAfb,EAAQc,IAAmBf,GAASC,EAAQc,OAG5B,iBAAhBd,EAAQe,KAAoBhB,EAAQC,EAAQe,gBAIzD,CDhIwB,IAAI7B,IAhBN,CACpB,cACA,UACA,OACA,aACA,OACA,aACA,YACA,WACA,QACA,UACA,iBACA,mBACA,WElCF,MAAM8B,EAAkB,IAAIC,QACtBC,EAAU,IAAID,QACdE,EAAa,IAAIF,QACjBG,EAAW,IAAIH,QACfI,EAAc,IAAIJ,QAClBK,EAAQ,IAAIL,QCwEZ,SAAUM,EAAaC,GAC3B,MAAO,CAAEC,IAAKD,EAAOvD,IAAKuD,EAC5B,CC1EA,SAASE,EACPC,EACAnE,GAKA,GFUuBoE,EEVVD,EFWRE,MAAMC,QAAQF,KAIfV,EAAQa,IAAIH,KAGZT,EAAWY,IAAIH,MAIfA,EAAOI,MAAMC,GAAQnC,EAASmC,KAChCjB,EAAgBkB,IAAIN,GACpBV,EAAQgB,IAAIN,GACL,IAETT,EAAWe,IAAIN,GACR,KE3BL,MAAO,CAAC,SAAUD,GFShB,IAAmBC,EEPvB,GFyGI,SAAmBO,GACvB,IAAKN,MAAMC,QAAQK,GACjB,OAAO,EAET,GAAIf,EAASW,IAAII,GACf,OAAO,EAET,GAAId,EAAYU,IAAII,GAClB,OAAO,EAET,GAAsB,IAAlBA,EAAOrE,OAET,OADAuD,EAAYa,IAAIC,IACT,EAET,IAAKN,MAAMC,QAAQK,EAAO,KAA4B,IAArBA,EAAO,GAAGrE,OAEzC,OADAuD,EAAYa,IAAIC,IACT,EAET,MAAMC,EAAYD,EAAO,GAAGrE,OAC5B,IAAK,MAAMuE,KAAOF,EAChB,IAAKN,MAAMC,QAAQO,IAAQA,EAAIvE,SAAWsE,GAAaC,EAAIC,KAAKC,IAASzC,EAASyC,IAEhF,OADAlB,EAAYa,IAAIC,IACT,EAMX,OAHAnB,EAAgBkB,IAAIC,GACpBb,EAAMY,IAAIC,GACVf,EAASc,IAAIC,IACN,CACT,CEtIMK,CAASb,GACX,MAAO,CAAC,SAAUA,GAEpB,IAAK7B,EAAS6B,GACZ,MAAM,IAAIzD,EAAW,kCAAkCyD,EAASnE,GAElE,MAAO,CAAC,SAAUmE,EACpB,CAEA,SAASc,EACPC,GAEA,MAAO,EAAEf,GAAQnE,KACf,MAAOmF,EAAWC,GAAWlB,EAAiCC,EAAOnE,GACrE,MAAkB,WAAdmF,EACKD,EAAGE,GAEW,WAAdD,EACAC,EAAQC,IAAIC,GAAOJ,EAAGI,IAGtBF,EAAQC,IAAIR,GAAOA,EAAIQ,IAAIC,GAAOJ,EAAGI,KAGlD,CAEA,MAucaC,EAA+B,CAC1C9D,KAAM,OACN+D,UAzcmD,CACnDC,IAAO,CACLC,SAAUT,EAAYK,GAAO9E,KAAKiF,IAAIH,IACtCtB,MAAOD,EAAa,GACpB4B,KAAM,CACJC,SAAU,OACVC,QAAS,CAAEC,KAAM,CAAC,SAAU,SAAU,WACtCC,KAAM,CACJC,EAAG,CAAEF,KAAM,CAAC,SAAU,SAAU,YAElCG,SAAU,CAAC,CAAEC,cAAe,CAAC,OAC7BC,YAAa,8OACbC,QAAS,CAAC,YAAa,YAAa,WAAY,WAAY,eAC5DC,SAAU,CACR,sCAEA,sCAEA,uCAEA,yCAEA,8CAEA,wDAKNC,KAAQ,CACNZ,SAAUT,EAAYK,GAAO9E,KAAK8F,KAAKhB,IACvCtB,MAAOD,EAAa,GACpB4B,KAAM,CACJC,SAAU,OACVC,QAAS,CAAEC,KAAM,CAAC,SAAU,SAAU,WACtCC,KAAM,CACJC,EAAG,CAAEF,KAAM,CAAC,SAAU,SAAU,YAElCG,SAAU,CAAC,CAAEC,cAAe,CAAC,OAC7BC,YAAa,oQACbC,QAAS,CAAC,WAAY,aAAc,YAAa,aACjDC,SAAU,CACR,wCAEA,wCAEA,2CAEA,gDAEA,0DAKNE,KAAQ,CACNb,SAAUT,EAAYK,GAAO9E,KAAK+F,KAAKjB,IACvCtB,MAAOD,EAAa,GACpB4B,KAAM,CACJC,SAAU,OACVC,QAAS,CAAEC,KAAM,CAAC,SAAU,SAAU,WACtCC,KAAM,CACJC,EAAG,CAAEF,KAAM,CAAC,SAAU,SAAU,YAElCG,SAAU,CAAC,CAAEC,cAAe,CAAC,OAC7BC,YAAa,0PACbC,QAAS,CAAC,aAAc,WAAY,YAAa,aACjDC,SAAU,CACR,wCAEA,wCAEA,2CAEA,sDAEA,kEAKNG,MAAS,CACPd,SAAUT,EAAYK,GAAO9E,KAAKgG,MAAMlB,IACxCtB,MAAOD,EAAa,GACpB4B,KAAM,CACJC,SAAU,OACVC,QAAS,CAAEC,KAAM,CAAC,SAAU,SAAU,WACtCC,KAAM,CACJC,EAAG,CAAEF,KAAM,CAAC,SAAU,SAAU,YAElCG,SAAU,CAAC,CAAEC,cAAe,CAAC,OAC7BC,YAAa,2QACbC,QAAS,CAAC,YAAa,YAAa,aAAc,cAClDC,SAAU,CACR,2CAEA,2CAEA,2CAEA,qDAEA,gEAKNI,IAAO,CACLf,SAAUT,EAAYK,GAAO9E,KAAKiG,IAAInB,IACtCtB,MAAOD,EAAa,GACpB4B,KAAM,CACJC,SAAU,OACVC,QAAS,CAAEC,KAAM,CAAC,SAAU,SAAU,WACtCC,KAAM,CACJC,EAAG,CAAEF,KAAM,CAAC,SAAU,SAAU,YAElCG,SAAU,CAAC,CAAEC,cAAe,CAAC,OAC7BC,YAAa,kPACbC,QAAS,CAAC,YAAa,YAAa,WAAY,WAAY,eAC5DC,SAAU,CACR,sCAEA,sCAEA,uCAEA,yCAEA,8CAEA,wDAKNK,KAAQ,CACNhB,SAAUT,EAAYK,GAAO9E,KAAKkG,KAAKpB,IACvCtB,MAAOD,EAAa,GACpB4B,KAAM,CACJC,SAAU,OACVC,QAAS,CAAEC,KAAM,CAAC,SAAU,SAAU,WACtCC,KAAM,CACJC,EAAG,CAAEF,KAAM,CAAC,SAAU,SAAU,YAElCG,SAAU,CAAC,CAAEC,cAAe,CAAC,OAC7BC,YAAa,0QACbC,QAAS,CAAC,WAAY,aAAc,YAAa,aACjDC,SAAU,CACR,wCAEA,wCAEA,2CAEA,sDAEA,kEAKNM,KAAQ,CACNjB,SAAUT,EAAYK,GAAO9E,KAAKmG,KAAKrB,IACvCtB,MAAOD,EAAa,GACpB4B,KAAM,CACJC,SAAU,OACVC,QAAS,CAAEC,KAAM,CAAC,SAAU,SAAU,WACtCC,KAAM,CACJC,EAAG,CAAEF,KAAM,CAAC,SAAU,SAAU,YAElCG,SAAU,CAAC,CAAEC,cAAe,CAAC,OAC7BC,YAAa,8PACbC,QAAS,CAAC,aAAc,WAAY,YAAa,aACjDC,SAAU,CACR,wCAEA,wCAEA,2CAEA,sDAEA,kEAKNO,MAAS,CACPlB,SAAUT,EAAYK,GAAO9E,KAAKoG,MAAMtB,IACxCtB,MAAOD,EAAa,GACpB4B,KAAM,CACJC,SAAU,OACVC,QAAS,CAAEC,KAAM,CAAC,SAAU,SAAU,WACtCC,KAAM,CACJC,EAAG,CAAEF,KAAM,CAAC,SAAU,SAAU,YAElCG,SAAU,CAAC,CAAEC,cAAe,CAAC,OAC7BC,YAAa,+QACbC,QAAS,CAAC,YAAa,YAAa,aAAc,cAClDC,SAAU,CACR,0CAEA,0CAEA,4CAEA,2CAEA,kDAEA,4DAKNQ,IAAO,CACLnB,SAAUT,EAAYK,GAAO9E,KAAKqG,IAAIvB,IACtCtB,MAAOD,EAAa,GACpB4B,KAAM,CACJC,SAAU,OACVC,QAAS,CAAEC,KAAM,CAAC,SAAU,SAAU,WACtCC,KAAM,CACJC,EAAG,CAAEF,KAAM,CAAC,SAAU,SAAU,YAElCG,SAAU,CAAC,CAAEC,cAAe,CAAC,OAC7BC,YAAa,oPACbC,QAAS,CAAC,YAAa,YAAa,WAAY,WAAY,eAC5DC,SAAU,CACR,sCAEA,sCAEA,uCAEA,yCAEA,8CAEA,wDAKNS,KAAQ,CACNpB,SAAUT,EAAYK,GAAO9E,KAAKsG,KAAKxB,IACvCtB,MAAOD,EAAa,GACpB4B,KAAM,CACJC,SAAU,OACVC,QAAS,CAAEC,KAAM,CAAC,SAAU,SAAU,WACtCC,KAAM,CACJC,EAAG,CAAEF,KAAM,CAAC,SAAU,SAAU,YAElCG,SAAU,CAAC,CAAEC,cAAe,CAAC,OAC7BC,YAAa,6QACbC,QAAS,CAAC,WAAY,aAAc,YAAa,aACjDC,SAAU,CACR,wCAEA,wCAEA,2CAEA,sDAEA,kEAKNU,KAAQ,CACNrB,SAAUT,EAAYK,GAAO9E,KAAKuG,KAAKzB,IACvCtB,MAAOD,EAAa,GACpB4B,KAAM,CACJC,SAAU,OACVC,QAAS,CAAEC,KAAM,CAAC,SAAU,SAAU,WACtCC,KAAM,CACJC,EAAG,CAAEF,KAAM,CAAC,SAAU,SAAU,YAElCG,SAAU,CAAC,CAAEC,cAAe,CAAC,OAC7BC,YAAa,gQACbC,QAAS,CAAC,aAAc,WAAY,YAAa,aACjDC,SAAU,CACR,wCAEA,wCAEA,2CAEA,4CAKNW,MAAS,CACPtB,SAAUT,EAAYK,GAAO9E,KAAKwG,MAAM1B,IACxCtB,MAAOD,EAAa,GACpB4B,KAAM,CACJC,SAAU,OACVC,QAAS,CAAEC,KAAM,CAAC,SAAU,SAAU,WACtCC,KAAM,CACJC,EAAG,CAAEF,KAAM,CAAC,SAAU,SAAU,YAElCG,SAAU,CAAC,CAAEC,cAAe,CAAC,OAC7BC,YAAa,iRACbC,QAAS,CAAC,YAAa,YAAa,aAAc,cAClDC,SAAU,CACR,0CAEA,4CAEA,6CAEA,wDAEA,oEAKNY,GAAM,CACJvB,SAAUT,EAAYK,GAAO9E,KAAK0G,IAAI5B,IACtCtB,MAAOD,EAAa,GACpB4B,KAAM,CACJC,SAAU,OACVC,QAAS,CAAEC,KAAM,CAAC,SAAU,SAAU,WACtCC,KAAM,CACJC,EAAG,CAAEF,KAAM,CAAC,SAAU,SAAU,YAElCG,SAAU,CAAC,CAAEC,cAAe,CAAC,OAC7BC,YAAa,uQACbC,QAAS,CAAC,YAAa,aAAc,KACrCC,SAAU,CACR,uCAEA,sCAEA,oCAEA,4CAEA,sDAKNc,KAAQ,CACNzB,SAAUT,EAAYK,GAAO9E,KAAK2G,KAAK7B,IACvCtB,MAAOD,EAAa,GACpB4B,KAAM,CACJC,SAAU,OACVC,QAAS,CAAEC,KAAM,CAAC,SAAU,SAAU,WACtCC,KAAM,CACJC,EAAG,CAAEF,KAAM,CAAC,SAAU,SAAU,YAElCG,SAAU,CAAC,CAAEC,cAAe,CAAC,OAC7BC,YAAa,8PACbC,QAAS,CAAC,UAAW,cACrBC,SAAU,CACR,2CAEA,6CAEA,0CAEA,gDAEA,0DAKNe,MAAS,CACP1B,SAAUT,EAAYK,GAAO9E,KAAK4G,MAAM9B,IACxCtB,MAAOD,EAAa,GACpB4B,KAAM,CACJC,SAAU,OACVC,QAAS,CAAEC,KAAM,CAAC,SAAU,SAAU,WACtCC,KAAM,CACJC,EAAG,CAAEF,KAAM,CAAC,SAAU,SAAU,YAElCG,SAAU,CAAC,CAAEC,cAAe,CAAC,OAC7BC,YAAa,iQACbC,QAAS,CAAC,UAAW,aACrBC,SAAU,CACR,6CAEA,gDAEA,4CAEA,kDAEA,4DAKN,SAAU,CACRX,SAAUT,EAAYK,GAAQA,EAAM9E,KAAK6G,GAAM,KAC/CrD,MAAOD,EAAa,GACpB4B,KAAM,CACJC,SAAU,OACVC,QAAS,CAAEC,KAAM,CAAC,SAAU,SAAU,WACtCC,KAAM,CACJC,EAAG,CAAEF,KAAM,CAAC,SAAU,SAAU,YAElCG,SAAU,CAAC,CAAEC,cAAe,CAAC,OAC7BC,YAAa,qOACbC,QAAS,CAAC,cAAe,WAAY,WAAY,YACjDC,SAAU,CACR,4CAEA,6CAEA,8CAEA,8CAEA,uDAEA,mEAKN,SAAU,CACRX,SAAUT,EAAYK,GAAc,IAANA,EAAa9E,KAAK6G,IAChDrD,MAAOD,EAAa,GACpB4B,KAAM,CACJC,SAAU,OACVC,QAAS,CAAEC,KAAM,CAAC,SAAU,SAAU,WACtCC,KAAM,CACJC,EAAG,CAAEF,KAAM,CAAC,SAAU,SAAU,YAElCG,SAAU,CAAC,CAAEC,cAAe,CAAC,OAC7BC,YAAa,qOACbC,QAAS,CAAC,eACVC,SAAU,CACR,4CAEA,6CAEA,iDAEA,qDAEA,0DAEA"}

package/dist/modules/matrix.esm.js

@@ -0,0 +1,2 @@
+function t(t){if(!t.position||!t.code)return"";const e=t.position.column-1,r=t.code.length-e-1;return`${" ".repeat(Math.max(e,0))}^${" ".repeat(Math.max(r,0))}`}class e extends Error{sourceCodeInfo;shortMessage;constructor(r,a){const n=r instanceof Error?r.message:`${r}`;super(function(e,r){if(!r)return e;const a=`${r.position.line}:${r.position.column}`;return`${e}${r.filePath?`\n${r.filePath}:${a}`:`\nLocation ${a}`}\n${r.code}\n${t(r)}`}(n,a)),this.shortMessage=n,this.sourceCodeInfo=a,Object.setPrototypeOf(this,e.prototype),this.name="DvalaError"}getCodeMarker(){return this.sourceCodeInfo&&t(this.sourceCodeInfo)}}const r={Number:1,String:2,NormalExpression:3,SpecialExpression:4,UserDefinedSymbol:5,NormalBuiltinSymbol:6,SpecialBuiltinSymbol:7,ReservedSymbol:8,Binding:9,Spread:10},a=new Set(Object.values(r));const n=new Set(["UserDefined","Partial","Comp","Constantly","Juxt","Complement","EveryPred","SomePred","Fnull","Builtin","SpecialBuiltin","NativeJsFunction","Module"]);function i(t){return null!==t&&"object"==typeof t&&("^^fn^^"in t&&"functionType"in t&&("string"==typeof(e=t.functionType)&&n.has(e)));var e}function o(t){return!(!Array.isArray(t)||t.length<2)&&("number"==typeof(e=t[0])&&a.has(e));var e}function m(t){return i(t)?`<function ${t.name||"λ"}>`:o(t)?`${e=t[0],Object.keys(r).find(t=>r[t]===e)}-node`:null===t?"null":"object"==typeof t&&t instanceof RegExp?`${t}`:"object"==typeof t&&t instanceof Error?t.toString():JSON.stringify(t);var e}function l(t,e={}){return"number"==typeof t&&(!Number.isNaN(t)&&(!(e.integer&&!Number.isInteger(t))&&(!(e.finite&&!Number.isFinite(t))&&((!e.zero||0===t)&&((!e.nonZero||0!==t)&&(!(e.positive&&t<=0)&&(!(e.negative&&t>=0)&&(!(e.nonPositive&&t>0)&&(!(e.nonNegative&&t<0)&&(!("number"==typeof e.gt&&t<=e.gt)&&(!("number"==typeof e.gte&&t<e.gte)&&(!("number"==typeof e.lt&&t>=e.lt)&&!("number"==typeof e.lte&&t>e.lte)))))))))))))}function s(t,r,a={}){if(!l(t,a))throw new e(`Expected ${function(t){if(t.zero)return"zero";const e=function(t){return t.positive?"positive":t.negative?"negative":t.nonNegative?"non negative":t.nonPositive?"non positive":t.nonZero?"non zero":""}(t),r=t.integer?"integer":"number",a=t.finite?"finite":"",n=function(t){return"number"!=typeof t.gt&&"number"!=typeof t.gte||"number"!=typeof t.lt&&"number"!=typeof t.lte?"number"==typeof t.gt||"number"==typeof t.gte?"number"==typeof t.gt?`n > ${t.gt}`:`n >= ${t.gte}`:"number"==typeof t.lt||"number"==typeof t.lte?"number"==typeof t.lt?`n < ${t.lt}`:`n <= ${t.lte}`:"":`${"number"==typeof t.gt?`${t.gt} < n `:`${t.gte} <= n `}${"number"==typeof t.lt?`< ${t.lt}`:`<= ${t.lte}`}`}(t);return[e,a,r,n].filter(t=>!!t).join(" ")}(a)}, got ${m(t)}.`,function(t,e){return t?.sourceCodeInfo??e}(t,r))}const u=new WeakSet,c=new WeakSet,x=new WeakSet,p=new WeakSet,f=new WeakSet,d=new WeakSet;function g(t,r){if(!function(t){return!(!Array.isArray(t)||!c.has(t)&&(x.has(t)||(t.every(t=>l(t))?(u.add(t),c.add(t),0):(x.add(t),1))))}(t))throw new e(`Expected a vector, but got ${t}`,r)}function h(t){if(!Array.isArray(t))return!1;if(p.has(t))return!0;if(f.has(t))return!1;if(0===t.length)return f.add(t),!1;if(!Array.isArray(t[0])||0===t[0].length)return f.add(t),!1;const e=t[0].length;for(const r of t)if(!Array.isArray(r)||r.length!==e||r.some(t=>!l(t)))return f.add(t),!1;return u.add(t),d.add(t),p.add(t),!0}function y(t,r){if(!h(t))throw new e(`Expected a matrix, but got ${t}`,r)}function b(t,r){if(!h(t))throw new e(`Expected a matrix, but got ${t}`,r);if(t.length!==t[0].length)throw new e(`Expected square matrix, but got ${t.length} and ${t[0].length}`,r)}function v(t,e,r=1e-10){if(t===e)return!0;const a=Math.abs(t-e);if(0===t||0===e||a<r)return a<r;return a/(Math.abs(t)+Math.abs(e))<r}function w(t){return Math.abs(t)<1e-10}function N(t){return{min:t,max:t}}const k={mul:{category:"matrix",description:"Multiplies two `matrices` using standard `matrix` multiplication based on **dot products** of rows and columns.",returns:{type:"matrix"},args:{a:{type:"matrix"},b:{type:"matrix"}},variants:[{argumentNames:["a","b"]}],examples:["let { mul } = import(matrix);\nmul([[1, 2], [3, 4]], [[5, 6], [7, 8]])","let { mul } = import(matrix);\nmul([[1, 2, 3], [4, 5, 6]], [[7, 8], [9, 10], [11, 12]])"],seeAlso:["matrix.det","matrix.inv"]},det:{category:"matrix",description:"Calculates the **determinant** of a square matrix.",returns:{type:"number"},args:{m:{type:"matrix",description:"The `matrix` to calculate the determinant of."}},variants:[{argumentNames:["m"]}],examples:["let { det } = import(matrix);\ndet([[1, 2], [3, 4]])","let { det } = import(matrix);\ndet([[1, 2, 3], [4, 5, 6], [7, 8, 9]])"],seeAlso:["matrix.inv","matrix.cofactor","matrix.adj","matrix.trace","matrix.rank","matrix.invertible?","matrix.mul","matrix.minor"]},inv:{category:"matrix",description:"Calculates the **inverse** of a square matrix.",returns:{type:"matrix"},args:{m:{type:"matrix",description:"The `matrix` to calculate the inverse of."}},variants:[{argumentNames:["m"]}],examples:["let { inv } = import(matrix);\ninv([[1, 2], [3, 4]])","let { inv } = import(matrix);\ninv([[1, 2, 3], [4, 5, 7], [7, 8, 10]])"],seeAlso:["matrix.det","matrix.adj","matrix.invertible?","linear-algebra.solve","matrix.mul","matrix.orthogonal-matrix?"]},adj:{category:"matrix",description:"Calculates the **adjugate** of a square matrix.",returns:{type:"matrix"},args:{m:{type:"matrix",description:"The `matrix` to calculate the adjugate of."}},variants:[{argumentNames:["m"]}],examples:["let { adj } = import(matrix);\nadj([[1, 2], [3, 4]])","let { adj } = import(matrix);\nadj([[1, 2, 3], [4, 5, 6], [7, 8, 9]])","let { adj } = import(matrix);\nadj([[1, 2, 3], [7, 8, 9], [4, 5, 6]])"],seeAlso:["matrix.cofactor","matrix.det","matrix.inv"]},cofactor:{category:"matrix",description:"Calculates the **cofactor** of a square matrix.",returns:{type:"matrix"},args:{m:{type:"matrix",description:"The `matrix` to calculate the cofactor of."}},variants:[{argumentNames:["m"]}],examples:["let { cofactor } = import(matrix);\ncofactor([[1, 2], [3, 4]])","let { cofactor } = import(matrix);\ncofactor([[1, 2, 3], [4, 5, 6], [7, 8, 9]])","let { cofactor } = import(matrix);\ncofactor([[1, 2, 3], [7, 8, 9], [4, 5, 6]])"],seeAlso:["matrix.adj","matrix.minor","matrix.det"]},minor:{category:"matrix",description:"Calculates the **minor** of a square matrix.",returns:{type:"matrix"},args:{m:{type:"matrix",description:"The `matrix` to calculate the minor of."},row:{type:"integer",description:"The row index of the element to calculate the minor for."},col:{type:"integer",description:"The column index of the element to calculate the minor for."}},variants:[{argumentNames:["m","row","col"]}],examples:["let { minor } = import(matrix);\nminor([[1, 2], [3, 4]], 0, 1)","let { minor } = import(matrix);\nminor([[1, 2, 3], [4, 5, 6], [7, 8, 9]], 1, 1)"],seeAlso:["matrix.cofactor","matrix.det"]},trace:{category:"matrix",description:"Calculates the **trace** of a square matrix.",returns:{type:"number"},args:{m:{type:"matrix",description:"The `matrix` to calculate the trace of."}},variants:[{argumentNames:["m"]}],examples:["let { trace } = import(matrix);\ntrace([[1, 2], [3, 4]])","let { trace } = import(matrix);\ntrace([[1, 2, 3], [4, 5, 6], [7, 8, 9]])"],seeAlso:["matrix.det","matrix.diagonal?"]},"symmetric?":{category:"matrix",description:"Checks if a `matrix` is **symmetric**.",returns:{type:"boolean"},args:{m:{type:"matrix",description:"The `matrix` to check for symmetry."}},variants:[{argumentNames:["m"]}],examples:["let { symmetric? } = import(matrix);\nsymmetric?([[1, 2], [2, 1]])","let { symmetric? } = import(matrix);\nsymmetric?([[1, 2, 3], [2, 1, 4], [3, 4, 1]])"],seeAlso:["matrix.orthogonal-matrix?","matrix.diagonal?","matrix.square?","matrix.hilbert"]},"triangular?":{category:"matrix",description:"Checks if a `matrix` is **triangular**.",returns:{type:"boolean"},args:{m:{type:"matrix",description:"The `matrix` to check for triangularity."}},variants:[{argumentNames:["m"]}],examples:["let { triangular? } = import(matrix);\ntriangular?([[2, 0], [0, 1]])","let { triangular? } = import(matrix);\ntriangular?([[1, 2, 3], [0, 4, 5], [0, 0, 6]])"],seeAlso:["matrix.upper-triangular?","matrix.lower-triangular?","matrix.diagonal?","matrix.banded?"]},"upper-triangular?":{category:"matrix",description:"Checks if a `matrix` is **upper triangular**.",returns:{type:"boolean"},args:{m:{type:"matrix",description:"The `matrix` to check for upper triangularity."}},variants:[{argumentNames:["m"]}],examples:["let { upper-triangular? } = import(matrix);\nupper-triangular?([[1, 2], [0, 3]])","let { upper-triangular? } = import(matrix);\nupper-triangular?([[1, 2, 3], [0, 4, 5], [0, 0, 6]])"],seeAlso:["matrix.lower-triangular?","matrix.triangular?","matrix.diagonal?"]},"lower-triangular?":{category:"matrix",description:"Checks if a `matrix` is **lower triangular**.",returns:{type:"boolean"},args:{m:{type:"matrix",description:"The `matrix` to check for lower triangularity."}},variants:[{argumentNames:["m"]}],examples:["let { lower-triangular? } = import(matrix);\nlower-triangular?([[1, 0], [2, 3]])","let { lower-triangular? } = import(matrix);\nlower-triangular?([[1, 0, 0], [2, 3, 0], [4, 5, 6]])"],seeAlso:["matrix.upper-triangular?","matrix.triangular?","matrix.diagonal?"]},"diagonal?":{category:"matrix",description:"Checks if a `matrix` is **diagonal**.",returns:{type:"boolean"},args:{m:{type:"matrix",description:"The `matrix` to check for diagonal property."}},variants:[{argumentNames:["m"]}],examples:["let { diagonal? } = import(matrix);\ndiagonal?([[1, 0], [0, 2]])","let { diagonal? } = import(matrix);\ndiagonal?([[1, 0, 0], [0, 2, 0], [0, 0, 3]])","let { diagonal? } = import(matrix);\ndiagonal?([[1, 0, 0], [2, 2, 2], [0, 0, 3]])"],seeAlso:["matrix.identity?","matrix.symmetric?","matrix.triangular?","matrix.trace","matrix.upper-triangular?","matrix.lower-triangular?","matrix.band","matrix.banded?"]},"square?":{category:"matrix",description:"Checks if a `matrix` is **square**.",returns:{type:"boolean"},args:{m:{type:"matrix",description:"The `matrix` to check for square property."}},variants:[{argumentNames:["m"]}],examples:["let { square? } = import(matrix);\nsquare?([[1, 2], [3, 4]])","let { square? } = import(matrix);\nsquare?([[1, 2, 3], [4, 5, 6], [7, 8, 9]])","let { square? } = import(matrix);\nsquare?([[1, 2, 3], [4, 5, 6]])"],seeAlso:["matrix.symmetric?","matrix.identity?","matrix.invertible?"]},"orthogonal-matrix?":{category:"matrix",description:"Checks if a `matrix` is **orthogonal**.",returns:{type:"boolean"},args:{m:{type:"matrix",description:"The `matrix` to check for **orthogonality**."}},variants:[{argumentNames:["m"]}],examples:["let { orthogonal-matrix? } = import(matrix);\northogonal-matrix?([[1, 0], [0, 1]])","let { orthogonal-matrix? } = import(matrix);\northogonal-matrix?([[1, 0], [0, -1]])","let { orthogonal-matrix? } = import(matrix);\northogonal-matrix?([[1, 2], [3, 4]])"],seeAlso:["matrix.symmetric?","matrix.inv","matrix.identity?","linear-algebra.orthogonal?"]},"identity?":{category:"matrix",description:"Checks if a `matrix` is an **identity matrix**.",returns:{type:"boolean"},args:{m:{type:"matrix",description:"The `matrix` to check for identity property."}},variants:[{argumentNames:["m"]}],examples:["let { identity? } = import(matrix);\nidentity?([[1, 0], [0, 1]])","let { identity? } = import(matrix);\nidentity?([[1, 0, 0], [0, 1, 0], [0, 0, 1]])","let { identity? } = import(matrix);\nidentity?([[1, 0, 0], [0, 1, 0], [0, 0, 0]])"],seeAlso:["matrix.diagonal?","matrix.square?","matrix.orthogonal-matrix?"]},"invertible?":{category:"matrix",description:"Checks if a `matrix` is **invertible**.",returns:{type:"boolean"},args:{m:{type:"matrix",description:"The `matrix` to check for invertibility."}},variants:[{argumentNames:["m"]}],examples:["let { invertible? } = import(matrix);\ninvertible?([[1, 2], [3, 4]])","let { invertible? } = import(matrix);\ninvertible?([[1, 2, 3], [4, 5, 6], [7, 8, 9]])","let { invertible? } = import(matrix);\ninvertible?([[1, 2], [2, 4]])"],seeAlso:["matrix.det","matrix.inv","matrix.rank","matrix.square?"]},hilbert:{category:"matrix",description:"Generates a **Hilbert matrix** of size `n`.",returns:{type:"matrix"},args:{n:{type:"integer",description:"The size of the Hilbert matrix."}},variants:[{argumentNames:["n"]}],examples:["let { hilbert } = import(matrix);\nhilbert(3)","let { hilbert } = import(matrix);\nhilbert(4)"],seeAlso:["matrix.vandermonde","matrix.symmetric?"]},vandermonde:{category:"matrix",description:"Generates a **Vandermonde matrix** from a vector.",returns:{type:"matrix"},args:{v:{type:"vector",description:"The vector to generate the Vandermonde matrix from."}},variants:[{argumentNames:["v"]}],examples:["let { vandermonde } = import(matrix);\nvandermonde([1, 2, 3])","let { vandermonde } = import(matrix);\nvandermonde([1, 0, 1])"],seeAlso:["matrix.hilbert","matrix.band"]},band:{category:"matrix",description:"Generates a **banded matrix** of size `n` with lower band index `lband` and upper band index `uband`.",returns:{type:"matrix"},args:{n:{type:"integer",description:"The size of the banded matrix."},lband:{type:"integer",description:"The lower band index."},uband:{type:"integer",description:"The upper band index."}},variants:[{argumentNames:["n","lband","uband"]}],examples:["let { band } = import(matrix);\nband(3, 1, 1)","let { band } = import(matrix);\nband(4, 1, 2)"],seeAlso:["matrix.banded?","matrix.diagonal?","matrix.vandermonde"]},"banded?":{category:"matrix",description:"Checks if a `matrix` is **banded** with lower band index `lband` and upper band index `uband`.",returns:{type:"boolean"},args:{m:{type:"matrix",description:"The `matrix` to check for **banded** property."},lband:{type:"integer",description:"The lower band index."},uband:{type:"integer",description:"The upper band index."}},variants:[{argumentNames:["m","lband","uband"]}],examples:["let { banded? } = import(matrix);\nbanded?([\n  [1, 1, 1, 0],\n  [1, 1, 1, 1],\n  [1, 1, 1, 1],\n  [0, 1, 1, 1],\n], 2, 2)","let { banded? } = import(matrix);\nbanded?([\n  [1, 1, 1, 0],\n  [1, 1, 1, 1],\n  [1, 1, 1, 1],\n  [0, 1, 1, 1],\n], 1, 1)"],seeAlso:["matrix.band","matrix.triangular?","matrix.diagonal?"]},rank:{category:"matrix",description:"Calculates the **rank** of a matrix using **Gaussian elimination**.",returns:{type:"number"},args:{m:{type:"matrix",description:"The `matrix` to calculate the rank of."}},variants:[{argumentNames:["m"]}],examples:["let { rank } = import(matrix);\nrank([[1, 0, 0], [0, 1, 0], [0, 0, 1]])","let { rank } = import(matrix);\nrank([[1, 2, 3], [4, 5, 6], [7, 8, 9]])","let { rank } = import(matrix);\nrank([[2, 4, 6], [3, 6, 9], [4, 8, 12]])"],seeAlso:["matrix.det","matrix.invertible?","linear-algebra.rref"]},"frobenius-norm":{category:"matrix",description:"Calculates the **Frobenius norm** of a matrix.",returns:{type:"number"},args:{m:{type:"matrix",description:"The `matrix` to calculate the Frobenius norm of."}},variants:[{argumentNames:["m"]}],examples:["let { frobenius-norm } = import(matrix);\nfrobenius-norm([[1, 2], [3, 4]])","let { frobenius-norm } = import(matrix);\nfrobenius-norm([[1, 2, 3], [4, 5, 6], [7, 8, 9]])"],seeAlso:["matrix.one-norm","matrix.inf-norm","matrix.max-norm"]},"one-norm":{category:"matrix",description:"Calculates the **one-norm** (column norm) of a matrix.",returns:{type:"number"},args:{m:{type:"matrix",description:"The `matrix` to calculate the one-norm of."}},variants:[{argumentNames:["m"]}],examples:["let { one-norm } = import(matrix);\none-norm([[1, 2], [3, 4]])","let { one-norm } = import(matrix);\none-norm([[1, 2, 3], [4, 5, 6], [7, 8, 9]])"],seeAlso:["matrix.frobenius-norm","matrix.inf-norm","matrix.max-norm"]},"inf-norm":{category:"matrix",description:"Calculates the **infinity norm** of a matrix.",returns:{type:"number"},args:{m:{type:"matrix",description:"The `matrix` to calculate the infinity norm of."}},variants:[{argumentNames:["m"]}],examples:["let { inf-norm } = import(matrix);\ninf-norm([[1, 2], [3, 4]])","let { inf-norm } = import(matrix);\ninf-norm([[1, 2, 3], [4, 5, 6], [7, 8, 9]])"],seeAlso:["matrix.frobenius-norm","matrix.one-norm","matrix.max-norm"]},"max-norm":{category:"matrix",description:"Calculates the **max norm** of a matrix.",returns:{type:"number"},args:{m:{type:"matrix",description:"The `matrix` to calculate the max norm of."}},variants:[{argumentNames:["m"]}],examples:["let { max-norm } = import(matrix);\nmax-norm([[1, 2], [3, 4]])","let { max-norm } = import(matrix);\nmax-norm([[1, 2, 3], [4, 5, 6], [7, 8, 9]])"],seeAlso:["matrix.frobenius-norm","matrix.one-norm","matrix.inf-norm"]}};function A(t){const e=t.length,r=[];for(let a=0;a<e;a++)r[a]=[...t[a]];if(1===e)return r[0][0];if(2===e)return r[0][0]*r[1][1]-r[0][1]*r[1][0];let a=1;for(let t=0;t<e-1;t+=1){let n=t;for(let a=t+1;a<e;a+=1)Math.abs(r[a][t])>Math.abs(r[n][t])&&(n=a);if(w(r[n][t]))return 0;n!==t&&([r[t],r[n]]=[r[n],r[t]],a=-a);for(let a=t+1;a<e;a+=1){const n=r[a][t]/r[t][t];for(let i=t;i<e;i++)r[a][i]-=n*r[t][i]}}let n=a;for(let t=0;t<e;t++)n*=r[t][t];return n}function T(t,e,r){const a=t.length,n=[];for(let i=0;i<a;i++)if(i!==e){const e=[];for(let n=0;n<a;n++)n!==r&&e.push(t[i][n]);n.push(e)}return n}function $(t){const e=t.length,r=[];for(let a=0;a<e;a++){r[a]=[];for(let n=0;n<e;n++){const e=(-1)**(a+n)*A(T(t,n,a));r[a][n]=e}}return r}function C(t){return t.length===t[0].length}function M(t){if(!C(t))return!1;const e=t.length;for(let r=0;r<e;r++)for(let a=0;a<e;a++)if(r===a){if(!v(t[r][a],1))return!1}else if(!w(t[r][a]))return!1;return!0}function q(t,e){if(0===t.length||0===e.length||t[0].length!==e.length)throw new Error("Matrix dimensions do not match for multiplication");const r=t.length,a=t[0].length,n=e[0].length,i=Array(r).fill(0).map(()=>Array(n).fill(0));for(let o=0;o<r;o++)for(let r=0;r<a;r++){const a=t[o][r];for(let t=0;t<n;t++)i[o][t]+=a*e[r][t]}return i}const j={mul:{evaluate:([t,r],a)=>{y(t,a),y(r,a);try{return q(t,r)}catch(n){throw new e(`The number of columns in the first matrix must be equal to the number of rows in the second matrix, but got ${t[0].length} and ${r.length}`,a)}},arity:N(2)},det:{evaluate:([t],e)=>(b(t,e),A(t)),arity:N(1)},inv:{evaluate:([t],r)=>{b(t,r);const a=function(t){const e=t.length;if(1===e){const e=t[0][0];return w(e)?null:[[1/e]]}const r=A(t);if(w(r))return null;const a=$(t),n=[];for(let t=0;t<e;t++){n[t]=[];for(let i=0;i<e;i++)n[t][i]=a[t][i]/r}return n}(t);if(null===a)throw new e("The matrix must be invertible",r);return a},arity:N(1)},adj:{evaluate:([t],e)=>(b(t,e),$(t)),arity:N(1)},cofactor:{evaluate:([t],e)=>(b(t,e),function(t){const e=t.length,r=[];for(let a=0;a<e;a++){r[a]=[];for(let n=0;n<e;n++){const e=T(t,a,n),i=(-1)**(a+n);r[a][n]=i*A(e)}}return r}(t)),arity:N(1)},minor:{evaluate:([t,e,r],a)=>(y(t,a),s(e,a,{integer:!0,nonNegative:!0,lte:t.length}),s(r,a,{integer:!0,nonNegative:!0,lte:t[0].length}),T(t,e,r)),arity:N(3)},trace:{evaluate:([t],e)=>(b(t,e),function(t){return t.reduce((t,e,r)=>t+e[r],0)}(t)),arity:N(1)},"symmetric?":{evaluate:([t],e)=>(y(t,e),function(t){const e=t.length;if(!C(t))return!1;for(let r=0;r<e;r+=1)for(let e=0;e<r;e+=1)if(t[r][e]!==t[e][r])return!1;return!0}(t)),arity:N(1)},"triangular?":{evaluate:([t],e)=>(y(t,e),function(t){if(!C(t))return!1;const e=t.length;let r=!0,a=!0;for(let n=0;n<e;n++)for(let i=0;i<e;i++){if(n>i&&0!==t[n][i]&&(r=!1,!a))return!1;n<i&&0!==t[n][i]&&(a=!1)}return r||a}(t)),arity:N(1)},"upper-triangular?":{evaluate:([t],e)=>(y(t,e),function(t){if(!C(t))return!1;const e=t.length;for(let r=0;r<e;r++)for(let e=0;e<r;e++)if(0!==t[r][e])return!1;return!0}(t)),arity:N(1)},"lower-triangular?":{evaluate:([t],e)=>(y(t,e),function(t){if(!C(t))return!1;const e=t.length;if(!t.every(t=>t.length===e))return!1;for(let r=0;r<e;r++)for(let a=r+1;a<e;a++)if(0!==t[r][a])return!1;return!0}(t)),arity:N(1)},"diagonal?":{evaluate:([t],e)=>(y(t,e),function(t){if(!C(t))return!1;const e=t.length;for(let r=0;r<e;r+=1)for(let a=0;a<e;a+=1)if(r!==a&&0!==t[r][a])return!1;return!0}(t)),arity:N(1)},"square?":{evaluate:([t],e)=>(y(t,e),C(t)),arity:N(1)},"orthogonal-matrix?":{evaluate:([t],e)=>(y(t,e),function(t){return!!C(t)&&M(q(t,function(t){const e=[];for(let r=0;r<t[0].length;r+=1){const a=[];for(let e=0;e<t.length;e+=1)a.push(t[e][r]);e.push(a)}return e}(t)))}(t)),arity:N(1)},"identity?":{evaluate:([t],e)=>(y(t,e),M(t)),arity:N(1)},"invertible?":{evaluate:([t],e)=>(y(t,e),!!function(t){return!!h(t)&&t.length===t[0].length}(t)&&!w(A(t))),arity:N(1)},hilbert:{evaluate:([t],e)=>{s(t,e,{integer:!0,positive:!0});const r=[];for(let e=0;e<t;e+=1){const a=[];for(let r=0;r<t;r+=1)a.push(1/(e+r+1));r.push(a)}return r},arity:N(1)},vandermonde:{evaluate:([t],e)=>{g(t,e);const r=[];for(let e=0;e<t.length;e+=1){const a=[];for(let r=0;r<t.length;r+=1)a.push(t[e]**r);r.push(a)}return r},arity:N(1)},band:{evaluate:([t,e,r],a)=>(s(t,a,{integer:!0,positive:!0}),s(e,a,{integer:!0,nonNegative:!0,lt:t}),s(r,a,{integer:!0,nonNegative:!0,lte:t}),function(t,e,r){const a=Array.from({length:t},()=>Array.from({length:t},()=>0));for(let n=0;n<t;n++)for(let i=Math.max(0,n-e);i<=Math.min(t-1,n+r);i++)a[n][i]=1;return a}(t,e,r)),arity:N(3)},"banded?":{evaluate:([t,e,r],a)=>{y(t,a);const n=Math.max(t.length,t[0].length);return s(e,a,{integer:!0,nonNegative:!0,lt:n}),s(r,a,{integer:!0,nonNegative:!0,lt:n}),function(t,e,r){const a=t.length,n=t[0].length;for(let i=0;i<a;i++)for(let a=0;a<n;a++)if(0!==t[i][a]&&(i-a>e||a-i>r))return!1;return!0}(t,e,r)},arity:N(3)},rank:{evaluate:([t],e)=>{y(t,e);const[,r]=function(t){const e=t.map(t=>[...t]),r=e.length,a=e[0].length;let n=0,i=0;for(let t=0;t<a;t++){let o=-1;for(let a=i;a<r;a++)if(!w(e[a][t])){o=a;break}if(-1===o)continue;n+=1,o!==i&&([e[o],e[i]]=[e[i],e[o]]);const m=e[i][t];for(let r=t;r<a;r++)e[i][r]/=m;for(let n=0;n<r;n++)if(n!==i&&!w(e[n][t])){const r=e[n][t];for(let o=t;o<a;o++)e[n][o]-=r*e[i][o]}if(i++,i===r)break}return[e,n]}(t);return r},arity:N(1)},"frobenius-norm":{evaluate:([t],e)=>(y(t,e),Math.sqrt(t.reduce((t,e)=>t+e.reduce((t,e)=>t+e*e,0),0))),arity:N(1)},"one-norm":{evaluate:([t],e)=>(y(t,e),function(t){const e=t.length,r=t[0].length;let a=0;for(let n=0;n<r;n+=1){let r=0;for(let a=0;a<e;a+=1)r+=Math.abs(t[a][n]);a=Math.max(a,r)}return a}(t)),arity:N(1)},"inf-norm":{evaluate:([t],e)=>(y(t,e),t.reduce((t,e)=>Math.max(t,e.reduce((t,e)=>t+Math.abs(e),0)),0)),arity:N(1)},"max-norm":{evaluate:([t],e)=>(y(t,e),t.reduce((t,e)=>{const r=e.reduce((t,e)=>Math.max(t,Math.abs(e)),0);return Math.max(t,r)},0)),arity:N(1)}};for(const[t,e]of Object.entries(k))j[t]&&(j[t].docs=e);const S={name:"matrix",functions:j};export{S as matrixModule};
+//# sourceMappingURL=matrix.esm.js.map

package/dist/modules/matrix.esm.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"matrix.esm.js","sources":["../../src/utils/debug/getCodeMarker.ts","../../src/errors.ts","../../src/constants/constants.ts","../../src/utils/debug/debugTools.ts","../../src/utils/symbols.ts","../../src/typeGuards/number.ts","../../src/utils/debug/getSourceCodeInfo.ts","../../src/typeGuards/annotatedArrays.ts","../../src/utils/index.ts","../../src/utils/arity.ts","../../src/builtin/modules/matrix/docs.ts","../../src/builtin/modules/matrix/helpers/determinant.ts","../../src/builtin/modules/matrix/helpers/minor.ts","../../src/builtin/modules/matrix/helpers/adjugate.ts","../../src/builtin/modules/matrix/helpers/isSquare.ts","../../src/builtin/modules/matrix/helpers/isIdentity.ts","../../src/builtin/modules/matrix/helpers/matrixMultiply.ts","../../src/builtin/modules/matrix/index.ts","../../src/builtin/modules/matrix/helpers/inverse.ts","../../src/builtin/modules/matrix/helpers/cofactor.ts","../../src/builtin/modules/matrix/helpers/trace.ts","../../src/builtin/modules/matrix/helpers/isSymetric.ts","../../src/builtin/modules/matrix/helpers/isTriangular.ts","../../src/builtin/modules/matrix/helpers/isDiagonal.ts","../../src/builtin/modules/matrix/helpers/isOrthogonal.ts","../../src/builtin/modules/grid/transpose.ts","../../src/builtin/modules/matrix/helpers/band.ts","../../src/builtin/modules/matrix/helpers/isBanded.ts","../../src/builtin/modules/linear-algebra/helpers/gaussJordanElimination.ts","../../src/builtin/modules/matrix/helpers/norm1.ts"],"sourcesContent":["import type { SourceCodeInfo } from '../../tokenizer/token'\n\nexport function getCodeMarker(sourceCodeInfo: SourceCodeInfo): string {\n  if (!sourceCodeInfo.position || !sourceCodeInfo.code)\n    return ''\n\n  const leftPadding = sourceCodeInfo.position.column - 1\n  const rightPadding = sourceCodeInfo.code.length - leftPadding - 1\n  return `${' '.repeat(Math.max(leftPadding, 0))}^${' '.repeat(Math.max(rightPadding, 0))}`\n}\n","import { getCodeMarker } from '../src/utils/debug/getCodeMarker'\nimport type { Arr } from './interface'\nimport type { SourceCodeInfo } from './tokenizer/token'\n\nfunction getDvalaErrorMessage(message: string, sourceCodeInfo?: SourceCodeInfo) {\n  if (!sourceCodeInfo) {\n    return message\n  }\n  const location = `${sourceCodeInfo.position.line}:${sourceCodeInfo.position.column}`\n  const filePathLine = sourceCodeInfo.filePath\n    ? `\\n${sourceCodeInfo.filePath}:${location}`\n    : `\\nLocation ${location}`\n  const codeLine = `\\n${sourceCodeInfo.code}`\n  const codeMarker = `\\n${getCodeMarker(sourceCodeInfo)}`\n  return `${message}${filePathLine}${codeLine}${codeMarker}`\n}\n\nexport class RecurSignal extends Error {\n  public params: Arr\n  constructor(params: Arr) {\n    super(`recur, params: ${params}`)\n    Object.setPrototypeOf(this, RecurSignal.prototype)\n    this.name = 'RecurSignal'\n    this.params = params\n  }\n}\n\nexport class DvalaError extends Error {\n  public readonly sourceCodeInfo?: SourceCodeInfo\n  public readonly shortMessage: string\n  constructor(err: unknown, sourceCodeInfo: SourceCodeInfo | undefined) {\n    const message = err instanceof Error\n      ? err.message\n      : `${err}`\n\n    super(getDvalaErrorMessage(message, sourceCodeInfo))\n    this.shortMessage = message\n    this.sourceCodeInfo = sourceCodeInfo\n    Object.setPrototypeOf(this, DvalaError.prototype)\n    this.name = 'DvalaError'\n  }\n\n  public getCodeMarker(): string | undefined {\n    return this.sourceCodeInfo && getCodeMarker(this.sourceCodeInfo)\n  }\n}\n\nexport class UserDefinedError extends DvalaError {\n  public userMessage: string\n  constructor(message: string, sourceCodeInfo?: SourceCodeInfo) {\n    super(message, sourceCodeInfo)\n    this.userMessage = message\n    Object.setPrototypeOf(this, UserDefinedError.prototype)\n    this.name = 'UserDefinedError'\n  }\n}\n\nexport class AssertionError extends DvalaError {\n  constructor(message: string | Error, sourceCodeInfo?: SourceCodeInfo) {\n    super(message, sourceCodeInfo)\n    Object.setPrototypeOf(this, AssertionError.prototype)\n    this.name = 'AssertionError'\n  }\n}\n\nexport class UndefinedSymbolError extends DvalaError {\n  public symbol: string\n  constructor(symbolName: string, sourceCodeInfo?: SourceCodeInfo) {\n    const message = `Undefined symbol '${symbolName}'.`\n    super(message, sourceCodeInfo)\n    this.symbol = symbolName\n    Object.setPrototypeOf(this, UndefinedSymbolError.prototype)\n    this.name = 'UndefinedSymbolError'\n  }\n}\n\nexport function isDvalaError(error: unknown): error is DvalaError {\n  return error instanceof DvalaError\n}\n","export const NodeTypes = {\n  Number: 1,\n  String: 2,\n  NormalExpression: 3,\n  SpecialExpression: 4,\n  UserDefinedSymbol: 5,\n  NormalBuiltinSymbol: 6,\n  SpecialBuiltinSymbol: 7,\n  ReservedSymbol: 8,\n  Binding: 9,\n  Spread: 10,\n} as const\n\nconst NodeTypesSet = new Set(Object.values(NodeTypes))\n\nexport type NodeType = typeof NodeTypes[keyof typeof NodeTypes]\n\nexport function getNodeTypeName(type: NodeType): keyof typeof NodeTypes {\n  return Object.keys(NodeTypes).find(key => NodeTypes[key as keyof typeof NodeTypes] === type) as keyof typeof NodeTypes\n}\n\n// TODO, is this needed?\nexport function isNodeType(type: unknown): type is NodeType {\n  return typeof type === 'number' && NodeTypesSet.has(type as NodeType)\n}\n\nconst functionTypes = [\n  'UserDefined',\n  'Partial',\n  'Comp',\n  'Constantly',\n  'Juxt',\n  'Complement',\n  'EveryPred',\n  'SomePred',\n  'Fnull',\n  'Builtin',\n  'SpecialBuiltin',\n  'NativeJsFunction',\n  'Module',\n] as const\n\nconst functionTypeSet = new Set(functionTypes)\n\nexport type FunctionType = typeof functionTypes[number]\n\nexport function isFunctionType(type: unknown): type is FunctionType {\n  return typeof type === 'string' && functionTypeSet.has(type as FunctionType)\n}\n","import { getNodeTypeName, isFunctionType, isNodeType } from '../../constants/constants'\nimport type { AstNode, DvalaFunction } from '../../parser/types'\nimport { FUNCTION_SYMBOL } from '../symbols'\n\nfunction isDvalaFunction(func: unknown): func is DvalaFunction {\n  if (func === null || typeof func !== 'object')\n    return false\n\n  return FUNCTION_SYMBOL in func && 'functionType' in func && isFunctionType(func.functionType)\n}\n\nfunction isNode(value: unknown): value is AstNode {\n  if (!Array.isArray(value) || value.length < 2)\n    return false\n  return isNodeType(value[0])\n}\n\nexport function valueToString(value: unknown): string {\n  if (isDvalaFunction(value))\n    // eslint-disable-next-line ts/no-unsafe-member-access\n    return `<function ${(value as any).name || '\\u03BB'}>`\n\n  if (isNode(value))\n    return `${getNodeTypeName(value[0])}-node`\n\n  if (value === null)\n    return 'null'\n\n  if (typeof value === 'object' && value instanceof RegExp)\n    return `${value}`\n\n  if (typeof value === 'object' && value instanceof Error)\n    return value.toString()\n\n  return JSON.stringify(value)\n}\n","export const FUNCTION_SYMBOL = '^^fn^^'\nexport const REGEXP_SYMBOL = '^^re^^'\nexport const EFFECT_SYMBOL = '^^ef^^'\n","import { DvalaError } from '../errors'\nimport type { SourceCodeInfo } from '../tokenizer/token'\nimport { valueToString } from '../utils/debug/debugTools'\nimport { getSourceCodeInfo } from '../utils/debug/getSourceCodeInfo'\n\ntype SignOptions =\n  | {\n    positive?: true\n    negative?: never\n    nonPositive?: never\n    nonNegative?: never\n    zero?: never\n    nonZero?: never\n  }\n  | {\n    positive?: never\n    negative?: true\n    nonPositive?: never\n    nonNegative?: never\n    zero?: never\n    nonZero?: never\n  }\n  | {\n    positive?: never\n    negative?: never\n    nonPositive?: true\n    nonNegative?: never\n    zero?: never\n    nonZero?: never\n  }\n  | {\n    positive?: never\n    negative?: never\n    nonPositive?: never\n    nonNegative?: true\n    zero?: never\n    nonZero?: never\n  }\n  | {\n    positive?: never\n    negative?: never\n    nonPositive?: never\n    nonNegative?: never\n    zero?: true\n    nonZero?: never\n  }\n  | {\n    positive?: never\n    negative?: never\n    nonPositive?: never\n    nonNegative?: never\n    zero?: never\n    nonZero?: true\n  }\n\ntype GtOptions =\n  | {\n    gt?: number\n    gte?: never\n  }\n  | {\n    gt?: never\n    gte?: number\n  }\n\ntype LtOptions =\n  | {\n    lt?: number\n    lte?: never\n  }\n  | {\n    lt?: never\n    lte?: number\n  }\n\ntype NumberOptions = {\n  integer?: true\n  finite?: true\n} & SignOptions &\nGtOptions &\nLtOptions\n\nfunction getRangeString(options: NumberOptions): string {\n  const hasUpperAndLowerBound\n    = (typeof options.gt === 'number' || typeof options.gte === 'number')\n      && (typeof options.lt === 'number' || typeof options.lte === 'number')\n  if (hasUpperAndLowerBound) {\n    return `${typeof options.gt === 'number' ? `${options.gt} < n ` : `${options.gte} <= n `}${\n      typeof options.lt === 'number' ? `< ${options.lt}` : `<= ${options.lte}`\n    }`\n  }\n  else if (typeof options.gt === 'number' || typeof options.gte === 'number') {\n    return `${typeof options.gt === 'number' ? `n > ${options.gt}` : `n >= ${options.gte}`}`\n  }\n  else if (typeof options.lt === 'number' || typeof options.lte === 'number') {\n    return `${typeof options.lt === 'number' ? `n < ${options.lt}` : `n <= ${options.lte}`}`\n  }\n  else { return '' }\n}\n\nfunction getSignString(options: NumberOptions): string {\n  return options.positive\n    ? 'positive'\n    : options.negative\n      ? 'negative'\n      : options.nonNegative\n        ? 'non negative'\n        : options.nonPositive\n          ? 'non positive'\n          : options.nonZero\n            ? 'non zero'\n            : ''\n}\n\nfunction getNumberTypeName(options: NumberOptions): string {\n  if (options.zero)\n    return 'zero'\n\n  const sign = getSignString(options)\n  const numberType = options.integer ? 'integer' : 'number'\n  const finite = options.finite ? 'finite' : ''\n  const range = getRangeString(options)\n\n  return [sign, finite, numberType, range].filter(x => !!x).join(' ')\n}\n\nexport function isNumber(value: unknown, options: NumberOptions = {}): value is number {\n  if (typeof value !== 'number')\n    return false\n\n  if (Number.isNaN(value))\n    return false\n\n  if (options.integer && !Number.isInteger(value))\n    return false\n\n  if (options.finite && !Number.isFinite(value))\n    return false\n\n  if (options.zero && value !== 0)\n    return false\n\n  if (options.nonZero && value === 0)\n    return false\n\n  if (options.positive && value <= 0)\n    return false\n\n  if (options.negative && value >= 0)\n    return false\n\n  if (options.nonPositive && value > 0)\n    return false\n\n  if (options.nonNegative && value < 0)\n    return false\n\n  if (typeof options.gt === 'number' && value <= options.gt)\n    return false\n\n  if (typeof options.gte === 'number' && value < options.gte)\n    return false\n\n  if (typeof options.lt === 'number' && value >= options.lt)\n    return false\n\n  if (typeof options.lte === 'number' && value > options.lte)\n    return false\n\n  return true\n}\n\nexport function assertNumber(\n  value: unknown,\n  sourceCodeInfo?: SourceCodeInfo,\n  options: NumberOptions = {},\n): asserts value is number {\n  if (!isNumber(value, options)) {\n    throw new DvalaError(\n      `Expected ${getNumberTypeName(options)}, got ${valueToString(value)}.`,\n      getSourceCodeInfo(value, sourceCodeInfo),\n    )\n  }\n}\n\nexport function asNumber(\n  value: unknown,\n  sourceCodeInfo: SourceCodeInfo | undefined,\n  options: NumberOptions = {},\n): number {\n  assertNumber(value, sourceCodeInfo, options)\n  return value\n}\n","import type { SourceCodeInfo } from '../../tokenizer/token'\n\nexport function getSourceCodeInfo(anyValue: any, sourceCodeInfo: SourceCodeInfo | undefined): SourceCodeInfo | undefined {\n  // eslint-disable-next-line ts/no-unsafe-return, ts/no-unsafe-member-access\n  return anyValue?.sourceCodeInfo ?? sourceCodeInfo\n}\n","import { DvalaError } from '../errors'\nimport type { Any } from '../interface'\nimport type { SourceCodeInfo } from '../tokenizer/token'\nimport { isNumber } from './number'\n\nconst annotatedArrays = new WeakSet<unknown[]>()\nconst vectors = new WeakSet<unknown[]>()\nconst notVectors = new WeakSet<unknown[]>()\nconst matrices = new WeakSet<unknown[]>()\nconst notMatrices = new WeakSet<unknown[]>()\nconst grids = new WeakSet<unknown[]>()\nconst notGrids = new WeakSet<unknown[]>()\n\nexport function annotate<T>(value: T): T {\n  if (!Array.isArray(value)) {\n    return value\n  }\n  if (annotatedArrays.has(value)) {\n    return value\n  }\n  isVector(value)\n  if (!isMatrix(value)) {\n    isGrid(value)\n  }\n\n  return value\n}\nexport function isVector(vector: unknown): vector is number[] {\n  if (!Array.isArray(vector)) {\n    return false\n  }\n\n  if (vectors.has(vector)) {\n    return true\n  }\n  if (notVectors.has(vector)) {\n    return false\n  }\n\n  if (vector.every(elem => isNumber(elem))) {\n    annotatedArrays.add(vector)\n    vectors.add(vector)\n    return true\n  }\n  notVectors.add(vector)\n  return false\n}\n\nexport function assertVector(vector: unknown, sourceCodeInfo: SourceCodeInfo | undefined): asserts vector is number[] {\n  if (!isVector(vector)) {\n    throw new DvalaError(`Expected a vector, but got ${vector}`, sourceCodeInfo)\n  }\n}\n\nexport function is2dVector(vector: unknown): vector is [number, number] {\n  if (!isVector(vector)) {\n    return false\n  }\n  return vector.length === 2\n}\nexport function assert2dVector(vector: unknown, sourceCodeInfo: SourceCodeInfo | undefined): asserts vector is [number, number] {\n  if (!is2dVector(vector)) {\n    throw new DvalaError(`Expected a 2d vector, but got ${vector}`, sourceCodeInfo)\n  }\n}\n\nexport function is3dVector(vector: unknown): vector is [number, number, number] {\n  if (!isVector(vector)) {\n    return false\n  }\n  return vector.length === 3\n}\nexport function assert3dVector(vector: unknown, sourceCodeInfo: SourceCodeInfo | undefined): asserts vector is [number, number, number] {\n  if (!is3dVector(vector)) {\n    throw new DvalaError(`Expected a 3d vector, but got ${vector}`, sourceCodeInfo)\n  }\n}\n\nexport function assertNonEmptyVector(vector: unknown, sourceCodeInfo: SourceCodeInfo | undefined): asserts vector is number[] {\n  assertVector(vector, sourceCodeInfo)\n  if (vector.length === 0) {\n    throw new DvalaError(`Expected a non empty vector, but got ${vector}`, sourceCodeInfo)\n  }\n}\n\nexport function isGrid(grid: unknown): grid is unknown[][] {\n  if (!Array.isArray(grid)) {\n    return false\n  }\n  if (grids.has(grid)) {\n    return true\n  }\n  if (notGrids.has(grid)) {\n    return false\n  }\n  if (grid.length === 0) {\n    notGrids.add(grid)\n    return false\n  }\n  if (!Array.isArray(grid[0])) {\n    notGrids.add(grid)\n    return false\n  }\n  const nbrOfCols = grid[0].length\n  for (const row of grid.slice(1)) {\n    if (!Array.isArray(row)) {\n      notGrids.add(grid)\n      return false\n    }\n    if (row.length !== nbrOfCols) {\n      notGrids.add(grid)\n      return false\n    }\n  }\n  annotatedArrays.add(grid)\n  grids.add(grid)\n  return true\n}\n\nexport function assertGrid(grid: unknown, sourceCodeInfo: SourceCodeInfo | undefined): asserts grid is Any[][] {\n  if (!isGrid(grid)) {\n    throw new DvalaError(`Expected a grid, but got ${grid}`, sourceCodeInfo)\n  }\n}\n\nexport function isMatrix(matrix: unknown): matrix is number[][] {\n  if (!Array.isArray(matrix)) {\n    return false\n  }\n  if (matrices.has(matrix)) {\n    return true\n  }\n  if (notMatrices.has(matrix)) {\n    return false\n  }\n  if (matrix.length === 0) {\n    notMatrices.add(matrix)\n    return false\n  }\n  if (!Array.isArray(matrix[0]) || matrix[0].length === 0) {\n    notMatrices.add(matrix)\n    return false\n  }\n  const nbrOfCols = matrix[0].length\n  for (const row of matrix) {\n    if (!Array.isArray(row) || row.length !== nbrOfCols || row.some(cell => !isNumber(cell))) {\n      notMatrices.add(matrix)\n      return false\n    }\n  }\n  annotatedArrays.add(matrix)\n  grids.add(matrix)\n  matrices.add(matrix)\n  return true\n}\n\nexport function assertMatrix(matrix: unknown, sourceCodeInfo: SourceCodeInfo | undefined): asserts matrix is number[][] {\n  if (!isMatrix(matrix)) {\n    throw new DvalaError(`Expected a matrix, but got ${matrix}`, sourceCodeInfo)\n  }\n}\n\nexport function assertSquareMatrix(matrix: unknown, sourceCodeInfo: SourceCodeInfo | undefined): asserts matrix is number[][] {\n  if (!isMatrix(matrix)) {\n    throw new DvalaError(`Expected a matrix, but got ${matrix}`, sourceCodeInfo)\n  }\n  if (matrix.length !== matrix[0]!.length) {\n    throw new DvalaError(`Expected square matrix, but got ${matrix.length} and ${matrix[0]!.length}`, sourceCodeInfo)\n  }\n}\n\nexport function isSquareMatrix(matrix: unknown): matrix is number[][] {\n  if (!isMatrix(matrix)) {\n    return false\n  }\n  if (matrix.length !== matrix[0]!.length) {\n    return false\n  }\n  return true\n}\n","import type { Any, Coll, Obj } from '../interface'\nimport type { SourceCodeInfo } from '../tokenizer/token'\nimport { asAny, isColl, isObj, isRegularExpression } from '../typeGuards/dvala'\nimport { isNumber } from '../typeGuards/number'\nimport { asString, assertStringOrNumber } from '../typeGuards/string'\nimport { isUnknownRecord } from '../typeGuards'\nimport { DvalaError } from '../errors'\n\nexport function collHasKey(coll: unknown, key: string | number): boolean {\n  if (!isColl(coll))\n    return false\n\n  if (typeof coll === 'string' || Array.isArray(coll)) {\n    if (!isNumber(key, { integer: true }))\n      return false\n\n    return key >= 0 && key < coll.length\n  }\n  return !!Object.getOwnPropertyDescriptor(coll, key)\n}\n\nexport function compare<T extends string | number>(a: T, b: T, sourceCodeInfo: SourceCodeInfo | undefined): number {\n  assertStringOrNumber(a, sourceCodeInfo)\n  assertStringOrNumber(b, sourceCodeInfo)\n\n  if (typeof a === 'string' && typeof b === 'string') {\n    return a < b ? -1 : a > b ? 1 : 0\n  }\n  if (typeof a === 'number' && typeof b === 'number') {\n    return Math.sign((a) - (b))\n  }\n  throw new DvalaError(`Cannot compare values of different types: ${typeof a} and ${typeof b}`, sourceCodeInfo)\n}\n\nexport function deepEqual(a: unknown, b: unknown, sourceCodeInfo?: SourceCodeInfo): boolean {\n  if (a === b)\n    return true\n\n  if (typeof a === 'number' && typeof b === 'number')\n    return approxEqual(a, b)\n\n  if (Array.isArray(a) && Array.isArray(b)) {\n    if (a.length !== b.length)\n      return false\n\n    for (let i = 0; i < a.length; i += 1) {\n      if (!deepEqual(asAny(a[i], sourceCodeInfo), asAny(b[i], sourceCodeInfo), sourceCodeInfo))\n        return false\n    }\n    return true\n  }\n  if (isRegularExpression(a) && isRegularExpression(b))\n    return a.s === b.s && a.f === b.f\n\n  if (isUnknownRecord(a) && isUnknownRecord(b)) {\n    const aKeys = Object.keys(a)\n    const bKeys = Object.keys(b)\n    if (aKeys.length !== bKeys.length)\n      return false\n\n    for (let i = 0; i < aKeys.length; i += 1) {\n      const key = asString(aKeys[i], sourceCodeInfo)\n      if (!deepEqual(a[key], b[key], sourceCodeInfo))\n        return false\n    }\n    return true\n  }\n  return false\n}\n\nexport function toNonNegativeInteger(num: number): number {\n  return Math.max(0, Math.ceil(num))\n}\n\nexport function toAny(value: unknown): Any {\n  return (value ?? null) as Any\n}\n\nfunction clone<T>(value: T): T {\n  if (isObj(value)) {\n    return Object.entries(value).reduce((result: Obj, entry) => {\n      const [key, val] = entry\n      result[key] = clone(val)\n      return result\n    }, {}) as T\n  }\n  if (Array.isArray(value))\n    // eslint-disable-next-line ts/no-unsafe-return\n    return value.map(item => clone(item)) as unknown as T\n\n  return value\n}\n\nexport function cloneColl<T extends Coll>(value: T): T {\n  return clone(value)\n}\n\nexport function joinSets<T>(...results: Set<T>[]): Set<T> {\n  const result = new Set<T>()\n  for (const symbols of results)\n    symbols.forEach(symbol => result.add(symbol))\n\n  return result\n}\n\nexport function addToSet<T>(target: Set<T>, source: Set<T>): void {\n  source.forEach(symbol => target.add(symbol))\n}\n\nexport const EPSILON = 1e-10\n\nexport function approxEqual(a: number, b: number, epsilon: number = EPSILON): boolean {\n  if (a === b) {\n    return true\n  }\n\n  const diff = Math.abs(a - b)\n\n  if (a === 0 || b === 0 || diff < epsilon) {\n    // Use absolute error for values near zero\n    return diff < epsilon\n  }\n  const absA = Math.abs(a)\n  const absB = Math.abs(b)\n\n  // Use relative error for larger values\n  return diff / (absA + absB) < epsilon\n}\n\nexport function approxZero(value: number): boolean {\n  return Math.abs(value) < EPSILON\n}\n\nexport function smartTrim(str: string, minIndent = 0): string {\n  const lines = str.split('\\n')\n  while (lines[0]?.match(/^\\s*$/)) {\n    lines.shift() // Remove leading empty lines\n  }\n  while (lines[lines.length - 1]?.match(/^\\s*$/)) {\n    lines.pop() // Remove trailing empty lines\n  }\n  const indent = lines.reduce((acc, line) => {\n    if (line.match(/^\\s*$/))\n      return acc // Skip empty lines\n    const lineIndent = line.match(/^\\s*/)![0].length\n    return Math.min(acc, lineIndent)\n  }, Infinity)\n  return lines.map(line => ' '.repeat(minIndent) + line.slice(indent)).join('\\n').trimEnd()\n}\n","import type { Arity } from '../builtin/interface'\nimport { DvalaError } from '../errors'\nimport type { FunctionLike } from '../parser/types'\nimport type { SourceCodeInfo } from '../tokenizer/token'\nimport { isColl } from '../typeGuards/dvala'\nimport { valueToString } from './debug/debugTools'\n\nexport function arityAccepts(arity: Arity, nbrOfParams: number): boolean {\n  const { min, max } = arity\n  if (typeof min === 'number' && nbrOfParams < min) {\n    return false\n  }\n  if (typeof max === 'number' && nbrOfParams > max) {\n    return false\n  }\n  return true\n}\n\nexport function arityAcceptsMin(arity: Arity, nbrOfParams: number): boolean {\n  const { min } = arity\n  if (typeof min === 'number' && nbrOfParams < min) {\n    return false\n  }\n  return true\n}\n\nexport function getCommonArityFromFunctions(params: FunctionLike[]): Arity | null {\n  return params.reduce((acc: Arity | null, param): Arity | null => {\n    if (acc === null) {\n      return null\n    }\n    const arity: Arity = (typeof param === 'number' || isColl(param)) ? toFixedArity(1) : param.arity\n    const { min: aMin, max: aMax } = arity\n    const { min: bMin, max: bMax } = acc\n    const min = typeof aMin === 'number' && typeof bMin === 'number'\n      ? Math.max(aMin, bMin)\n      : typeof aMin === 'number' ? aMin : typeof bMin === 'number' ? bMin : undefined\n    const max = typeof aMax === 'number' && typeof bMax === 'number'\n      ? Math.min(aMax, bMax)\n      : typeof aMax === 'number' ? aMax : typeof bMax === 'number' ? bMax : undefined\n\n    if (typeof min === 'number' && typeof max === 'number' && min > max) {\n      return null\n    }\n\n    return { min, max }\n  }, {})\n}\n\nexport function getArityFromFunction(param: FunctionLike): Arity {\n  return (typeof param === 'number' || isColl(param)) ? toFixedArity(1) : param.arity\n}\n\nexport function assertNumberOfParams(arity: Arity, length: number, sourceCodeInfo: SourceCodeInfo | undefined): void {\n  const { min, max } = arity\n  if (typeof min === 'number' && length < min) {\n    throw new DvalaError(\n      `Wrong number of arguments, expected at least ${min}, got ${valueToString(length)}.`,\n      sourceCodeInfo,\n    )\n  }\n\n  if (typeof max === 'number' && length > max) {\n    throw new DvalaError(\n      `Wrong number of arguments, expected at most ${max}, got ${valueToString(length)}.`,\n      sourceCodeInfo,\n    )\n  }\n}\n\nexport function canBeOperator(count: Arity): boolean {\n  if (typeof count.max === 'number' && count.max < 2) {\n    return false\n  }\n\n  if (typeof count.min === 'number' && count.min > 2) {\n    return false\n  }\n\n  return true\n}\n\nexport function toFixedArity(arity: number): Arity {\n  return { min: arity, max: arity }\n}\n","import type { FunctionDocs } from '../../interface'\n\nexport const moduleDocs: Record<string, FunctionDocs> = {\n  'mul': {\n    category: 'matrix',\n    description: 'Multiplies two `matrices` using standard `matrix` multiplication based on **dot products** of rows and columns.',\n    returns: {\n      type: 'matrix',\n    },\n    args: {\n      a: {\n        type: 'matrix',\n      },\n      b: {\n        type: 'matrix',\n      },\n    },\n    variants: [\n      {\n        argumentNames: [\n          'a',\n          'b',\n        ],\n      },\n    ],\n    examples: [\n      'let { mul } = import(matrix);\\nmul([[1, 2], [3, 4]], [[5, 6], [7, 8]])',\n      'let { mul } = import(matrix);\\nmul([[1, 2, 3], [4, 5, 6]], [[7, 8], [9, 10], [11, 12]])',\n    ],\n    seeAlso: ['matrix.det', 'matrix.inv'],\n  },\n  'det': {\n    category: 'matrix',\n    description: 'Calculates the **determinant** of a square matrix.',\n    returns: {\n      type: 'number',\n    },\n    args: {\n      m: {\n        type: 'matrix',\n        description: 'The `matrix` to calculate the determinant of.',\n      },\n    },\n    variants: [\n      {\n        argumentNames: [\n          'm',\n        ],\n      },\n    ],\n    examples: [\n      'let { det } = import(matrix);\\ndet([[1, 2], [3, 4]])',\n      'let { det } = import(matrix);\\ndet([[1, 2, 3], [4, 5, 6], [7, 8, 9]])',\n    ],\n    seeAlso: ['matrix.inv', 'matrix.cofactor', 'matrix.adj', 'matrix.trace', 'matrix.rank', 'matrix.invertible?', 'matrix.mul', 'matrix.minor'],\n  },\n  'inv': {\n    category: 'matrix',\n    description: 'Calculates the **inverse** of a square matrix.',\n    returns: {\n      type: 'matrix',\n    },\n    args: {\n      m: {\n        type: 'matrix',\n        description: 'The `matrix` to calculate the inverse of.',\n      },\n    },\n    variants: [\n      {\n        argumentNames: [\n          'm',\n        ],\n      },\n    ],\n    examples: [\n      'let { inv } = import(matrix);\\ninv([[1, 2], [3, 4]])',\n      'let { inv } = import(matrix);\\ninv([[1, 2, 3], [4, 5, 7], [7, 8, 10]])',\n    ],\n    seeAlso: ['matrix.det', 'matrix.adj', 'matrix.invertible?', 'linear-algebra.solve', 'matrix.mul', 'matrix.orthogonal-matrix?'],\n  },\n  'adj': {\n    category: 'matrix',\n    description: 'Calculates the **adjugate** of a square matrix.',\n    returns: {\n      type: 'matrix',\n    },\n    args: {\n      m: {\n        type: 'matrix',\n        description: 'The `matrix` to calculate the adjugate of.',\n      },\n    },\n    variants: [\n      {\n        argumentNames: [\n          'm',\n        ],\n      },\n    ],\n    examples: [\n      'let { adj } = import(matrix);\\nadj([[1, 2], [3, 4]])',\n      'let { adj } = import(matrix);\\nadj([[1, 2, 3], [4, 5, 6], [7, 8, 9]])',\n      'let { adj } = import(matrix);\\nadj([[1, 2, 3], [7, 8, 9], [4, 5, 6]])',\n    ],\n    seeAlso: ['matrix.cofactor', 'matrix.det', 'matrix.inv'],\n  },\n  'cofactor': {\n    category: 'matrix',\n    description: 'Calculates the **cofactor** of a square matrix.',\n    returns: {\n      type: 'matrix',\n    },\n    args: {\n      m: {\n        type: 'matrix',\n        description: 'The `matrix` to calculate the cofactor of.',\n      },\n    },\n    variants: [\n      {\n        argumentNames: [\n          'm',\n        ],\n      },\n    ],\n    examples: [\n      'let { cofactor } = import(matrix);\\ncofactor([[1, 2], [3, 4]])',\n      'let { cofactor } = import(matrix);\\ncofactor([[1, 2, 3], [4, 5, 6], [7, 8, 9]])',\n      'let { cofactor } = import(matrix);\\ncofactor([[1, 2, 3], [7, 8, 9], [4, 5, 6]])',\n    ],\n    seeAlso: ['matrix.adj', 'matrix.minor', 'matrix.det'],\n  },\n  'minor': {\n    category: 'matrix',\n    description: 'Calculates the **minor** of a square matrix.',\n    returns: {\n      type: 'matrix',\n    },\n    args: {\n      m: {\n        type: 'matrix',\n        description: 'The `matrix` to calculate the minor of.',\n      },\n      row: {\n        type: 'integer',\n        description: 'The row index of the element to calculate the minor for.',\n      },\n      col: {\n        type: 'integer',\n        description: 'The column index of the element to calculate the minor for.',\n      },\n    },\n    variants: [\n      {\n        argumentNames: [\n          'm',\n          'row',\n          'col',\n        ],\n      },\n    ],\n    examples: [\n      'let { minor } = import(matrix);\\nminor([[1, 2], [3, 4]], 0, 1)',\n      'let { minor } = import(matrix);\\nminor([[1, 2, 3], [4, 5, 6], [7, 8, 9]], 1, 1)',\n    ],\n    seeAlso: ['matrix.cofactor', 'matrix.det'],\n  },\n  'trace': {\n    category: 'matrix',\n    description: 'Calculates the **trace** of a square matrix.',\n    returns: {\n      type: 'number',\n    },\n    args: {\n      m: {\n        type: 'matrix',\n        description: 'The `matrix` to calculate the trace of.',\n      },\n    },\n    variants: [\n      {\n        argumentNames: [\n          'm',\n        ],\n      },\n    ],\n    examples: [\n      'let { trace } = import(matrix);\\ntrace([[1, 2], [3, 4]])',\n      'let { trace } = import(matrix);\\ntrace([[1, 2, 3], [4, 5, 6], [7, 8, 9]])',\n    ],\n    seeAlso: ['matrix.det', 'matrix.diagonal?'],\n  },\n  'symmetric?': {\n    category: 'matrix',\n    description: 'Checks if a `matrix` is **symmetric**.',\n    returns: {\n      type: 'boolean',\n    },\n    args: {\n      m: {\n        type: 'matrix',\n        description: 'The `matrix` to check for symmetry.',\n      },\n    },\n    variants: [\n      {\n        argumentNames: [\n          'm',\n        ],\n      },\n    ],\n    examples: [\n      'let { symmetric? } = import(matrix);\\nsymmetric?([[1, 2], [2, 1]])',\n      'let { symmetric? } = import(matrix);\\nsymmetric?([[1, 2, 3], [2, 1, 4], [3, 4, 1]])',\n    ],\n    seeAlso: ['matrix.orthogonal-matrix?', 'matrix.diagonal?', 'matrix.square?', 'matrix.hilbert'],\n  },\n  'triangular?': {\n    category: 'matrix',\n    description: 'Checks if a `matrix` is **triangular**.',\n    returns: {\n      type: 'boolean',\n    },\n    args: {\n      m: {\n        type: 'matrix',\n        description: 'The `matrix` to check for triangularity.',\n      },\n    },\n    variants: [\n      {\n        argumentNames: [\n          'm',\n        ],\n      },\n    ],\n    examples: [\n      'let { triangular? } = import(matrix);\\ntriangular?([[2, 0], [0, 1]])',\n      'let { triangular? } = import(matrix);\\ntriangular?([[1, 2, 3], [0, 4, 5], [0, 0, 6]])',\n    ],\n    seeAlso: ['matrix.upper-triangular?', 'matrix.lower-triangular?', 'matrix.diagonal?', 'matrix.banded?'],\n  },\n  'upper-triangular?': {\n    category: 'matrix',\n    description: 'Checks if a `matrix` is **upper triangular**.',\n    returns: {\n      type: 'boolean',\n    },\n    args: {\n      m: {\n        type: 'matrix',\n        description: 'The `matrix` to check for upper triangularity.',\n      },\n    },\n    variants: [\n      {\n        argumentNames: [\n          'm',\n        ],\n      },\n    ],\n    examples: [\n      'let { upper-triangular? } = import(matrix);\\nupper-triangular?([[1, 2], [0, 3]])',\n      'let { upper-triangular? } = import(matrix);\\nupper-triangular?([[1, 2, 3], [0, 4, 5], [0, 0, 6]])',\n    ],\n    seeAlso: ['matrix.lower-triangular?', 'matrix.triangular?', 'matrix.diagonal?'],\n  },\n  'lower-triangular?': {\n    category: 'matrix',\n    description: 'Checks if a `matrix` is **lower triangular**.',\n    returns: {\n      type: 'boolean',\n    },\n    args: {\n      m: {\n        type: 'matrix',\n        description: 'The `matrix` to check for lower triangularity.',\n      },\n    },\n    variants: [\n      {\n        argumentNames: [\n          'm',\n        ],\n      },\n    ],\n    examples: [\n      'let { lower-triangular? } = import(matrix);\\nlower-triangular?([[1, 0], [2, 3]])',\n      'let { lower-triangular? } = import(matrix);\\nlower-triangular?([[1, 0, 0], [2, 3, 0], [4, 5, 6]])',\n    ],\n    seeAlso: ['matrix.upper-triangular?', 'matrix.triangular?', 'matrix.diagonal?'],\n  },\n  'diagonal?': {\n    category: 'matrix',\n    description: 'Checks if a `matrix` is **diagonal**.',\n    returns: {\n      type: 'boolean',\n    },\n    args: {\n      m: {\n        type: 'matrix',\n        description: 'The `matrix` to check for diagonal property.',\n      },\n    },\n    variants: [\n      {\n        argumentNames: [\n          'm',\n        ],\n      },\n    ],\n    examples: [\n      'let { diagonal? } = import(matrix);\\ndiagonal?([[1, 0], [0, 2]])',\n      'let { diagonal? } = import(matrix);\\ndiagonal?([[1, 0, 0], [0, 2, 0], [0, 0, 3]])',\n      'let { diagonal? } = import(matrix);\\ndiagonal?([[1, 0, 0], [2, 2, 2], [0, 0, 3]])',\n    ],\n    seeAlso: ['matrix.identity?', 'matrix.symmetric?', 'matrix.triangular?', 'matrix.trace', 'matrix.upper-triangular?', 'matrix.lower-triangular?', 'matrix.band', 'matrix.banded?'],\n  },\n  'square?': {\n    category: 'matrix',\n    description: 'Checks if a `matrix` is **square**.',\n    returns: {\n      type: 'boolean',\n    },\n    args: {\n      m: {\n        type: 'matrix',\n        description: 'The `matrix` to check for square property.',\n      },\n    },\n    variants: [\n      {\n        argumentNames: [\n          'm',\n        ],\n      },\n    ],\n    examples: [\n      'let { square? } = import(matrix);\\nsquare?([[1, 2], [3, 4]])',\n      'let { square? } = import(matrix);\\nsquare?([[1, 2, 3], [4, 5, 6], [7, 8, 9]])',\n      'let { square? } = import(matrix);\\nsquare?([[1, 2, 3], [4, 5, 6]])',\n    ],\n    seeAlso: ['matrix.symmetric?', 'matrix.identity?', 'matrix.invertible?'],\n  },\n  'orthogonal-matrix?': {\n    category: 'matrix',\n    description: 'Checks if a `matrix` is **orthogonal**.',\n    returns: {\n      type: 'boolean',\n    },\n    args: {\n      m: {\n        type: 'matrix',\n        description: 'The `matrix` to check for **orthogonality**.',\n      },\n    },\n    variants: [\n      {\n        argumentNames: [\n          'm',\n        ],\n      },\n    ],\n    examples: [\n      'let { orthogonal-matrix? } = import(matrix);\\northogonal-matrix?([[1, 0], [0, 1]])',\n      'let { orthogonal-matrix? } = import(matrix);\\northogonal-matrix?([[1, 0], [0, -1]])',\n      'let { orthogonal-matrix? } = import(matrix);\\northogonal-matrix?([[1, 2], [3, 4]])',\n    ],\n    seeAlso: ['matrix.symmetric?', 'matrix.inv', 'matrix.identity?', 'linear-algebra.orthogonal?'],\n  },\n  'identity?': {\n    category: 'matrix',\n    description: 'Checks if a `matrix` is an **identity matrix**.',\n    returns: {\n      type: 'boolean',\n    },\n    args: {\n      m: {\n        type: 'matrix',\n        description: 'The `matrix` to check for identity property.',\n      },\n    },\n    variants: [\n      {\n        argumentNames: [\n          'm',\n        ],\n      },\n    ],\n    examples: [\n      'let { identity? } = import(matrix);\\nidentity?([[1, 0], [0, 1]])',\n      'let { identity? } = import(matrix);\\nidentity?([[1, 0, 0], [0, 1, 0], [0, 0, 1]])',\n      'let { identity? } = import(matrix);\\nidentity?([[1, 0, 0], [0, 1, 0], [0, 0, 0]])',\n    ],\n    seeAlso: ['matrix.diagonal?', 'matrix.square?', 'matrix.orthogonal-matrix?'],\n  },\n  'invertible?': {\n    category: 'matrix',\n    description: 'Checks if a `matrix` is **invertible**.',\n    returns: {\n      type: 'boolean',\n    },\n    args: {\n      m: {\n        type: 'matrix',\n        description: 'The `matrix` to check for invertibility.',\n      },\n    },\n    variants: [\n      {\n        argumentNames: [\n          'm',\n        ],\n      },\n    ],\n    examples: [\n      'let { invertible? } = import(matrix);\\ninvertible?([[1, 2], [3, 4]])',\n      'let { invertible? } = import(matrix);\\ninvertible?([[1, 2, 3], [4, 5, 6], [7, 8, 9]])',\n      'let { invertible? } = import(matrix);\\ninvertible?([[1, 2], [2, 4]])',\n    ],\n    seeAlso: ['matrix.det', 'matrix.inv', 'matrix.rank', 'matrix.square?'],\n  },\n  'hilbert': {\n    category: 'matrix',\n    description: 'Generates a **Hilbert matrix** of size `n`.',\n    returns: {\n      type: 'matrix',\n    },\n    args: {\n      n: {\n        type: 'integer',\n        description: 'The size of the Hilbert matrix.',\n      },\n    },\n    variants: [\n      {\n        argumentNames: [\n          'n',\n        ],\n      },\n    ],\n    examples: [\n      'let { hilbert } = import(matrix);\\nhilbert(3)',\n      'let { hilbert } = import(matrix);\\nhilbert(4)',\n    ],\n    seeAlso: ['matrix.vandermonde', 'matrix.symmetric?'],\n  },\n  'vandermonde': {\n    category: 'matrix',\n    description: 'Generates a **Vandermonde matrix** from a vector.',\n    returns: {\n      type: 'matrix',\n    },\n    args: {\n      v: {\n        type: 'vector',\n        description: 'The vector to generate the Vandermonde matrix from.',\n      },\n    },\n    variants: [\n      {\n        argumentNames: [\n          'v',\n        ],\n      },\n    ],\n    examples: [\n      'let { vandermonde } = import(matrix);\\nvandermonde([1, 2, 3])',\n      'let { vandermonde } = import(matrix);\\nvandermonde([1, 0, 1])',\n    ],\n    seeAlso: ['matrix.hilbert', 'matrix.band'],\n  },\n  'band': {\n    category: 'matrix',\n    description: 'Generates a **banded matrix** of size `n` with lower band index `lband` and upper band index `uband`.',\n    returns: {\n      type: 'matrix',\n    },\n    args: {\n      n: {\n        type: 'integer',\n        description: 'The size of the banded matrix.',\n      },\n      lband: {\n        type: 'integer',\n        description: 'The lower band index.',\n      },\n      uband: {\n        type: 'integer',\n        description: 'The upper band index.',\n      },\n    },\n    variants: [\n      {\n        argumentNames: [\n          'n',\n          'lband',\n          'uband',\n        ],\n      },\n    ],\n    examples: [\n      'let { band } = import(matrix);\\nband(3, 1, 1)',\n      'let { band } = import(matrix);\\nband(4, 1, 2)',\n    ],\n    seeAlso: ['matrix.banded?', 'matrix.diagonal?', 'matrix.vandermonde'],\n  },\n  'banded?': {\n    category: 'matrix',\n    description: 'Checks if a `matrix` is **banded** with lower band index `lband` and upper band index `uband`.',\n    returns: {\n      type: 'boolean',\n    },\n    args: {\n      m: {\n        type: 'matrix',\n        description: 'The `matrix` to check for **banded** property.',\n      },\n      lband: {\n        type: 'integer',\n        description: 'The lower band index.',\n      },\n      uband: {\n        type: 'integer',\n        description: 'The upper band index.',\n      },\n    },\n    variants: [\n      {\n        argumentNames: [\n          'm',\n          'lband',\n          'uband',\n        ],\n      },\n    ],\n    examples: [\n      'let { banded? } = import(matrix);\\nbanded?([\\n  [1, 1, 1, 0],\\n  [1, 1, 1, 1],\\n  [1, 1, 1, 1],\\n  [0, 1, 1, 1],\\n], 2, 2)',\n      'let { banded? } = import(matrix);\\nbanded?([\\n  [1, 1, 1, 0],\\n  [1, 1, 1, 1],\\n  [1, 1, 1, 1],\\n  [0, 1, 1, 1],\\n], 1, 1)',\n    ],\n    seeAlso: ['matrix.band', 'matrix.triangular?', 'matrix.diagonal?'],\n  },\n  'rank': {\n    category: 'matrix',\n    description: 'Calculates the **rank** of a matrix using **Gaussian elimination**.',\n    returns: {\n      type: 'number',\n    },\n    args: {\n      m: {\n        type: 'matrix',\n        description: 'The `matrix` to calculate the rank of.',\n      },\n    },\n    variants: [\n      {\n        argumentNames: [\n          'm',\n        ],\n      },\n    ],\n    examples: [\n      'let { rank } = import(matrix);\\nrank([[1, 0, 0], [0, 1, 0], [0, 0, 1]])',\n      'let { rank } = import(matrix);\\nrank([[1, 2, 3], [4, 5, 6], [7, 8, 9]])',\n      'let { rank } = import(matrix);\\nrank([[2, 4, 6], [3, 6, 9], [4, 8, 12]])',\n    ],\n    seeAlso: ['matrix.det', 'matrix.invertible?', 'linear-algebra.rref'],\n  },\n  'frobenius-norm': {\n    category: 'matrix',\n    description: 'Calculates the **Frobenius norm** of a matrix.',\n    returns: {\n      type: 'number',\n    },\n    args: {\n      m: {\n        type: 'matrix',\n        description: 'The `matrix` to calculate the Frobenius norm of.',\n      },\n    },\n    variants: [\n      {\n        argumentNames: [\n          'm',\n        ],\n      },\n    ],\n    examples: [\n      'let { frobenius-norm } = import(matrix);\\nfrobenius-norm([[1, 2], [3, 4]])',\n      'let { frobenius-norm } = import(matrix);\\nfrobenius-norm([[1, 2, 3], [4, 5, 6], [7, 8, 9]])',\n    ],\n    seeAlso: ['matrix.one-norm', 'matrix.inf-norm', 'matrix.max-norm'],\n  },\n  'one-norm': {\n    category: 'matrix',\n    description: 'Calculates the **one-norm** (column norm) of a matrix.',\n    returns: {\n      type: 'number',\n    },\n    args: {\n      m: {\n        type: 'matrix',\n        description: 'The `matrix` to calculate the one-norm of.',\n      },\n    },\n    variants: [\n      {\n        argumentNames: [\n          'm',\n        ],\n      },\n    ],\n    examples: [\n      'let { one-norm } = import(matrix);\\none-norm([[1, 2], [3, 4]])',\n      'let { one-norm } = import(matrix);\\none-norm([[1, 2, 3], [4, 5, 6], [7, 8, 9]])',\n    ],\n    seeAlso: ['matrix.frobenius-norm', 'matrix.inf-norm', 'matrix.max-norm'],\n  },\n  'inf-norm': {\n    category: 'matrix',\n    description: 'Calculates the **infinity norm** of a matrix.',\n    returns: {\n      type: 'number',\n    },\n    args: {\n      m: {\n        type: 'matrix',\n        description: 'The `matrix` to calculate the infinity norm of.',\n      },\n    },\n    variants: [\n      {\n        argumentNames: [\n          'm',\n        ],\n      },\n    ],\n    examples: [\n      'let { inf-norm } = import(matrix);\\ninf-norm([[1, 2], [3, 4]])',\n      'let { inf-norm } = import(matrix);\\ninf-norm([[1, 2, 3], [4, 5, 6], [7, 8, 9]])',\n    ],\n    seeAlso: ['matrix.frobenius-norm', 'matrix.one-norm', 'matrix.max-norm'],\n  },\n  'max-norm': {\n    category: 'matrix',\n    description: 'Calculates the **max norm** of a matrix.',\n    returns: {\n      type: 'number',\n    },\n    args: {\n      m: {\n        type: 'matrix',\n        description: 'The `matrix` to calculate the max norm of.',\n      },\n    },\n    variants: [\n      {\n        argumentNames: [\n          'm',\n        ],\n      },\n    ],\n    examples: [\n      'let { max-norm } = import(matrix);\\nmax-norm([[1, 2], [3, 4]])',\n      'let { max-norm } = import(matrix);\\nmax-norm([[1, 2, 3], [4, 5, 6], [7, 8, 9]])',\n    ],\n    seeAlso: ['matrix.frobenius-norm', 'matrix.one-norm', 'matrix.inf-norm'],\n  },\n}\n","import { approxZero } from '../../../../utils'\n\n/**\n * Calculates the determinant of a matrix using Gaussian Elimination\n * @param matrix A square matrix represented as a 2D array\n * @returns The determinant of the matrix\n */\nexport function determinant(matrix: number[][]): number {\n  // First, make a deep copy of the matrix to avoid modifying the original\n  const n = matrix.length\n  const A: number[][] = []\n\n  for (let i = 0; i < n; i++) {\n    A[i] = [...matrix[i]!]\n  }\n\n  // Handle special cases for small matrices\n  if (n === 1) {\n    return A[0]![0]!\n  }\n\n  if (n === 2) {\n    return A[0]![0]! * A[1]![1]! - A[0]![1]! * A[1]![0]!\n  }\n\n  // For larger matrices, use Gaussian elimination\n  let sign = 1 // Track sign changes from row swaps\n\n  // Perform Gaussian elimination to get an upper triangular matrix\n  for (let i = 0; i < n - 1; i += 1) {\n    // Find pivot (maximum element in current column)\n    let maxRow = i\n\n    for (let j = i + 1; j < n; j += 1) {\n      if (Math.abs(A[j]![i]!) > Math.abs(A[maxRow]![i]!)) {\n        maxRow = j\n      }\n    }\n\n    // If the pivot is zero, the determinant is zero\n    if (approxZero(A[maxRow]![i]!)) {\n      return 0\n    }\n\n    // Swap rows if necessary\n    if (maxRow !== i) {\n      [A[i], A[maxRow]] = [A[maxRow]!, A[i]!] // ES6 array destructuring for swap\n      sign = -sign // Each row swap changes the sign\n    }\n\n    // Eliminate entries below the pivot\n    for (let j = i + 1; j < n; j += 1) {\n      const factor = A[j]![i]! / A[i]![i]!\n\n      // Subtract (factor * pivot row) from current row\n      for (let k = i; k < n; k++) {\n        A[j]![k]! -= factor * A[i]![k]!\n      }\n    }\n  }\n\n  // Calculate determinant as the product of diagonal elements\n  let det = sign\n  for (let i = 0; i < n; i++) {\n    det *= A[i]![i]!\n  }\n\n  return det\n}\n","export function minor(matrix: number[][], row: number, col: number): number[][] {\n  const n = matrix.length\n  const result: number[][] = []\n\n  for (let i = 0; i < n; i++) {\n    if (i !== row) {\n      const minorRow: number[] = []\n      for (let j = 0; j < n; j++) {\n        if (j !== col) {\n          minorRow.push(matrix[i]![j]!)\n        }\n      }\n      result.push(minorRow)\n    }\n  }\n\n  return result\n}\n","import { determinant } from './determinant'\nimport { minor } from './minor'\n\nexport function adjugate(matrix: number[][]): number[][] {\n  const n = matrix.length\n  const adj: number[][] = []\n  for (let i = 0; i < n; i++) {\n    adj[i] = []\n    for (let j = 0; j < n; j++) {\n      const min = minor(matrix, j, i)\n      const sign = (-1) ** (i + j)\n      const cofactor = sign * determinant(min)\n      adj[i]![j] = cofactor\n    }\n  }\n  return adj\n}\n","export function isSquare(matrix: number[][]): boolean {\n  return matrix.length === matrix[0]!.length\n}\n","import { approxEqual, approxZero } from '../../../../utils'\nimport { isSquare } from './isSquare'\n\nexport function isIdentity(matrix: number[][]): boolean {\n  if (!isSquare(matrix)) {\n    return false\n  }\n  const n = matrix.length\n\n  for (let i = 0; i < n; i++) {\n    for (let j = 0; j < n; j++) {\n      if (i === j) {\n        if (!approxEqual(matrix[i]![j]!, 1)) {\n          return false\n        }\n      }\n      else {\n        if (!approxZero(matrix[i]![j]!)) {\n          return false\n        }\n      }\n    }\n  }\n\n  return true\n}\n","/**\n * Performs cache-optimized matrix multiplication.\n * @param A The first input matrix (m x n)\n * @param B The second input matrix (n x p)\n * @returns The result matrix C (m x p) where C = A × B\n */\nexport function matrixMultiply(A: number[][], B: number[][]): number[][] {\n  // Check if matrices can be multiplied\n  if (A.length === 0 || B.length === 0 || A[0]!.length !== B.length) {\n    throw new Error('Matrix dimensions do not match for multiplication')\n  }\n\n  const m = A.length // Number of rows in A\n  const n = A[0]!.length // Number of columns in A / Number of rows in B\n  const p = B[0]!.length // Number of columns in B\n\n  // Initialize result matrix C with zeros\n  const C: number[][] = (Array(m).fill(0) as number[]).map(() => Array(p).fill(0) as number[])\n\n  // Perform multiplication with cache-optimized loop order (i-k-j)\n  for (let i = 0; i < m; i++) {\n    for (let k = 0; k < n; k++) {\n      const aik = A[i]![k]! // Cache this value to avoid repeated lookups\n      for (let j = 0; j < p; j++) {\n        C[i]![j]! += aik * B[k]![j]!\n      }\n    }\n  }\n\n  return C\n}\n","import { DvalaError } from '../../../errors'\nimport { assertMatrix, assertSquareMatrix, assertVector, isSquareMatrix } from '../../../typeGuards/annotatedArrays'\nimport { assertNumber } from '../../../typeGuards/number'\nimport { approxZero } from '../../../utils'\nimport { toFixedArity } from '../../../utils/arity'\nimport type { BuiltinNormalExpressions } from '../../../builtin/interface'\nimport type { DvalaModule } from '../interface'\nimport { gaussJordanElimination } from '../linear-algebra/helpers/gaussJordanElimination'\nimport { moduleDocs } from './docs'\nimport { adjugate } from './helpers/adjugate'\nimport { band } from './helpers/band'\nimport { cofactor } from './helpers/cofactor'\nimport { determinant } from './helpers/determinant'\nimport { inverse } from './helpers/inverse'\nimport { isBanded } from './helpers/isBanded'\nimport { isDiagonal } from './helpers/isDiagonal'\nimport { isIdentity } from './helpers/isIdentity'\nimport { isOrthogonal } from './helpers/isOrthogonal'\nimport { isSquare } from './helpers/isSquare'\nimport { isSymetric } from './helpers/isSymetric'\nimport { isTriangular, isTriangularLower, isTriangularUpper } from './helpers/isTriangular'\nimport { matrixMultiply } from './helpers/matrixMultiply'\nimport { minor } from './helpers/minor'\nimport { norm1 } from './helpers/norm1'\nimport { trace } from './helpers/trace'\n\nexport const matrixNormalExpression: BuiltinNormalExpressions = {\n  'mul': {\n    evaluate: ([matrix1, matrix2], sourceCodeInfo): number[][] => {\n      assertMatrix(matrix1, sourceCodeInfo)\n      assertMatrix(matrix2, sourceCodeInfo)\n      try {\n        return matrixMultiply(matrix1, matrix2)\n      }\n      catch (error) {\n        throw new DvalaError(`The number of columns in the first matrix must be equal to the number of rows in the second matrix, but got ${matrix1[0]!.length} and ${matrix2.length}`, sourceCodeInfo)\n      }\n    },\n    arity: toFixedArity(2),\n  },\n  'det': {\n    evaluate: ([matrix], sourceCodeInfo): number => {\n      assertSquareMatrix(matrix, sourceCodeInfo)\n      return determinant(matrix)\n    },\n    arity: toFixedArity(1),\n  },\n  'inv': {\n    evaluate: ([matrix], sourceCodeInfo): number[][] => {\n      assertSquareMatrix(matrix, sourceCodeInfo)\n      const result = inverse(matrix)\n      if (result === null) {\n        throw new DvalaError('The matrix must be invertible', sourceCodeInfo)\n      }\n      return result\n    },\n    arity: toFixedArity(1),\n  },\n  'adj': {\n    evaluate: ([matrix], sourceCodeInfo): number[][] => {\n      assertSquareMatrix(matrix, sourceCodeInfo)\n      return adjugate(matrix)\n    },\n    arity: toFixedArity(1),\n  },\n  'cofactor': {\n    evaluate: ([matrix], sourceCodeInfo): number[][] => {\n      assertSquareMatrix(matrix, sourceCodeInfo)\n      return cofactor(matrix)\n    },\n    arity: toFixedArity(1),\n  },\n  'minor': {\n    evaluate: ([matrix, row, col], sourceCodeInfo): number[][] => {\n      assertMatrix(matrix, sourceCodeInfo)\n      assertNumber(row, sourceCodeInfo, { integer: true, nonNegative: true, lte: matrix.length })\n      assertNumber(col, sourceCodeInfo, { integer: true, nonNegative: true, lte: matrix[0]!.length })\n\n      return minor(matrix, row, col)\n    },\n    arity: toFixedArity(3),\n  },\n  'trace': {\n    evaluate: ([matrix], sourceCodeInfo): number => {\n      assertSquareMatrix(matrix, sourceCodeInfo)\n      return trace(matrix)\n    },\n    arity: toFixedArity(1),\n  },\n  'symmetric?': {\n    evaluate: ([matrix], sourceCodeInfo): boolean => {\n      assertMatrix(matrix, sourceCodeInfo)\n      return isSymetric(matrix)\n    },\n    arity: toFixedArity(1),\n  },\n  'triangular?': {\n    evaluate: ([matrix], sourceCodeInfo): boolean => {\n      assertMatrix(matrix, sourceCodeInfo)\n      return isTriangular(matrix)\n    },\n    arity: toFixedArity(1),\n  },\n  'upper-triangular?': {\n    evaluate: ([matrix], sourceCodeInfo): boolean => {\n      assertMatrix(matrix, sourceCodeInfo)\n      return isTriangularUpper(matrix)\n    },\n    arity: toFixedArity(1),\n  },\n  'lower-triangular?': {\n    evaluate: ([matrix], sourceCodeInfo): boolean => {\n      assertMatrix(matrix, sourceCodeInfo)\n      return isTriangularLower(matrix)\n    },\n    arity: toFixedArity(1),\n  },\n  'diagonal?': {\n    evaluate: ([matrix], sourceCodeInfo): boolean => {\n      assertMatrix(matrix, sourceCodeInfo)\n      return isDiagonal(matrix)\n    },\n    arity: toFixedArity(1),\n  },\n  'square?': {\n    evaluate: ([matrix], sourceCodeInfo): boolean => {\n      assertMatrix(matrix, sourceCodeInfo)\n      return isSquare(matrix)\n    },\n    arity: toFixedArity(1),\n  },\n  'orthogonal-matrix?': {\n    evaluate: ([matrix], sourceCodeInfo): boolean => {\n      assertMatrix(matrix, sourceCodeInfo)\n      return isOrthogonal(matrix)\n    },\n    arity: toFixedArity(1),\n  },\n  'identity?': {\n    evaluate: ([matrix], sourceCodeInfo): boolean => {\n      assertMatrix(matrix, sourceCodeInfo)\n      return isIdentity(matrix)\n    },\n    arity: toFixedArity(1),\n  },\n  'invertible?': {\n    evaluate: ([matrix], sourceCodeInfo): boolean => {\n      assertMatrix(matrix, sourceCodeInfo)\n      if (!isSquareMatrix(matrix)) {\n        return false\n      }\n      return !approxZero(determinant(matrix))\n    },\n    arity: toFixedArity(1),\n  },\n  'hilbert': {\n    evaluate: ([size], sourceCodeInfo): number[][] => {\n      assertNumber(size, sourceCodeInfo, { integer: true, positive: true })\n      const result: number[][] = []\n      for (let i = 0; i < size; i += 1) {\n        const row: number[] = []\n        for (let j = 0; j < size; j += 1) {\n          row.push(1 / (i + j + 1))\n        }\n        result.push(row)\n      }\n      return result\n    },\n    arity: toFixedArity(1),\n  },\n  'vandermonde': {\n    evaluate: ([vector], sourceCodeInfo): number[][] => {\n      assertVector(vector, sourceCodeInfo)\n      const result: number[][] = []\n      for (let i = 0; i < vector.length; i += 1) {\n        const row: number[] = []\n        for (let j = 0; j < vector.length; j += 1) {\n          row.push((vector[i]!) ** j)\n        }\n        result.push(row)\n      }\n      return result\n    },\n    arity: toFixedArity(1),\n  },\n  'band': {\n    evaluate: ([n, lband, uband], sourceCodeInfo): number[][] => {\n      assertNumber(n, sourceCodeInfo, { integer: true, positive: true })\n      assertNumber(lband, sourceCodeInfo, { integer: true, nonNegative: true, lt: n })\n      assertNumber(uband, sourceCodeInfo, { integer: true, nonNegative: true, lte: n })\n      return band(n, lband, uband)\n    },\n    arity: toFixedArity(3),\n  },\n  'banded?': {\n    evaluate: ([matrix, lband, uband], sourceCodeInfo): boolean => {\n      assertMatrix(matrix, sourceCodeInfo)\n      const maxBand = Math.max(matrix.length, matrix[0]!.length)\n      assertNumber(lband, sourceCodeInfo, { integer: true, nonNegative: true, lt: maxBand })\n      assertNumber(uband, sourceCodeInfo, { integer: true, nonNegative: true, lt: maxBand })\n      return isBanded(matrix, lband, uband)\n    },\n    arity: toFixedArity(3),\n  },\n  'rank': {\n    evaluate: ([matrix], sourceCodeInfo): number => {\n      assertMatrix(matrix, sourceCodeInfo)\n      const [, result] = gaussJordanElimination(matrix)\n      return result\n    },\n    arity: toFixedArity(1),\n  },\n  // Frobenius norm\n  'frobenius-norm': {\n    evaluate: ([matrix], sourceCodeInfo): number => {\n      assertMatrix(matrix, sourceCodeInfo)\n      return Math.sqrt(matrix.reduce((sum, row) => sum + row.reduce((rowSum, cell) => rowSum + cell * cell, 0), 0))\n    },\n    arity: toFixedArity(1),\n  },\n  // one-norm (column norm)\n  'one-norm': {\n    evaluate: ([matrix], sourceCodeInfo): number => {\n      assertMatrix(matrix, sourceCodeInfo)\n      return norm1(matrix)\n    },\n    arity: toFixedArity(1),\n  },\n  // Infinity norm\n  'inf-norm': {\n    evaluate: ([matrix], sourceCodeInfo): number => {\n      assertMatrix(matrix, sourceCodeInfo)\n      return matrix.reduce((max, row) => Math.max(max, row.reduce((sum, cell) => sum + Math.abs(cell), 0)), 0)\n    },\n    arity: toFixedArity(1),\n  },\n  // Max norm\n  'max-norm': {\n    evaluate: ([matrix], sourceCodeInfo): number => {\n      assertMatrix(matrix, sourceCodeInfo)\n      return matrix.reduce((maxVal, row) => {\n        const rowMax = row.reduce((max, val) => Math.max(max, Math.abs(val)), 0)\n        return Math.max(maxVal, rowMax)\n      }, 0)\n    },\n    arity: toFixedArity(1),\n  },\n}\n\nfor (const [key, docs] of Object.entries(moduleDocs)) {\n  if (matrixNormalExpression[key])\n    matrixNormalExpression[key].docs = docs\n}\n\nexport const matrixModule: DvalaModule = {\n  name: 'matrix',\n  functions: matrixNormalExpression,\n}\n","/**\n * Calculate the inverse of a matrix using the adjugate method\n * @param matrix The input matrix\n * @returns The inverse matrix or null if the matrix is not invertible\n */\n\nimport { approxZero } from '../../../../utils'\nimport { adjugate } from './adjugate'\nimport { determinant } from './determinant'\n\nexport function inverse(matrix: number[][]): number[][] | null {\n  const n = matrix.length\n\n  // Special case for 1x1 matrix - handle it directly\n  if (n === 1) {\n    const element = matrix[0]![0]!\n    if (approxZero(element)) {\n      return null // Not invertible\n    }\n    return [[1 / element]]\n  }\n\n  // Calculate determinant\n  const det = determinant(matrix)\n\n  // Check if matrix is invertible\n  if (approxZero(det)) {\n    return null // Matrix is not invertible\n  }\n\n  // Get the adjugate matrix\n  const adj = adjugate(matrix)\n\n  // Calculate the inverse: inverse = adjugate / determinant\n  const inverseMatrix: number[][] = []\n  for (let i = 0; i < n; i++) {\n    inverseMatrix[i] = []\n    for (let j = 0; j < n; j++) {\n      inverseMatrix[i]![j] = adj[i]![j]! / det\n    }\n  }\n\n  return inverseMatrix\n}\n","import { determinant } from './determinant'\nimport { minor } from './minor'\n\nexport function cofactor(matrix: number[][]): number[][] {\n  const n = matrix.length\n  const cofactors: number[][] = []\n\n  // Create a new matrix to store cofactors\n  for (let i = 0; i < n; i++) {\n    cofactors[i] = []\n    for (let j = 0; j < n; j++) {\n      // Get the minor by removing row i and column j\n      const min = minor(matrix, i, j)\n      const sign = (-1) ** (i + j)\n      cofactors[i]![j] = sign * determinant(min)\n    }\n  }\n\n  return cofactors\n}\n","/**\n * Calculates the trace of a square matrix.\n * The trace is defined as the sum of the elements on the main diagonal.\n *\n * @param matrix - A 2D array representing a square matrix.\n * @returns The trace of the matrix.\n */\nexport function trace(matrix: number[][]): number {\n  return matrix.reduce((sum, row, i) => sum + row[i]!, 0)\n}\n","import { isSquare } from './isSquare'\n\n/**\n * Checks if a given matrix is symmetric.\n * A matrix is symmetric if it is square and its transpose is equal to itself.\n *\n * @param matrix - A 2D array representing the matrix.\n * @returns `true` if the matrix is symmetric, otherwise `false`.\n */\nexport function isSymetric(matrix: number[][]): boolean {\n  const rows = matrix.length\n\n  // Check if the matrix is square\n  if (!isSquare(matrix)) {\n    return false\n  }\n\n  // Check symmetry\n  for (let i = 0; i < rows; i += 1) {\n    for (let j = 0; j < i; j += 1) {\n      if (matrix[i]![j]! !== matrix[j]![i]!) {\n        return false\n      }\n    }\n  }\n\n  return true\n}\n","import { isSquare } from './isSquare'\n\n/**\n * Determines whether a given matrix is triangular.\n *\n * A triangular matrix is a square matrix where all elements\n * below or above the main diagonal are zero. This function\n * checks if the matrix is square and symmetric.\n *\n * @param matrix - A two-dimensional array of numbers representing the matrix.\n * @returns `true` if the matrix is triangular, otherwise `false`.\n */\nexport function isTriangular(matrix: number[][]): boolean {\n  if (!isSquare(matrix)) {\n    return false\n  }\n\n  const rows = matrix.length\n\n  let isUpperTriangular = true\n  let isLowerTriangular = true\n\n  for (let i = 0; i < rows; i++) {\n    for (let j = 0; j < rows; j++) {\n      if (i > j && matrix[i]![j] !== 0) {\n        isUpperTriangular = false\n        if (!isLowerTriangular) {\n          return false\n        }\n      }\n      if (i < j && matrix[i]![j] !== 0) {\n        isLowerTriangular = false\n      }\n    }\n  }\n\n  return isUpperTriangular || isLowerTriangular\n}\n\nexport function isTriangularUpper(matrix: number[][]): boolean {\n  if (!isSquare(matrix)) {\n    return false\n  }\n  const rows = matrix.length\n\n  for (let i = 0; i < rows; i++) {\n    for (let j = 0; j < i; j++) {\n      if (matrix[i]![j] !== 0) {\n        return false\n      }\n    }\n  }\n\n  return true\n}\n\nexport function isTriangularLower(matrix: number[][]): boolean {\n  if (!isSquare(matrix)) {\n    return false\n  }\n  const rows = matrix.length\n\n  // Check if the matrix is square\n  if (!matrix.every(row => row.length === rows)) {\n    return false\n  }\n\n  for (let i = 0; i < rows; i++) {\n    for (let j = i + 1; j < rows; j++) {\n      if (matrix[i]![j] !== 0) {\n        return false\n      }\n    }\n  }\n\n  return true\n}\n","import { isSquare } from './isSquare'\n\n/**\n * Checks if a given matrix is diagonal.\n *\n * A matrix is considered diagonal if it is square (i.e., the number of rows equals the number of columns)\n * and all elements outside the main diagonal are zero.\n *\n * @param matrix - A two-dimensional array of numbers representing the matrix to check.\n * @returns `true` if the matrix is diagonal, otherwise `false`.\n */\nexport function isDiagonal(matrix: number[][]): boolean {\n  if (!isSquare(matrix)) {\n    return false\n  }\n  const rows = matrix.length\n\n  for (let i = 0; i < rows; i += 1) {\n    for (let j = 0; j < rows; j += 1) {\n      if (i !== j && matrix[i]![j] !== 0) {\n        return false\n      }\n    }\n  }\n\n  return true\n}\n","import { transpose } from '../../../modules/grid/transpose'\nimport { matrixMultiply } from './matrixMultiply'\nimport { isIdentity } from './isIdentity'\nimport { isSquare } from './isSquare'\n\nexport function isOrthogonal(matrix: number[][]): boolean {\n  if (!isSquare(matrix)) {\n    return false\n  }\n\n  // Calculate matrix transpose\n  const transposed = transpose(matrix)\n\n  // Check if matrix * transpose = Identity\n  const product = matrixMultiply(matrix, transposed)\n\n  // Check if the product is an identity matrix\n  return isIdentity(product)\n}\n","import type { Any } from '../../../interface'\n\nexport function transpose<T extends Any>(grid: T[][]): T[][] {\n  const result: T[][] = []\n\n  for (let i = 0; i < grid[0]!.length; i += 1) {\n    const row: T[] = []\n    for (let j = 0; j < grid.length; j += 1) {\n      row.push(grid[j]![i]!)\n    }\n    result.push(row)\n  }\n  return result\n}\n","/**\n * Creates a band matrix with specified lower and upper bandwidths\n *\n * @param n Size of the square matrix\n * @param lband Lower bandwidth (number of non-zero diagonals below main diagonal)\n * @param uband Upper bandwidth (number of non-zero diagonals above main diagonal)\n * @returns A 2D array representing the band matrix with 1s in the band and 0s elsewhere\n */\nexport function band(n: number, lband: number, uband: number): number[][] {\n  // Create an n×n matrix filled with zeros\n  const matrix: number[][] = Array.from({ length: n }, () => Array.from({ length: n }, () => 0))\n\n  // Fill the band with 1s\n  for (let i = 0; i < n; i++) {\n    for (let j = Math.max(0, i - lband); j <= Math.min(n - 1, i + uband); j++) {\n      matrix[i]![j] = 1\n    }\n  }\n\n  return matrix\n}\n","/**\n * Checks if a matrix is banded with the given lower and upper bandwidth.\n * A matrix is banded if all non-zero elements are within 'lower' diagonals\n * below the main diagonal and 'upper' diagonals above the main diagonal.\n *\n * @param matrix - The matrix to check, represented as a 2D array of numbers\n * @param lower - Number of non-zero diagonals below the main diagonal\n * @param upper - Number of non-zero diagonals above the main diagonal\n * @returns true if the matrix is banded with the given parameters, false otherwise\n */\nexport function isBanded(matrix: number[][], lower: number, upper: number): boolean {\n  const rows = matrix.length\n  const cols = matrix[0]!.length\n\n  // Check each element in the matrix\n  for (let i = 0; i < rows; i++) {\n    for (let j = 0; j < cols; j++) {\n      // If we find a non-zero element outside the band, return false\n      if (matrix[i]![j] !== 0 && (i - j > lower || j - i > upper)) {\n        return false\n      }\n    }\n  }\n\n  // All elements outside the band are zero\n  return true\n}\n","import { approxZero } from '../../../../utils'\n\n/**\n * Performs Gauss-Jordan elimination on a matrix, transforming it to reduced row echelon form\n *\n * @param matrix - The input matrix\n * @returns A tuple containing the reduced row echelon form matrix and the rank\n */\nexport function gaussJordanElimination(matrix: number[][]): [number[][], number] {\n  // Create a copy of the matrix to avoid modifying the original\n  const m = matrix.map(row => [...row])\n  const rows = m.length\n  const cols = m[0]!.length\n\n  let rank = 0\n  let rowsProcessed = 0\n\n  // Row reduction to reduced row echelon form\n  for (let col = 0; col < cols; col++) {\n    // Find the pivot\n    let pivotRow = -1\n\n    for (let row = rowsProcessed; row < rows; row++) {\n      if (!approxZero(m[row]![col]!)) {\n        pivotRow = row\n        break\n      }\n    }\n\n    if (pivotRow === -1)\n      continue // No pivot in this column\n\n    // Increase rank\n    rank += 1\n\n    // Swap rows\n    if (pivotRow !== rowsProcessed) {\n      [m[pivotRow], m[rowsProcessed]] = [m[rowsProcessed]!, m[pivotRow]!]\n    }\n\n    // Get the pivot value\n    const pivotValue = m[rowsProcessed]![col]!\n\n    // Normalize the pivot row (always, for RREF)\n    for (let j = col; j < cols; j++) {\n      m[rowsProcessed]![j]! /= pivotValue\n    }\n\n    // Eliminate above and below (full Gauss-Jordan)\n    for (let row = 0; row < rows; row++) {\n      if (row !== rowsProcessed && !approxZero(m[row]![col]!)) {\n        const factor = m[row]![col]!\n        for (let j = col; j < cols; j++) {\n          m[row]![j]! -= factor * m[rowsProcessed]![j]!\n        }\n      }\n    }\n\n    rowsProcessed++\n    if (rowsProcessed === rows)\n      break\n  }\n\n  return [m, rank]\n}\n","// Assuming a matrix is represented as a 2D array\nexport function norm1(matrix: number[][]): number {\n  const numRows = matrix.length\n  const numCols = matrix[0]!.length\n\n  let maxColSum = 0\n\n  // Iterate through each column\n  for (let j = 0; j < numCols; j += 1) {\n    let colSum = 0\n\n    // Sum the absolute values of all elements in this column\n    for (let i = 0; i < numRows; i += 1) {\n      colSum += Math.abs(matrix[i]![j]!)\n    }\n\n    // Update the maximum column sum if necessary\n    maxColSum = Math.max(maxColSum, colSum)\n  }\n\n  return maxColSum\n}\n"],"names":["getCodeMarker","sourceCodeInfo","position","code","leftPadding","column","rightPadding","length","repeat","Math","max","DvalaError","Error","shortMessage","constructor","err","message","super","location","line","filePath","getDvalaErrorMessage","this","Object","setPrototypeOf","prototype","name","NodeTypes","Number","String","NormalExpression","SpecialExpression","UserDefinedSymbol","NormalBuiltinSymbol","SpecialBuiltinSymbol","ReservedSymbol","Binding","Spread","NodeTypesSet","Set","values","functionTypeSet","isDvalaFunction","func","type","functionType","has","isNode","value","Array","isArray","valueToString","keys","find","key","RegExp","toString","JSON","stringify","isNumber","options","isNaN","integer","isInteger","finite","isFinite","zero","nonZero","positive","negative","nonPositive","nonNegative","gt","gte","lt","lte","assertNumber","sign","getSignString","numberType","range","getRangeString","filter","x","join","getNumberTypeName","anyValue","getSourceCodeInfo","annotatedArrays","WeakSet","vectors","notVectors","matrices","notMatrices","grids","assertVector","vector","every","elem","add","isVector","isMatrix","matrix","nbrOfCols","row","some","cell","assertMatrix","assertSquareMatrix","approxEqual","a","b","epsilon","diff","abs","approxZero","toFixedArity","arity","min","moduleDocs","mul","category","description","returns","args","variants","argumentNames","examples","seeAlso","det","m","inv","adj","cofactor","minor","col","trace","hilbert","n","vandermonde","v","band","lband","uband","rank","determinant","A","i","maxRow","j","factor","k","result","minorRow","push","adjugate","isSquare","isIdentity","matrixMultiply","B","p","C","fill","map","aik","matrixNormalExpression","evaluate","matrix1","matrix2","error","element","inverseMatrix","inverse","cofactors","reduce","sum","rows","isSymetric","isUpperTriangular","isLowerTriangular","isTriangular","isTriangularUpper","isTriangularLower","isDiagonal","grid","transpose","isOrthogonal","isSquareMatrix","size","from","maxBand","lower","upper","cols","isBanded","rowsProcessed","pivotRow","pivotValue","gaussJordanElimination","sqrt","rowSum","numRows","numCols","maxColSum","colSum","norm1","maxVal","rowMax","val","docs","entries","matrixModule","functions"],"mappings":"AAEM,SAAUA,EAAcC,GAC5B,IAAKA,EAAeC,WAAaD,EAAeE,KAC9C,MAAO,GAET,MAAMC,EAAcH,EAAeC,SAASG,OAAS,EAC/CC,EAAeL,EAAeE,KAAKI,OAASH,EAAc,EAChE,MAAO,GAAG,IAAII,OAAOC,KAAKC,IAAIN,EAAa,OAAO,IAAII,OAAOC,KAAKC,IAAIJ,EAAc,KACtF,CCkBM,MAAOK,UAAmBC,MACdX,eACAY,aAChB,WAAAC,CAAYC,EAAcd,GACxB,MAAMe,EAAUD,aAAeH,MAC3BG,EAAIC,QACJ,GAAGD,IAEPE,MA/BJ,SAA8BD,EAAiBf,GAC7C,IAAKA,EACH,OAAOe,EAET,MAAME,EAAW,GAAGjB,EAAeC,SAASiB,QAAQlB,EAAeC,SAASG,SAM5E,MAAO,GAAGW,IALWf,EAAemB,SAChC,KAAKnB,EAAemB,YAAYF,IAChC,cAAcA,QACIjB,EAAeE,SACbH,EAAcC,IAExC,CAoBUoB,CAAqBL,EAASf,IACpCqB,KAAKT,aAAeG,EACpBM,KAAKrB,eAAiBA,EACtBsB,OAAOC,eAAeF,KAAMX,EAAWc,WACvCH,KAAKI,KAAO,YACb,CAEM,aAAA1B,GACL,OAAOsB,KAAKrB,gBAAkBD,EAAcsB,KAAKrB,eAClD,EC5CI,MAAM0B,EAAY,CACvBC,OAAQ,EACRC,OAAQ,EACRC,iBAAkB,EAClBC,kBAAmB,EACnBC,kBAAmB,EACnBC,oBAAqB,EACrBC,qBAAsB,EACtBC,eAAgB,EAChBC,QAAS,EACTC,OAAQ,IAGJC,EAAe,IAAIC,IAAIhB,OAAOiB,OAAOb,IAa3C,MAgBMc,EAAkB,IAAIF,IAhBN,CACpB,cACA,UACA,OACA,aACA,OACA,aACA,YACA,WACA,QACA,UACA,iBACA,mBACA,WCnCF,SAASG,EAAgBC,GACvB,OAAa,OAATA,GAAiC,iBAATA,ICLC,WDQHA,GAAQ,iBAAkBA,IDuC7B,iBADMC,ECtC8CD,EAAKE,eDuC7CJ,EAAgBK,IAAIF,KADnD,IAAyBA,CCrC/B,CAEA,SAASG,EAAOC,GACd,SAAKC,MAAMC,QAAQF,IAAUA,EAAMzC,OAAS,KDWrB,iBADEqC,ECRPI,EAAM,KDSWV,EAAaQ,IAAIF,IADhD,IAAqBA,CCP3B,CAEM,SAAUO,EAAcH,GAC5B,OAAIN,EAAgBM,GAEX,aAAcA,EAActB,MAAQ,OAEzCqB,EAAOC,GACF,GDNqBJ,ECMFI,EAAM,GDL3BzB,OAAO6B,KAAKzB,GAAW0B,KAAKC,GAAO3B,EAAU2B,KAAmCV,UCOzE,OAAVI,EACK,OAEY,iBAAVA,GAAsBA,aAAiBO,OACzC,GAAGP,IAES,iBAAVA,GAAsBA,aAAiBpC,MACzCoC,EAAMQ,WAERC,KAAKC,UAAUV,GDjBlB,IAA0BJ,CCkBhC,UE2FgBe,EAASX,EAAgBY,EAAyB,IAChE,MAAqB,iBAAVZ,KAGPpB,OAAOiC,MAAMb,OAGbY,EAAQE,UAAYlC,OAAOmC,UAAUf,QAGrCY,EAAQI,SAAWpC,OAAOqC,SAASjB,QAGnCY,EAAQM,MAAkB,IAAVlB,OAGhBY,EAAQO,SAAqB,IAAVnB,OAGnBY,EAAQQ,UAAYpB,GAAS,OAG7BY,EAAQS,UAAYrB,GAAS,OAG7BY,EAAQU,aAAetB,EAAQ,OAG/BY,EAAQW,aAAevB,EAAQ,OAGT,iBAAfY,EAAQY,IAAmBxB,GAASY,EAAQY,QAG5B,iBAAhBZ,EAAQa,KAAoBzB,EAAQY,EAAQa,SAG7B,iBAAfb,EAAQc,IAAmB1B,GAASY,EAAQc,OAG5B,iBAAhBd,EAAQe,KAAoB3B,EAAQY,EAAQe,gBAIzD,CAEM,SAAUC,EACd5B,EACA/C,EACA2D,EAAyB,CAAA,GAEzB,IAAKD,EAASX,EAAOY,GACnB,MAAM,IAAIjD,EACR,YAjEN,SAA2BiD,GACzB,GAAIA,EAAQM,KACV,MAAO,OAET,MAAMW,EAlBR,SAAuBjB,GACrB,OAAOA,EAAQQ,SACX,WACAR,EAAQS,SACN,WACAT,EAAQW,YACN,eACAX,EAAQU,YACN,eACAV,EAAQO,QACN,WACA,EACd,CAMeW,CAAclB,GACrBmB,EAAanB,EAAQE,QAAU,UAAY,SAC3CE,EAASJ,EAAQI,OAAS,SAAW,GACrCgB,EAvCR,SAAwBpB,GAItB,MAF2B,iBAAfA,EAAQY,IAA0C,iBAAhBZ,EAAQa,KACxB,iBAAfb,EAAQc,IAA0C,iBAAhBd,EAAQe,IAM1B,iBAAff,EAAQY,IAA0C,iBAAhBZ,EAAQa,IACxB,iBAAfb,EAAQY,GAAkB,OAAOZ,EAAQY,KAAO,QAAQZ,EAAQa,MAEpD,iBAAfb,EAAQc,IAA0C,iBAAhBd,EAAQe,IACxB,iBAAff,EAAQc,GAAkB,OAAOd,EAAQc,KAAO,QAAQd,EAAQe,MAErE,GAVL,GAAyB,iBAAff,EAAQY,GAAkB,GAAGZ,EAAQY,UAAY,GAAGZ,EAAQa,cACrD,iBAAfb,EAAQc,GAAkB,KAAKd,EAAQc,KAAO,MAAMd,EAAQe,OAUzE,CAuBgBM,CAAerB,GAE7B,MAAO,CAACiB,EAAMb,EAAQe,EAAYC,GAAOE,OAAOC,KAAOA,GAAGC,KAAK,IACjE,CAuDkBC,CAAkBzB,WAAiBT,EAAcH,MCjLnD,SAAkBsC,EAAerF,GAE/C,OAAOqF,GAAUrF,gBAAkBA,CACrC,CD+KMsF,CAAkBvC,EAAO/C,GAG/B,CElLA,MAAMuF,EAAkB,IAAIC,QACtBC,EAAU,IAAID,QACdE,EAAa,IAAIF,QACjBG,EAAW,IAAIH,QACfI,EAAc,IAAIJ,QAClBK,EAAQ,IAAIL,QAsCF,SAAAM,EAAaC,EAAiB/F,GAC5C,IAtBI,SAAmB+F,GACvB,SAAK/C,MAAMC,QAAQ8C,KAIfN,EAAQ5C,IAAIkD,KAGZL,EAAW7C,IAAIkD,KAIfA,EAAOC,MAAMC,GAAQvC,EAASuC,KAChCV,EAAgBW,IAAIH,GACpBN,EAAQS,IAAIH,GACL,IAETL,EAAWQ,IAAIH,GACR,KACT,CAGOI,CAASJ,GACZ,MAAM,IAAIrF,EAAW,8BAA8BqF,IAAU/F,EAEjE,CAyEM,SAAUoG,EAASC,GACvB,IAAKrD,MAAMC,QAAQoD,GACjB,OAAO,EAET,GAAIV,EAAS9C,IAAIwD,GACf,OAAO,EAET,GAAIT,EAAY/C,IAAIwD,GAClB,OAAO,EAET,GAAsB,IAAlBA,EAAO/F,OAET,OADAsF,EAAYM,IAAIG,IACT,EAET,IAAKrD,MAAMC,QAAQoD,EAAO,KAA4B,IAArBA,EAAO,GAAG/F,OAEzC,OADAsF,EAAYM,IAAIG,IACT,EAET,MAAMC,EAAYD,EAAO,GAAG/F,OAC5B,IAAK,MAAMiG,KAAOF,EAChB,IAAKrD,MAAMC,QAAQsD,IAAQA,EAAIjG,SAAWgG,GAAaC,EAAIC,KAAKC,IAAS/C,EAAS+C,IAEhF,OADAb,EAAYM,IAAIG,IACT,EAMX,OAHAd,EAAgBW,IAAIG,GACpBR,EAAMK,IAAIG,GACVV,EAASO,IAAIG,IACN,CACT,CAEgB,SAAAK,EAAaL,EAAiBrG,GAC5C,IAAKoG,EAASC,GACZ,MAAM,IAAI3F,EAAW,8BAA8B2F,IAAUrG,EAEjE,CAEgB,SAAA2G,EAAmBN,EAAiBrG,GAClD,IAAKoG,EAASC,GACZ,MAAM,IAAI3F,EAAW,8BAA8B2F,IAAUrG,GAE/D,GAAIqG,EAAO/F,SAAW+F,EAAO,GAAI/F,OAC/B,MAAM,IAAII,EAAW,mCAAmC2F,EAAO/F,cAAc+F,EAAO,GAAI/F,SAAUN,EAEtG,CC1DM,SAAU4G,EAAYC,EAAWC,EAAWC,EAF3B,OAGrB,GAAIF,IAAMC,EACR,OAAO,EAGT,MAAME,EAAOxG,KAAKyG,IAAIJ,EAAIC,GAE1B,GAAU,IAAND,GAAiB,IAANC,GAAWE,EAAOD,EAE/B,OAAOC,EAAOD,EAMhB,OAAOC,GAJMxG,KAAKyG,IAAIJ,GACTrG,KAAKyG,IAAIH,IAGQC,CAChC,CAEM,SAAUG,EAAWnE,GACzB,OAAOvC,KAAKyG,IAAIlE,GArBK,KAsBvB,CCjDM,SAAUoE,EAAaC,GAC3B,MAAO,CAAEC,IAAKD,EAAO3G,IAAK2G,EAC5B,CClFO,MAAME,EAA2C,CACtDC,IAAO,CACLC,SAAU,SACVC,YAAa,kHACbC,QAAS,CACP/E,KAAM,UAERgF,KAAM,CACJd,EAAG,CACDlE,KAAM,UAERmE,EAAG,CACDnE,KAAM,WAGViF,SAAU,CACR,CACEC,cAAe,CACb,IACA,OAINC,SAAU,CACR,yEACA,2FAEFC,QAAS,CAAC,aAAc,eAE1BC,IAAO,CACLR,SAAU,SACVC,YAAa,qDACbC,QAAS,CACP/E,KAAM,UAERgF,KAAM,CACJM,EAAG,CACDtF,KAAM,SACN8E,YAAa,kDAGjBG,SAAU,CACR,CACEC,cAAe,CACb,OAINC,SAAU,CACR,uDACA,yEAEFC,QAAS,CAAC,aAAc,kBAAmB,aAAc,eAAgB,cAAe,qBAAsB,aAAc,iBAE9HG,IAAO,CACLV,SAAU,SACVC,YAAa,iDACbC,QAAS,CACP/E,KAAM,UAERgF,KAAM,CACJM,EAAG,CACDtF,KAAM,SACN8E,YAAa,8CAGjBG,SAAU,CACR,CACEC,cAAe,CACb,OAINC,SAAU,CACR,uDACA,0EAEFC,QAAS,CAAC,aAAc,aAAc,qBAAsB,uBAAwB,aAAc,8BAEpGI,IAAO,CACLX,SAAU,SACVC,YAAa,kDACbC,QAAS,CACP/E,KAAM,UAERgF,KAAM,CACJM,EAAG,CACDtF,KAAM,SACN8E,YAAa,+CAGjBG,SAAU,CACR,CACEC,cAAe,CACb,OAINC,SAAU,CACR,uDACA,wEACA,yEAEFC,QAAS,CAAC,kBAAmB,aAAc,eAE7CK,SAAY,CACVZ,SAAU,SACVC,YAAa,kDACbC,QAAS,CACP/E,KAAM,UAERgF,KAAM,CACJM,EAAG,CACDtF,KAAM,SACN8E,YAAa,+CAGjBG,SAAU,CACR,CACEC,cAAe,CACb,OAINC,SAAU,CACR,iEACA,kFACA,mFAEFC,QAAS,CAAC,aAAc,eAAgB,eAE1CM,MAAS,CACPb,SAAU,SACVC,YAAa,+CACbC,QAAS,CACP/E,KAAM,UAERgF,KAAM,CACJM,EAAG,CACDtF,KAAM,SACN8E,YAAa,2CAEflB,IAAK,CACH5D,KAAM,UACN8E,YAAa,4DAEfa,IAAK,CACH3F,KAAM,UACN8E,YAAa,gEAGjBG,SAAU,CACR,CACEC,cAAe,CACb,IACA,MACA,SAINC,SAAU,CACR,iEACA,mFAEFC,QAAS,CAAC,kBAAmB,eAE/BQ,MAAS,CACPf,SAAU,SACVC,YAAa,+CACbC,QAAS,CACP/E,KAAM,UAERgF,KAAM,CACJM,EAAG,CACDtF,KAAM,SACN8E,YAAa,4CAGjBG,SAAU,CACR,CACEC,cAAe,CACb,OAINC,SAAU,CACR,2DACA,6EAEFC,QAAS,CAAC,aAAc,qBAE1B,aAAc,CACZP,SAAU,SACVC,YAAa,yCACbC,QAAS,CACP/E,KAAM,WAERgF,KAAM,CACJM,EAAG,CACDtF,KAAM,SACN8E,YAAa,wCAGjBG,SAAU,CACR,CACEC,cAAe,CACb,OAINC,SAAU,CACR,qEACA,uFAEFC,QAAS,CAAC,4BAA6B,mBAAoB,iBAAkB,mBAE/E,cAAe,CACbP,SAAU,SACVC,YAAa,0CACbC,QAAS,CACP/E,KAAM,WAERgF,KAAM,CACJM,EAAG,CACDtF,KAAM,SACN8E,YAAa,6CAGjBG,SAAU,CACR,CACEC,cAAe,CACb,OAINC,SAAU,CACR,uEACA,yFAEFC,QAAS,CAAC,2BAA4B,2BAA4B,mBAAoB,mBAExF,oBAAqB,CACnBP,SAAU,SACVC,YAAa,gDACbC,QAAS,CACP/E,KAAM,WAERgF,KAAM,CACJM,EAAG,CACDtF,KAAM,SACN8E,YAAa,mDAGjBG,SAAU,CACR,CACEC,cAAe,CACb,OAINC,SAAU,CACR,mFACA,qGAEFC,QAAS,CAAC,2BAA4B,qBAAsB,qBAE9D,oBAAqB,CACnBP,SAAU,SACVC,YAAa,gDACbC,QAAS,CACP/E,KAAM,WAERgF,KAAM,CACJM,EAAG,CACDtF,KAAM,SACN8E,YAAa,mDAGjBG,SAAU,CACR,CACEC,cAAe,CACb,OAINC,SAAU,CACR,mFACA,qGAEFC,QAAS,CAAC,2BAA4B,qBAAsB,qBAE9D,YAAa,CACXP,SAAU,SACVC,YAAa,wCACbC,QAAS,CACP/E,KAAM,WAERgF,KAAM,CACJM,EAAG,CACDtF,KAAM,SACN8E,YAAa,iDAGjBG,SAAU,CACR,CACEC,cAAe,CACb,OAINC,SAAU,CACR,mEACA,oFACA,qFAEFC,QAAS,CAAC,mBAAoB,oBAAqB,qBAAsB,eAAgB,2BAA4B,2BAA4B,cAAe,mBAElK,UAAW,CACTP,SAAU,SACVC,YAAa,sCACbC,QAAS,CACP/E,KAAM,WAERgF,KAAM,CACJM,EAAG,CACDtF,KAAM,SACN8E,YAAa,+CAGjBG,SAAU,CACR,CACEC,cAAe,CACb,OAINC,SAAU,CACR,+DACA,gFACA,sEAEFC,QAAS,CAAC,oBAAqB,mBAAoB,uBAErD,qBAAsB,CACpBP,SAAU,SACVC,YAAa,0CACbC,QAAS,CACP/E,KAAM,WAERgF,KAAM,CACJM,EAAG,CACDtF,KAAM,SACN8E,YAAa,iDAGjBG,SAAU,CACR,CACEC,cAAe,CACb,OAINC,SAAU,CACR,qFACA,sFACA,sFAEFC,QAAS,CAAC,oBAAqB,aAAc,mBAAoB,+BAEnE,YAAa,CACXP,SAAU,SACVC,YAAa,kDACbC,QAAS,CACP/E,KAAM,WAERgF,KAAM,CACJM,EAAG,CACDtF,KAAM,SACN8E,YAAa,iDAGjBG,SAAU,CACR,CACEC,cAAe,CACb,OAINC,SAAU,CACR,mEACA,oFACA,qFAEFC,QAAS,CAAC,mBAAoB,iBAAkB,8BAElD,cAAe,CACbP,SAAU,SACVC,YAAa,0CACbC,QAAS,CACP/E,KAAM,WAERgF,KAAM,CACJM,EAAG,CACDtF,KAAM,SACN8E,YAAa,6CAGjBG,SAAU,CACR,CACEC,cAAe,CACb,OAINC,SAAU,CACR,uEACA,wFACA,wEAEFC,QAAS,CAAC,aAAc,aAAc,cAAe,mBAEvDS,QAAW,CACThB,SAAU,SACVC,YAAa,8CACbC,QAAS,CACP/E,KAAM,UAERgF,KAAM,CACJc,EAAG,CACD9F,KAAM,UACN8E,YAAa,oCAGjBG,SAAU,CACR,CACEC,cAAe,CACb,OAINC,SAAU,CACR,gDACA,iDAEFC,QAAS,CAAC,qBAAsB,sBAElCW,YAAe,CACblB,SAAU,SACVC,YAAa,oDACbC,QAAS,CACP/E,KAAM,UAERgF,KAAM,CACJgB,EAAG,CACDhG,KAAM,SACN8E,YAAa,wDAGjBG,SAAU,CACR,CACEC,cAAe,CACb,OAINC,SAAU,CACR,gEACA,iEAEFC,QAAS,CAAC,iBAAkB,gBAE9Ba,KAAQ,CACNpB,SAAU,SACVC,YAAa,wGACbC,QAAS,CACP/E,KAAM,UAERgF,KAAM,CACJc,EAAG,CACD9F,KAAM,UACN8E,YAAa,kCAEfoB,MAAO,CACLlG,KAAM,UACN8E,YAAa,yBAEfqB,MAAO,CACLnG,KAAM,UACN8E,YAAa,0BAGjBG,SAAU,CACR,CACEC,cAAe,CACb,IACA,QACA,WAINC,SAAU,CACR,gDACA,iDAEFC,QAAS,CAAC,iBAAkB,mBAAoB,uBAElD,UAAW,CACTP,SAAU,SACVC,YAAa,iGACbC,QAAS,CACP/E,KAAM,WAERgF,KAAM,CACJM,EAAG,CACDtF,KAAM,SACN8E,YAAa,kDAEfoB,MAAO,CACLlG,KAAM,UACN8E,YAAa,yBAEfqB,MAAO,CACLnG,KAAM,UACN8E,YAAa,0BAGjBG,SAAU,CACR,CACEC,cAAe,CACb,IACA,QACA,WAINC,SAAU,CACR,6HACA,8HAEFC,QAAS,CAAC,cAAe,qBAAsB,qBAEjDgB,KAAQ,CACNvB,SAAU,SACVC,YAAa,sEACbC,QAAS,CACP/E,KAAM,UAERgF,KAAM,CACJM,EAAG,CACDtF,KAAM,SACN8E,YAAa,2CAGjBG,SAAU,CACR,CACEC,cAAe,CACb,OAINC,SAAU,CACR,0EACA,0EACA,4EAEFC,QAAS,CAAC,aAAc,qBAAsB,wBAEhD,iBAAkB,CAChBP,SAAU,SACVC,YAAa,iDACbC,QAAS,CACP/E,KAAM,UAERgF,KAAM,CACJM,EAAG,CACDtF,KAAM,SACN8E,YAAa,qDAGjBG,SAAU,CACR,CACEC,cAAe,CACb,OAINC,SAAU,CACR,6EACA,+FAEFC,QAAS,CAAC,kBAAmB,kBAAmB,oBAElD,WAAY,CACVP,SAAU,SACVC,YAAa,yDACbC,QAAS,CACP/E,KAAM,UAERgF,KAAM,CACJM,EAAG,CACDtF,KAAM,SACN8E,YAAa,+CAGjBG,SAAU,CACR,CACEC,cAAe,CACb,OAINC,SAAU,CACR,iEACA,mFAEFC,QAAS,CAAC,wBAAyB,kBAAmB,oBAExD,WAAY,CACVP,SAAU,SACVC,YAAa,gDACbC,QAAS,CACP/E,KAAM,UAERgF,KAAM,CACJM,EAAG,CACDtF,KAAM,SACN8E,YAAa,oDAGjBG,SAAU,CACR,CACEC,cAAe,CACb,OAINC,SAAU,CACR,iEACA,mFAEFC,QAAS,CAAC,wBAAyB,kBAAmB,oBAExD,WAAY,CACVP,SAAU,SACVC,YAAa,2CACbC,QAAS,CACP/E,KAAM,UAERgF,KAAM,CACJM,EAAG,CACDtF,KAAM,SACN8E,YAAa,+CAGjBG,SAAU,CACR,CACEC,cAAe,CACb,OAINC,SAAU,CACR,iEACA,mFAEFC,QAAS,CAAC,wBAAyB,kBAAmB,qBCppBpD,SAAUiB,EAAY3C,GAE1B,MAAMoC,EAAIpC,EAAO/F,OACX2I,EAAgB,GAEtB,IAAK,IAAIC,EAAI,EAAGA,EAAIT,EAAGS,IACrBD,EAAEC,GAAK,IAAI7C,EAAO6C,IAIpB,GAAU,IAANT,EACF,OAAOQ,EAAE,GAAI,GAGf,GAAU,IAANR,EACF,OAAOQ,EAAE,GAAI,GAAMA,EAAE,GAAI,GAAMA,EAAE,GAAI,GAAMA,EAAE,GAAI,GAInD,IAAIrE,EAAO,EAGX,IAAK,IAAIsE,EAAI,EAAGA,EAAIT,EAAI,EAAGS,GAAK,EAAG,CAEjC,IAAIC,EAASD,EAEb,IAAK,IAAIE,EAAIF,EAAI,EAAGE,EAAIX,EAAGW,GAAK,EAC1B5I,KAAKyG,IAAIgC,EAAEG,GAAIF,IAAO1I,KAAKyG,IAAIgC,EAAEE,GAASD,MAC5CC,EAASC,GAKb,GAAIlC,EAAW+B,EAAEE,GAASD,IACxB,OAAO,EAILC,IAAWD,KACZD,EAAEC,GAAID,EAAEE,IAAW,CAACF,EAAEE,GAAUF,EAAEC,IACnCtE,GAAQA,GAIV,IAAK,IAAIwE,EAAIF,EAAI,EAAGE,EAAIX,EAAGW,GAAK,EAAG,CACjC,MAAMC,EAASJ,EAAEG,GAAIF,GAAMD,EAAEC,GAAIA,GAGjC,IAAK,IAAII,EAAIJ,EAAGI,EAAIb,EAAGa,IACrBL,EAAEG,GAAIE,IAAOD,EAASJ,EAAEC,GAAII,EAE/B,CACF,CAGD,IAAItB,EAAMpD,EACV,IAAK,IAAIsE,EAAI,EAAGA,EAAIT,EAAGS,IACrBlB,GAAOiB,EAAEC,GAAIA,GAGf,OAAOlB,CACT,UCpEgBK,EAAMhC,EAAoBE,EAAa+B,GACrD,MAAMG,EAAIpC,EAAO/F,OACXiJ,EAAqB,GAE3B,IAAK,IAAIL,EAAI,EAAGA,EAAIT,EAAGS,IACrB,GAAIA,IAAM3C,EAAK,CACb,MAAMiD,EAAqB,GAC3B,IAAK,IAAIJ,EAAI,EAAGA,EAAIX,EAAGW,IACjBA,IAAMd,GACRkB,EAASC,KAAKpD,EAAO6C,GAAIE,IAG7BG,EAAOE,KAAKD,EACb,CAGH,OAAOD,CACT,CCdM,SAAUG,EAASrD,GACvB,MAAMoC,EAAIpC,EAAO/F,OACX6H,EAAkB,GACxB,IAAK,IAAIe,EAAI,EAAGA,EAAIT,EAAGS,IAAK,CAC1Bf,EAAIe,GAAK,GACT,IAAK,IAAIE,EAAI,EAAGA,EAAIX,EAAGW,IAAK,CAC1B,MAEMhB,IADS,KAAOc,EAAIE,GACFJ,EAFZX,EAAMhC,EAAQ+C,EAAGF,IAG7Bf,EAAIe,GAAIE,GAAKhB,CACd,CACF,CACD,OAAOD,CACT,CChBM,SAAUwB,EAAStD,GACvB,OAAOA,EAAO/F,SAAW+F,EAAO,GAAI/F,MACtC,CCCM,SAAUsJ,EAAWvD,GACzB,IAAKsD,EAAStD,GACZ,OAAO,EAET,MAAMoC,EAAIpC,EAAO/F,OAEjB,IAAK,IAAI4I,EAAI,EAAGA,EAAIT,EAAGS,IACrB,IAAK,IAAIE,EAAI,EAAGA,EAAIX,EAAGW,IACrB,GAAIF,IAAME,GACR,IAAKxC,EAAYP,EAAO6C,GAAIE,GAAK,GAC/B,OAAO,OAIT,IAAKlC,EAAWb,EAAO6C,GAAIE,IACzB,OAAO,EAMf,OAAO,CACT,CCnBgB,SAAAS,EAAeZ,EAAea,GAE5C,GAAiB,IAAbb,EAAE3I,QAA6B,IAAbwJ,EAAExJ,QAAgB2I,EAAE,GAAI3I,SAAWwJ,EAAExJ,OACzD,MAAM,IAAIK,MAAM,qDAGlB,MAAMsH,EAAIgB,EAAE3I,OACNmI,EAAIQ,EAAE,GAAI3I,OACVyJ,EAAID,EAAE,GAAIxJ,OAGV0J,EAAiBhH,MAAMiF,GAAGgC,KAAK,GAAgBC,IAAI,IAAMlH,MAAM+G,GAAGE,KAAK,IAG7E,IAAK,IAAIf,EAAI,EAAGA,EAAIjB,EAAGiB,IACrB,IAAK,IAAII,EAAI,EAAGA,EAAIb,EAAGa,IAAK,CAC1B,MAAMa,EAAMlB,EAAEC,GAAII,GAClB,IAAK,IAAIF,EAAI,EAAGA,EAAIW,EAAGX,IACrBY,EAAEd,GAAIE,IAAOe,EAAML,EAAER,GAAIF,EAE5B,CAGH,OAAOY,CACT,CCJO,MAAMI,EAAmD,CAC9D7C,IAAO,CACL8C,SAAU,EAAEC,EAASC,GAAUvK,KAC7B0G,EAAa4D,EAAStK,GACtB0G,EAAa6D,EAASvK,GACtB,IACE,OAAO6J,EAAeS,EAASC,EAChC,CACD,MAAOC,GACL,MAAM,IAAI9J,EAAW,+GAA+G4J,EAAQ,GAAIhK,cAAciK,EAAQjK,SAAUN,EACjL,GAEHoH,MAAOD,EAAa,IAEtBa,IAAO,CACLqC,SAAU,EAAEhE,GAASrG,KACnB2G,EAAmBN,EAAQrG,GACpBgJ,EAAY3C,IAErBe,MAAOD,EAAa,IAEtBe,IAAO,CACLmC,SAAU,EAAEhE,GAASrG,KACnB2G,EAAmBN,EAAQrG,GAC3B,MAAMuJ,ECxCN,SAAkBlD,GACtB,MAAMoC,EAAIpC,EAAO/F,OAGjB,GAAU,IAANmI,EAAS,CACX,MAAMgC,EAAUpE,EAAO,GAAI,GAC3B,OAAIa,EAAWuD,GACN,KAEF,CAAC,CAAC,EAAIA,GACd,CAGD,MAAMzC,EAAMgB,EAAY3C,GAGxB,GAAIa,EAAWc,GACb,OAAO,KAIT,MAAMG,EAAMuB,EAASrD,GAGfqE,EAA4B,GAClC,IAAK,IAAIxB,EAAI,EAAGA,EAAIT,EAAGS,IAAK,CAC1BwB,EAAcxB,GAAK,GACnB,IAAK,IAAIE,EAAI,EAAGA,EAAIX,EAAGW,IACrBsB,EAAcxB,GAAIE,GAAKjB,EAAIe,GAAIE,GAAMpB,CAExC,CAED,OAAO0C,CACT,CDOqBC,CAAQtE,GACvB,GAAe,OAAXkD,EACF,MAAM,IAAI7I,EAAW,gCAAiCV,GAExD,OAAOuJ,GAETnC,MAAOD,EAAa,IAEtBgB,IAAO,CACLkC,SAAU,EAAEhE,GAASrG,KACnB2G,EAAmBN,EAAQrG,GACpB0J,EAASrD,IAElBe,MAAOD,EAAa,IAEtBiB,SAAY,CACViC,SAAU,EAAEhE,GAASrG,KACnB2G,EAAmBN,EAAQrG,GEhE3B,SAAmBqG,GACvB,MAAMoC,EAAIpC,EAAO/F,OACXsK,EAAwB,GAG9B,IAAK,IAAI1B,EAAI,EAAGA,EAAIT,EAAGS,IAAK,CAC1B0B,EAAU1B,GAAK,GACf,IAAK,IAAIE,EAAI,EAAGA,EAAIX,EAAGW,IAAK,CAE1B,MAAM/B,EAAMgB,EAAMhC,EAAQ6C,EAAGE,GACvBxE,IAAS,KAAOsE,EAAIE,GAC1BwB,EAAU1B,GAAIE,GAAKxE,EAAOoE,EAAY3B,EACvC,CACF,CAED,OAAOuD,CACT,CFiDaxC,CAAS/B,IAElBe,MAAOD,EAAa,IAEtBkB,MAAS,CACPgC,SAAU,EAAEhE,EAAQE,EAAK+B,GAAMtI,KAC7B0G,EAAaL,EAAQrG,GACrB2E,EAAa4B,EAAKvG,EAAgB,CAAE6D,SAAS,EAAMS,aAAa,EAAMI,IAAK2B,EAAO/F,SAClFqE,EAAa2D,EAAKtI,EAAgB,CAAE6D,SAAS,EAAMS,aAAa,EAAMI,IAAK2B,EAAO,GAAI/F,SAE/E+H,EAAMhC,EAAQE,EAAK+B,IAE5BlB,MAAOD,EAAa,IAEtBoB,MAAS,CACP8B,SAAU,EAAEhE,GAASrG,KACnB2G,EAAmBN,EAAQrG,GG7E3B,SAAgBqG,GACpB,OAAOA,EAAOwE,OAAO,CAACC,EAAKvE,EAAK2C,IAAM4B,EAAMvE,EAAI2C,GAAK,EACvD,CH4EaX,CAAMlC,IAEfe,MAAOD,EAAa,IAEtB,aAAc,CACZkD,SAAU,EAAEhE,GAASrG,KACnB0G,EAAaL,EAAQrG,GIlFrB,SAAqBqG,GACzB,MAAM0E,EAAO1E,EAAO/F,OAGpB,IAAKqJ,EAAStD,GACZ,OAAO,EAIT,IAAK,IAAI6C,EAAI,EAAGA,EAAI6B,EAAM7B,GAAK,EAC7B,IAAK,IAAIE,EAAI,EAAGA,EAAIF,EAAGE,GAAK,EAC1B,GAAI/C,EAAO6C,GAAIE,KAAQ/C,EAAO+C,GAAIF,GAChC,OAAO,EAKb,OAAO,CACT,CJiEa8B,CAAW3E,IAEpBe,MAAOD,EAAa,IAEtB,cAAe,CACbkD,SAAU,EAAEhE,GAASrG,KACnB0G,EAAaL,EAAQrG,GKtFrB,SAAuBqG,GAC3B,IAAKsD,EAAStD,GACZ,OAAO,EAGT,MAAM0E,EAAO1E,EAAO/F,OAEpB,IAAI2K,GAAoB,EACpBC,GAAoB,EAExB,IAAK,IAAIhC,EAAI,EAAGA,EAAI6B,EAAM7B,IACxB,IAAK,IAAIE,EAAI,EAAGA,EAAI2B,EAAM3B,IAAK,CAC7B,GAAIF,EAAIE,GAAuB,IAAlB/C,EAAO6C,GAAIE,KACtB6B,GAAoB,GACfC,GACH,OAAO,EAGPhC,EAAIE,GAAuB,IAAlB/C,EAAO6C,GAAIE,KACtB8B,GAAoB,EAEvB,CAGH,OAAOD,GAAqBC,CAC9B,CL8DaC,CAAa9E,IAEtBe,MAAOD,EAAa,IAEtB,oBAAqB,CACnBkD,SAAU,EAAEhE,GAASrG,KACnB0G,EAAaL,EAAQrG,GKlErB,SAA4BqG,GAChC,IAAKsD,EAAStD,GACZ,OAAO,EAET,MAAM0E,EAAO1E,EAAO/F,OAEpB,IAAK,IAAI4I,EAAI,EAAGA,EAAI6B,EAAM7B,IACxB,IAAK,IAAIE,EAAI,EAAGA,EAAIF,EAAGE,IACrB,GAAsB,IAAlB/C,EAAO6C,GAAIE,GACb,OAAO,EAKb,OAAO,CACT,CLoDagC,CAAkB/E,IAE3Be,MAAOD,EAAa,IAEtB,oBAAqB,CACnBkD,SAAU,EAAEhE,GAASrG,KACnB0G,EAAaL,EAAQrG,GKxDrB,SAA4BqG,GAChC,IAAKsD,EAAStD,GACZ,OAAO,EAET,MAAM0E,EAAO1E,EAAO/F,OAGpB,IAAK+F,EAAOL,MAAMO,GAAOA,EAAIjG,SAAWyK,GACtC,OAAO,EAGT,IAAK,IAAI7B,EAAI,EAAGA,EAAI6B,EAAM7B,IACxB,IAAK,IAAIE,EAAIF,EAAI,EAAGE,EAAI2B,EAAM3B,IAC5B,GAAsB,IAAlB/C,EAAO6C,GAAIE,GACb,OAAO,EAKb,OAAO,CACT,CLqCaiC,CAAkBhF,IAE3Be,MAAOD,EAAa,IAEtB,YAAa,CACXkD,SAAU,EAAEhE,GAASrG,KACnB0G,EAAaL,EAAQrG,GM5GrB,SAAqBqG,GACzB,IAAKsD,EAAStD,GACZ,OAAO,EAET,MAAM0E,EAAO1E,EAAO/F,OAEpB,IAAK,IAAI4I,EAAI,EAAGA,EAAI6B,EAAM7B,GAAK,EAC7B,IAAK,IAAIE,EAAI,EAAGA,EAAI2B,EAAM3B,GAAK,EAC7B,GAAIF,IAAME,GAAuB,IAAlB/C,EAAO6C,GAAIE,GACxB,OAAO,EAKb,OAAO,CACT,CN8FakC,CAAWjF,IAEpBe,MAAOD,EAAa,IAEtB,UAAW,CACTkD,SAAU,EAAEhE,GAASrG,KACnB0G,EAAaL,EAAQrG,GACd2J,EAAStD,IAElBe,MAAOD,EAAa,IAEtB,qBAAsB,CACpBkD,SAAU,EAAEhE,GAASrG,KACnB0G,EAAaL,EAAQrG,GOhIrB,SAAuBqG,GAC3B,QAAKsD,EAAStD,IAWPuD,EAHSC,EAAexD,ECZ3B,SAAmCkF,GACvC,MAAMhC,EAAgB,GAEtB,IAAK,IAAIL,EAAI,EAAGA,EAAIqC,EAAK,GAAIjL,OAAQ4I,GAAK,EAAG,CAC3C,MAAM3C,EAAW,GACjB,IAAK,IAAI6C,EAAI,EAAGA,EAAImC,EAAKjL,OAAQ8I,GAAK,EACpC7C,EAAIkD,KAAK8B,EAAKnC,GAAIF,IAEpBK,EAAOE,KAAKlD,EACb,CACD,OAAOgD,CACT,CDFqBiC,CAAUnF,IAO/B,CPoHaoF,CAAapF,IAEtBe,MAAOD,EAAa,IAEtB,YAAa,CACXkD,SAAU,EAAEhE,GAASrG,KACnB0G,EAAaL,EAAQrG,GACd4J,EAAWvD,IAEpBe,MAAOD,EAAa,IAEtB,cAAe,CACbkD,SAAU,EAAEhE,GAASrG,KACnB0G,EAAaL,EAAQrG,KVwBrB,SAAyBqG,GAC7B,QAAKD,EAASC,IAGVA,EAAO/F,SAAW+F,EAAO,GAAI/F,MAInC,CU/BWoL,CAAerF,KAGZa,EAAW8B,EAAY3C,KAEjCe,MAAOD,EAAa,IAEtBqB,QAAW,CACT6B,SAAU,EAAEsB,GAAO3L,KACjB2E,EAAagH,EAAM3L,EAAgB,CAAE6D,SAAS,EAAMM,UAAU,IAC9D,MAAMoF,EAAqB,GAC3B,IAAK,IAAIL,EAAI,EAAGA,EAAIyC,EAAMzC,GAAK,EAAG,CAChC,MAAM3C,EAAgB,GACtB,IAAK,IAAI6C,EAAI,EAAGA,EAAIuC,EAAMvC,GAAK,EAC7B7C,EAAIkD,KAAK,GAAKP,EAAIE,EAAI,IAExBG,EAAOE,KAAKlD,EACb,CACD,OAAOgD,GAETnC,MAAOD,EAAa,IAEtBuB,YAAe,CACb2B,SAAU,EAAEtE,GAAS/F,KACnB8F,EAAaC,EAAQ/F,GACrB,MAAMuJ,EAAqB,GAC3B,IAAK,IAAIL,EAAI,EAAGA,EAAInD,EAAOzF,OAAQ4I,GAAK,EAAG,CACzC,MAAM3C,EAAgB,GACtB,IAAK,IAAI6C,EAAI,EAAGA,EAAIrD,EAAOzF,OAAQ8I,GAAK,EACtC7C,EAAIkD,KAAM1D,EAAOmD,IAAQE,GAE3BG,EAAOE,KAAKlD,EACb,CACD,OAAOgD,GAETnC,MAAOD,EAAa,IAEtByB,KAAQ,CACNyB,SAAU,EAAE5B,EAAGI,EAAOC,GAAQ9I,KAC5B2E,EAAa8D,EAAGzI,EAAgB,CAAE6D,SAAS,EAAMM,UAAU,IAC3DQ,EAAakE,EAAO7I,EAAgB,CAAE6D,SAAS,EAAMS,aAAa,EAAMG,GAAIgE,IAC5E9D,EAAamE,EAAO9I,EAAgB,CAAE6D,SAAS,EAAMS,aAAa,EAAMI,IAAK+D,aSrL9DA,EAAWI,EAAeC,GAE7C,MAAMzC,EAAqBrD,MAAM4I,KAAK,CAAEtL,OAAQmI,GAAK,IAAMzF,MAAM4I,KAAK,CAAEtL,OAAQmI,GAAK,IAAM,IAG3F,IAAK,IAAIS,EAAI,EAAGA,EAAIT,EAAGS,IACrB,IAAK,IAAIE,EAAI5I,KAAKC,IAAI,EAAGyI,EAAIL,GAAQO,GAAK5I,KAAK6G,IAAIoB,EAAI,EAAGS,EAAIJ,GAAQM,IACpE/C,EAAO6C,GAAIE,GAAK,EAIpB,OAAO/C,CACT,CT0KauC,CAAKH,EAAGI,EAAOC,IAExB1B,MAAOD,EAAa,IAEtB,UAAW,CACTkD,SAAU,EAAEhE,EAAQwC,EAAOC,GAAQ9I,KACjC0G,EAAaL,EAAQrG,GACrB,MAAM6L,EAAUrL,KAAKC,IAAI4F,EAAO/F,OAAQ+F,EAAO,GAAI/F,QAGnD,OAFAqE,EAAakE,EAAO7I,EAAgB,CAAE6D,SAAS,EAAMS,aAAa,EAAMG,GAAIoH,IAC5ElH,EAAamE,EAAO9I,EAAgB,CAAE6D,SAAS,EAAMS,aAAa,EAAMG,GAAIoH,aU7LzDxF,EAAoByF,EAAeC,GAC1D,MAAMhB,EAAO1E,EAAO/F,OACd0L,EAAO3F,EAAO,GAAI/F,OAGxB,IAAK,IAAI4I,EAAI,EAAGA,EAAI6B,EAAM7B,IACxB,IAAK,IAAIE,EAAI,EAAGA,EAAI4C,EAAM5C,IAExB,GAAsB,IAAlB/C,EAAO6C,GAAIE,KAAaF,EAAIE,EAAI0C,GAAS1C,EAAIF,EAAI6C,GACnD,OAAO,EAMb,OAAO,CACT,CV8KaE,CAAS5F,EAAQwC,EAAOC,IAEjC1B,MAAOD,EAAa,IAEtB4B,KAAQ,CACNsB,SAAU,EAAEhE,GAASrG,KACnB0G,EAAaL,EAAQrG,GACrB,OAASuJ,GWvMT,SAAiClD,GAErC,MAAM4B,EAAI5B,EAAO6D,IAAI3D,GAAO,IAAIA,IAC1BwE,EAAO9C,EAAE3H,OACT0L,EAAO/D,EAAE,GAAI3H,OAEnB,IAAIyI,EAAO,EACPmD,EAAgB,EAGpB,IAAK,IAAI5D,EAAM,EAAGA,EAAM0D,EAAM1D,IAAO,CAEnC,IAAI6D,GAAY,EAEhB,IAAK,IAAI5F,EAAM2F,EAAe3F,EAAMwE,EAAMxE,IACxC,IAAKW,EAAWe,EAAE1B,GAAM+B,IAAQ,CAC9B6D,EAAW5F,EACX,KACD,CAGH,IAAkB,IAAd4F,EACF,SAGFpD,GAAQ,EAGJoD,IAAaD,KACdjE,EAAEkE,GAAWlE,EAAEiE,IAAkB,CAACjE,EAAEiE,GAAiBjE,EAAEkE,KAI1D,MAAMC,EAAanE,EAAEiE,GAAgB5D,GAGrC,IAAK,IAAIc,EAAId,EAAKc,EAAI4C,EAAM5C,IAC1BnB,EAAEiE,GAAgB9C,IAAOgD,EAI3B,IAAK,IAAI7F,EAAM,EAAGA,EAAMwE,EAAMxE,IAC5B,GAAIA,IAAQ2F,IAAkBhF,EAAWe,EAAE1B,GAAM+B,IAAQ,CACvD,MAAMe,EAASpB,EAAE1B,GAAM+B,GACvB,IAAK,IAAIc,EAAId,EAAKc,EAAI4C,EAAM5C,IAC1BnB,EAAE1B,GAAM6C,IAAOC,EAASpB,EAAEiE,GAAgB9C,EAE7C,CAIH,GADA8C,IACIA,IAAkBnB,EACpB,KACH,CAED,MAAO,CAAC9C,EAAGc,EACb,CX+IyBsD,CAAuBhG,GAC1C,OAAOkD,GAETnC,MAAOD,EAAa,IAGtB,iBAAkB,CAChBkD,SAAU,EAAEhE,GAASrG,KACnB0G,EAAaL,EAAQrG,GACdQ,KAAK8L,KAAKjG,EAAOwE,OAAO,CAACC,EAAKvE,IAAQuE,EAAMvE,EAAIsE,OAAO,CAAC0B,EAAQ9F,IAAS8F,EAAS9F,EAAOA,EAAM,GAAI,KAE5GW,MAAOD,EAAa,IAGtB,WAAY,CACVkD,SAAU,EAAEhE,GAASrG,KACnB0G,EAAaL,EAAQrG,GY9NrB,SAAgBqG,GACpB,MAAMmG,EAAUnG,EAAO/F,OACjBmM,EAAUpG,EAAO,GAAI/F,OAE3B,IAAIoM,EAAY,EAGhB,IAAK,IAAItD,EAAI,EAAGA,EAAIqD,EAASrD,GAAK,EAAG,CACnC,IAAIuD,EAAS,EAGb,IAAK,IAAIzD,EAAI,EAAGA,EAAIsD,EAAStD,GAAK,EAChCyD,GAAUnM,KAAKyG,IAAIZ,EAAO6C,GAAIE,IAIhCsD,EAAYlM,KAAKC,IAAIiM,EAAWC,EACjC,CAED,OAAOD,CACT,CZ2MaE,CAAMvG,IAEfe,MAAOD,EAAa,IAGtB,WAAY,CACVkD,SAAU,EAAEhE,GAASrG,KACnB0G,EAAaL,EAAQrG,GACdqG,EAAOwE,OAAO,CAACpK,EAAK8F,IAAQ/F,KAAKC,IAAIA,EAAK8F,EAAIsE,OAAO,CAACC,EAAKrE,IAASqE,EAAMtK,KAAKyG,IAAIR,GAAO,IAAK,IAExGW,MAAOD,EAAa,IAGtB,WAAY,CACVkD,SAAU,EAAEhE,GAASrG,KACnB0G,EAAaL,EAAQrG,GACdqG,EAAOwE,OAAO,CAACgC,EAAQtG,KAC5B,MAAMuG,EAASvG,EAAIsE,OAAO,CAACpK,EAAKsM,IAAQvM,KAAKC,IAAIA,EAAKD,KAAKyG,IAAI8F,IAAO,GACtE,OAAOvM,KAAKC,IAAIoM,EAAQC,IACvB,IAEL1F,MAAOD,EAAa,KAIxB,IAAK,MAAO9D,EAAK2J,KAAS1L,OAAO2L,QAAQ3F,GACnC8C,EAAuB/G,KACzB+G,EAAuB/G,GAAK2J,KAAOA,GAG1B,MAAAE,EAA4B,CACvCzL,KAAM,SACN0L,UAAW/C"}

package/dist/modules/matrix.js

@@ -0,0 +1,2 @@
+"use strict";function t(t){if(!t.position||!t.code)return"";const e=t.position.column-1,r=t.code.length-e-1;return`${" ".repeat(Math.max(e,0))}^${" ".repeat(Math.max(r,0))}`}class e extends Error{sourceCodeInfo;shortMessage;constructor(r,a){const n=r instanceof Error?r.message:`${r}`;super(function(e,r){if(!r)return e;const a=`${r.position.line}:${r.position.column}`;return`${e}${r.filePath?`\n${r.filePath}:${a}`:`\nLocation ${a}`}\n${r.code}\n${t(r)}`}(n,a)),this.shortMessage=n,this.sourceCodeInfo=a,Object.setPrototypeOf(this,e.prototype),this.name="DvalaError"}getCodeMarker(){return this.sourceCodeInfo&&t(this.sourceCodeInfo)}}const r={Number:1,String:2,NormalExpression:3,SpecialExpression:4,UserDefinedSymbol:5,NormalBuiltinSymbol:6,SpecialBuiltinSymbol:7,ReservedSymbol:8,Binding:9,Spread:10},a=new Set(Object.values(r));const n=new Set(["UserDefined","Partial","Comp","Constantly","Juxt","Complement","EveryPred","SomePred","Fnull","Builtin","SpecialBuiltin","NativeJsFunction","Module"]);function i(t){return null!==t&&"object"==typeof t&&("^^fn^^"in t&&"functionType"in t&&("string"==typeof(e=t.functionType)&&n.has(e)));var e}function o(t){return!(!Array.isArray(t)||t.length<2)&&("number"==typeof(e=t[0])&&a.has(e));var e}function m(t){return i(t)?`<function ${t.name||"λ"}>`:o(t)?`${e=t[0],Object.keys(r).find(t=>r[t]===e)}-node`:null===t?"null":"object"==typeof t&&t instanceof RegExp?`${t}`:"object"==typeof t&&t instanceof Error?t.toString():JSON.stringify(t);var e}function l(t,e={}){return"number"==typeof t&&(!Number.isNaN(t)&&(!(e.integer&&!Number.isInteger(t))&&(!(e.finite&&!Number.isFinite(t))&&((!e.zero||0===t)&&((!e.nonZero||0!==t)&&(!(e.positive&&t<=0)&&(!(e.negative&&t>=0)&&(!(e.nonPositive&&t>0)&&(!(e.nonNegative&&t<0)&&(!("number"==typeof e.gt&&t<=e.gt)&&(!("number"==typeof e.gte&&t<e.gte)&&(!("number"==typeof e.lt&&t>=e.lt)&&!("number"==typeof e.lte&&t>e.lte)))))))))))))}function s(t,r,a={}){if(!l(t,a))throw new e(`Expected ${function(t){if(t.zero)return"zero";const e=function(t){return t.positive?"positive":t.negative?"negative":t.nonNegative?"non negative":t.nonPositive?"non positive":t.nonZero?"non zero":""}(t),r=t.integer?"integer":"number",a=t.finite?"finite":"",n=function(t){return"number"!=typeof t.gt&&"number"!=typeof t.gte||"number"!=typeof t.lt&&"number"!=typeof t.lte?"number"==typeof t.gt||"number"==typeof t.gte?"number"==typeof t.gt?`n > ${t.gt}`:`n >= ${t.gte}`:"number"==typeof t.lt||"number"==typeof t.lte?"number"==typeof t.lt?`n < ${t.lt}`:`n <= ${t.lte}`:"":`${"number"==typeof t.gt?`${t.gt} < n `:`${t.gte} <= n `}${"number"==typeof t.lt?`< ${t.lt}`:`<= ${t.lte}`}`}(t);return[e,a,r,n].filter(t=>!!t).join(" ")}(a)}, got ${m(t)}.`,function(t,e){return t?.sourceCodeInfo??e}(t,r))}const u=new WeakSet,c=new WeakSet,x=new WeakSet,p=new WeakSet,f=new WeakSet,d=new WeakSet;function g(t,r){if(!function(t){return!(!Array.isArray(t)||!c.has(t)&&(x.has(t)||(t.every(t=>l(t))?(u.add(t),c.add(t),0):(x.add(t),1))))}(t))throw new e(`Expected a vector, but got ${t}`,r)}function h(t){if(!Array.isArray(t))return!1;if(p.has(t))return!0;if(f.has(t))return!1;if(0===t.length)return f.add(t),!1;if(!Array.isArray(t[0])||0===t[0].length)return f.add(t),!1;const e=t[0].length;for(const r of t)if(!Array.isArray(r)||r.length!==e||r.some(t=>!l(t)))return f.add(t),!1;return u.add(t),d.add(t),p.add(t),!0}function y(t,r){if(!h(t))throw new e(`Expected a matrix, but got ${t}`,r)}function b(t,r){if(!h(t))throw new e(`Expected a matrix, but got ${t}`,r);if(t.length!==t[0].length)throw new e(`Expected square matrix, but got ${t.length} and ${t[0].length}`,r)}function v(t,e,r=1e-10){if(t===e)return!0;const a=Math.abs(t-e);if(0===t||0===e||a<r)return a<r;return a/(Math.abs(t)+Math.abs(e))<r}function w(t){return Math.abs(t)<1e-10}function N(t){return{min:t,max:t}}const k={mul:{category:"matrix",description:"Multiplies two `matrices` using standard `matrix` multiplication based on **dot products** of rows and columns.",returns:{type:"matrix"},args:{a:{type:"matrix"},b:{type:"matrix"}},variants:[{argumentNames:["a","b"]}],examples:["let { mul } = import(matrix);\nmul([[1, 2], [3, 4]], [[5, 6], [7, 8]])","let { mul } = import(matrix);\nmul([[1, 2, 3], [4, 5, 6]], [[7, 8], [9, 10], [11, 12]])"],seeAlso:["matrix.det","matrix.inv"]},det:{category:"matrix",description:"Calculates the **determinant** of a square matrix.",returns:{type:"number"},args:{m:{type:"matrix",description:"The `matrix` to calculate the determinant of."}},variants:[{argumentNames:["m"]}],examples:["let { det } = import(matrix);\ndet([[1, 2], [3, 4]])","let { det } = import(matrix);\ndet([[1, 2, 3], [4, 5, 6], [7, 8, 9]])"],seeAlso:["matrix.inv","matrix.cofactor","matrix.adj","matrix.trace","matrix.rank","matrix.invertible?","matrix.mul","matrix.minor"]},inv:{category:"matrix",description:"Calculates the **inverse** of a square matrix.",returns:{type:"matrix"},args:{m:{type:"matrix",description:"The `matrix` to calculate the inverse of."}},variants:[{argumentNames:["m"]}],examples:["let { inv } = import(matrix);\ninv([[1, 2], [3, 4]])","let { inv } = import(matrix);\ninv([[1, 2, 3], [4, 5, 7], [7, 8, 10]])"],seeAlso:["matrix.det","matrix.adj","matrix.invertible?","linear-algebra.solve","matrix.mul","matrix.orthogonal-matrix?"]},adj:{category:"matrix",description:"Calculates the **adjugate** of a square matrix.",returns:{type:"matrix"},args:{m:{type:"matrix",description:"The `matrix` to calculate the adjugate of."}},variants:[{argumentNames:["m"]}],examples:["let { adj } = import(matrix);\nadj([[1, 2], [3, 4]])","let { adj } = import(matrix);\nadj([[1, 2, 3], [4, 5, 6], [7, 8, 9]])","let { adj } = import(matrix);\nadj([[1, 2, 3], [7, 8, 9], [4, 5, 6]])"],seeAlso:["matrix.cofactor","matrix.det","matrix.inv"]},cofactor:{category:"matrix",description:"Calculates the **cofactor** of a square matrix.",returns:{type:"matrix"},args:{m:{type:"matrix",description:"The `matrix` to calculate the cofactor of."}},variants:[{argumentNames:["m"]}],examples:["let { cofactor } = import(matrix);\ncofactor([[1, 2], [3, 4]])","let { cofactor } = import(matrix);\ncofactor([[1, 2, 3], [4, 5, 6], [7, 8, 9]])","let { cofactor } = import(matrix);\ncofactor([[1, 2, 3], [7, 8, 9], [4, 5, 6]])"],seeAlso:["matrix.adj","matrix.minor","matrix.det"]},minor:{category:"matrix",description:"Calculates the **minor** of a square matrix.",returns:{type:"matrix"},args:{m:{type:"matrix",description:"The `matrix` to calculate the minor of."},row:{type:"integer",description:"The row index of the element to calculate the minor for."},col:{type:"integer",description:"The column index of the element to calculate the minor for."}},variants:[{argumentNames:["m","row","col"]}],examples:["let { minor } = import(matrix);\nminor([[1, 2], [3, 4]], 0, 1)","let { minor } = import(matrix);\nminor([[1, 2, 3], [4, 5, 6], [7, 8, 9]], 1, 1)"],seeAlso:["matrix.cofactor","matrix.det"]},trace:{category:"matrix",description:"Calculates the **trace** of a square matrix.",returns:{type:"number"},args:{m:{type:"matrix",description:"The `matrix` to calculate the trace of."}},variants:[{argumentNames:["m"]}],examples:["let { trace } = import(matrix);\ntrace([[1, 2], [3, 4]])","let { trace } = import(matrix);\ntrace([[1, 2, 3], [4, 5, 6], [7, 8, 9]])"],seeAlso:["matrix.det","matrix.diagonal?"]},"symmetric?":{category:"matrix",description:"Checks if a `matrix` is **symmetric**.",returns:{type:"boolean"},args:{m:{type:"matrix",description:"The `matrix` to check for symmetry."}},variants:[{argumentNames:["m"]}],examples:["let { symmetric? } = import(matrix);\nsymmetric?([[1, 2], [2, 1]])","let { symmetric? } = import(matrix);\nsymmetric?([[1, 2, 3], [2, 1, 4], [3, 4, 1]])"],seeAlso:["matrix.orthogonal-matrix?","matrix.diagonal?","matrix.square?","matrix.hilbert"]},"triangular?":{category:"matrix",description:"Checks if a `matrix` is **triangular**.",returns:{type:"boolean"},args:{m:{type:"matrix",description:"The `matrix` to check for triangularity."}},variants:[{argumentNames:["m"]}],examples:["let { triangular? } = import(matrix);\ntriangular?([[2, 0], [0, 1]])","let { triangular? } = import(matrix);\ntriangular?([[1, 2, 3], [0, 4, 5], [0, 0, 6]])"],seeAlso:["matrix.upper-triangular?","matrix.lower-triangular?","matrix.diagonal?","matrix.banded?"]},"upper-triangular?":{category:"matrix",description:"Checks if a `matrix` is **upper triangular**.",returns:{type:"boolean"},args:{m:{type:"matrix",description:"The `matrix` to check for upper triangularity."}},variants:[{argumentNames:["m"]}],examples:["let { upper-triangular? } = import(matrix);\nupper-triangular?([[1, 2], [0, 3]])","let { upper-triangular? } = import(matrix);\nupper-triangular?([[1, 2, 3], [0, 4, 5], [0, 0, 6]])"],seeAlso:["matrix.lower-triangular?","matrix.triangular?","matrix.diagonal?"]},"lower-triangular?":{category:"matrix",description:"Checks if a `matrix` is **lower triangular**.",returns:{type:"boolean"},args:{m:{type:"matrix",description:"The `matrix` to check for lower triangularity."}},variants:[{argumentNames:["m"]}],examples:["let { lower-triangular? } = import(matrix);\nlower-triangular?([[1, 0], [2, 3]])","let { lower-triangular? } = import(matrix);\nlower-triangular?([[1, 0, 0], [2, 3, 0], [4, 5, 6]])"],seeAlso:["matrix.upper-triangular?","matrix.triangular?","matrix.diagonal?"]},"diagonal?":{category:"matrix",description:"Checks if a `matrix` is **diagonal**.",returns:{type:"boolean"},args:{m:{type:"matrix",description:"The `matrix` to check for diagonal property."}},variants:[{argumentNames:["m"]}],examples:["let { diagonal? } = import(matrix);\ndiagonal?([[1, 0], [0, 2]])","let { diagonal? } = import(matrix);\ndiagonal?([[1, 0, 0], [0, 2, 0], [0, 0, 3]])","let { diagonal? } = import(matrix);\ndiagonal?([[1, 0, 0], [2, 2, 2], [0, 0, 3]])"],seeAlso:["matrix.identity?","matrix.symmetric?","matrix.triangular?","matrix.trace","matrix.upper-triangular?","matrix.lower-triangular?","matrix.band","matrix.banded?"]},"square?":{category:"matrix",description:"Checks if a `matrix` is **square**.",returns:{type:"boolean"},args:{m:{type:"matrix",description:"The `matrix` to check for square property."}},variants:[{argumentNames:["m"]}],examples:["let { square? } = import(matrix);\nsquare?([[1, 2], [3, 4]])","let { square? } = import(matrix);\nsquare?([[1, 2, 3], [4, 5, 6], [7, 8, 9]])","let { square? } = import(matrix);\nsquare?([[1, 2, 3], [4, 5, 6]])"],seeAlso:["matrix.symmetric?","matrix.identity?","matrix.invertible?"]},"orthogonal-matrix?":{category:"matrix",description:"Checks if a `matrix` is **orthogonal**.",returns:{type:"boolean"},args:{m:{type:"matrix",description:"The `matrix` to check for **orthogonality**."}},variants:[{argumentNames:["m"]}],examples:["let { orthogonal-matrix? } = import(matrix);\northogonal-matrix?([[1, 0], [0, 1]])","let { orthogonal-matrix? } = import(matrix);\northogonal-matrix?([[1, 0], [0, -1]])","let { orthogonal-matrix? } = import(matrix);\northogonal-matrix?([[1, 2], [3, 4]])"],seeAlso:["matrix.symmetric?","matrix.inv","matrix.identity?","linear-algebra.orthogonal?"]},"identity?":{category:"matrix",description:"Checks if a `matrix` is an **identity matrix**.",returns:{type:"boolean"},args:{m:{type:"matrix",description:"The `matrix` to check for identity property."}},variants:[{argumentNames:["m"]}],examples:["let { identity? } = import(matrix);\nidentity?([[1, 0], [0, 1]])","let { identity? } = import(matrix);\nidentity?([[1, 0, 0], [0, 1, 0], [0, 0, 1]])","let { identity? } = import(matrix);\nidentity?([[1, 0, 0], [0, 1, 0], [0, 0, 0]])"],seeAlso:["matrix.diagonal?","matrix.square?","matrix.orthogonal-matrix?"]},"invertible?":{category:"matrix",description:"Checks if a `matrix` is **invertible**.",returns:{type:"boolean"},args:{m:{type:"matrix",description:"The `matrix` to check for invertibility."}},variants:[{argumentNames:["m"]}],examples:["let { invertible? } = import(matrix);\ninvertible?([[1, 2], [3, 4]])","let { invertible? } = import(matrix);\ninvertible?([[1, 2, 3], [4, 5, 6], [7, 8, 9]])","let { invertible? } = import(matrix);\ninvertible?([[1, 2], [2, 4]])"],seeAlso:["matrix.det","matrix.inv","matrix.rank","matrix.square?"]},hilbert:{category:"matrix",description:"Generates a **Hilbert matrix** of size `n`.",returns:{type:"matrix"},args:{n:{type:"integer",description:"The size of the Hilbert matrix."}},variants:[{argumentNames:["n"]}],examples:["let { hilbert } = import(matrix);\nhilbert(3)","let { hilbert } = import(matrix);\nhilbert(4)"],seeAlso:["matrix.vandermonde","matrix.symmetric?"]},vandermonde:{category:"matrix",description:"Generates a **Vandermonde matrix** from a vector.",returns:{type:"matrix"},args:{v:{type:"vector",description:"The vector to generate the Vandermonde matrix from."}},variants:[{argumentNames:["v"]}],examples:["let { vandermonde } = import(matrix);\nvandermonde([1, 2, 3])","let { vandermonde } = import(matrix);\nvandermonde([1, 0, 1])"],seeAlso:["matrix.hilbert","matrix.band"]},band:{category:"matrix",description:"Generates a **banded matrix** of size `n` with lower band index `lband` and upper band index `uband`.",returns:{type:"matrix"},args:{n:{type:"integer",description:"The size of the banded matrix."},lband:{type:"integer",description:"The lower band index."},uband:{type:"integer",description:"The upper band index."}},variants:[{argumentNames:["n","lband","uband"]}],examples:["let { band } = import(matrix);\nband(3, 1, 1)","let { band } = import(matrix);\nband(4, 1, 2)"],seeAlso:["matrix.banded?","matrix.diagonal?","matrix.vandermonde"]},"banded?":{category:"matrix",description:"Checks if a `matrix` is **banded** with lower band index `lband` and upper band index `uband`.",returns:{type:"boolean"},args:{m:{type:"matrix",description:"The `matrix` to check for **banded** property."},lband:{type:"integer",description:"The lower band index."},uband:{type:"integer",description:"The upper band index."}},variants:[{argumentNames:["m","lband","uband"]}],examples:["let { banded? } = import(matrix);\nbanded?([\n  [1, 1, 1, 0],\n  [1, 1, 1, 1],\n  [1, 1, 1, 1],\n  [0, 1, 1, 1],\n], 2, 2)","let { banded? } = import(matrix);\nbanded?([\n  [1, 1, 1, 0],\n  [1, 1, 1, 1],\n  [1, 1, 1, 1],\n  [0, 1, 1, 1],\n], 1, 1)"],seeAlso:["matrix.band","matrix.triangular?","matrix.diagonal?"]},rank:{category:"matrix",description:"Calculates the **rank** of a matrix using **Gaussian elimination**.",returns:{type:"number"},args:{m:{type:"matrix",description:"The `matrix` to calculate the rank of."}},variants:[{argumentNames:["m"]}],examples:["let { rank } = import(matrix);\nrank([[1, 0, 0], [0, 1, 0], [0, 0, 1]])","let { rank } = import(matrix);\nrank([[1, 2, 3], [4, 5, 6], [7, 8, 9]])","let { rank } = import(matrix);\nrank([[2, 4, 6], [3, 6, 9], [4, 8, 12]])"],seeAlso:["matrix.det","matrix.invertible?","linear-algebra.rref"]},"frobenius-norm":{category:"matrix",description:"Calculates the **Frobenius norm** of a matrix.",returns:{type:"number"},args:{m:{type:"matrix",description:"The `matrix` to calculate the Frobenius norm of."}},variants:[{argumentNames:["m"]}],examples:["let { frobenius-norm } = import(matrix);\nfrobenius-norm([[1, 2], [3, 4]])","let { frobenius-norm } = import(matrix);\nfrobenius-norm([[1, 2, 3], [4, 5, 6], [7, 8, 9]])"],seeAlso:["matrix.one-norm","matrix.inf-norm","matrix.max-norm"]},"one-norm":{category:"matrix",description:"Calculates the **one-norm** (column norm) of a matrix.",returns:{type:"number"},args:{m:{type:"matrix",description:"The `matrix` to calculate the one-norm of."}},variants:[{argumentNames:["m"]}],examples:["let { one-norm } = import(matrix);\none-norm([[1, 2], [3, 4]])","let { one-norm } = import(matrix);\none-norm([[1, 2, 3], [4, 5, 6], [7, 8, 9]])"],seeAlso:["matrix.frobenius-norm","matrix.inf-norm","matrix.max-norm"]},"inf-norm":{category:"matrix",description:"Calculates the **infinity norm** of a matrix.",returns:{type:"number"},args:{m:{type:"matrix",description:"The `matrix` to calculate the infinity norm of."}},variants:[{argumentNames:["m"]}],examples:["let { inf-norm } = import(matrix);\ninf-norm([[1, 2], [3, 4]])","let { inf-norm } = import(matrix);\ninf-norm([[1, 2, 3], [4, 5, 6], [7, 8, 9]])"],seeAlso:["matrix.frobenius-norm","matrix.one-norm","matrix.max-norm"]},"max-norm":{category:"matrix",description:"Calculates the **max norm** of a matrix.",returns:{type:"number"},args:{m:{type:"matrix",description:"The `matrix` to calculate the max norm of."}},variants:[{argumentNames:["m"]}],examples:["let { max-norm } = import(matrix);\nmax-norm([[1, 2], [3, 4]])","let { max-norm } = import(matrix);\nmax-norm([[1, 2, 3], [4, 5, 6], [7, 8, 9]])"],seeAlso:["matrix.frobenius-norm","matrix.one-norm","matrix.inf-norm"]}};function A(t){const e=t.length,r=[];for(let a=0;a<e;a++)r[a]=[...t[a]];if(1===e)return r[0][0];if(2===e)return r[0][0]*r[1][1]-r[0][1]*r[1][0];let a=1;for(let t=0;t<e-1;t+=1){let n=t;for(let a=t+1;a<e;a+=1)Math.abs(r[a][t])>Math.abs(r[n][t])&&(n=a);if(w(r[n][t]))return 0;n!==t&&([r[t],r[n]]=[r[n],r[t]],a=-a);for(let a=t+1;a<e;a+=1){const n=r[a][t]/r[t][t];for(let i=t;i<e;i++)r[a][i]-=n*r[t][i]}}let n=a;for(let t=0;t<e;t++)n*=r[t][t];return n}function T(t,e,r){const a=t.length,n=[];for(let i=0;i<a;i++)if(i!==e){const e=[];for(let n=0;n<a;n++)n!==r&&e.push(t[i][n]);n.push(e)}return n}function $(t){const e=t.length,r=[];for(let a=0;a<e;a++){r[a]=[];for(let n=0;n<e;n++){const e=(-1)**(a+n)*A(T(t,n,a));r[a][n]=e}}return r}function C(t){return t.length===t[0].length}function M(t){if(!C(t))return!1;const e=t.length;for(let r=0;r<e;r++)for(let a=0;a<e;a++)if(r===a){if(!v(t[r][a],1))return!1}else if(!w(t[r][a]))return!1;return!0}function q(t,e){if(0===t.length||0===e.length||t[0].length!==e.length)throw new Error("Matrix dimensions do not match for multiplication");const r=t.length,a=t[0].length,n=e[0].length,i=Array(r).fill(0).map(()=>Array(n).fill(0));for(let o=0;o<r;o++)for(let r=0;r<a;r++){const a=t[o][r];for(let t=0;t<n;t++)i[o][t]+=a*e[r][t]}return i}const j={mul:{evaluate:([t,r],a)=>{y(t,a),y(r,a);try{return q(t,r)}catch(n){throw new e(`The number of columns in the first matrix must be equal to the number of rows in the second matrix, but got ${t[0].length} and ${r.length}`,a)}},arity:N(2)},det:{evaluate:([t],e)=>(b(t,e),A(t)),arity:N(1)},inv:{evaluate:([t],r)=>{b(t,r);const a=function(t){const e=t.length;if(1===e){const e=t[0][0];return w(e)?null:[[1/e]]}const r=A(t);if(w(r))return null;const a=$(t),n=[];for(let t=0;t<e;t++){n[t]=[];for(let i=0;i<e;i++)n[t][i]=a[t][i]/r}return n}(t);if(null===a)throw new e("The matrix must be invertible",r);return a},arity:N(1)},adj:{evaluate:([t],e)=>(b(t,e),$(t)),arity:N(1)},cofactor:{evaluate:([t],e)=>(b(t,e),function(t){const e=t.length,r=[];for(let a=0;a<e;a++){r[a]=[];for(let n=0;n<e;n++){const e=T(t,a,n),i=(-1)**(a+n);r[a][n]=i*A(e)}}return r}(t)),arity:N(1)},minor:{evaluate:([t,e,r],a)=>(y(t,a),s(e,a,{integer:!0,nonNegative:!0,lte:t.length}),s(r,a,{integer:!0,nonNegative:!0,lte:t[0].length}),T(t,e,r)),arity:N(3)},trace:{evaluate:([t],e)=>(b(t,e),function(t){return t.reduce((t,e,r)=>t+e[r],0)}(t)),arity:N(1)},"symmetric?":{evaluate:([t],e)=>(y(t,e),function(t){const e=t.length;if(!C(t))return!1;for(let r=0;r<e;r+=1)for(let e=0;e<r;e+=1)if(t[r][e]!==t[e][r])return!1;return!0}(t)),arity:N(1)},"triangular?":{evaluate:([t],e)=>(y(t,e),function(t){if(!C(t))return!1;const e=t.length;let r=!0,a=!0;for(let n=0;n<e;n++)for(let i=0;i<e;i++){if(n>i&&0!==t[n][i]&&(r=!1,!a))return!1;n<i&&0!==t[n][i]&&(a=!1)}return r||a}(t)),arity:N(1)},"upper-triangular?":{evaluate:([t],e)=>(y(t,e),function(t){if(!C(t))return!1;const e=t.length;for(let r=0;r<e;r++)for(let e=0;e<r;e++)if(0!==t[r][e])return!1;return!0}(t)),arity:N(1)},"lower-triangular?":{evaluate:([t],e)=>(y(t,e),function(t){if(!C(t))return!1;const e=t.length;if(!t.every(t=>t.length===e))return!1;for(let r=0;r<e;r++)for(let a=r+1;a<e;a++)if(0!==t[r][a])return!1;return!0}(t)),arity:N(1)},"diagonal?":{evaluate:([t],e)=>(y(t,e),function(t){if(!C(t))return!1;const e=t.length;for(let r=0;r<e;r+=1)for(let a=0;a<e;a+=1)if(r!==a&&0!==t[r][a])return!1;return!0}(t)),arity:N(1)},"square?":{evaluate:([t],e)=>(y(t,e),C(t)),arity:N(1)},"orthogonal-matrix?":{evaluate:([t],e)=>(y(t,e),function(t){return!!C(t)&&M(q(t,function(t){const e=[];for(let r=0;r<t[0].length;r+=1){const a=[];for(let e=0;e<t.length;e+=1)a.push(t[e][r]);e.push(a)}return e}(t)))}(t)),arity:N(1)},"identity?":{evaluate:([t],e)=>(y(t,e),M(t)),arity:N(1)},"invertible?":{evaluate:([t],e)=>(y(t,e),!!function(t){return!!h(t)&&t.length===t[0].length}(t)&&!w(A(t))),arity:N(1)},hilbert:{evaluate:([t],e)=>{s(t,e,{integer:!0,positive:!0});const r=[];for(let e=0;e<t;e+=1){const a=[];for(let r=0;r<t;r+=1)a.push(1/(e+r+1));r.push(a)}return r},arity:N(1)},vandermonde:{evaluate:([t],e)=>{g(t,e);const r=[];for(let e=0;e<t.length;e+=1){const a=[];for(let r=0;r<t.length;r+=1)a.push(t[e]**r);r.push(a)}return r},arity:N(1)},band:{evaluate:([t,e,r],a)=>(s(t,a,{integer:!0,positive:!0}),s(e,a,{integer:!0,nonNegative:!0,lt:t}),s(r,a,{integer:!0,nonNegative:!0,lte:t}),function(t,e,r){const a=Array.from({length:t},()=>Array.from({length:t},()=>0));for(let n=0;n<t;n++)for(let i=Math.max(0,n-e);i<=Math.min(t-1,n+r);i++)a[n][i]=1;return a}(t,e,r)),arity:N(3)},"banded?":{evaluate:([t,e,r],a)=>{y(t,a);const n=Math.max(t.length,t[0].length);return s(e,a,{integer:!0,nonNegative:!0,lt:n}),s(r,a,{integer:!0,nonNegative:!0,lt:n}),function(t,e,r){const a=t.length,n=t[0].length;for(let i=0;i<a;i++)for(let a=0;a<n;a++)if(0!==t[i][a]&&(i-a>e||a-i>r))return!1;return!0}(t,e,r)},arity:N(3)},rank:{evaluate:([t],e)=>{y(t,e);const[,r]=function(t){const e=t.map(t=>[...t]),r=e.length,a=e[0].length;let n=0,i=0;for(let t=0;t<a;t++){let o=-1;for(let a=i;a<r;a++)if(!w(e[a][t])){o=a;break}if(-1===o)continue;n+=1,o!==i&&([e[o],e[i]]=[e[i],e[o]]);const m=e[i][t];for(let r=t;r<a;r++)e[i][r]/=m;for(let n=0;n<r;n++)if(n!==i&&!w(e[n][t])){const r=e[n][t];for(let o=t;o<a;o++)e[n][o]-=r*e[i][o]}if(i++,i===r)break}return[e,n]}(t);return r},arity:N(1)},"frobenius-norm":{evaluate:([t],e)=>(y(t,e),Math.sqrt(t.reduce((t,e)=>t+e.reduce((t,e)=>t+e*e,0),0))),arity:N(1)},"one-norm":{evaluate:([t],e)=>(y(t,e),function(t){const e=t.length,r=t[0].length;let a=0;for(let n=0;n<r;n+=1){let r=0;for(let a=0;a<e;a+=1)r+=Math.abs(t[a][n]);a=Math.max(a,r)}return a}(t)),arity:N(1)},"inf-norm":{evaluate:([t],e)=>(y(t,e),t.reduce((t,e)=>Math.max(t,e.reduce((t,e)=>t+Math.abs(e),0)),0)),arity:N(1)},"max-norm":{evaluate:([t],e)=>(y(t,e),t.reduce((t,e)=>{const r=e.reduce((t,e)=>Math.max(t,Math.abs(e)),0);return Math.max(t,r)},0)),arity:N(1)}};for(const[t,e]of Object.entries(k))j[t]&&(j[t].docs=e);const S={name:"matrix",functions:j};exports.matrixModule=S;
+//# sourceMappingURL=matrix.js.map