@visactor/vutils 0.18.1 → 0.18.2

package/cjs/math.d.ts

@@ -61,6 +61,7 @@ export declare function crossProductPoint(dir1: {
     x: number;
     y: number;
 }): number;
+export declare function dotProduct(a: number[] | Float32Array, b: number[] | Float32Array): number;
 export declare function fuzzyEqualNumber(a: number, b: number): boolean;
 export declare function fuzzyEqualVec(a: vec2, b: vec2): boolean;
 export declare function fixPrecision(num: number, precision?: number): number;

package/cjs/math.js

@@ -30,6 +30,12 @@ function crossProductPoint(dir1, dir2) {
     return dir1.x * dir2.y - dir1.y * dir2.x;
 }
 
+function dotProduct(a, b) {
+    let ret = 0;
+    for (let i = 0; i < a.length; ++i) ret += a[i] * b[i];
+    return ret;
+}
+
 function fuzzyEqualNumber(a, b) {
     return (0, exports.abs)(a - b) < exports.epsilon;
 }
@@ -57,7 +63,7 @@ function precisionSub(a, b) {
 
 Object.defineProperty(exports, "__esModule", {
     value: !0
-}), exports.precisionSub = exports.precisionAdd = exports.getDecimalPlaces = exports.fixPrecision = exports.fuzzyEqualVec = exports.fuzzyEqualNumber = exports.crossProductPoint = exports.crossProduct = exports.lengthFromPointToLine = exports.pointAt = exports.asin = exports.acos = exports.pow = exports.sqrt = exports.sin = exports.min = exports.max = exports.cos = exports.atan2 = exports.abs = exports.pi2 = exports.SUBDIVISION_MAX_ITERATIONS = exports.SUBDIVISION_PRECISION = exports.NEWTON_MIN_SLOPE = exports.NEWTON_ITERATIONS = exports.tau = exports.halfPi = exports.pi = exports.epsilon = void 0, 
+}), exports.precisionSub = exports.precisionAdd = exports.getDecimalPlaces = exports.fixPrecision = exports.fuzzyEqualVec = exports.fuzzyEqualNumber = exports.dotProduct = exports.crossProductPoint = exports.crossProduct = exports.lengthFromPointToLine = exports.pointAt = exports.asin = exports.acos = exports.pow = exports.sqrt = exports.sin = exports.min = exports.max = exports.cos = exports.atan2 = exports.abs = exports.pi2 = exports.SUBDIVISION_MAX_ITERATIONS = exports.SUBDIVISION_PRECISION = exports.NEWTON_MIN_SLOPE = exports.NEWTON_ITERATIONS = exports.tau = exports.halfPi = exports.pi = exports.epsilon = void 0, 
 exports.epsilon = 1e-12, exports.pi = Math.PI, exports.halfPi = exports.pi / 2, 
 exports.tau = 2 * exports.pi, exports.NEWTON_ITERATIONS = 4, exports.NEWTON_MIN_SLOPE = .001, 
 exports.SUBDIVISION_PRECISION = 1e-7, exports.SUBDIVISION_MAX_ITERATIONS = 10, exports.pi2 = 2 * Math.PI, 
@@ -65,7 +71,7 @@ exports.abs = Math.abs, exports.atan2 = Math.atan2, exports.cos = Math.cos, expo
 exports.min = Math.min, exports.sin = Math.sin, exports.sqrt = Math.sqrt, exports.pow = Math.pow, 
 exports.acos = acos, exports.asin = asin, exports.pointAt = pointAt, exports.lengthFromPointToLine = lengthFromPointToLine, 
 exports.crossProduct = crossProduct, exports.crossProductPoint = crossProductPoint, 
-exports.fuzzyEqualNumber = fuzzyEqualNumber, exports.fuzzyEqualVec = fuzzyEqualVec, 
+exports.dotProduct = dotProduct, exports.fuzzyEqualNumber = fuzzyEqualNumber, exports.fuzzyEqualVec = fuzzyEqualVec, 
 exports.fixPrecision = fixPrecision, exports.getDecimalPlaces = getDecimalPlaces, 
 exports.precisionAdd = precisionAdd, exports.precisionSub = precisionSub;
 //# sourceMappingURL=math.js.map

package/cjs/math.js.map

@@ -1 +1 @@
-{"version":3,"sources":["../src/math.ts"],"names":[],"mappings":";;;AAAa,QAAA,OAAO,GAAG,KAAK,CAAC;AAChB,QAAA,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC;AACb,QAAA,MAAM,GAAG,UAAE,GAAG,CAAC,CAAC;AAChB,QAAA,GAAG,GAAG,CAAC,GAAG,UAAE,CAAC;AACb,QAAA,iBAAiB,GAAG,CAAC,CAAC;AACtB,QAAA,gBAAgB,GAAG,KAAK,CAAC;AACzB,QAAA,qBAAqB,GAAG,SAAS,CAAC;AAClC,QAAA,0BAA0B,GAAG,EAAE,CAAC;AAChC,QAAA,GAAG,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC;AA2BlB,QAAA,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC;AACf,QAAA,KAAK,GAAG,IAAI,CAAC,KAAK,CAAC;AACnB,QAAA,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC;AACf,QAAA,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC;AACf,QAAA,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC;AACf,QAAA,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC;AACf,QAAA,IAAI,GAAG,IAAI,CAAC,IAAI,CAAC;AACjB,QAAA,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC;AAE5B,SAAgB,IAAI,CAAC,CAAS;IAC5B,OAAO,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,UAAE,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;AAChD,CAAC;AAFD,oBAEC;AAED,SAAgB,IAAI,CAAC,CAAS;IAC5B,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,cAAM,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,cAAM,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;AAC5D,CAAC;AAFD,oBAEC;AAWD,SAAgB,OAAO,CACrB,EAA6B,EAC7B,EAA6B,EAC7B,EAA6B,EAC7B,EAA6B,EAC7B,CAAS;IAKT,IAAI,CAAqB,CAAC;IAC1B,IAAI,CAAqB,CAAC;IAC1B,IAAI,OAAO,EAAE,KAAK,QAAQ,IAAI,OAAO,EAAE,KAAK,QAAQ,EAAE;QACpD,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,EAAE,CAAC;KAC3B;IACD,IAAI,OAAO,EAAE,KAAK,QAAQ,IAAI,OAAO,EAAE,KAAK,QAAQ,EAAE;QACpD,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,EAAE,CAAC;KAC3B;IACD,OAAO;QACL,CAAC;QACD,CAAC;KACF,CAAC;AACJ,CAAC;AAtBD,0BAsBC;AAED,SAAgB,qBAAqB,CACnC,KAA+B,EAC/B,MAAgC,EAChC,MAAgC;IAGhC,MAAM,KAAK,GAAG,MAAM,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC;IAClC,MAAM,KAAK,GAAG,MAAM,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC;IAClC,MAAM,KAAK,GAAG,KAAK,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC;IACjC,MAAM,KAAK,GAAG,KAAK,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC;IACjC,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,GAAG,KAAK,GAAG,KAAK,GAAG,KAAK,CAAC,CAAC;IAClD,MAAM,MAAM,GAAG,IAAI,CAAC,IAAI,CAAC,KAAK,GAAG,KAAK,GAAG,KAAK,GAAG,KAAK,CAAC,CAAC;IACxD,OAAO,CAAC,GAAG,MAAM,CAAC;AACpB,CAAC;AAbD,sDAaC;AAED,SAAgB,YAAY,CAAC,IAAU,EAAE,IAAU;IACjD,OAAO,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC;AAC/C,CAAC;AAFD,oCAEC;AAED,SAAgB,iBAAiB,CAAC,IAA8B,EAAE,IAA8B;IAC9F,OAAO,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC;AAC3C,CAAC;AAFD,8CAEC;AAED,SAAgB,gBAAgB,CAAC,CAAS,EAAE,CAAS;IACnD,OAAO,IAAA,WAAG,EAAC,CAAC,GAAG,CAAC,CAAC,GAAG,eAAO,CAAC;AAC9B,CAAC;AAFD,4CAEC;AAED,SAAgB,aAAa,CAAC,CAAO,EAAE,CAAO;IAC5C,OAAO,IAAA,WAAG,EAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,IAAA,WAAG,EAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,eAAO,CAAC;AACvD,CAAC;AAFD,sCAEC;AAED,SAAgB,YAAY,CAAC,GAAW,EAAE,SAAS,GAAG,EAAE;IACtD,OAAO,IAAI,CAAC,KAAK,CAAC,GAAG,GAAG,SAAS,CAAC,GAAG,SAAS,CAAC;AACjD,CAAC;AAFD,oCAEC;AAED,SAAgB,gBAAgB,CAAC,CAAS;IACxC,MAAM,IAAI,GAAG,CAAC,CAAC,QAAQ,EAAE,CAAC,KAAK,CAAC,MAAM,CAAC,CAAC;IACxC,MAAM,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,IAAI,EAAE,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC;IACjE,OAAO,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;AACvB,CAAC;AAJD,4CAIC;AAED,SAAgB,YAAY,CAAC,CAAS,EAAE,CAAS;IAC/C,OAAO,YAAY,CAAC,CAAC,GAAG,CAAC,EAAE,EAAE,IAAI,IAAI,CAAC,GAAG,CAAC,gBAAgB,CAAC,CAAC,CAAC,EAAE,gBAAgB,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;AACvF,CAAC;AAFD,oCAEC;AAED,SAAgB,YAAY,CAAC,CAAS,EAAE,CAAS;IAC/C,OAAO,YAAY,CAAC,CAAC,GAAG,CAAC,EAAE,EAAE,IAAI,IAAI,CAAC,GAAG,CAAC,gBAAgB,CAAC,CAAC,CAAC,EAAE,gBAAgB,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;AACvF,CAAC;AAFD,oCAEC","file":"math.js","sourcesContent":["export const epsilon = 1e-12;\nexport const pi = Math.PI;\nexport const halfPi = pi / 2;\nexport const tau = 2 * pi;\nexport const NEWTON_ITERATIONS = 4;\nexport const NEWTON_MIN_SLOPE = 0.001;\nexport const SUBDIVISION_PRECISION = 0.0000001;\nexport const SUBDIVISION_MAX_ITERATIONS = 10;\nexport const pi2 = Math.PI * 2;\n\nexport type vec2 = [number, number] | Float32Array;\nexport type vec3 = [number, number, number] | Float32Array;\nexport type vec4 = [number, number, number, number] | Float32Array;\nexport type vec8 = [number, number, number, number, number, number, number, number] | Float32Array;\nexport type mat4 =\n  | [\n      number,\n      number,\n      number,\n      number,\n      number,\n      number,\n      number,\n      number,\n      number,\n      number,\n      number,\n      number,\n      number,\n      number,\n      number,\n      number\n    ]\n  | Float32Array;\n\nexport const abs = Math.abs;\nexport const atan2 = Math.atan2;\nexport const cos = Math.cos;\nexport const max = Math.max;\nexport const min = Math.min;\nexport const sin = Math.sin;\nexport const sqrt = Math.sqrt;\nexport const pow = Math.pow;\n\nexport function acos(x: number) {\n  return x > 1 ? 0 : x < -1 ? pi : Math.acos(x);\n}\n\nexport function asin(x: number) {\n  return x >= 1 ? halfPi : x <= -1 ? -halfPi : Math.asin(x);\n}\n\n/**\n * 根据比例获直线上的点\n * @param {number} x1 起始点 x\n * @param {number} y1 起始点 y\n * @param {number} x2 结束点 x\n * @param {number} y2 结束点 y\n * @param {number} t 指定比例\n * @return {object} 包含 x, y 的点\n */\nexport function pointAt(\n  x1: number | null | undefined,\n  y1: number | null | undefined,\n  x2: number | null | undefined,\n  y2: number | null | undefined,\n  t: number\n): {\n  x: undefined | number;\n  y: undefined | number;\n} {\n  let x: number | undefined;\n  let y: number | undefined;\n  if (typeof x1 === 'number' && typeof x2 === 'number') {\n    x = (1 - t) * x1 + t * x2;\n  }\n  if (typeof y1 === 'number' && typeof y2 === 'number') {\n    y = (1 - t) * y1 + t * y2;\n  }\n  return {\n    x,\n    y\n  };\n}\n\nexport function lengthFromPointToLine(\n  point: { x: number; y: number },\n  point1: { x: number; y: number },\n  point2: { x: number; y: number }\n): number {\n  // 面积\n  const dir1X = point2.x - point1.x;\n  const dir1Y = point2.y - point1.y;\n  const dir2X = point.x - point1.x;\n  const dir2Y = point.y - point1.y;\n  const s = Math.abs(dir1X * dir2Y - dir2X * dir1Y);\n  const length = Math.sqrt(dir1X * dir1X + dir1Y * dir1Y);\n  return s / length;\n}\n\nexport function crossProduct(dir1: vec2, dir2: vec2): number {\n  return dir1[0] * dir2[1] - dir1[1] * dir2[0];\n}\n\nexport function crossProductPoint(dir1: { x: number; y: number }, dir2: { x: number; y: number }): number {\n  return dir1.x * dir2.y - dir1.y * dir2.x;\n}\n\nexport function fuzzyEqualNumber(a: number, b: number): boolean {\n  return abs(a - b) < epsilon;\n}\n\nexport function fuzzyEqualVec(a: vec2, b: vec2): boolean {\n  return abs(a[0] - b[0]) + abs(a[1] - b[1]) < epsilon;\n}\n\nexport function fixPrecision(num: number, precision = 10) {\n  return Math.round(num * precision) / precision;\n}\n\nexport function getDecimalPlaces(n: number): number {\n  const dStr = n.toString().split(/[eE]/);\n  const s = (dStr[0].split('.')[1] || '').length - (+dStr[1] || 0);\n  return s > 0 ? s : 0;\n}\n\nexport function precisionAdd(a: number, b: number) {\n  return fixPrecision(a + b, 10 ** Math.max(getDecimalPlaces(a), getDecimalPlaces(b)));\n}\n\nexport function precisionSub(a: number, b: number) {\n  return fixPrecision(a - b, 10 ** Math.max(getDecimalPlaces(a), getDecimalPlaces(b)));\n}\n"]}
+{"version":3,"sources":["../src/math.ts"],"names":[],"mappings":";;;AAAa,QAAA,OAAO,GAAG,KAAK,CAAC;AAChB,QAAA,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC;AACb,QAAA,MAAM,GAAG,UAAE,GAAG,CAAC,CAAC;AAChB,QAAA,GAAG,GAAG,CAAC,GAAG,UAAE,CAAC;AACb,QAAA,iBAAiB,GAAG,CAAC,CAAC;AACtB,QAAA,gBAAgB,GAAG,KAAK,CAAC;AACzB,QAAA,qBAAqB,GAAG,SAAS,CAAC;AAClC,QAAA,0BAA0B,GAAG,EAAE,CAAC;AAChC,QAAA,GAAG,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC;AA2BlB,QAAA,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC;AACf,QAAA,KAAK,GAAG,IAAI,CAAC,KAAK,CAAC;AACnB,QAAA,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC;AACf,QAAA,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC;AACf,QAAA,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC;AACf,QAAA,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC;AACf,QAAA,IAAI,GAAG,IAAI,CAAC,IAAI,CAAC;AACjB,QAAA,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC;AAE5B,SAAgB,IAAI,CAAC,CAAS;IAC5B,OAAO,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,UAAE,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;AAChD,CAAC;AAFD,oBAEC;AAED,SAAgB,IAAI,CAAC,CAAS;IAC5B,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,cAAM,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,cAAM,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;AAC5D,CAAC;AAFD,oBAEC;AAWD,SAAgB,OAAO,CACrB,EAA6B,EAC7B,EAA6B,EAC7B,EAA6B,EAC7B,EAA6B,EAC7B,CAAS;IAKT,IAAI,CAAqB,CAAC;IAC1B,IAAI,CAAqB,CAAC;IAC1B,IAAI,OAAO,EAAE,KAAK,QAAQ,IAAI,OAAO,EAAE,KAAK,QAAQ,EAAE;QACpD,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,EAAE,CAAC;KAC3B;IACD,IAAI,OAAO,EAAE,KAAK,QAAQ,IAAI,OAAO,EAAE,KAAK,QAAQ,EAAE;QACpD,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,EAAE,CAAC;KAC3B;IACD,OAAO;QACL,CAAC;QACD,CAAC;KACF,CAAC;AACJ,CAAC;AAtBD,0BAsBC;AAED,SAAgB,qBAAqB,CACnC,KAA+B,EAC/B,MAAgC,EAChC,MAAgC;IAGhC,MAAM,KAAK,GAAG,MAAM,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC;IAClC,MAAM,KAAK,GAAG,MAAM,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC;IAClC,MAAM,KAAK,GAAG,KAAK,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC;IACjC,MAAM,KAAK,GAAG,KAAK,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC;IACjC,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,GAAG,KAAK,GAAG,KAAK,GAAG,KAAK,CAAC,CAAC;IAClD,MAAM,MAAM,GAAG,IAAI,CAAC,IAAI,CAAC,KAAK,GAAG,KAAK,GAAG,KAAK,GAAG,KAAK,CAAC,CAAC;IACxD,OAAO,CAAC,GAAG,MAAM,CAAC;AACpB,CAAC;AAbD,sDAaC;AAED,SAAgB,YAAY,CAAC,IAAU,EAAE,IAAU;IACjD,OAAO,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC;AAC/C,CAAC;AAFD,oCAEC;AAED,SAAgB,iBAAiB,CAAC,IAA8B,EAAE,IAA8B;IAC9F,OAAO,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC;AAC3C,CAAC;AAFD,8CAEC;AAED,SAAgB,UAAU,CAAC,CAA0B,EAAE,CAA0B;IAC/E,IAAI,GAAG,GAAG,CAAC,CAAC;IACZ,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,MAAM,EAAE,EAAE,CAAC,EAAE;QACjC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;KACpB;IACD,OAAO,GAAG,CAAC;AACb,CAAC;AAND,gCAMC;AAED,SAAgB,gBAAgB,CAAC,CAAS,EAAE,CAAS;IACnD,OAAO,IAAA,WAAG,EAAC,CAAC,GAAG,CAAC,CAAC,GAAG,eAAO,CAAC;AAC9B,CAAC;AAFD,4CAEC;AAED,SAAgB,aAAa,CAAC,CAAO,EAAE,CAAO;IAC5C,OAAO,IAAA,WAAG,EAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,IAAA,WAAG,EAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,eAAO,CAAC;AACvD,CAAC;AAFD,sCAEC;AAED,SAAgB,YAAY,CAAC,GAAW,EAAE,SAAS,GAAG,EAAE;IACtD,OAAO,IAAI,CAAC,KAAK,CAAC,GAAG,GAAG,SAAS,CAAC,GAAG,SAAS,CAAC;AACjD,CAAC;AAFD,oCAEC;AAED,SAAgB,gBAAgB,CAAC,CAAS;IACxC,MAAM,IAAI,GAAG,CAAC,CAAC,QAAQ,EAAE,CAAC,KAAK,CAAC,MAAM,CAAC,CAAC;IACxC,MAAM,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,IAAI,EAAE,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC;IACjE,OAAO,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;AACvB,CAAC;AAJD,4CAIC;AAED,SAAgB,YAAY,CAAC,CAAS,EAAE,CAAS;IAC/C,OAAO,YAAY,CAAC,CAAC,GAAG,CAAC,EAAE,EAAE,IAAI,IAAI,CAAC,GAAG,CAAC,gBAAgB,CAAC,CAAC,CAAC,EAAE,gBAAgB,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;AACvF,CAAC;AAFD,oCAEC;AAED,SAAgB,YAAY,CAAC,CAAS,EAAE,CAAS;IAC/C,OAAO,YAAY,CAAC,CAAC,GAAG,CAAC,EAAE,EAAE,IAAI,IAAI,CAAC,GAAG,CAAC,gBAAgB,CAAC,CAAC,CAAC,EAAE,gBAAgB,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;AACvF,CAAC;AAFD,oCAEC","file":"math.js","sourcesContent":["export const epsilon = 1e-12;\nexport const pi = Math.PI;\nexport const halfPi = pi / 2;\nexport const tau = 2 * pi;\nexport const NEWTON_ITERATIONS = 4;\nexport const NEWTON_MIN_SLOPE = 0.001;\nexport const SUBDIVISION_PRECISION = 0.0000001;\nexport const SUBDIVISION_MAX_ITERATIONS = 10;\nexport const pi2 = Math.PI * 2;\n\nexport type vec2 = [number, number] | Float32Array;\nexport type vec3 = [number, number, number] | Float32Array;\nexport type vec4 = [number, number, number, number] | Float32Array;\nexport type vec8 = [number, number, number, number, number, number, number, number] | Float32Array;\nexport type mat4 =\n  | [\n      number,\n      number,\n      number,\n      number,\n      number,\n      number,\n      number,\n      number,\n      number,\n      number,\n      number,\n      number,\n      number,\n      number,\n      number,\n      number\n    ]\n  | Float32Array;\n\nexport const abs = Math.abs;\nexport const atan2 = Math.atan2;\nexport const cos = Math.cos;\nexport const max = Math.max;\nexport const min = Math.min;\nexport const sin = Math.sin;\nexport const sqrt = Math.sqrt;\nexport const pow = Math.pow;\n\nexport function acos(x: number) {\n  return x > 1 ? 0 : x < -1 ? pi : Math.acos(x);\n}\n\nexport function asin(x: number) {\n  return x >= 1 ? halfPi : x <= -1 ? -halfPi : Math.asin(x);\n}\n\n/**\n * 根据比例获直线上的点\n * @param {number} x1 起始点 x\n * @param {number} y1 起始点 y\n * @param {number} x2 结束点 x\n * @param {number} y2 结束点 y\n * @param {number} t 指定比例\n * @return {object} 包含 x, y 的点\n */\nexport function pointAt(\n  x1: number | null | undefined,\n  y1: number | null | undefined,\n  x2: number | null | undefined,\n  y2: number | null | undefined,\n  t: number\n): {\n  x: undefined | number;\n  y: undefined | number;\n} {\n  let x: number | undefined;\n  let y: number | undefined;\n  if (typeof x1 === 'number' && typeof x2 === 'number') {\n    x = (1 - t) * x1 + t * x2;\n  }\n  if (typeof y1 === 'number' && typeof y2 === 'number') {\n    y = (1 - t) * y1 + t * y2;\n  }\n  return {\n    x,\n    y\n  };\n}\n\nexport function lengthFromPointToLine(\n  point: { x: number; y: number },\n  point1: { x: number; y: number },\n  point2: { x: number; y: number }\n): number {\n  // 面积\n  const dir1X = point2.x - point1.x;\n  const dir1Y = point2.y - point1.y;\n  const dir2X = point.x - point1.x;\n  const dir2Y = point.y - point1.y;\n  const s = Math.abs(dir1X * dir2Y - dir2X * dir1Y);\n  const length = Math.sqrt(dir1X * dir1X + dir1Y * dir1Y);\n  return s / length;\n}\n\nexport function crossProduct(dir1: vec2, dir2: vec2): number {\n  return dir1[0] * dir2[1] - dir1[1] * dir2[0];\n}\n\nexport function crossProductPoint(dir1: { x: number; y: number }, dir2: { x: number; y: number }): number {\n  return dir1.x * dir2.y - dir1.y * dir2.x;\n}\n\nexport function dotProduct(a: number[] | Float32Array, b: number[] | Float32Array): number {\n  let ret = 0;\n  for (let i = 0; i < a.length; ++i) {\n    ret += a[i] * b[i];\n  }\n  return ret;\n}\n\nexport function fuzzyEqualNumber(a: number, b: number): boolean {\n  return abs(a - b) < epsilon;\n}\n\nexport function fuzzyEqualVec(a: vec2, b: vec2): boolean {\n  return abs(a[0] - b[0]) + abs(a[1] - b[1]) < epsilon;\n}\n\nexport function fixPrecision(num: number, precision = 10) {\n  return Math.round(num * precision) / precision;\n}\n\nexport function getDecimalPlaces(n: number): number {\n  const dStr = n.toString().split(/[eE]/);\n  const s = (dStr[0].split('.')[1] || '').length - (+dStr[1] || 0);\n  return s > 0 ? s : 0;\n}\n\nexport function precisionAdd(a: number, b: number) {\n  return fixPrecision(a + b, 10 ** Math.max(getDecimalPlaces(a), getDecimalPlaces(b)));\n}\n\nexport function precisionSub(a: number, b: number) {\n  return fixPrecision(a - b, 10 ** Math.max(getDecimalPlaces(a), getDecimalPlaces(b)));\n}\n"]}

package/dist/index.js

@@ -1664,6 +1664,13 @@
 	function crossProductPoint(dir1, dir2) {
 	    return dir1.x * dir2.y - dir1.y * dir2.x;
 	}
+	function dotProduct(a, b) {
+	    let ret = 0;
+	    for (let i = 0; i < a.length; ++i) {
+	        ret += a[i] * b[i];
+	    }
+	    return ret;
+	}
 	function fuzzyEqualNumber(a, b) {
 	    return abs(a - b) < epsilon;
 	}
@@ -4784,6 +4791,163 @@
 	    return { x: lng, y: lat };
 	}
 
+	const SMALL = 1e-10;
+
+	function intersectionArea(circles, stats) {
+	    const intersectionPoints = getIntersectionPoints(circles);
+	    const innerPoints = intersectionPoints.filter(function (p) {
+	        return containedInCircles(p, circles);
+	    });
+	    let arcArea = 0;
+	    let polygonArea = 0;
+	    const arcs = [];
+	    if (innerPoints.length > 1) {
+	        const center = getCenter(innerPoints);
+	        for (let i = 0; i < innerPoints.length; ++i) {
+	            const p = innerPoints[i];
+	            p.angle = Math.atan2(p.x - center.x, p.y - center.y);
+	        }
+	        innerPoints.sort(function (a, b) {
+	            return b.angle - a.angle;
+	        });
+	        let p2 = innerPoints[innerPoints.length - 1];
+	        for (let i = 0; i < innerPoints.length; ++i) {
+	            const p1 = innerPoints[i];
+	            polygonArea += (p2.x + p1.x) * (p1.y - p2.y);
+	            const midPoint = { x: (p1.x + p2.x) / 2, y: (p1.y + p2.y) / 2 };
+	            let arc = null;
+	            for (let j = 0; j < p1.parentIndex.length; ++j) {
+	                if (p2.parentIndex.indexOf(p1.parentIndex[j]) > -1) {
+	                    const circle = circles[p1.parentIndex[j]];
+	                    const a1 = Math.atan2(p1.x - circle.x, p1.y - circle.y);
+	                    const a2 = Math.atan2(p2.x - circle.x, p2.y - circle.y);
+	                    let angleDiff = a2 - a1;
+	                    if (angleDiff < 0) {
+	                        angleDiff += 2 * Math.PI;
+	                    }
+	                    const a = a2 - angleDiff / 2;
+	                    let width = PointService.distancePP(midPoint, {
+	                        x: circle.x + circle.radius * Math.sin(a),
+	                        y: circle.y + circle.radius * Math.cos(a)
+	                    });
+	                    if (width > circle.radius * 2) {
+	                        width = circle.radius * 2;
+	                    }
+	                    if (arc === null || arc.width > width) {
+	                        arc = { circle: circle, width: width, p1: p1, p2: p2 };
+	                    }
+	                }
+	            }
+	            if (arc !== null) {
+	                arcs.push(arc);
+	                arcArea += circleArea(arc.circle.radius, arc.width);
+	                p2 = p1;
+	            }
+	        }
+	    }
+	    else {
+	        let smallest = circles[0];
+	        for (let i = 1; i < circles.length; ++i) {
+	            if (circles[i].radius < smallest.radius) {
+	                smallest = circles[i];
+	            }
+	        }
+	        let disjoint = false;
+	        for (let i = 0; i < circles.length; ++i) {
+	            if (PointService.distancePP(circles[i], smallest) > Math.abs(smallest.radius - circles[i].radius)) {
+	                disjoint = true;
+	                break;
+	            }
+	        }
+	        if (disjoint) {
+	            arcArea = polygonArea = 0;
+	        }
+	        else {
+	            arcArea = smallest.radius * smallest.radius * Math.PI;
+	            arcs.push({
+	                circle: smallest,
+	                p1: { x: smallest.x, y: smallest.y + smallest.radius },
+	                p2: { x: smallest.x - SMALL, y: smallest.y + smallest.radius },
+	                width: smallest.radius * 2
+	            });
+	        }
+	    }
+	    polygonArea /= 2;
+	    if (stats) {
+	        stats.area = arcArea + polygonArea;
+	        stats.arcArea = arcArea;
+	        stats.polygonArea = polygonArea;
+	        stats.arcs = arcs;
+	        stats.innerPoints = innerPoints;
+	        stats.intersectionPoints = intersectionPoints;
+	    }
+	    return arcArea + polygonArea;
+	}
+	function containedInCircles(point, circles) {
+	    for (let i = 0; i < circles.length; ++i) {
+	        if (PointService.distancePP(point, circles[i]) > circles[i].radius + SMALL) {
+	            return false;
+	        }
+	    }
+	    return true;
+	}
+	function getIntersectionPoints(circles) {
+	    const ret = [];
+	    for (let i = 0; i < circles.length; ++i) {
+	        for (let j = i + 1; j < circles.length; ++j) {
+	            const intersect = circleCircleIntersection(circles[i], circles[j]);
+	            for (let k = 0; k < intersect.length; ++k) {
+	                const p = intersect[k];
+	                p.parentIndex = [i, j];
+	                ret.push(p);
+	            }
+	        }
+	    }
+	    return ret;
+	}
+	function circleArea(r, width) {
+	    return r * r * Math.acos(1 - width / r) - (r - width) * Math.sqrt(width * (2 * r - width));
+	}
+	function circleOverlap(r1, r2, d) {
+	    if (d >= r1 + r2) {
+	        return 0;
+	    }
+	    if (d <= Math.abs(r1 - r2)) {
+	        return Math.PI * Math.min(r1, r2) * Math.min(r1, r2);
+	    }
+	    const w1 = r1 - (d * d - r2 * r2 + r1 * r1) / (2 * d);
+	    const w2 = r2 - (d * d - r1 * r1 + r2 * r2) / (2 * d);
+	    return circleArea(r1, w1) + circleArea(r2, w2);
+	}
+	function circleCircleIntersection(p1, p2) {
+	    const d = PointService.distancePP(p1, p2);
+	    const r1 = p1.radius;
+	    const r2 = p2.radius;
+	    if (d >= r1 + r2 || d <= Math.abs(r1 - r2)) {
+	        return [];
+	    }
+	    const a = (r1 * r1 - r2 * r2 + d * d) / (2 * d);
+	    const h = Math.sqrt(r1 * r1 - a * a);
+	    const x0 = p1.x + (a * (p2.x - p1.x)) / d;
+	    const y0 = p1.y + (a * (p2.y - p1.y)) / d;
+	    const rx = -(p2.y - p1.y) * (h / d);
+	    const ry = -(p2.x - p1.x) * (h / d);
+	    return [
+	        { x: x0 + rx, y: y0 - ry },
+	        { x: x0 - rx, y: y0 + ry }
+	    ];
+	}
+	function getCenter(points) {
+	    const center = { x: 0, y: 0 };
+	    for (let i = 0; i < points.length; ++i) {
+	        center.x += points[i].x;
+	        center.y += points[i].y;
+	    }
+	    center.x /= points.length;
+	    center.y /= points.length;
+	    return center;
+	}
+
 	class TimeUtil {
 	    static getInstance() {
 	        if (!TimeUtil.instance) {
@@ -5462,6 +5626,277 @@
 	    };
 	}
 
+	function zeros(x) {
+	    const r = new Array(x);
+	    for (let i = 0; i < x; ++i) {
+	        r[i] = 0;
+	    }
+	    return r;
+	}
+	function zerosM(x, y) {
+	    return zeros(x).map(function () {
+	        return zeros(y);
+	    });
+	}
+	function norm2(a) {
+	    return Math.sqrt(dotProduct(a, a));
+	}
+	function scale(ret, value, c) {
+	    for (let i = 0; i < value.length; ++i) {
+	        ret[i] = value[i] * c;
+	    }
+	}
+	function weightedSum(ret, w1, v1, w2, v2) {
+	    for (let j = 0; j < ret.length; ++j) {
+	        ret[j] = w1 * v1[j] + w2 * v2[j];
+	    }
+	}
+	function gemv(output, A, x) {
+	    for (let i = 0; i < output.length; ++i) {
+	        output[i] = dotProduct(A[i], x);
+	    }
+	}
+
+	function nelderMead(f, x0, parameters) {
+	    parameters = parameters || {};
+	    const maxIterations = parameters.maxIterations || x0.length * 200;
+	    const nonZeroDelta = parameters.nonZeroDelta || 1.05;
+	    const zeroDelta = parameters.zeroDelta || 0.001;
+	    const minErrorDelta = parameters.minErrorDelta || 1e-6;
+	    const minTolerance = parameters.minErrorDelta || 1e-5;
+	    const rho = parameters.rho !== undefined ? parameters.rho : 1;
+	    const chi = parameters.chi !== undefined ? parameters.chi : 2;
+	    const psi = parameters.psi !== undefined ? parameters.psi : -0.5;
+	    const sigma = parameters.sigma !== undefined ? parameters.sigma : 0.5;
+	    let maxDiff;
+	    const N = x0.length;
+	    const simplex = new Array(N + 1);
+	    simplex[0] = x0;
+	    simplex[0].fx = f(x0);
+	    simplex[0].id = 0;
+	    for (let i = 0; i < N; ++i) {
+	        const point = x0.slice();
+	        point[i] = point[i] ? point[i] * nonZeroDelta : zeroDelta;
+	        simplex[i + 1] = point;
+	        simplex[i + 1].fx = f(point);
+	        simplex[i + 1].id = i + 1;
+	    }
+	    function updateSimplex(value) {
+	        for (let i = 0; i < value.length; i++) {
+	            simplex[N][i] = value[i];
+	        }
+	        simplex[N].fx = value.fx;
+	    }
+	    const sortOrder = function (a, b) {
+	        return a.fx - b.fx;
+	    };
+	    const centroid = x0.slice();
+	    const reflected = x0.slice();
+	    const contracted = x0.slice();
+	    const expanded = x0.slice();
+	    for (let iteration = 0; iteration < maxIterations; ++iteration) {
+	        simplex.sort(sortOrder);
+	        if (parameters.history) {
+	            const sortedSimplex = simplex.map(function (x) {
+	                const state = x.slice();
+	                state.fx = x.fx;
+	                state.id = x.id;
+	                return state;
+	            });
+	            sortedSimplex.sort(function (a, b) {
+	                return a.id - b.id;
+	            });
+	            parameters.history.push({ x: simplex[0].slice(), fx: simplex[0].fx, simplex: sortedSimplex });
+	        }
+	        maxDiff = 0;
+	        for (let i = 0; i < N; ++i) {
+	            maxDiff = Math.max(maxDiff, Math.abs(simplex[0][i] - simplex[1][i]));
+	        }
+	        if (Math.abs(simplex[0].fx - simplex[N].fx) < minErrorDelta && maxDiff < minTolerance) {
+	            break;
+	        }
+	        for (let i = 0; i < N; ++i) {
+	            centroid[i] = 0;
+	            for (let j = 0; j < N; ++j) {
+	                centroid[i] += simplex[j][i];
+	            }
+	            centroid[i] /= N;
+	        }
+	        const worst = simplex[N];
+	        weightedSum(reflected, 1 + rho, centroid, -rho, worst);
+	        reflected.fx = f(reflected);
+	        if (reflected.fx < simplex[0].fx) {
+	            weightedSum(expanded, 1 + chi, centroid, -chi, worst);
+	            expanded.fx = f(expanded);
+	            if (expanded.fx < reflected.fx) {
+	                updateSimplex(expanded);
+	            }
+	            else {
+	                updateSimplex(reflected);
+	            }
+	        }
+	        else if (reflected.fx >= simplex[N - 1].fx) {
+	            let shouldReduce = false;
+	            if (reflected.fx > worst.fx) {
+	                weightedSum(contracted, 1 + psi, centroid, -psi, worst);
+	                contracted.fx = f(contracted);
+	                if (contracted.fx < worst.fx) {
+	                    updateSimplex(contracted);
+	                }
+	                else {
+	                    shouldReduce = true;
+	                }
+	            }
+	            else {
+	                weightedSum(contracted, 1 - psi * rho, centroid, psi * rho, worst);
+	                contracted.fx = f(contracted);
+	                if (contracted.fx < reflected.fx) {
+	                    updateSimplex(contracted);
+	                }
+	                else {
+	                    shouldReduce = true;
+	                }
+	            }
+	            if (shouldReduce) {
+	                if (sigma >= 1) {
+	                    break;
+	                }
+	                for (let i = 1; i < simplex.length; ++i) {
+	                    weightedSum(simplex[i], 1 - sigma, simplex[0], sigma, simplex[i]);
+	                    simplex[i].fx = f(simplex[i]);
+	                }
+	            }
+	        }
+	        else {
+	            updateSimplex(reflected);
+	        }
+	    }
+	    simplex.sort(sortOrder);
+	    return { fx: simplex[0].fx, x: simplex[0] };
+	}
+
+	function wolfeLineSearch(f, pk, current, next, a, c1, c2) {
+	    const phi0 = current.fx;
+	    const phiPrime0 = dotProduct(current.fxprime, pk);
+	    let phi = phi0;
+	    let phi_old = phi0;
+	    let phiPrime = phiPrime0;
+	    let a0 = 0;
+	    a = a || 1;
+	    c1 = c1 || 1e-6;
+	    c2 = c2 || 0.1;
+	    function zoom(a_lo, a_high, phi_lo) {
+	        for (let iteration = 0; iteration < 16; ++iteration) {
+	            a = (a_lo + a_high) / 2;
+	            weightedSum(next.x, 1.0, current.x, a, pk);
+	            phi = next.fx = f(next.x, next.fxprime);
+	            phiPrime = dotProduct(next.fxprime, pk);
+	            if (phi > phi0 + c1 * a * phiPrime0 || phi >= phi_lo) {
+	                a_high = a;
+	            }
+	            else {
+	                if (Math.abs(phiPrime) <= -c2 * phiPrime0) {
+	                    return a;
+	                }
+	                if (phiPrime * (a_high - a_lo) >= 0) {
+	                    a_high = a_lo;
+	                }
+	                a_lo = a;
+	                phi_lo = phi;
+	            }
+	        }
+	        return 0;
+	    }
+	    for (let iteration = 0; iteration < 10; ++iteration) {
+	        weightedSum(next.x, 1.0, current.x, a, pk);
+	        phi = next.fx = f(next.x, next.fxprime);
+	        phiPrime = dotProduct(next.fxprime, pk);
+	        if (phi > phi0 + c1 * a * phiPrime0 || (iteration && phi >= phi_old)) {
+	            return zoom(a0, a, phi_old);
+	        }
+	        if (Math.abs(phiPrime) <= -c2 * phiPrime0) {
+	            return a;
+	        }
+	        if (phiPrime >= 0) {
+	            return zoom(a, a0, phi);
+	        }
+	        phi_old = phi;
+	        a0 = a;
+	        a *= 2;
+	    }
+	    return a;
+	}
+
+	function conjugateGradient(f, initial, params) {
+	    let current = { x: initial.slice(), fx: 0, fxprime: initial.slice() };
+	    let next = { x: initial.slice(), fx: 0, fxprime: initial.slice() };
+	    const yk = initial.slice();
+	    let temp;
+	    let a = 1;
+	    params = params || {};
+	    const maxIterations = params.maxIterations || initial.length * 20;
+	    current.fx = f(current.x, current.fxprime);
+	    const pk = current.fxprime.slice();
+	    scale(pk, current.fxprime, -1);
+	    for (let i = 0; i < maxIterations; ++i) {
+	        a = wolfeLineSearch(f, pk, current, next, a);
+	        if (params.history) {
+	            params.history.push({ x: current.x.slice(), fx: current.fx, fxprime: current.fxprime.slice(), alpha: a });
+	        }
+	        if (!a) {
+	            scale(pk, current.fxprime, -1);
+	        }
+	        else {
+	            weightedSum(yk, 1, next.fxprime, -1, current.fxprime);
+	            const delta_k = dotProduct(current.fxprime, current.fxprime);
+	            const beta_k = Math.max(0, dotProduct(yk, next.fxprime) / delta_k);
+	            weightedSum(pk, beta_k, pk, -1, next.fxprime);
+	            temp = current;
+	            current = next;
+	            next = temp;
+	        }
+	        if (norm2(current.fxprime) <= 1e-5) {
+	            break;
+	        }
+	    }
+	    if (params.history) {
+	        params.history.push({ x: current.x.slice(), fx: current.fx, fxprime: current.fxprime.slice(), alpha: a });
+	    }
+	    return current;
+	}
+
+	function findZeroOfFunction(f, a, b, parameters) {
+	    parameters = parameters || {};
+	    const maxIterations = parameters.maxIterations || 100;
+	    const tolerance = parameters.tolerance || 1e-10;
+	    const fA = f(a);
+	    const fB = f(b);
+	    let delta = b - a;
+	    if (fA * fB > 0) {
+	        const logger = Logger.getInstance();
+	        logger.error('Initial bisect points must have opposite signs');
+	        return NaN;
+	    }
+	    if (fA === 0) {
+	        return a;
+	    }
+	    if (fB === 0) {
+	        return b;
+	    }
+	    for (let i = 0; i < maxIterations; ++i) {
+	        delta /= 2;
+	        const mid = a + delta;
+	        const fMid = f(mid);
+	        if (fMid * fA >= 0) {
+	            a = mid;
+	        }
+	        if (Math.abs(delta) < tolerance || fMid === 0) {
+	            return mid;
+	        }
+	    }
+	    return a + delta;
+	}
+
 	exports.AABBBounds = AABBBounds;
 	exports.Bounds = Bounds;
 	exports.Color = Color;
@@ -5488,6 +5923,7 @@
 	exports.PolarPoint = PolarPoint;
 	exports.RGB = RGB;
 	exports.SECOND = SECOND;
+	exports.SMALL = SMALL;
 	exports.SUBDIVISION_MAX_ITERATIONS = SUBDIVISION_MAX_ITERATIONS;
 	exports.SUBDIVISION_PRECISION = SUBDIVISION_PRECISION;
 	exports.TextMeasure = TextMeasure;
@@ -5501,6 +5937,9 @@
 	exports.asin = asin;
 	exports.atan2 = atan2;
 	exports.bisect = bisect;
+	exports.circleArea = circleArea;
+	exports.circleCircleIntersection = circleCircleIntersection;
+	exports.circleOverlap = circleOverlap;
 	exports.clamp = clamp;
 	exports.clampAngleByDegree = clampAngleByDegree;
 	exports.clampAngleByRadian = clampAngleByRadian;
@@ -5510,7 +5949,9 @@
 	exports.clamper = clamper;
 	exports.clone = clone;
 	exports.cloneDeep = cloneDeep;
+	exports.conjugateGradient = conjugateGradient;
 	exports.constant = constant;
+	exports.containedInCircles = containedInCircles;
 	exports.cos = cos;
 	exports.crossProduct = crossProduct;
 	exports.crossProductPoint = crossProductPoint;
@@ -5525,9 +5966,11 @@
 	exports.degreeToRadian = degreeToRadian;
 	exports.destination = destination;
 	exports.deviation = deviation;
+	exports.dotProduct = dotProduct;
 	exports.eastAsianCharacterInfo = eastAsianCharacterInfo;
 	exports.epsilon = epsilon;
 	exports.exponent = exponent;
+	exports.findZeroOfFunction = findZeroOfFunction;
 	exports.fixPrecision = fixPrecision;
 	exports.flattenArray = flattenArray;
 	exports.formatNumerals = formatNumerals;
@@ -5537,12 +5980,14 @@
 	exports.fullYearSetterName = fullYearSetterName;
 	exports.fuzzyEqualNumber = fuzzyEqualNumber;
 	exports.fuzzyEqualVec = fuzzyEqualVec;
+	exports.gemv = gemv;
 	exports.generateCeil = generateCeil;
 	exports.generateCount = generateCount;
 	exports.generateStepInterval = generateStepInterval;
 	exports.get = get;
 	exports.getAABBFromPoints = getAABBFromPoints;
 	exports.getAngleByPoint = getAngleByPoint;
+	exports.getCenter = getCenter;
 	exports.getContainerSize = getContainerSize;
 	exports.getContextFont = getContextFont;
 	exports.getDecimalPlaces = getDecimalPlaces;
@@ -5575,6 +6020,7 @@
 	exports.interpolateNumberRound = interpolateNumberRound;
 	exports.interpolateRgb = interpolateRgb;
 	exports.interpolateString = interpolateString;
+	exports.intersectionArea = intersectionArea;
 	exports.isArray = isArray;
 	exports.isArrayLike = isArrayLike;
 	exports.isBase64 = isBase64;
@@ -5637,6 +6083,8 @@
 	exports.monthGetterName = monthGetterName;
 	exports.monthOffset = monthOffset;
 	exports.monthSetterName = monthSetterName;
+	exports.nelderMead = nelderMead;
+	exports.norm2 = norm2;
 	exports.normalTransform = normalTransform;
 	exports.normalizePadding = normalizePadding;
 	exports.numberSpecifierReg = numberSpecifierReg;
@@ -5664,6 +6112,7 @@
 	exports.rectInsideAnotherRect = rectInsideAnotherRect;
 	exports.rgbToHex = rgbToHex;
 	exports.rgbToHsl = rgbToHsl;
+	exports.scale = scale;
 	exports.secondCount = secondCount;
 	exports.secondField = secondField;
 	exports.secondFloor = secondFloor;
@@ -5711,9 +6160,12 @@
 	exports.utcYearOffset = utcYearOffset;
 	exports.uuid = uuid;
 	exports.variance = variance;
+	exports.weightedSum = weightedSum;
 	exports.yearCount = yearCount;
 	exports.yearField = yearField;
 	exports.yearFloor = yearFloor;
 	exports.yearOffset = yearOffset;
+	exports.zeros = zeros;
+	exports.zerosM = zerosM;
 
 }));