@performant-software/geospatial 3.1.25-beta.1 → 3.1.25-beta.2

package/dist/index.cjs548.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"index.cjs548.js","sources":["../../../node_modules/proj4/lib/datumUtils.js"],"sourcesContent":["'use strict';\nimport { PJD_3PARAM, PJD_7PARAM, HALF_PI } from './constants/values';\nexport function compareDatums(source, dest) {\n  if (source.datum_type !== dest.datum_type) {\n    return false; // false, datums are not equal\n  } else if (source.a !== dest.a || Math.abs(source.es - dest.es) > 0.000000000050) {\n    // the tolerance for es is to ensure that GRS80 and WGS84\n    // are considered identical\n    return false;\n  } else if (source.datum_type === PJD_3PARAM) {\n    return (source.datum_params[0] === dest.datum_params[0] && source.datum_params[1] === dest.datum_params[1] && source.datum_params[2] === dest.datum_params[2]);\n  } else if (source.datum_type === PJD_7PARAM) {\n    return (source.datum_params[0] === dest.datum_params[0] && source.datum_params[1] === dest.datum_params[1] && source.datum_params[2] === dest.datum_params[2] && source.datum_params[3] === dest.datum_params[3] && source.datum_params[4] === dest.datum_params[4] && source.datum_params[5] === dest.datum_params[5] && source.datum_params[6] === dest.datum_params[6]);\n  } else {\n    return true; // datums are equal\n  }\n} // cs_compare_datums()\n\n/*\n * The function Convert_Geodetic_To_Geocentric converts geodetic coordinates\n * (latitude, longitude, and height) to geocentric coordinates (X, Y, Z),\n * according to the current ellipsoid parameters.\n *\n *    Latitude  : Geodetic latitude in radians                     (input)\n *    Longitude : Geodetic longitude in radians                    (input)\n *    Height    : Geodetic height, in meters                       (input)\n *    X         : Calculated Geocentric X coordinate, in meters    (output)\n *    Y         : Calculated Geocentric Y coordinate, in meters    (output)\n *    Z         : Calculated Geocentric Z coordinate, in meters    (output)\n *\n */\nexport function geodeticToGeocentric(p, es, a) {\n  var Longitude = p.x;\n  var Latitude = p.y;\n  var Height = p.z ? p.z : 0; // Z value not always supplied\n\n  var Rn; /*  Earth radius at location  */\n  var Sin_Lat; /*  Math.sin(Latitude)  */\n  var Sin2_Lat; /*  Square of Math.sin(Latitude)  */\n  var Cos_Lat; /*  Math.cos(Latitude)  */\n\n  /*\n   ** Don't blow up if Latitude is just a little out of the value\n   ** range as it may just be a rounding issue.  Also removed longitude\n   ** test, it should be wrapped by Math.cos() and Math.sin().  NFW for PROJ.4, Sep/2001.\n   */\n  if (Latitude < -HALF_PI && Latitude > -1.001 * HALF_PI) {\n    Latitude = -HALF_PI;\n  } else if (Latitude > HALF_PI && Latitude < 1.001 * HALF_PI) {\n    Latitude = HALF_PI;\n  } else if (Latitude < -HALF_PI) {\n    /* Latitude out of range */\n    // ..reportError('geocent:lat out of range:' + Latitude);\n    return { x: -Infinity, y: -Infinity, z: p.z };\n  } else if (Latitude > HALF_PI) {\n    /* Latitude out of range */\n    return { x: Infinity, y: Infinity, z: p.z };\n  }\n\n  if (Longitude > Math.PI) {\n    Longitude -= (2 * Math.PI);\n  }\n  Sin_Lat = Math.sin(Latitude);\n  Cos_Lat = Math.cos(Latitude);\n  Sin2_Lat = Sin_Lat * Sin_Lat;\n  Rn = a / (Math.sqrt(1.0e0 - es * Sin2_Lat));\n  return {\n    x: (Rn + Height) * Cos_Lat * Math.cos(Longitude),\n    y: (Rn + Height) * Cos_Lat * Math.sin(Longitude),\n    z: ((Rn * (1 - es)) + Height) * Sin_Lat\n  };\n} // cs_geodetic_to_geocentric()\n\nexport function geocentricToGeodetic(p, es, a, b) {\n  /* local defintions and variables */\n  /* end-criterium of loop, accuracy of sin(Latitude) */\n  var genau = 1e-12;\n  var genau2 = (genau * genau);\n  var maxiter = 30;\n\n  var P; /* distance between semi-minor axis and location */\n  var RR; /* distance between center and location */\n  var CT; /* sin of geocentric latitude */\n  var ST; /* cos of geocentric latitude */\n  var RX;\n  var RK;\n  var RN; /* Earth radius at location */\n  var CPHI0; /* cos of start or old geodetic latitude in iterations */\n  var SPHI0; /* sin of start or old geodetic latitude in iterations */\n  var CPHI; /* cos of searched geodetic latitude */\n  var SPHI; /* sin of searched geodetic latitude */\n  var SDPHI; /* end-criterium: addition-theorem of sin(Latitude(iter)-Latitude(iter-1)) */\n  var iter; /* # of continous iteration, max. 30 is always enough (s.a.) */\n\n  var X = p.x;\n  var Y = p.y;\n  var Z = p.z ? p.z : 0.0; // Z value not always supplied\n  var Longitude;\n  var Latitude;\n  var Height;\n\n  P = Math.sqrt(X * X + Y * Y);\n  RR = Math.sqrt(X * X + Y * Y + Z * Z);\n\n  /*      special cases for latitude and longitude */\n  if (P / a < genau) {\n    /*  special case, if P=0. (X=0., Y=0.) */\n    Longitude = 0.0;\n\n    /*  if (X,Y,Z)=(0.,0.,0.) then Height becomes semi-minor axis\n     *  of ellipsoid (=center of mass), Latitude becomes PI/2 */\n    if (RR / a < genau) {\n      Latitude = HALF_PI;\n      Height = -b;\n      return {\n        x: p.x,\n        y: p.y,\n        z: p.z\n      };\n    }\n  } else {\n    /*  ellipsoidal (geodetic) longitude\n     *  interval: -PI < Longitude <= +PI */\n    Longitude = Math.atan2(Y, X);\n  }\n\n  /* --------------------------------------------------------------\n   * Following iterative algorithm was developped by\n   * \"Institut for Erdmessung\", University of Hannover, July 1988.\n   * Internet: www.ife.uni-hannover.de\n   * Iterative computation of CPHI,SPHI and Height.\n   * Iteration of CPHI and SPHI to 10**-12 radian resp.\n   * 2*10**-7 arcsec.\n   * --------------------------------------------------------------\n   */\n  CT = Z / RR;\n  ST = P / RR;\n  RX = 1.0 / Math.sqrt(1.0 - es * (2.0 - es) * ST * ST);\n  CPHI0 = ST * (1.0 - es) * RX;\n  SPHI0 = CT * RX;\n  iter = 0;\n\n  /* loop to find sin(Latitude) resp. Latitude\n   * until |sin(Latitude(iter)-Latitude(iter-1))| < genau */\n  do {\n    iter++;\n    RN = a / Math.sqrt(1.0 - es * SPHI0 * SPHI0);\n\n    /*  ellipsoidal (geodetic) height */\n    Height = P * CPHI0 + Z * SPHI0 - RN * (1.0 - es * SPHI0 * SPHI0);\n\n    RK = es * RN / (RN + Height);\n    RX = 1.0 / Math.sqrt(1.0 - RK * (2.0 - RK) * ST * ST);\n    CPHI = ST * (1.0 - RK) * RX;\n    SPHI = CT * RX;\n    SDPHI = SPHI * CPHI0 - CPHI * SPHI0;\n    CPHI0 = CPHI;\n    SPHI0 = SPHI;\n  }\n  while (SDPHI * SDPHI > genau2 && iter < maxiter);\n\n  /*      ellipsoidal (geodetic) latitude */\n  Latitude = Math.atan(SPHI / Math.abs(CPHI));\n  return {\n    x: Longitude,\n    y: Latitude,\n    z: Height\n  };\n} // cs_geocentric_to_geodetic()\n\n/****************************************************************/\n// pj_geocentic_to_wgs84( p )\n//  p = point to transform in geocentric coordinates (x,y,z)\n\n/** point object, nothing fancy, just allows values to be\n    passed back and forth by reference rather than by value.\n    Other point classes may be used as long as they have\n    x and y properties, which will get modified in the transform method.\n*/\nexport function geocentricToWgs84(p, datum_type, datum_params) {\n  if (datum_type === PJD_3PARAM) {\n    // if( x[io] === HUGE_VAL )\n    //    continue;\n    return {\n      x: p.x + datum_params[0],\n      y: p.y + datum_params[1],\n      z: p.z + datum_params[2]\n    };\n  } else if (datum_type === PJD_7PARAM) {\n    var Dx_BF = datum_params[0];\n    var Dy_BF = datum_params[1];\n    var Dz_BF = datum_params[2];\n    var Rx_BF = datum_params[3];\n    var Ry_BF = datum_params[4];\n    var Rz_BF = datum_params[5];\n    var M_BF = datum_params[6];\n    // if( x[io] === HUGE_VAL )\n    //    continue;\n    return {\n      x: M_BF * (p.x - Rz_BF * p.y + Ry_BF * p.z) + Dx_BF,\n      y: M_BF * (Rz_BF * p.x + p.y - Rx_BF * p.z) + Dy_BF,\n      z: M_BF * (-Ry_BF * p.x + Rx_BF * p.y + p.z) + Dz_BF\n    };\n  }\n} // cs_geocentric_to_wgs84\n\n/****************************************************************/\n// pj_geocentic_from_wgs84()\n//  coordinate system definition,\n//  point to transform in geocentric coordinates (x,y,z)\nexport function geocentricFromWgs84(p, datum_type, datum_params) {\n  if (datum_type === PJD_3PARAM) {\n    // if( x[io] === HUGE_VAL )\n    //    continue;\n    return {\n      x: p.x - datum_params[0],\n      y: p.y - datum_params[1],\n      z: p.z - datum_params[2]\n    };\n  } else if (datum_type === PJD_7PARAM) {\n    var Dx_BF = datum_params[0];\n    var Dy_BF = datum_params[1];\n    var Dz_BF = datum_params[2];\n    var Rx_BF = datum_params[3];\n    var Ry_BF = datum_params[4];\n    var Rz_BF = datum_params[5];\n    var M_BF = datum_params[6];\n    var x_tmp = (p.x - Dx_BF) / M_BF;\n    var y_tmp = (p.y - Dy_BF) / M_BF;\n    var z_tmp = (p.z - Dz_BF) / M_BF;\n    // if( x[io] === HUGE_VAL )\n    //    continue;\n\n    return {\n      x: x_tmp + Rz_BF * y_tmp - Ry_BF * z_tmp,\n      y: -Rz_BF * x_tmp + y_tmp + Rx_BF * z_tmp,\n      z: Ry_BF * x_tmp - Rx_BF * y_tmp + z_tmp\n    };\n  } // cs_geocentric_from_wgs84()\n}\n"],"names":["compareDatums","source","dest","PJD_3PARAM","PJD_7PARAM","geodeticToGeocentric","p","es","a","Longitude","Latitude","Height","Rn","Sin_Lat","Sin2_Lat","Cos_Lat","HALF_PI","geocentricToGeodetic","b","genau","genau2","maxiter","P","RR","CT","ST","RX","RK","RN","CPHI0","SPHI0","CPHI","SPHI","SDPHI","iter","X","Y","Z","geocentricToWgs84","datum_type","datum_params","Dx_BF","Dy_BF","Dz_BF","Rx_BF","Ry_BF","Rz_BF","M_BF","geocentricFromWgs84","x_tmp","y_tmp","z_tmp"],"mappings":"qHAEO,SAASA,EAAcC,EAAQC,EAAM,CAC1C,OAAID,EAAO,aAAeC,EAAK,YAEpBD,EAAO,IAAMC,EAAK,GAAK,KAAK,IAAID,EAAO,GAAKC,EAAK,EAAE,EAAI,MADzD,GAKED,EAAO,aAAeE,aACvBF,EAAO,aAAa,CAAC,IAAMC,EAAK,aAAa,CAAC,GAAKD,EAAO,aAAa,CAAC,IAAMC,EAAK,aAAa,CAAC,GAAKD,EAAO,aAAa,CAAC,IAAMC,EAAK,aAAa,CAAC,EACnJD,EAAO,aAAeG,aACvBH,EAAO,aAAa,CAAC,IAAMC,EAAK,aAAa,CAAC,GAAKD,EAAO,aAAa,CAAC,IAAMC,EAAK,aAAa,CAAC,GAAKD,EAAO,aAAa,CAAC,IAAMC,EAAK,aAAa,CAAC,GAAKD,EAAO,aAAa,CAAC,IAAMC,EAAK,aAAa,CAAC,GAAKD,EAAO,aAAa,CAAC,IAAMC,EAAK,aAAa,CAAC,GAAKD,EAAO,aAAa,CAAC,IAAMC,EAAK,aAAa,CAAC,GAAKD,EAAO,aAAa,CAAC,IAAMC,EAAK,aAAa,CAAC,EAEjW,EAEX,CAeO,SAASG,EAAqBC,EAAGC,EAAIC,EAAG,CAC7C,IAAIC,EAAYH,EAAE,EACdI,EAAWJ,EAAE,EACbK,EAASL,EAAE,EAAIA,EAAE,EAAI,EAErBM,EACAC,EACAC,EACAC,EAOJ,GAAIL,EAAW,CAACM,EAAAA,SAAWN,EAAW,OAASM,EAAAA,QAC7CN,EAAW,CAACM,EAAAA,gBACHN,EAAWM,EAAAA,SAAWN,EAAW,MAAQM,EAAAA,QAClDN,EAAWM,EAAAA,YACb,IAAWN,EAAW,CAACM,UAGrB,MAAO,CAAE,EAAG,KAAW,EAAG,KAAW,EAAGV,EAAE,CAAA,EAC5C,GAAWI,EAAWM,UAEpB,MAAO,CAAE,EAAG,IAAU,EAAG,IAAU,EAAGV,EAAE,CAAA,EAG1C,OAAIG,EAAY,KAAK,KACnBA,GAAc,EAAI,KAAK,IAEzBI,EAAU,KAAK,IAAIH,CAAQ,EAC3BK,EAAU,KAAK,IAAIL,CAAQ,EAC3BI,EAAWD,EAAUA,EACrBD,EAAKJ,EAAK,KAAK,KAAK,EAAQD,EAAKO,CAAQ,EAClC,CACL,GAAIF,EAAKD,GAAUI,EAAU,KAAK,IAAIN,CAAS,EAC/C,GAAIG,EAAKD,GAAUI,EAAU,KAAK,IAAIN,CAAS,EAC/C,GAAKG,GAAM,EAAIL,GAAOI,GAAUE,CAAA,CAEpC,CAEO,SAASI,EAAqBX,EAAGC,EAAIC,EAAGU,EAAG,CAGhD,IAAIC,EAAQ,MACRC,EAAUD,EAAQA,EAClBE,EAAU,GAEVC,EACAC,EACAC,EACAC,EACAC,EACAC,EACAC,EACAC,EACAC,EACAC,EACAC,EACAC,EACAC,EAEAC,EAAI7B,EAAE,EACN8B,EAAI9B,EAAE,EACN+B,EAAI/B,EAAE,EAAIA,EAAE,EAAI,EAChBG,EACAC,EACAC,EAMJ,GAJAW,EAAI,KAAK,KAAKa,EAAIA,EAAIC,EAAIA,CAAC,EAC3Bb,EAAK,KAAK,KAAKY,EAAIA,EAAIC,EAAIA,EAAIC,EAAIA,CAAC,EAGhCf,EAAId,EAAIW,GAMV,GAJAV,EAAY,EAIRc,EAAKf,EAAIW,EACX,OAAAT,EAAWM,EAAAA,QACXL,EAAS,CAACO,EACH,CACL,EAAGZ,EAAE,EACL,EAAGA,EAAE,EACL,EAAGA,EAAE,CAAA,OAMTG,EAAY,KAAK,MAAM2B,EAAGD,CAAC,EAY7BX,EAAKa,EAAId,EACTE,EAAKH,EAAIC,EACTG,EAAK,EAAM,KAAK,KAAK,EAAMnB,GAAM,EAAMA,GAAMkB,EAAKA,CAAE,EACpDI,EAAQJ,GAAM,EAAMlB,GAAMmB,EAC1BI,EAAQN,EAAKE,EACbQ,EAAO,EAIP,GACEA,IACAN,EAAKpB,EAAI,KAAK,KAAK,EAAMD,EAAKuB,EAAQA,CAAK,EAG3CnB,EAASW,EAAIO,EAAQQ,EAAIP,EAAQF,GAAM,EAAMrB,EAAKuB,EAAQA,GAE1DH,EAAKpB,EAAKqB,GAAMA,EAAKjB,GACrBe,EAAK,EAAM,KAAK,KAAK,EAAMC,GAAM,EAAMA,GAAMF,EAAKA,CAAE,EACpDM,EAAON,GAAM,EAAME,GAAMD,EACzBM,EAAOR,EAAKE,EACZO,EAAQD,EAAOH,EAAQE,EAAOD,EAC9BD,EAAQE,EACRD,EAAQE,QAEHC,EAAQA,EAAQb,GAAUc,EAAOb,GAGxC,OAAAX,EAAW,KAAK,KAAKsB,EAAO,KAAK,IAAID,CAAI,CAAC,EACnC,CACL,EAAGtB,EACH,EAAGC,EACH,EAAGC,CAAA,CAEP,CAWO,SAAS2B,EAAkBhC,EAAGiC,EAAYC,EAAc,CAC7D,GAAID,IAAepC,EAAAA,WAGjB,MAAO,CACL,EAAGG,EAAE,EAAIkC,EAAa,CAAC,EACvB,EAAGlC,EAAE,EAAIkC,EAAa,CAAC,EACvB,EAAGlC,EAAE,EAAIkC,EAAa,CAAC,CAAA,EAE3B,GAAWD,IAAenC,aAAY,CACpC,IAAIqC,EAAQD,EAAa,CAAC,EACtBE,EAAQF,EAAa,CAAC,EACtBG,EAAQH,EAAa,CAAC,EACtBI,EAAQJ,EAAa,CAAC,EACtBK,EAAQL,EAAa,CAAC,EACtBM,EAAQN,EAAa,CAAC,EACtBO,EAAOP,EAAa,CAAC,EAGzB,MAAO,CACL,EAAGO,GAAQzC,EAAE,EAAIwC,EAAQxC,EAAE,EAAIuC,EAAQvC,EAAE,GAAKmC,EAC9C,EAAGM,GAAQD,EAAQxC,EAAE,EAAIA,EAAE,EAAIsC,EAAQtC,EAAE,GAAKoC,EAC9C,EAAGK,GAAQ,CAACF,EAAQvC,EAAE,EAAIsC,EAAQtC,EAAE,EAAIA,EAAE,GAAKqC,CAAA,CAEnD,CACF,CAMO,SAASK,EAAoB1C,EAAGiC,EAAYC,EAAc,CAC/D,GAAID,IAAepC,EAAAA,WAGjB,MAAO,CACL,EAAGG,EAAE,EAAIkC,EAAa,CAAC,EACvB,EAAGlC,EAAE,EAAIkC,EAAa,CAAC,EACvB,EAAGlC,EAAE,EAAIkC,EAAa,CAAC,CAAA,EAE3B,GAAWD,IAAenC,aAAY,CACpC,IAAIqC,EAAQD,EAAa,CAAC,EACtBE,EAAQF,EAAa,CAAC,EACtBG,EAAQH,EAAa,CAAC,EACtBI,EAAQJ,EAAa,CAAC,EACtBK,EAAQL,EAAa,CAAC,EACtBM,EAAQN,EAAa,CAAC,EACtBO,EAAOP,EAAa,CAAC,EACrBS,GAAS3C,EAAE,EAAImC,GAASM,EACxBG,GAAS5C,EAAE,EAAIoC,GAASK,EACxBI,GAAS7C,EAAE,EAAIqC,GAASI,EAI5B,MAAO,CACL,EAAGE,EAAQH,EAAQI,EAAQL,EAAQM,EACnC,EAAG,CAACL,EAAQG,EAAQC,EAAQN,EAAQO,EACpC,EAAGN,EAAQI,EAAQL,EAAQM,EAAQC,CAAA,CAEvC,CACF","x_google_ignoreList":[0]}

package/dist/index.cjs549.js

@@ -0,0 +1,2 @@
+"use strict";Object.defineProperties(exports,{__esModule:{value:!0},[Symbol.toStringTag]:{value:"Module"}});function t(){}function e(i){return i}var n=["longlat","identity"];const r={init:t,forward:e,inverse:e,names:n};exports.default=r;exports.forward=e;exports.init=t;exports.inverse=e;exports.names=n;
+//# sourceMappingURL=index.cjs549.js.map

package/dist/index.cjs549.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"index.cjs549.js","sources":["../../../node_modules/proj4/lib/projections/longlat.js"],"sourcesContent":["export function init() {\n  // no-op for longlat\n}\n\nfunction identity(pt) {\n  return pt;\n}\nexport { identity as forward };\nexport { identity as inverse };\nexport var names = ['longlat', 'identity'];\nexport default {\n  init: init,\n  forward: identity,\n  inverse: identity,\n  names: names\n};\n"],"names":["init","identity","pt","names","longlat"],"mappings":"4GAAO,SAASA,GAAO,CAEvB,CAEA,SAASC,EAASC,EAAI,CACpB,OAAOA,CACT,CAGO,IAAIC,EAAQ,CAAC,UAAW,UAAU,EACzC,MAAAC,EAAe,CACb,KAAAJ,EACA,QAASC,EACT,QAASA,EACT,MAAAE,CACF","x_google_ignoreList":[0]}

package/dist/index.cjs550.js

@@ -0,0 +1,2 @@
+"use strict";Object.defineProperties(exports,{__esModule:{value:!0},[Symbol.toStringTag]:{value:"Module"}});const M=require("./index.cjs531.js"),n=require("./index.cjs521.js"),v=require("./index.cjs532.js"),k=require("./index.cjs533.js"),e=require("./index.cjs386.js");function o(){var t=this.b/this.a;this.es=1-t*t,"x0"in this||(this.x0=0),"y0"in this||(this.y0=0),this.e=Math.sqrt(this.es),this.lat_ts?this.sphere?this.k0=Math.cos(this.lat_ts):this.k0=M(this.e,Math.sin(this.lat_ts),Math.cos(this.lat_ts)):this.k0||(this.k?this.k0=this.k:this.k0=1)}function l(t){var h=t.x,s=t.y;if(s*e.R2D>90&&s*e.R2D<-90&&h*e.R2D>180&&h*e.R2D<-180)return null;var r,i;if(Math.abs(Math.abs(s)-e.HALF_PI)<=e.EPSLN)return null;if(this.sphere)r=this.x0+this.a*this.k0*n(h-this.long0,this.over),i=this.y0+this.a*this.k0*Math.log(Math.tan(e.FORTPI+.5*s));else{var a=Math.sin(s),f=v(this.e,s,a);r=this.x0+this.a*this.k0*n(h-this.long0,this.over),i=this.y0-this.a*this.k0*Math.log(f)}return t.x=r,t.y=i,t}function u(t){var h=t.x-this.x0,s=t.y-this.y0,r,i;if(this.sphere)i=e.HALF_PI-2*Math.atan(Math.exp(-s/(this.a*this.k0)));else{var a=Math.exp(-s/(this.a*this.k0));if(i=k(this.e,a),i===-9999)return null}return r=n(this.long0+h/(this.a*this.k0),this.over),t.x=r,t.y=i,t}var c=["Mercator","Popular Visualisation Pseudo Mercator","Mercator_1SP","Mercator_Auxiliary_Sphere","Mercator_Variant_A","merc"];const x={init:o,forward:l,inverse:u,names:c};exports.default=x;exports.forward=l;exports.init=o;exports.inverse=u;exports.names=c;
+//# sourceMappingURL=index.cjs550.js.map

package/dist/index.cjs550.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"index.cjs550.js","sources":["../../../node_modules/proj4/lib/projections/merc.js"],"sourcesContent":["import msfnz from '../common/msfnz';\n\nimport adjust_lon from '../common/adjust_lon';\nimport tsfnz from '../common/tsfnz';\nimport phi2z from '../common/phi2z';\nimport { FORTPI, R2D, EPSLN, HALF_PI } from '../constants/values';\n\n/**\n * @typedef {Object} LocalThis\n * @property {number} es\n * @property {number} e\n * @property {number} k\n */\n\n/** @this {import('../defs.js').ProjectionDefinition & LocalThis} */\nexport function init() {\n  var con = this.b / this.a;\n  this.es = 1 - con * con;\n  if (!('x0' in this)) {\n    this.x0 = 0;\n  }\n  if (!('y0' in this)) {\n    this.y0 = 0;\n  }\n  this.e = Math.sqrt(this.es);\n  if (this.lat_ts) {\n    if (this.sphere) {\n      this.k0 = Math.cos(this.lat_ts);\n    } else {\n      this.k0 = msfnz(this.e, Math.sin(this.lat_ts), Math.cos(this.lat_ts));\n    }\n  } else {\n    if (!this.k0) {\n      if (this.k) {\n        this.k0 = this.k;\n      } else {\n        this.k0 = 1;\n      }\n    }\n  }\n}\n\n/* Mercator forward equations--mapping lat,long to x,y\n  -------------------------------------------------- */\n\nexport function forward(p) {\n  var lon = p.x;\n  var lat = p.y;\n  // convert to radians\n  if (lat * R2D > 90 && lat * R2D < -90 && lon * R2D > 180 && lon * R2D < -180) {\n    return null;\n  }\n\n  var x, y;\n  if (Math.abs(Math.abs(lat) - HALF_PI) <= EPSLN) {\n    return null;\n  } else {\n    if (this.sphere) {\n      x = this.x0 + this.a * this.k0 * adjust_lon(lon - this.long0, this.over);\n      y = this.y0 + this.a * this.k0 * Math.log(Math.tan(FORTPI + 0.5 * lat));\n    } else {\n      var sinphi = Math.sin(lat);\n      var ts = tsfnz(this.e, lat, sinphi);\n      x = this.x0 + this.a * this.k0 * adjust_lon(lon - this.long0, this.over);\n      y = this.y0 - this.a * this.k0 * Math.log(ts);\n    }\n    p.x = x;\n    p.y = y;\n    return p;\n  }\n}\n\n/* Mercator inverse equations--mapping x,y to lat/long\n  -------------------------------------------------- */\nexport function inverse(p) {\n  var x = p.x - this.x0;\n  var y = p.y - this.y0;\n  var lon, lat;\n\n  if (this.sphere) {\n    lat = HALF_PI - 2 * Math.atan(Math.exp(-y / (this.a * this.k0)));\n  } else {\n    var ts = Math.exp(-y / (this.a * this.k0));\n    lat = phi2z(this.e, ts);\n    if (lat === -9999) {\n      return null;\n    }\n  }\n  lon = adjust_lon(this.long0 + x / (this.a * this.k0), this.over);\n\n  p.x = lon;\n  p.y = lat;\n  return p;\n}\n\nexport var names = ['Mercator', 'Popular Visualisation Pseudo Mercator', 'Mercator_1SP', 'Mercator_Auxiliary_Sphere', 'Mercator_Variant_A', 'merc'];\nexport default {\n  init: init,\n  forward: forward,\n  inverse: inverse,\n  names: names\n};\n"],"names":["init","con","msfnz","forward","p","lon","lat","R2D","x","y","HALF_PI","EPSLN","adjust_lon","FORTPI","sinphi","ts","tsfnz","inverse","phi2z","names","merc"],"mappings":"6QAeO,SAASA,GAAO,CACrB,IAAIC,EAAM,KAAK,EAAI,KAAK,EACxB,KAAK,GAAK,EAAIA,EAAMA,EACd,OAAQ,OACZ,KAAK,GAAK,GAEN,OAAQ,OACZ,KAAK,GAAK,GAEZ,KAAK,EAAI,KAAK,KAAK,KAAK,EAAE,EACtB,KAAK,OACH,KAAK,OACP,KAAK,GAAK,KAAK,IAAI,KAAK,MAAM,EAE9B,KAAK,GAAKC,EAAM,KAAK,EAAG,KAAK,IAAI,KAAK,MAAM,EAAG,KAAK,IAAI,KAAK,MAAM,CAAC,EAGjE,KAAK,KACJ,KAAK,EACP,KAAK,GAAK,KAAK,EAEf,KAAK,GAAK,EAIlB,CAKO,SAASC,EAAQC,EAAG,CACzB,IAAIC,EAAMD,EAAE,EACRE,EAAMF,EAAE,EAEZ,GAAIE,EAAMC,EAAAA,IAAM,IAAMD,EAAMC,EAAAA,IAAM,KAAOF,EAAME,EAAAA,IAAM,KAAOF,EAAME,EAAAA,IAAM,KACtE,OAAO,KAGT,IAAIC,EAAGC,EACP,GAAI,KAAK,IAAI,KAAK,IAAIH,CAAG,EAAII,EAAAA,OAAO,GAAKC,QACvC,OAAO,KAEP,GAAI,KAAK,OACPH,EAAI,KAAK,GAAK,KAAK,EAAI,KAAK,GAAKI,EAAWP,EAAM,KAAK,MAAO,KAAK,IAAI,EACvEI,EAAI,KAAK,GAAK,KAAK,EAAI,KAAK,GAAK,KAAK,IAAI,KAAK,IAAII,EAAAA,OAAS,GAAMP,CAAG,CAAC,MACjE,CACL,IAAIQ,EAAS,KAAK,IAAIR,CAAG,EACrBS,EAAKC,EAAM,KAAK,EAAGV,EAAKQ,CAAM,EAClCN,EAAI,KAAK,GAAK,KAAK,EAAI,KAAK,GAAKI,EAAWP,EAAM,KAAK,MAAO,KAAK,IAAI,EACvEI,EAAI,KAAK,GAAK,KAAK,EAAI,KAAK,GAAK,KAAK,IAAIM,CAAE,CAC9C,CACA,OAAAX,EAAE,EAAII,EACNJ,EAAE,EAAIK,EACCL,CAEX,CAIO,SAASa,EAAQb,EAAG,CACzB,IAAII,EAAIJ,EAAE,EAAI,KAAK,GACfK,EAAIL,EAAE,EAAI,KAAK,GACfC,EAAKC,EAET,GAAI,KAAK,OACPA,EAAMI,EAAAA,QAAU,EAAI,KAAK,KAAK,KAAK,IAAI,CAACD,GAAK,KAAK,EAAI,KAAK,GAAG,CAAC,MAC1D,CACL,IAAIM,EAAK,KAAK,IAAI,CAACN,GAAK,KAAK,EAAI,KAAK,GAAG,EAEzC,GADAH,EAAMY,EAAM,KAAK,EAAGH,CAAE,EAClBT,IAAQ,MACV,OAAO,IAEX,CACA,OAAAD,EAAMO,EAAW,KAAK,MAAQJ,GAAK,KAAK,EAAI,KAAK,IAAK,KAAK,IAAI,EAE/DJ,EAAE,EAAIC,EACND,EAAE,EAAIE,EACCF,CACT,CAEO,IAAIe,EAAQ,CAAC,WAAY,wCAAyC,eAAgB,4BAA6B,qBAAsB,MAAM,EAClJ,MAAAC,EAAe,CACb,KAAApB,EACA,QAAAG,EACA,QAAAc,EACA,MAAAE,CACF","x_google_ignoreList":[0]}

package/dist/index.cjs551.js

@@ -0,0 +1,2 @@
+"use strict";var e={MERIT:{a:6378137,rf:298.257,ellipseName:"MERIT 1983"},SGS85:{a:6378136,rf:298.257,ellipseName:"Soviet Geodetic System 85"},GRS80:{a:6378137,rf:298.257222101,ellipseName:"GRS 1980(IUGG, 1980)"},IAU76:{a:6378140,rf:298.257,ellipseName:"IAU 1976"},airy:{a:6377563396e-3,b:635625691e-2,ellipseName:"Airy 1830"},APL4:{a:6378137,rf:298.25,ellipseName:"Appl. Physics. 1965"},NWL9D:{a:6378145,rf:298.25,ellipseName:"Naval Weapons Lab., 1965"},mod_airy:{a:6377340189e-3,b:6356034446e-3,ellipseName:"Modified Airy"},andrae:{a:637710443e-2,rf:300,ellipseName:"Andrae 1876 (Den., Iclnd.)"},aust_SA:{a:6378160,rf:298.25,ellipseName:"Australian Natl & S. Amer. 1969"},GRS67:{a:6378160,rf:298.247167427,ellipseName:"GRS 67(IUGG 1967)"},bessel:{a:6377397155e-3,rf:299.1528128,ellipseName:"Bessel 1841"},bess_nam:{a:6377483865e-3,rf:299.1528128,ellipseName:"Bessel 1841 (Namibia)"},clrk66:{a:63782064e-1,b:63565838e-1,ellipseName:"Clarke 1866"},clrk80:{a:6378249145e-3,rf:293.4663,ellipseName:"Clarke 1880 mod."},clrk80ign:{a:63782492e-1,b:6356515,rf:293.4660213,ellipseName:"Clarke 1880 (IGN)"},clrk58:{a:6378293645208759e-9,rf:294.2606763692654,ellipseName:"Clarke 1858"},CPM:{a:63757387e-1,rf:334.29,ellipseName:"Comm. des Poids et Mesures 1799"},delmbr:{a:6376428,rf:311.5,ellipseName:"Delambre 1810 (Belgium)"},engelis:{a:637813605e-2,rf:298.2566,ellipseName:"Engelis 1985"},evrst30:{a:6377276345e-3,rf:300.8017,ellipseName:"Everest 1830"},evrst48:{a:6377304063e-3,rf:300.8017,ellipseName:"Everest 1948"},evrst56:{a:6377301243e-3,rf:300.8017,ellipseName:"Everest 1956"},evrst69:{a:6377295664e-3,rf:300.8017,ellipseName:"Everest 1969"},evrstSS:{a:6377298556e-3,rf:300.8017,ellipseName:"Everest (Sabah & Sarawak)"},fschr60:{a:6378166,rf:298.3,ellipseName:"Fischer (Mercury Datum) 1960"},fschr60m:{a:6378155,rf:298.3,ellipseName:"Fischer 1960"},fschr68:{a:6378150,rf:298.3,ellipseName:"Fischer 1968"},helmert:{a:6378200,rf:298.3,ellipseName:"Helmert 1906"},hough:{a:6378270,rf:297,ellipseName:"Hough"},intl:{a:6378388,rf:297,ellipseName:"International 1909 (Hayford)"},kaula:{a:6378163,rf:298.24,ellipseName:"Kaula 1961"},lerch:{a:6378139,rf:298.257,ellipseName:"Lerch 1979"},mprts:{a:6397300,rf:191,ellipseName:"Maupertius 1738"},new_intl:{a:63781575e-1,b:63567722e-1,ellipseName:"New International 1967"},plessis:{a:6376523,rf:6355863,ellipseName:"Plessis 1817 (France)"},krass:{a:6378245,rf:298.3,ellipseName:"Krassovsky, 1942"},SEasia:{a:6378155,b:63567733205e-4,ellipseName:"Southeast Asia"},walbeck:{a:6376896,b:63558348467e-4,ellipseName:"Walbeck"},WGS60:{a:6378165,rf:298.3,ellipseName:"WGS 60"},WGS66:{a:6378145,rf:298.25,ellipseName:"WGS 66"},WGS7:{a:6378135,rf:298.26,ellipseName:"WGS 72"},WGS84:{a:6378137,rf:298.257223563,ellipseName:"WGS 84"},sphere:{a:6370997,b:6370997,ellipseName:"Normal Sphere (r=6370997)"}};module.exports=e;
+//# sourceMappingURL=index.cjs551.js.map

package/dist/index.cjs551.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"index.cjs551.js","sources":["../../../node_modules/proj4/lib/constants/Ellipsoid.js"],"sourcesContent":["var ellipsoids = {\n  MERIT: {\n    a: 6378137,\n    rf: 298.257,\n    ellipseName: 'MERIT 1983'\n  },\n  SGS85: {\n    a: 6378136,\n    rf: 298.257,\n    ellipseName: 'Soviet Geodetic System 85'\n  },\n  GRS80: {\n    a: 6378137,\n    rf: 298.257222101,\n    ellipseName: 'GRS 1980(IUGG, 1980)'\n  },\n  IAU76: {\n    a: 6378140,\n    rf: 298.257,\n    ellipseName: 'IAU 1976'\n  },\n  airy: {\n    a: 6377563.396,\n    b: 6356256.91,\n    ellipseName: 'Airy 1830'\n  },\n  APL4: {\n    a: 6378137,\n    rf: 298.25,\n    ellipseName: 'Appl. Physics. 1965'\n  },\n  NWL9D: {\n    a: 6378145,\n    rf: 298.25,\n    ellipseName: 'Naval Weapons Lab., 1965'\n  },\n  mod_airy: {\n    a: 6377340.189,\n    b: 6356034.446,\n    ellipseName: 'Modified Airy'\n  },\n  andrae: {\n    a: 6377104.43,\n    rf: 300,\n    ellipseName: 'Andrae 1876 (Den., Iclnd.)'\n  },\n  aust_SA: {\n    a: 6378160,\n    rf: 298.25,\n    ellipseName: 'Australian Natl & S. Amer. 1969'\n  },\n  GRS67: {\n    a: 6378160,\n    rf: 298.247167427,\n    ellipseName: 'GRS 67(IUGG 1967)'\n  },\n  bessel: {\n    a: 6377397.155,\n    rf: 299.1528128,\n    ellipseName: 'Bessel 1841'\n  },\n  bess_nam: {\n    a: 6377483.865,\n    rf: 299.1528128,\n    ellipseName: 'Bessel 1841 (Namibia)'\n  },\n  clrk66: {\n    a: 6378206.4,\n    b: 6356583.8,\n    ellipseName: 'Clarke 1866'\n  },\n  clrk80: {\n    a: 6378249.145,\n    rf: 293.4663,\n    ellipseName: 'Clarke 1880 mod.'\n  },\n  clrk80ign: {\n    a: 6378249.2,\n    b: 6356515,\n    rf: 293.4660213,\n    ellipseName: 'Clarke 1880 (IGN)'\n  },\n  clrk58: {\n    a: 6378293.645208759,\n    rf: 294.2606763692654,\n    ellipseName: 'Clarke 1858'\n  },\n  CPM: {\n    a: 6375738.7,\n    rf: 334.29,\n    ellipseName: 'Comm. des Poids et Mesures 1799'\n  },\n  delmbr: {\n    a: 6376428,\n    rf: 311.5,\n    ellipseName: 'Delambre 1810 (Belgium)'\n  },\n  engelis: {\n    a: 6378136.05,\n    rf: 298.2566,\n    ellipseName: 'Engelis 1985'\n  },\n  evrst30: {\n    a: 6377276.345,\n    rf: 300.8017,\n    ellipseName: 'Everest 1830'\n  },\n  evrst48: {\n    a: 6377304.063,\n    rf: 300.8017,\n    ellipseName: 'Everest 1948'\n  },\n  evrst56: {\n    a: 6377301.243,\n    rf: 300.8017,\n    ellipseName: 'Everest 1956'\n  },\n  evrst69: {\n    a: 6377295.664,\n    rf: 300.8017,\n    ellipseName: 'Everest 1969'\n  },\n  evrstSS: {\n    a: 6377298.556,\n    rf: 300.8017,\n    ellipseName: 'Everest (Sabah & Sarawak)'\n  },\n  fschr60: {\n    a: 6378166,\n    rf: 298.3,\n    ellipseName: 'Fischer (Mercury Datum) 1960'\n  },\n  fschr60m: {\n    a: 6378155,\n    rf: 298.3,\n    ellipseName: 'Fischer 1960'\n  },\n  fschr68: {\n    a: 6378150,\n    rf: 298.3,\n    ellipseName: 'Fischer 1968'\n  },\n  helmert: {\n    a: 6378200,\n    rf: 298.3,\n    ellipseName: 'Helmert 1906'\n  },\n  hough: {\n    a: 6378270,\n    rf: 297,\n    ellipseName: 'Hough'\n  },\n  intl: {\n    a: 6378388,\n    rf: 297,\n    ellipseName: 'International 1909 (Hayford)'\n  },\n  kaula: {\n    a: 6378163,\n    rf: 298.24,\n    ellipseName: 'Kaula 1961'\n  },\n  lerch: {\n    a: 6378139,\n    rf: 298.257,\n    ellipseName: 'Lerch 1979'\n  },\n  mprts: {\n    a: 6397300,\n    rf: 191,\n    ellipseName: 'Maupertius 1738'\n  },\n  new_intl: {\n    a: 6378157.5,\n    b: 6356772.2,\n    ellipseName: 'New International 1967'\n  },\n  plessis: {\n    a: 6376523,\n    rf: 6355863,\n    ellipseName: 'Plessis 1817 (France)'\n  },\n  krass: {\n    a: 6378245,\n    rf: 298.3,\n    ellipseName: 'Krassovsky, 1942'\n  },\n  SEasia: {\n    a: 6378155,\n    b: 6356773.3205,\n    ellipseName: 'Southeast Asia'\n  },\n  walbeck: {\n    a: 6376896,\n    b: 6355834.8467,\n    ellipseName: 'Walbeck'\n  },\n  WGS60: {\n    a: 6378165,\n    rf: 298.3,\n    ellipseName: 'WGS 60'\n  },\n  WGS66: {\n    a: 6378145,\n    rf: 298.25,\n    ellipseName: 'WGS 66'\n  },\n  WGS7: {\n    a: 6378135,\n    rf: 298.26,\n    ellipseName: 'WGS 72'\n  },\n  WGS84: {\n    a: 6378137,\n    rf: 298.257223563,\n    ellipseName: 'WGS 84'\n  },\n  sphere: {\n    a: 6370997,\n    b: 6370997,\n    ellipseName: 'Normal Sphere (r=6370997)'\n  }\n};\n\nexport default ellipsoids;\n"],"names":["ellipsoids"],"mappings":"aAAA,IAAIA,EAAa,CACf,MAAO,CACL,EAAG,QACH,GAAI,QACJ,YAAa,YAAA,EAEf,MAAO,CACL,EAAG,QACH,GAAI,QACJ,YAAa,2BAAA,EAEf,MAAO,CACL,EAAG,QACH,GAAI,cACJ,YAAa,sBAAA,EAEf,MAAO,CACL,EAAG,QACH,GAAI,QACJ,YAAa,UAAA,EAEf,KAAM,CACJ,EAAG,cACH,EAAG,aACH,YAAa,WAAA,EAEf,KAAM,CACJ,EAAG,QACH,GAAI,OACJ,YAAa,qBAAA,EAEf,MAAO,CACL,EAAG,QACH,GAAI,OACJ,YAAa,0BAAA,EAEf,SAAU,CACR,EAAG,cACH,EAAG,cACH,YAAa,eAAA,EAEf,OAAQ,CACN,EAAG,aACH,GAAI,IACJ,YAAa,4BAAA,EAEf,QAAS,CACP,EAAG,QACH,GAAI,OACJ,YAAa,iCAAA,EAEf,MAAO,CACL,EAAG,QACH,GAAI,cACJ,YAAa,mBAAA,EAEf,OAAQ,CACN,EAAG,cACH,GAAI,YACJ,YAAa,aAAA,EAEf,SAAU,CACR,EAAG,cACH,GAAI,YACJ,YAAa,uBAAA,EAEf,OAAQ,CACN,EAAG,YACH,EAAG,YACH,YAAa,aAAA,EAEf,OAAQ,CACN,EAAG,cACH,GAAI,SACJ,YAAa,kBAAA,EAEf,UAAW,CACT,EAAG,YACH,EAAG,QACH,GAAI,YACJ,YAAa,mBAAA,EAEf,OAAQ,CACN,EAAG,oBACH,GAAI,kBACJ,YAAa,aAAA,EAEf,IAAK,CACH,EAAG,YACH,GAAI,OACJ,YAAa,iCAAA,EAEf,OAAQ,CACN,EAAG,QACH,GAAI,MACJ,YAAa,yBAAA,EAEf,QAAS,CACP,EAAG,aACH,GAAI,SACJ,YAAa,cAAA,EAEf,QAAS,CACP,EAAG,cACH,GAAI,SACJ,YAAa,cAAA,EAEf,QAAS,CACP,EAAG,cACH,GAAI,SACJ,YAAa,cAAA,EAEf,QAAS,CACP,EAAG,cACH,GAAI,SACJ,YAAa,cAAA,EAEf,QAAS,CACP,EAAG,cACH,GAAI,SACJ,YAAa,cAAA,EAEf,QAAS,CACP,EAAG,cACH,GAAI,SACJ,YAAa,2BAAA,EAEf,QAAS,CACP,EAAG,QACH,GAAI,MACJ,YAAa,8BAAA,EAEf,SAAU,CACR,EAAG,QACH,GAAI,MACJ,YAAa,cAAA,EAEf,QAAS,CACP,EAAG,QACH,GAAI,MACJ,YAAa,cAAA,EAEf,QAAS,CACP,EAAG,QACH,GAAI,MACJ,YAAa,cAAA,EAEf,MAAO,CACL,EAAG,QACH,GAAI,IACJ,YAAa,OAAA,EAEf,KAAM,CACJ,EAAG,QACH,GAAI,IACJ,YAAa,8BAAA,EAEf,MAAO,CACL,EAAG,QACH,GAAI,OACJ,YAAa,YAAA,EAEf,MAAO,CACL,EAAG,QACH,GAAI,QACJ,YAAa,YAAA,EAEf,MAAO,CACL,EAAG,QACH,GAAI,IACJ,YAAa,iBAAA,EAEf,SAAU,CACR,EAAG,YACH,EAAG,YACH,YAAa,wBAAA,EAEf,QAAS,CACP,EAAG,QACH,GAAI,QACJ,YAAa,uBAAA,EAEf,MAAO,CACL,EAAG,QACH,GAAI,MACJ,YAAa,kBAAA,EAEf,OAAQ,CACN,EAAG,QACH,EAAG,eACH,YAAa,gBAAA,EAEf,QAAS,CACP,EAAG,QACH,EAAG,eACH,YAAa,SAAA,EAEf,MAAO,CACL,EAAG,QACH,GAAI,MACJ,YAAa,QAAA,EAEf,MAAO,CACL,EAAG,QACH,GAAI,OACJ,YAAa,QAAA,EAEf,KAAM,CACJ,EAAG,QACH,GAAI,OACJ,YAAa,QAAA,EAEf,MAAO,CACL,EAAG,QACH,GAAI,cACJ,YAAa,QAAA,EAEf,OAAQ,CACN,EAAG,QACH,EAAG,QACH,YAAa,2BAAA,CAEjB","x_google_ignoreList":[0]}

package/dist/index.cjs552.js

@@ -0,0 +1,2 @@
+"use strict";const p=require("./index.cjs553.js"),c=require("./index.cjs554.js"),h=require("./index.cjs555.js"),s=require("./index.cjs556.js"),l=require("./index.cjs557.js");function e(t){var a=-1,o=t==null?0:t.length;for(this.clear();++a<o;){var r=t[a];this.set(r[0],r[1])}}e.prototype.clear=p;e.prototype.delete=c;e.prototype.get=h;e.prototype.has=s;e.prototype.set=l;module.exports=e;
+//# sourceMappingURL=index.cjs552.js.map

package/dist/index.cjs552.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"index.cjs552.js","sources":["../../../node_modules/lodash-es/_MapCache.js"],"sourcesContent":["import mapCacheClear from './_mapCacheClear.js';\nimport mapCacheDelete from './_mapCacheDelete.js';\nimport mapCacheGet from './_mapCacheGet.js';\nimport mapCacheHas from './_mapCacheHas.js';\nimport mapCacheSet from './_mapCacheSet.js';\n\n/**\n * Creates a map cache object to store key-value pairs.\n *\n * @private\n * @constructor\n * @param {Array} [entries] The key-value pairs to cache.\n */\nfunction MapCache(entries) {\n  var index = -1,\n      length = entries == null ? 0 : entries.length;\n\n  this.clear();\n  while (++index < length) {\n    var entry = entries[index];\n    this.set(entry[0], entry[1]);\n  }\n}\n\n// Add methods to `MapCache`.\nMapCache.prototype.clear = mapCacheClear;\nMapCache.prototype['delete'] = mapCacheDelete;\nMapCache.prototype.get = mapCacheGet;\nMapCache.prototype.has = mapCacheHas;\nMapCache.prototype.set = mapCacheSet;\n\nexport default MapCache;\n"],"names":["MapCache","entries","index","length","entry","mapCacheClear","mapCacheDelete","mapCacheGet","mapCacheHas","mapCacheSet"],"mappings":"8KAaA,SAASA,EAASC,EAAS,CACzB,IAAIC,EAAQ,GACRC,EAASF,GAAW,KAAO,EAAIA,EAAQ,OAG3C,IADA,KAAK,MAAA,EACE,EAAEC,EAAQC,GAAQ,CACvB,IAAIC,EAAQH,EAAQC,CAAK,EACzB,KAAK,IAAIE,EAAM,CAAC,EAAGA,EAAM,CAAC,CAAC,CAC7B,CACF,CAGAJ,EAAS,UAAU,MAAQK,EAC3BL,EAAS,UAAU,OAAYM,EAC/BN,EAAS,UAAU,IAAMO,EACzBP,EAAS,UAAU,IAAMQ,EACzBR,EAAS,UAAU,IAAMS","x_google_ignoreList":[0]}

package/dist/index.cjs553.js

@@ -0,0 +1,2 @@
+"use strict";const e=require("./index.cjs594.js"),s=require("./index.cjs573.js"),t=require("./index.cjs513.js");function a(){this.size=0,this.__data__={hash:new e,map:new(t||s),string:new e}}module.exports=a;
+//# sourceMappingURL=index.cjs553.js.map

package/dist/index.cjs553.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"index.cjs553.js","sources":["../../../node_modules/lodash-es/_mapCacheClear.js"],"sourcesContent":["import Hash from './_Hash.js';\nimport ListCache from './_ListCache.js';\nimport Map from './_Map.js';\n\n/**\n * Removes all key-value entries from the map.\n *\n * @private\n * @name clear\n * @memberOf MapCache\n */\nfunction mapCacheClear() {\n  this.size = 0;\n  this.__data__ = {\n    'hash': new Hash,\n    'map': new (Map || ListCache),\n    'string': new Hash\n  };\n}\n\nexport default mapCacheClear;\n"],"names":["mapCacheClear","Hash","Map","ListCache"],"mappings":"gHAWA,SAASA,GAAgB,CACvB,KAAK,KAAO,EACZ,KAAK,SAAW,CACd,KAAQ,IAAIC,EACZ,IAAO,IAAKC,GAAOC,GACnB,OAAU,IAAIF,CAAA,CAElB","x_google_ignoreList":[0]}

package/dist/index.cjs554.js

@@ -0,0 +1,2 @@
+"use strict";const r=require("./index.cjs595.js");function s(e){var t=r(this,e).delete(e);return this.size-=t?1:0,t}module.exports=s;
+//# sourceMappingURL=index.cjs554.js.map

package/dist/index.cjs554.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"index.cjs554.js","sources":["../../../node_modules/lodash-es/_mapCacheDelete.js"],"sourcesContent":["import getMapData from './_getMapData.js';\n\n/**\n * Removes `key` and its value from the map.\n *\n * @private\n * @name delete\n * @memberOf MapCache\n * @param {string} key The key of the value to remove.\n * @returns {boolean} Returns `true` if the entry was removed, else `false`.\n */\nfunction mapCacheDelete(key) {\n  var result = getMapData(this, key)['delete'](key);\n  this.size -= result ? 1 : 0;\n  return result;\n}\n\nexport default mapCacheDelete;\n"],"names":["mapCacheDelete","key","result","getMapData"],"mappings":"kDAWA,SAASA,EAAeC,EAAK,CAC3B,IAAIC,EAASC,EAAW,KAAMF,CAAG,EAAE,OAAUA,CAAG,EAChD,YAAK,MAAQC,EAAS,EAAI,EACnBA,CACT","x_google_ignoreList":[0]}

package/dist/index.cjs555.js

@@ -0,0 +1,2 @@
+"use strict";const e=require("./index.cjs595.js");function r(t){return e(this,t).get(t)}module.exports=r;
+//# sourceMappingURL=index.cjs555.js.map

package/dist/index.cjs555.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"index.cjs555.js","sources":["../../../node_modules/lodash-es/_mapCacheGet.js"],"sourcesContent":["import getMapData from './_getMapData.js';\n\n/**\n * Gets the map value for `key`.\n *\n * @private\n * @name get\n * @memberOf MapCache\n * @param {string} key The key of the value to get.\n * @returns {*} Returns the entry value.\n */\nfunction mapCacheGet(key) {\n  return getMapData(this, key).get(key);\n}\n\nexport default mapCacheGet;\n"],"names":["mapCacheGet","key","getMapData"],"mappings":"kDAWA,SAASA,EAAYC,EAAK,CACxB,OAAOC,EAAW,KAAMD,CAAG,EAAE,IAAIA,CAAG,CACtC","x_google_ignoreList":[0]}

package/dist/index.cjs556.js

@@ -0,0 +1,2 @@
+"use strict";const a=require("./index.cjs595.js");function s(t){return a(this,t).has(t)}module.exports=s;
+//# sourceMappingURL=index.cjs556.js.map

package/dist/index.cjs556.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"index.cjs556.js","sources":["../../../node_modules/lodash-es/_mapCacheHas.js"],"sourcesContent":["import getMapData from './_getMapData.js';\n\n/**\n * Checks if a map value for `key` exists.\n *\n * @private\n * @name has\n * @memberOf MapCache\n * @param {string} key The key of the entry to check.\n * @returns {boolean} Returns `true` if an entry for `key` exists, else `false`.\n */\nfunction mapCacheHas(key) {\n  return getMapData(this, key).has(key);\n}\n\nexport default mapCacheHas;\n"],"names":["mapCacheHas","key","getMapData"],"mappings":"kDAWA,SAASA,EAAYC,EAAK,CACxB,OAAOC,EAAW,KAAMD,CAAG,EAAE,IAAIA,CAAG,CACtC","x_google_ignoreList":[0]}

package/dist/index.cjs557.js

@@ -0,0 +1,2 @@
+"use strict";const a=require("./index.cjs595.js");function r(e,s){var t=a(this,e),i=t.size;return t.set(e,s),this.size+=t.size==i?0:1,this}module.exports=r;
+//# sourceMappingURL=index.cjs557.js.map

package/dist/index.cjs557.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"index.cjs557.js","sources":["../../../node_modules/lodash-es/_mapCacheSet.js"],"sourcesContent":["import getMapData from './_getMapData.js';\n\n/**\n * Sets the map `key` to `value`.\n *\n * @private\n * @name set\n * @memberOf MapCache\n * @param {string} key The key of the value to set.\n * @param {*} value The value to set.\n * @returns {Object} Returns the map cache instance.\n */\nfunction mapCacheSet(key, value) {\n  var data = getMapData(this, key),\n      size = data.size;\n\n  data.set(key, value);\n  this.size += data.size == size ? 0 : 1;\n  return this;\n}\n\nexport default mapCacheSet;\n"],"names":["mapCacheSet","key","value","data","getMapData","size"],"mappings":"kDAYA,SAASA,EAAYC,EAAKC,EAAO,CAC/B,IAAIC,EAAOC,EAAW,KAAMH,CAAG,EAC3BI,EAAOF,EAAK,KAEhB,OAAAA,EAAK,IAAIF,EAAKC,CAAK,EACnB,KAAK,MAAQC,EAAK,MAAQE,EAAO,EAAI,EAC9B,IACT","x_google_ignoreList":[0]}

package/dist/index.cjs558.js

@@ -0,0 +1,2 @@
+"use strict";function o(r){var e=-1,t=Array(r.size);return r.forEach(function(n,i){t[++e]=[i,n]}),t}module.exports=o;
+//# sourceMappingURL=index.cjs558.js.map

package/dist/index.cjs558.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"index.cjs558.js","sources":["../../../node_modules/lodash-es/_mapToArray.js"],"sourcesContent":["/**\n * Converts `map` to its key-value pairs.\n *\n * @private\n * @param {Object} map The map to convert.\n * @returns {Array} Returns the key-value pairs.\n */\nfunction mapToArray(map) {\n  var index = -1,\n      result = Array(map.size);\n\n  map.forEach(function(value, key) {\n    result[++index] = [key, value];\n  });\n  return result;\n}\n\nexport default mapToArray;\n"],"names":["mapToArray","map","index","result","value","key"],"mappings":"aAOA,SAASA,EAAWC,EAAK,CACvB,IAAIC,EAAQ,GACRC,EAAS,MAAMF,EAAI,IAAI,EAE3B,OAAAA,EAAI,QAAQ,SAASG,EAAOC,EAAK,CAC/BF,EAAO,EAAED,CAAK,EAAI,CAACG,EAAKD,CAAK,CAC/B,CAAC,EACMD,CACT","x_google_ignoreList":[0]}

package/dist/index.cjs559.js

@@ -0,0 +1,2 @@
+"use strict";function i(r){var e=-1,t=Array(r.size);return r.forEach(function(n){t[++e]=n}),t}module.exports=i;
+//# sourceMappingURL=index.cjs559.js.map

package/dist/index.cjs559.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"index.cjs559.js","sources":["../../../node_modules/lodash-es/_setToArray.js"],"sourcesContent":["/**\n * Converts `set` to an array of its values.\n *\n * @private\n * @param {Object} set The set to convert.\n * @returns {Array} Returns the values.\n */\nfunction setToArray(set) {\n  var index = -1,\n      result = Array(set.size);\n\n  set.forEach(function(value) {\n    result[++index] = value;\n  });\n  return result;\n}\n\nexport default setToArray;\n"],"names":["setToArray","set","index","result","value"],"mappings":"aAOA,SAASA,EAAWC,EAAK,CACvB,IAAIC,EAAQ,GACRC,EAAS,MAAMF,EAAI,IAAI,EAE3B,OAAAA,EAAI,QAAQ,SAASG,EAAO,CAC1BD,EAAO,EAAED,CAAK,EAAIE,CACpB,CAAC,EACMD,CACT","x_google_ignoreList":[0]}

package/dist/index.cjs56.js

@@ -1,2 +1,2 @@
-"use strict";const t=require("./index.cjs94.js"),c=require("./index.cjs38.js"),l=require("./index.cjs161.js"),f=require("./index.cjs229.js"),q=require("./index.cjs228.js"),v=t(function(s,e,o){var i,u;return c(e)?u=e:(e=q(e),i=e.slice(0,-1),e=e[e.length-1]),l(s,function(r){var n=u;if(!n){if(i&&i.length&&(r=f(r,i)),r==null)return;n=r[e]}return n==null?n:n.apply(r,o)})});module.exports=v;
+"use strict";const t=require("./index.cjs94.js"),c=require("./index.cjs38.js"),l=require("./index.cjs161.js"),f=require("./index.cjs246.js"),q=require("./index.cjs245.js"),v=t(function(s,e,o){var i,u;return c(e)?u=e:(e=q(e),i=e.slice(0,-1),e=e[e.length-1]),l(s,function(r){var n=u;if(!n){if(i&&i.length&&(r=f(r,i)),r==null)return;n=r[e]}return n==null?n:n.apply(r,o)})});module.exports=v;
 //# sourceMappingURL=index.cjs56.js.map

package/dist/index.cjs560.js

@@ -0,0 +1,2 @@
+"use strict";function u(e,n){for(var t=-1,l=e==null?0:e.length;++t<l;)if(n(e[t],t,e))return!0;return!1}module.exports=u;
+//# sourceMappingURL=index.cjs560.js.map

package/dist/index.cjs560.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"index.cjs560.js","sources":["../../../node_modules/lodash-es/_arraySome.js"],"sourcesContent":["/**\n * A specialized version of `_.some` for arrays without support for iteratee\n * shorthands.\n *\n * @private\n * @param {Array} [array] The array to iterate over.\n * @param {Function} predicate The function invoked per iteration.\n * @returns {boolean} Returns `true` if any element passes the predicate check,\n *  else `false`.\n */\nfunction arraySome(array, predicate) {\n  var index = -1,\n      length = array == null ? 0 : array.length;\n\n  while (++index < length) {\n    if (predicate(array[index], index, array)) {\n      return true;\n    }\n  }\n  return false;\n}\n\nexport default arraySome;\n"],"names":["arraySome","array","predicate","index","length"],"mappings":"aAUA,SAASA,EAAUC,EAAOC,EAAW,CAInC,QAHIC,EAAQ,GACRC,EAASH,GAAS,KAAO,EAAIA,EAAM,OAEhC,EAAEE,EAAQC,GACf,GAAIF,EAAUD,EAAME,CAAK,EAAGA,EAAOF,CAAK,EACtC,MAAO,GAGX,MAAO,EACT","x_google_ignoreList":[0]}

package/dist/index.cjs561.js

@@ -0,0 +1,2 @@
+"use strict";var r=/\s/;function n(t){for(var e=t.length;e--&&r.test(t.charAt(e)););return e}module.exports=n;
+//# sourceMappingURL=index.cjs561.js.map

package/dist/index.cjs561.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"index.cjs561.js","sources":["../../../node_modules/lodash-es/_trimmedEndIndex.js"],"sourcesContent":["/** Used to match a single whitespace character. */\nvar reWhitespace = /\\s/;\n\n/**\n * Used by `_.trim` and `_.trimEnd` to get the index of the last non-whitespace\n * character of `string`.\n *\n * @private\n * @param {string} string The string to inspect.\n * @returns {number} Returns the index of the last non-whitespace character.\n */\nfunction trimmedEndIndex(string) {\n  var index = string.length;\n\n  while (index-- && reWhitespace.test(string.charAt(index))) {}\n  return index;\n}\n\nexport default trimmedEndIndex;\n"],"names":["reWhitespace","trimmedEndIndex","string","index"],"mappings":"aACA,IAAIA,EAAe,KAUnB,SAASC,EAAgBC,EAAQ,CAG/B,QAFIC,EAAQD,EAAO,OAEZC,KAAWH,EAAa,KAAKE,EAAO,OAAOC,CAAK,CAAC,GAAG,CAC3D,OAAOA,CACT","x_google_ignoreList":[0]}

package/dist/index.cjs562.js

@@ -0,0 +1,2 @@
+"use strict";const e=require("./index.cjs438.js"),o=require("./index.cjs563.js"),n=require("./index.cjs429.js"),s=require("./index.cjs517.js");var c=/[\\^$.*+?()[\]{}|]/g,a=/^\[object .+?Constructor\]$/,i=Function.prototype,p=Object.prototype,u=i.toString,f=p.hasOwnProperty,g=RegExp("^"+u.call(f).replace(c,"\\$&").replace(/hasOwnProperty|(function).*?(?=\\\()| for .+?(?=\\\])/g,"$1.*?")+"$");function v(r){if(!n(r)||o(r))return!1;var t=e(r)?g:a;return t.test(s(r))}module.exports=v;
+//# sourceMappingURL=index.cjs562.js.map

package/dist/index.cjs562.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"index.cjs562.js","sources":["../../../node_modules/lodash-es/_baseIsNative.js"],"sourcesContent":["import isFunction from './isFunction.js';\nimport isMasked from './_isMasked.js';\nimport isObject from './isObject.js';\nimport toSource from './_toSource.js';\n\n/**\n * Used to match `RegExp`\n * [syntax characters](http://ecma-international.org/ecma-262/7.0/#sec-patterns).\n */\nvar reRegExpChar = /[\\\\^$.*+?()[\\]{}|]/g;\n\n/** Used to detect host constructors (Safari). */\nvar reIsHostCtor = /^\\[object .+?Constructor\\]$/;\n\n/** Used for built-in method references. */\nvar funcProto = Function.prototype,\n    objectProto = Object.prototype;\n\n/** Used to resolve the decompiled source of functions. */\nvar funcToString = funcProto.toString;\n\n/** Used to check objects for own properties. */\nvar hasOwnProperty = objectProto.hasOwnProperty;\n\n/** Used to detect if a method is native. */\nvar reIsNative = RegExp('^' +\n  funcToString.call(hasOwnProperty).replace(reRegExpChar, '\\\\$&')\n  .replace(/hasOwnProperty|(function).*?(?=\\\\\\()| for .+?(?=\\\\\\])/g, '$1.*?') + '$'\n);\n\n/**\n * The base implementation of `_.isNative` without bad shim checks.\n *\n * @private\n * @param {*} value The value to check.\n * @returns {boolean} Returns `true` if `value` is a native function,\n *  else `false`.\n */\nfunction baseIsNative(value) {\n  if (!isObject(value) || isMasked(value)) {\n    return false;\n  }\n  var pattern = isFunction(value) ? reIsNative : reIsHostCtor;\n  return pattern.test(toSource(value));\n}\n\nexport default baseIsNative;\n"],"names":["reRegExpChar","reIsHostCtor","funcProto","objectProto","funcToString","hasOwnProperty","reIsNative","baseIsNative","value","isObject","isMasked","pattern","isFunction","toSource"],"mappings":"+IASA,IAAIA,EAAe,sBAGfC,EAAe,8BAGfC,EAAY,SAAS,UACrBC,EAAc,OAAO,UAGrBC,EAAeF,EAAU,SAGzBG,EAAiBF,EAAY,eAG7BG,EAAa,OAAO,IACtBF,EAAa,KAAKC,CAAc,EAAE,QAAQL,EAAc,MAAM,EAC7D,QAAQ,yDAA0D,OAAO,EAAI,GAChF,EAUA,SAASO,EAAaC,EAAO,CAC3B,GAAI,CAACC,EAASD,CAAK,GAAKE,EAASF,CAAK,EACpC,MAAO,GAET,IAAIG,EAAUC,EAAWJ,CAAK,EAAIF,EAAaL,EAC/C,OAAOU,EAAQ,KAAKE,EAASL,CAAK,CAAC,CACrC","x_google_ignoreList":[0]}

package/dist/index.cjs563.js

@@ -0,0 +1,2 @@
+"use strict";const r=require("./index.cjs597.js");var s=function(){var e=/[^.]+$/.exec(r&&r.keys&&r.keys.IE_PROTO||"");return e?"Symbol(src)_1."+e:""}();function t(e){return!!s&&s in e}module.exports=t;
+//# sourceMappingURL=index.cjs563.js.map

package/dist/index.cjs563.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"index.cjs563.js","sources":["../../../node_modules/lodash-es/_isMasked.js"],"sourcesContent":["import coreJsData from './_coreJsData.js';\n\n/** Used to detect methods masquerading as native. */\nvar maskSrcKey = (function() {\n  var uid = /[^.]+$/.exec(coreJsData && coreJsData.keys && coreJsData.keys.IE_PROTO || '');\n  return uid ? ('Symbol(src)_1.' + uid) : '';\n}());\n\n/**\n * Checks if `func` has its source masked.\n *\n * @private\n * @param {Function} func The function to check.\n * @returns {boolean} Returns `true` if `func` is masked, else `false`.\n */\nfunction isMasked(func) {\n  return !!maskSrcKey && (maskSrcKey in func);\n}\n\nexport default isMasked;\n"],"names":["maskSrcKey","uid","coreJsData","isMasked","func"],"mappings":"kDAGA,IAAIA,EAAc,UAAW,CAC3B,IAAIC,EAAM,SAAS,KAAKC,GAAcA,EAAW,MAAQA,EAAW,KAAK,UAAY,EAAE,EACvF,OAAOD,EAAO,iBAAmBA,EAAO,EAC1C,EAAA,EASA,SAASE,EAASC,EAAM,CACtB,MAAO,CAAC,CAACJ,GAAeA,KAAcI,CACxC","x_google_ignoreList":[0]}

package/dist/index.cjs564.js

@@ -0,0 +1,9 @@
+"use strict";Object.defineProperty(exports,Symbol.toStringTag,{value:"Module"});const q=require("./index.cjs565.js");var It;function Te(){if(It)return q.__exports;It=1,Object.defineProperty(q.__exports,"__esModule",{value:!0});const Nt=Object.prototype.toString;function B(i){const e=Nt.call(i);return e.endsWith("Array]")&&!e.includes("Big")}function Tt(i,e={}){if(!B(i))throw new TypeError("input must be an array");if(i.length===0)throw new TypeError("input must not be empty");const{fromIndex:s=0,toIndex:t=i.length}=e;if(s<0||s>=i.length||!Number.isInteger(s))throw new Error("fromIndex must be a positive integer smaller than length");if(t<=s||t>i.length||!Number.isInteger(t))throw new Error("toIndex must be an integer greater than fromIndex and at most equal to length");let o=i[s];for(let r=s+1;r<t;r++)i[r]>o&&(o=i[r]);return o}function qt(i,e={}){if(!B(i))throw new TypeError("input must be an array");if(i.length===0)throw new TypeError("input must not be empty");const{fromIndex:s=0,toIndex:t=i.length}=e;if(s<0||s>=i.length||!Number.isInteger(s))throw new Error("fromIndex must be a positive integer smaller than length");if(t<=s||t>i.length||!Number.isInteger(t))throw new Error("toIndex must be an integer greater than fromIndex and at most equal to length");let o=i[s];for(let r=s+1;r<t;r++)i[r]<o&&(o=i[r]);return o}function dt(i,e={}){if(B(i)){if(i.length===0)throw new TypeError("input must not be empty")}else throw new TypeError("input must be an array");let s;if(e.output!==void 0){if(!B(e.output))throw new TypeError("output option must be an array if specified");s=e.output}else s=new Array(i.length);const t=qt(i),o=Tt(i);if(t===o)throw new RangeError("minimum and maximum input values are equal. Cannot rescale a constant array");const{min:r=e.autoMinMax?t:0,max:l=e.autoMinMax?o:1}=e;if(r>=l)throw new RangeError("min option must be smaller than max option");const n=(l-r)/(o-t);for(let h=0;h<i.length;h++)s[h]=(i[h]-t)*n+r;return s}const rt=" ".repeat(2),yt=" ".repeat(4);function _t(){return Mt(this)}function Mt(i,e={}){const{maxRows:s=15,maxColumns:t=10,maxNumSize:o=8,padMinus:r="auto"}=e;return`${i.constructor.name} {
+${rt}[
+${yt}${Dt(i,s,t,o,r)}
+${rt}]
+${rt}rows: ${i.rows}
+${rt}columns: ${i.columns}
+}`}function Dt(i,e,s,t,o){const{rows:r,columns:l}=i,n=Math.min(r,e),h=Math.min(l,s),u=[];if(o==="auto"){o=!1;t:for(let w=0;w<n;w++)for(let f=0;f<h;f++)if(i.get(w,f)<0){o=!0;break t}}for(let w=0;w<n;w++){let f=[];for(let a=0;a<h;a++)f.push(Ft(i.get(w,a),t,o));u.push(`${f.join(" ")}`)}return h!==l&&(u[u.length-1]+=` ... ${l-s} more columns`),n!==r&&u.push(`... ${r-e} more rows`),u.join(`
+${yt}`)}function Ft(i,e,s){return(i>=0&&s?` ${jt(i,e-1)}`:jt(i,e)).padEnd(e)}function jt(i,e){let s=i.toString();if(s.length<=e)return s;let t=i.toFixed(e);if(t.length>e&&(t=i.toFixed(Math.max(0,e-(t.length-e)))),t.length<=e&&!t.startsWith("0.000")&&!t.startsWith("-0.000"))return t;let o=i.toExponential(e);return o.length>e&&(o=i.toExponential(Math.max(0,e-(o.length-e)))),o.slice(0)}function zt(i,e){i.prototype.add=function(t){return typeof t=="number"?this.addS(t):this.addM(t)},i.prototype.addS=function(t){for(let o=0;o<this.rows;o++)for(let r=0;r<this.columns;r++)this.set(o,r,this.get(o,r)+t);return this},i.prototype.addM=function(t){if(t=e.checkMatrix(t),this.rows!==t.rows||this.columns!==t.columns)throw new RangeError("Matrices dimensions must be equal");for(let o=0;o<this.rows;o++)for(let r=0;r<this.columns;r++)this.set(o,r,this.get(o,r)+t.get(o,r));return this},i.add=function(t,o){return new e(t).add(o)},i.prototype.sub=function(t){return typeof t=="number"?this.subS(t):this.subM(t)},i.prototype.subS=function(t){for(let o=0;o<this.rows;o++)for(let r=0;r<this.columns;r++)this.set(o,r,this.get(o,r)-t);return this},i.prototype.subM=function(t){if(t=e.checkMatrix(t),this.rows!==t.rows||this.columns!==t.columns)throw new RangeError("Matrices dimensions must be equal");for(let o=0;o<this.rows;o++)for(let r=0;r<this.columns;r++)this.set(o,r,this.get(o,r)-t.get(o,r));return this},i.sub=function(t,o){return new e(t).sub(o)},i.prototype.subtract=i.prototype.sub,i.prototype.subtractS=i.prototype.subS,i.prototype.subtractM=i.prototype.subM,i.subtract=i.sub,i.prototype.mul=function(t){return typeof t=="number"?this.mulS(t):this.mulM(t)},i.prototype.mulS=function(t){for(let o=0;o<this.rows;o++)for(let r=0;r<this.columns;r++)this.set(o,r,this.get(o,r)*t);return this},i.prototype.mulM=function(t){if(t=e.checkMatrix(t),this.rows!==t.rows||this.columns!==t.columns)throw new RangeError("Matrices dimensions must be equal");for(let o=0;o<this.rows;o++)for(let r=0;r<this.columns;r++)this.set(o,r,this.get(o,r)*t.get(o,r));return this},i.mul=function(t,o){return new e(t).mul(o)},i.prototype.multiply=i.prototype.mul,i.prototype.multiplyS=i.prototype.mulS,i.prototype.multiplyM=i.prototype.mulM,i.multiply=i.mul,i.prototype.div=function(t){return typeof t=="number"?this.divS(t):this.divM(t)},i.prototype.divS=function(t){for(let o=0;o<this.rows;o++)for(let r=0;r<this.columns;r++)this.set(o,r,this.get(o,r)/t);return this},i.prototype.divM=function(t){if(t=e.checkMatrix(t),this.rows!==t.rows||this.columns!==t.columns)throw new RangeError("Matrices dimensions must be equal");for(let o=0;o<this.rows;o++)for(let r=0;r<this.columns;r++)this.set(o,r,this.get(o,r)/t.get(o,r));return this},i.div=function(t,o){return new e(t).div(o)},i.prototype.divide=i.prototype.div,i.prototype.divideS=i.prototype.divS,i.prototype.divideM=i.prototype.divM,i.divide=i.div,i.prototype.mod=function(t){return typeof t=="number"?this.modS(t):this.modM(t)},i.prototype.modS=function(t){for(let o=0;o<this.rows;o++)for(let r=0;r<this.columns;r++)this.set(o,r,this.get(o,r)%t);return this},i.prototype.modM=function(t){if(t=e.checkMatrix(t),this.rows!==t.rows||this.columns!==t.columns)throw new RangeError("Matrices dimensions must be equal");for(let o=0;o<this.rows;o++)for(let r=0;r<this.columns;r++)this.set(o,r,this.get(o,r)%t.get(o,r));return this},i.mod=function(t,o){return new e(t).mod(o)},i.prototype.modulus=i.prototype.mod,i.prototype.modulusS=i.prototype.modS,i.prototype.modulusM=i.prototype.modM,i.modulus=i.mod,i.prototype.and=function(t){return typeof t=="number"?this.andS(t):this.andM(t)},i.prototype.andS=function(t){for(let o=0;o<this.rows;o++)for(let r=0;r<this.columns;r++)this.set(o,r,this.get(o,r)&t);return this},i.prototype.andM=function(t){if(t=e.checkMatrix(t),this.rows!==t.rows||this.columns!==t.columns)throw new RangeError("Matrices dimensions must be equal");for(let o=0;o<this.rows;o++)for(let r=0;r<this.columns;r++)this.set(o,r,this.get(o,r)&t.get(o,r));return this},i.and=function(t,o){return new e(t).and(o)},i.prototype.or=function(t){return typeof t=="number"?this.orS(t):this.orM(t)},i.prototype.orS=function(t){for(let o=0;o<this.rows;o++)for(let r=0;r<this.columns;r++)this.set(o,r,this.get(o,r)|t);return this},i.prototype.orM=function(t){if(t=e.checkMatrix(t),this.rows!==t.rows||this.columns!==t.columns)throw new RangeError("Matrices dimensions must be equal");for(let o=0;o<this.rows;o++)for(let r=0;r<this.columns;r++)this.set(o,r,this.get(o,r)|t.get(o,r));return this},i.or=function(t,o){return new e(t).or(o)},i.prototype.xor=function(t){return typeof t=="number"?this.xorS(t):this.xorM(t)},i.prototype.xorS=function(t){for(let o=0;o<this.rows;o++)for(let r=0;r<this.columns;r++)this.set(o,r,this.get(o,r)^t);return this},i.prototype.xorM=function(t){if(t=e.checkMatrix(t),this.rows!==t.rows||this.columns!==t.columns)throw new RangeError("Matrices dimensions must be equal");for(let o=0;o<this.rows;o++)for(let r=0;r<this.columns;r++)this.set(o,r,this.get(o,r)^t.get(o,r));return this},i.xor=function(t,o){return new e(t).xor(o)},i.prototype.leftShift=function(t){return typeof t=="number"?this.leftShiftS(t):this.leftShiftM(t)},i.prototype.leftShiftS=function(t){for(let o=0;o<this.rows;o++)for(let r=0;r<this.columns;r++)this.set(o,r,this.get(o,r)<<t);return this},i.prototype.leftShiftM=function(t){if(t=e.checkMatrix(t),this.rows!==t.rows||this.columns!==t.columns)throw new RangeError("Matrices dimensions must be equal");for(let o=0;o<this.rows;o++)for(let r=0;r<this.columns;r++)this.set(o,r,this.get(o,r)<<t.get(o,r));return this},i.leftShift=function(t,o){return new e(t).leftShift(o)},i.prototype.signPropagatingRightShift=function(t){return typeof t=="number"?this.signPropagatingRightShiftS(t):this.signPropagatingRightShiftM(t)},i.prototype.signPropagatingRightShiftS=function(t){for(let o=0;o<this.rows;o++)for(let r=0;r<this.columns;r++)this.set(o,r,this.get(o,r)>>t);return this},i.prototype.signPropagatingRightShiftM=function(t){if(t=e.checkMatrix(t),this.rows!==t.rows||this.columns!==t.columns)throw new RangeError("Matrices dimensions must be equal");for(let o=0;o<this.rows;o++)for(let r=0;r<this.columns;r++)this.set(o,r,this.get(o,r)>>t.get(o,r));return this},i.signPropagatingRightShift=function(t,o){return new e(t).signPropagatingRightShift(o)},i.prototype.rightShift=function(t){return typeof t=="number"?this.rightShiftS(t):this.rightShiftM(t)},i.prototype.rightShiftS=function(t){for(let o=0;o<this.rows;o++)for(let r=0;r<this.columns;r++)this.set(o,r,this.get(o,r)>>>t);return this},i.prototype.rightShiftM=function(t){if(t=e.checkMatrix(t),this.rows!==t.rows||this.columns!==t.columns)throw new RangeError("Matrices dimensions must be equal");for(let o=0;o<this.rows;o++)for(let r=0;r<this.columns;r++)this.set(o,r,this.get(o,r)>>>t.get(o,r));return this},i.rightShift=function(t,o){return new e(t).rightShift(o)},i.prototype.zeroFillRightShift=i.prototype.rightShift,i.prototype.zeroFillRightShiftS=i.prototype.rightShiftS,i.prototype.zeroFillRightShiftM=i.prototype.rightShiftM,i.zeroFillRightShift=i.rightShift,i.prototype.not=function(){for(let t=0;t<this.rows;t++)for(let o=0;o<this.columns;o++)this.set(t,o,~this.get(t,o));return this},i.not=function(t){return new e(t).not()},i.prototype.abs=function(){for(let t=0;t<this.rows;t++)for(let o=0;o<this.columns;o++)this.set(t,o,Math.abs(this.get(t,o)));return this},i.abs=function(t){return new e(t).abs()},i.prototype.acos=function(){for(let t=0;t<this.rows;t++)for(let o=0;o<this.columns;o++)this.set(t,o,Math.acos(this.get(t,o)));return this},i.acos=function(t){return new e(t).acos()},i.prototype.acosh=function(){for(let t=0;t<this.rows;t++)for(let o=0;o<this.columns;o++)this.set(t,o,Math.acosh(this.get(t,o)));return this},i.acosh=function(t){return new e(t).acosh()},i.prototype.asin=function(){for(let t=0;t<this.rows;t++)for(let o=0;o<this.columns;o++)this.set(t,o,Math.asin(this.get(t,o)));return this},i.asin=function(t){return new e(t).asin()},i.prototype.asinh=function(){for(let t=0;t<this.rows;t++)for(let o=0;o<this.columns;o++)this.set(t,o,Math.asinh(this.get(t,o)));return this},i.asinh=function(t){return new e(t).asinh()},i.prototype.atan=function(){for(let t=0;t<this.rows;t++)for(let o=0;o<this.columns;o++)this.set(t,o,Math.atan(this.get(t,o)));return this},i.atan=function(t){return new e(t).atan()},i.prototype.atanh=function(){for(let t=0;t<this.rows;t++)for(let o=0;o<this.columns;o++)this.set(t,o,Math.atanh(this.get(t,o)));return this},i.atanh=function(t){return new e(t).atanh()},i.prototype.cbrt=function(){for(let t=0;t<this.rows;t++)for(let o=0;o<this.columns;o++)this.set(t,o,Math.cbrt(this.get(t,o)));return this},i.cbrt=function(t){return new e(t).cbrt()},i.prototype.ceil=function(){for(let t=0;t<this.rows;t++)for(let o=0;o<this.columns;o++)this.set(t,o,Math.ceil(this.get(t,o)));return this},i.ceil=function(t){return new e(t).ceil()},i.prototype.clz32=function(){for(let t=0;t<this.rows;t++)for(let o=0;o<this.columns;o++)this.set(t,o,Math.clz32(this.get(t,o)));return this},i.clz32=function(t){return new e(t).clz32()},i.prototype.cos=function(){for(let t=0;t<this.rows;t++)for(let o=0;o<this.columns;o++)this.set(t,o,Math.cos(this.get(t,o)));return this},i.cos=function(t){return new e(t).cos()},i.prototype.cosh=function(){for(let t=0;t<this.rows;t++)for(let o=0;o<this.columns;o++)this.set(t,o,Math.cosh(this.get(t,o)));return this},i.cosh=function(t){return new e(t).cosh()},i.prototype.exp=function(){for(let t=0;t<this.rows;t++)for(let o=0;o<this.columns;o++)this.set(t,o,Math.exp(this.get(t,o)));return this},i.exp=function(t){return new e(t).exp()},i.prototype.expm1=function(){for(let t=0;t<this.rows;t++)for(let o=0;o<this.columns;o++)this.set(t,o,Math.expm1(this.get(t,o)));return this},i.expm1=function(t){return new e(t).expm1()},i.prototype.floor=function(){for(let t=0;t<this.rows;t++)for(let o=0;o<this.columns;o++)this.set(t,o,Math.floor(this.get(t,o)));return this},i.floor=function(t){return new e(t).floor()},i.prototype.fround=function(){for(let t=0;t<this.rows;t++)for(let o=0;o<this.columns;o++)this.set(t,o,Math.fround(this.get(t,o)));return this},i.fround=function(t){return new e(t).fround()},i.prototype.log=function(){for(let t=0;t<this.rows;t++)for(let o=0;o<this.columns;o++)this.set(t,o,Math.log(this.get(t,o)));return this},i.log=function(t){return new e(t).log()},i.prototype.log1p=function(){for(let t=0;t<this.rows;t++)for(let o=0;o<this.columns;o++)this.set(t,o,Math.log1p(this.get(t,o)));return this},i.log1p=function(t){return new e(t).log1p()},i.prototype.log10=function(){for(let t=0;t<this.rows;t++)for(let o=0;o<this.columns;o++)this.set(t,o,Math.log10(this.get(t,o)));return this},i.log10=function(t){return new e(t).log10()},i.prototype.log2=function(){for(let t=0;t<this.rows;t++)for(let o=0;o<this.columns;o++)this.set(t,o,Math.log2(this.get(t,o)));return this},i.log2=function(t){return new e(t).log2()},i.prototype.round=function(){for(let t=0;t<this.rows;t++)for(let o=0;o<this.columns;o++)this.set(t,o,Math.round(this.get(t,o)));return this},i.round=function(t){return new e(t).round()},i.prototype.sign=function(){for(let t=0;t<this.rows;t++)for(let o=0;o<this.columns;o++)this.set(t,o,Math.sign(this.get(t,o)));return this},i.sign=function(t){return new e(t).sign()},i.prototype.sin=function(){for(let t=0;t<this.rows;t++)for(let o=0;o<this.columns;o++)this.set(t,o,Math.sin(this.get(t,o)));return this},i.sin=function(t){return new e(t).sin()},i.prototype.sinh=function(){for(let t=0;t<this.rows;t++)for(let o=0;o<this.columns;o++)this.set(t,o,Math.sinh(this.get(t,o)));return this},i.sinh=function(t){return new e(t).sinh()},i.prototype.sqrt=function(){for(let t=0;t<this.rows;t++)for(let o=0;o<this.columns;o++)this.set(t,o,Math.sqrt(this.get(t,o)));return this},i.sqrt=function(t){return new e(t).sqrt()},i.prototype.tan=function(){for(let t=0;t<this.rows;t++)for(let o=0;o<this.columns;o++)this.set(t,o,Math.tan(this.get(t,o)));return this},i.tan=function(t){return new e(t).tan()},i.prototype.tanh=function(){for(let t=0;t<this.rows;t++)for(let o=0;o<this.columns;o++)this.set(t,o,Math.tanh(this.get(t,o)));return this},i.tanh=function(t){return new e(t).tanh()},i.prototype.trunc=function(){for(let t=0;t<this.rows;t++)for(let o=0;o<this.columns;o++)this.set(t,o,Math.trunc(this.get(t,o)));return this},i.trunc=function(t){return new e(t).trunc()},i.pow=function(t,o){return new e(t).pow(o)},i.prototype.pow=function(t){return typeof t=="number"?this.powS(t):this.powM(t)},i.prototype.powS=function(t){for(let o=0;o<this.rows;o++)for(let r=0;r<this.columns;r++)this.set(o,r,this.get(o,r)**t);return this},i.prototype.powM=function(t){if(t=e.checkMatrix(t),this.rows!==t.rows||this.columns!==t.columns)throw new RangeError("Matrices dimensions must be equal");for(let o=0;o<this.rows;o++)for(let r=0;r<this.columns;r++)this.set(o,r,this.get(o,r)**t.get(o,r));return this}}function J(i,e,s){let t=s?i.rows:i.rows-1;if(e<0||e>t)throw new RangeError("Row index out of range")}function G(i,e,s){let t=s?i.columns:i.columns-1;if(e<0||e>t)throw new RangeError("Column index out of range")}function tt(i,e){if(e.to1DArray&&(e=e.to1DArray()),e.length!==i.columns)throw new RangeError("vector size must be the same as the number of columns");return e}function et(i,e){if(e.to1DArray&&(e=e.to1DArray()),e.length!==i.rows)throw new RangeError("vector size must be the same as the number of rows");return e}function gt(i,e){if(!B(e))throw new TypeError("row indices must be an array");for(let s=0;s<e.length;s++)if(e[s]<0||e[s]>=i.rows)throw new RangeError("row indices are out of range")}function wt(i,e){if(!B(e))throw new TypeError("column indices must be an array");for(let s=0;s<e.length;s++)if(e[s]<0||e[s]>=i.columns)throw new RangeError("column indices are out of range")}function mt(i,e,s,t,o){if(arguments.length!==5)throw new RangeError("expected 4 arguments");if(nt("startRow",e),nt("endRow",s),nt("startColumn",t),nt("endColumn",o),e>s||t>o||e<0||e>=i.rows||s<0||s>=i.rows||t<0||t>=i.columns||o<0||o>=i.columns)throw new RangeError("Submatrix indices are out of range")}function it(i,e=0){let s=[];for(let t=0;t<i;t++)s.push(e);return s}function nt(i,e){if(typeof e!="number")throw new TypeError(`${i} must be a number`)}function st(i){if(i.isEmpty())throw new Error("Empty matrix has no elements to index")}function Vt(i){let e=it(i.rows);for(let s=0;s<i.rows;++s)for(let t=0;t<i.columns;++t)e[s]+=i.get(s,t);return e}function Ct(i){let e=it(i.columns);for(let s=0;s<i.rows;++s)for(let t=0;t<i.columns;++t)e[t]+=i.get(s,t);return e}function Pt(i){let e=0;for(let s=0;s<i.rows;s++)for(let t=0;t<i.columns;t++)e+=i.get(s,t);return e}function $t(i){let e=it(i.rows,1);for(let s=0;s<i.rows;++s)for(let t=0;t<i.columns;++t)e[s]*=i.get(s,t);return e}function Lt(i){let e=it(i.columns,1);for(let s=0;s<i.rows;++s)for(let t=0;t<i.columns;++t)e[t]*=i.get(s,t);return e}function Bt(i){let e=1;for(let s=0;s<i.rows;s++)for(let t=0;t<i.columns;t++)e*=i.get(s,t);return e}function Ut(i,e,s){const t=i.rows,o=i.columns,r=[];for(let l=0;l<t;l++){let n=0,h=0,u=0;for(let w=0;w<o;w++)u=i.get(l,w)-s[l],n+=u,h+=u*u;e?r.push((h-n*n/o)/(o-1)):r.push((h-n*n/o)/o)}return r}function Ot(i,e,s){const t=i.rows,o=i.columns,r=[];for(let l=0;l<o;l++){let n=0,h=0,u=0;for(let w=0;w<t;w++)u=i.get(w,l)-s[l],n+=u,h+=u*u;e?r.push((h-n*n/t)/(t-1)):r.push((h-n*n/t)/t)}return r}function Qt(i,e,s){const t=i.rows,o=i.columns,r=t*o;let l=0,n=0,h=0;for(let u=0;u<t;u++)for(let w=0;w<o;w++)h=i.get(u,w)-s,l+=h,n+=h*h;return e?(n-l*l/r)/(r-1):(n-l*l/r)/r}function Wt(i,e){for(let s=0;s<i.rows;s++)for(let t=0;t<i.columns;t++)i.set(s,t,i.get(s,t)-e[s])}function Yt(i,e){for(let s=0;s<i.rows;s++)for(let t=0;t<i.columns;t++)i.set(s,t,i.get(s,t)-e[t])}function Jt(i,e){for(let s=0;s<i.rows;s++)for(let t=0;t<i.columns;t++)i.set(s,t,i.get(s,t)-e)}function Gt(i){const e=[];for(let s=0;s<i.rows;s++){let t=0;for(let o=0;o<i.columns;o++)t+=i.get(s,o)**2/(i.columns-1);e.push(Math.sqrt(t))}return e}function Kt(i,e){for(let s=0;s<i.rows;s++)for(let t=0;t<i.columns;t++)i.set(s,t,i.get(s,t)/e[s])}function Xt(i){const e=[];for(let s=0;s<i.columns;s++){let t=0;for(let o=0;o<i.rows;o++)t+=i.get(o,s)**2/(i.rows-1);e.push(Math.sqrt(t))}return e}function Zt(i,e){for(let s=0;s<i.rows;s++)for(let t=0;t<i.columns;t++)i.set(s,t,i.get(s,t)/e[t])}function xt(i){const e=i.size-1;let s=0;for(let t=0;t<i.columns;t++)for(let o=0;o<i.rows;o++)s+=i.get(o,t)**2/e;return Math.sqrt(s)}function At(i,e){for(let s=0;s<i.rows;s++)for(let t=0;t<i.columns;t++)i.set(s,t,i.get(s,t)/e)}class T{static from1DArray(e,s,t){if(e*s!==t.length)throw new RangeError("data length does not match given dimensions");let r=new S(e,s);for(let l=0;l<e;l++)for(let n=0;n<s;n++)r.set(l,n,t[l*s+n]);return r}static rowVector(e){let s=new S(1,e.length);for(let t=0;t<e.length;t++)s.set(0,t,e[t]);return s}static columnVector(e){let s=new S(e.length,1);for(let t=0;t<e.length;t++)s.set(t,0,e[t]);return s}static zeros(e,s){return new S(e,s)}static ones(e,s){return new S(e,s).fill(1)}static rand(e,s,t={}){if(typeof t!="object")throw new TypeError("options must be an object");const{random:o=Math.random}=t;let r=new S(e,s);for(let l=0;l<e;l++)for(let n=0;n<s;n++)r.set(l,n,o());return r}static randInt(e,s,t={}){if(typeof t!="object")throw new TypeError("options must be an object");const{min:o=0,max:r=1e3,random:l=Math.random}=t;if(!Number.isInteger(o))throw new TypeError("min must be an integer");if(!Number.isInteger(r))throw new TypeError("max must be an integer");if(o>=r)throw new RangeError("min must be smaller than max");let n=r-o,h=new S(e,s);for(let u=0;u<e;u++)for(let w=0;w<s;w++){let f=o+Math.round(l()*n);h.set(u,w,f)}return h}static eye(e,s,t){s===void 0&&(s=e),t===void 0&&(t=1);let o=Math.min(e,s),r=this.zeros(e,s);for(let l=0;l<o;l++)r.set(l,l,t);return r}static diag(e,s,t){let o=e.length;s===void 0&&(s=o),t===void 0&&(t=s);let r=Math.min(o,s,t),l=this.zeros(s,t);for(let n=0;n<r;n++)l.set(n,n,e[n]);return l}static min(e,s){e=this.checkMatrix(e),s=this.checkMatrix(s);let t=e.rows,o=e.columns,r=new S(t,o);for(let l=0;l<t;l++)for(let n=0;n<o;n++)r.set(l,n,Math.min(e.get(l,n),s.get(l,n)));return r}static max(e,s){e=this.checkMatrix(e),s=this.checkMatrix(s);let t=e.rows,o=e.columns,r=new this(t,o);for(let l=0;l<t;l++)for(let n=0;n<o;n++)r.set(l,n,Math.max(e.get(l,n),s.get(l,n)));return r}static checkMatrix(e){return T.isMatrix(e)?e:new S(e)}static isMatrix(e){return e!=null&&e.klass==="Matrix"}get size(){return this.rows*this.columns}apply(e){if(typeof e!="function")throw new TypeError("callback must be a function");for(let s=0;s<this.rows;s++)for(let t=0;t<this.columns;t++)e.call(this,s,t);return this}to1DArray(){let e=[];for(let s=0;s<this.rows;s++)for(let t=0;t<this.columns;t++)e.push(this.get(s,t));return e}to2DArray(){let e=[];for(let s=0;s<this.rows;s++){e.push([]);for(let t=0;t<this.columns;t++)e[s].push(this.get(s,t))}return e}toJSON(){return this.to2DArray()}isRowVector(){return this.rows===1}isColumnVector(){return this.columns===1}isVector(){return this.rows===1||this.columns===1}isSquare(){return this.rows===this.columns}isEmpty(){return this.rows===0||this.columns===0}isSymmetric(){if(this.isSquare()){for(let e=0;e<this.rows;e++)for(let s=0;s<=e;s++)if(this.get(e,s)!==this.get(s,e))return!1;return!0}return!1}isDistance(){if(!this.isSymmetric())return!1;for(let e=0;e<this.rows;e++)if(this.get(e,e)!==0)return!1;return!0}isEchelonForm(){let e=0,s=0,t=-1,o=!0,r=!1;for(;e<this.rows&&o;){for(s=0,r=!1;s<this.columns&&r===!1;)this.get(e,s)===0?s++:this.get(e,s)===1&&s>t?(r=!0,t=s):(o=!1,r=!0);e++}return o}isReducedEchelonForm(){let e=0,s=0,t=-1,o=!0,r=!1;for(;e<this.rows&&o;){for(s=0,r=!1;s<this.columns&&r===!1;)this.get(e,s)===0?s++:this.get(e,s)===1&&s>t?(r=!0,t=s):(o=!1,r=!0);for(let l=s+1;l<this.rows;l++)this.get(e,l)!==0&&(o=!1);e++}return o}echelonForm(){let e=this.clone(),s=0,t=0;for(;s<e.rows&&t<e.columns;){let o=s;for(let r=s;r<e.rows;r++)e.get(r,t)>e.get(o,t)&&(o=r);if(e.get(o,t)===0)t++;else{e.swapRows(s,o);let r=e.get(s,t);for(let l=t;l<e.columns;l++)e.set(s,l,e.get(s,l)/r);for(let l=s+1;l<e.rows;l++){let n=e.get(l,t)/e.get(s,t);e.set(l,t,0);for(let h=t+1;h<e.columns;h++)e.set(l,h,e.get(l,h)-e.get(s,h)*n)}s++,t++}}return e}reducedEchelonForm(){let e=this.echelonForm(),s=e.columns,t=e.rows,o=t-1;for(;o>=0;)if(e.maxRow(o)===0)o--;else{let r=0,l=!1;for(;r<t&&l===!1;)e.get(o,r)===1?l=!0:r++;for(let n=0;n<o;n++){let h=e.get(n,r);for(let u=r;u<s;u++){let w=e.get(n,u)-h*e.get(o,u);e.set(n,u,w)}}o--}return e}set(){throw new Error("set method is unimplemented")}get(){throw new Error("get method is unimplemented")}repeat(e={}){if(typeof e!="object")throw new TypeError("options must be an object");const{rows:s=1,columns:t=1}=e;if(!Number.isInteger(s)||s<=0)throw new TypeError("rows must be a positive integer");if(!Number.isInteger(t)||t<=0)throw new TypeError("columns must be a positive integer");let o=new S(this.rows*s,this.columns*t);for(let r=0;r<s;r++)for(let l=0;l<t;l++)o.setSubMatrix(this,this.rows*r,this.columns*l);return o}fill(e){for(let s=0;s<this.rows;s++)for(let t=0;t<this.columns;t++)this.set(s,t,e);return this}neg(){return this.mulS(-1)}getRow(e){J(this,e);let s=[];for(let t=0;t<this.columns;t++)s.push(this.get(e,t));return s}getRowVector(e){return S.rowVector(this.getRow(e))}setRow(e,s){J(this,e),s=tt(this,s);for(let t=0;t<this.columns;t++)this.set(e,t,s[t]);return this}swapRows(e,s){J(this,e),J(this,s);for(let t=0;t<this.columns;t++){let o=this.get(e,t);this.set(e,t,this.get(s,t)),this.set(s,t,o)}return this}getColumn(e){G(this,e);let s=[];for(let t=0;t<this.rows;t++)s.push(this.get(t,e));return s}getColumnVector(e){return S.columnVector(this.getColumn(e))}setColumn(e,s){G(this,e),s=et(this,s);for(let t=0;t<this.rows;t++)this.set(t,e,s[t]);return this}swapColumns(e,s){G(this,e),G(this,s);for(let t=0;t<this.rows;t++){let o=this.get(t,e);this.set(t,e,this.get(t,s)),this.set(t,s,o)}return this}addRowVector(e){e=tt(this,e);for(let s=0;s<this.rows;s++)for(let t=0;t<this.columns;t++)this.set(s,t,this.get(s,t)+e[t]);return this}subRowVector(e){e=tt(this,e);for(let s=0;s<this.rows;s++)for(let t=0;t<this.columns;t++)this.set(s,t,this.get(s,t)-e[t]);return this}mulRowVector(e){e=tt(this,e);for(let s=0;s<this.rows;s++)for(let t=0;t<this.columns;t++)this.set(s,t,this.get(s,t)*e[t]);return this}divRowVector(e){e=tt(this,e);for(let s=0;s<this.rows;s++)for(let t=0;t<this.columns;t++)this.set(s,t,this.get(s,t)/e[t]);return this}addColumnVector(e){e=et(this,e);for(let s=0;s<this.rows;s++)for(let t=0;t<this.columns;t++)this.set(s,t,this.get(s,t)+e[s]);return this}subColumnVector(e){e=et(this,e);for(let s=0;s<this.rows;s++)for(let t=0;t<this.columns;t++)this.set(s,t,this.get(s,t)-e[s]);return this}mulColumnVector(e){e=et(this,e);for(let s=0;s<this.rows;s++)for(let t=0;t<this.columns;t++)this.set(s,t,this.get(s,t)*e[s]);return this}divColumnVector(e){e=et(this,e);for(let s=0;s<this.rows;s++)for(let t=0;t<this.columns;t++)this.set(s,t,this.get(s,t)/e[s]);return this}mulRow(e,s){J(this,e);for(let t=0;t<this.columns;t++)this.set(e,t,this.get(e,t)*s);return this}mulColumn(e,s){G(this,e);for(let t=0;t<this.rows;t++)this.set(t,e,this.get(t,e)*s);return this}max(e){if(this.isEmpty())return NaN;switch(e){case"row":{const s=new Array(this.rows).fill(Number.NEGATIVE_INFINITY);for(let t=0;t<this.rows;t++)for(let o=0;o<this.columns;o++)this.get(t,o)>s[t]&&(s[t]=this.get(t,o));return s}case"column":{const s=new Array(this.columns).fill(Number.NEGATIVE_INFINITY);for(let t=0;t<this.rows;t++)for(let o=0;o<this.columns;o++)this.get(t,o)>s[o]&&(s[o]=this.get(t,o));return s}case void 0:{let s=this.get(0,0);for(let t=0;t<this.rows;t++)for(let o=0;o<this.columns;o++)this.get(t,o)>s&&(s=this.get(t,o));return s}default:throw new Error(`invalid option: ${e}`)}}maxIndex(){st(this);let e=this.get(0,0),s=[0,0];for(let t=0;t<this.rows;t++)for(let o=0;o<this.columns;o++)this.get(t,o)>e&&(e=this.get(t,o),s[0]=t,s[1]=o);return s}min(e){if(this.isEmpty())return NaN;switch(e){case"row":{const s=new Array(this.rows).fill(Number.POSITIVE_INFINITY);for(let t=0;t<this.rows;t++)for(let o=0;o<this.columns;o++)this.get(t,o)<s[t]&&(s[t]=this.get(t,o));return s}case"column":{const s=new Array(this.columns).fill(Number.POSITIVE_INFINITY);for(let t=0;t<this.rows;t++)for(let o=0;o<this.columns;o++)this.get(t,o)<s[o]&&(s[o]=this.get(t,o));return s}case void 0:{let s=this.get(0,0);for(let t=0;t<this.rows;t++)for(let o=0;o<this.columns;o++)this.get(t,o)<s&&(s=this.get(t,o));return s}default:throw new Error(`invalid option: ${e}`)}}minIndex(){st(this);let e=this.get(0,0),s=[0,0];for(let t=0;t<this.rows;t++)for(let o=0;o<this.columns;o++)this.get(t,o)<e&&(e=this.get(t,o),s[0]=t,s[1]=o);return s}maxRow(e){if(J(this,e),this.isEmpty())return NaN;let s=this.get(e,0);for(let t=1;t<this.columns;t++)this.get(e,t)>s&&(s=this.get(e,t));return s}maxRowIndex(e){J(this,e),st(this);let s=this.get(e,0),t=[e,0];for(let o=1;o<this.columns;o++)this.get(e,o)>s&&(s=this.get(e,o),t[1]=o);return t}minRow(e){if(J(this,e),this.isEmpty())return NaN;let s=this.get(e,0);for(let t=1;t<this.columns;t++)this.get(e,t)<s&&(s=this.get(e,t));return s}minRowIndex(e){J(this,e),st(this);let s=this.get(e,0),t=[e,0];for(let o=1;o<this.columns;o++)this.get(e,o)<s&&(s=this.get(e,o),t[1]=o);return t}maxColumn(e){if(G(this,e),this.isEmpty())return NaN;let s=this.get(0,e);for(let t=1;t<this.rows;t++)this.get(t,e)>s&&(s=this.get(t,e));return s}maxColumnIndex(e){G(this,e),st(this);let s=this.get(0,e),t=[0,e];for(let o=1;o<this.rows;o++)this.get(o,e)>s&&(s=this.get(o,e),t[0]=o);return t}minColumn(e){if(G(this,e),this.isEmpty())return NaN;let s=this.get(0,e);for(let t=1;t<this.rows;t++)this.get(t,e)<s&&(s=this.get(t,e));return s}minColumnIndex(e){G(this,e),st(this);let s=this.get(0,e),t=[0,e];for(let o=1;o<this.rows;o++)this.get(o,e)<s&&(s=this.get(o,e),t[0]=o);return t}diag(){let e=Math.min(this.rows,this.columns),s=[];for(let t=0;t<e;t++)s.push(this.get(t,t));return s}norm(e="frobenius"){switch(e){case"max":return this.max();case"frobenius":return Math.sqrt(this.dot(this));default:throw new RangeError(`unknown norm type: ${e}`)}}cumulativeSum(){let e=0;for(let s=0;s<this.rows;s++)for(let t=0;t<this.columns;t++)e+=this.get(s,t),this.set(s,t,e);return this}dot(e){T.isMatrix(e)&&(e=e.to1DArray());let s=this.to1DArray();if(s.length!==e.length)throw new RangeError("vectors do not have the same size");let t=0;for(let o=0;o<s.length;o++)t+=s[o]*e[o];return t}mmul(e){e=S.checkMatrix(e);let s=this.rows,t=this.columns,o=e.columns,r=new S(s,o),l=new Float64Array(t);for(let n=0;n<o;n++){for(let h=0;h<t;h++)l[h]=e.get(h,n);for(let h=0;h<s;h++){let u=0;for(let w=0;w<t;w++)u+=this.get(h,w)*l[w];r.set(h,n,u)}}return r}gram(){const e=this.rows,s=this.columns,t=new Float64Array(s*s);for(let r=0;r<e;r++)for(let l=0;l<s;l++){const n=this.get(r,l);if(n===0)continue;const h=l*s;for(let u=l;u<s;u++)t[h+u]+=n*this.get(r,u)}const o=new S(s,s);for(let r=0;r<s;r++){const l=r*s;for(let n=r;n<s;n++){const h=t[l+n];o.set(r,n,h),o.set(n,r,h)}}return o}mmulByTranspose(e){let s=this.rows,t=this.columns;if(e!==void 0&&e.length!==t)throw new RangeError("scale must have one value per column");let o=new S(s,s),r=new Float64Array(t);for(let l=0;l<s;l++){if(e===void 0)for(let n=0;n<t;n++)r[n]=this.get(l,n);else for(let n=0;n<t;n++)r[n]=e[n]*this.get(l,n);for(let n=l;n<s;n++){let h=0;for(let u=0;u<t;u++)h+=this.get(n,u)*r[u];o.set(n,l,h),o.set(l,n,h)}}return o}mpow(e){if(!this.isSquare())throw new RangeError("Matrix must be square");if(!Number.isInteger(e)||e<0)throw new RangeError("Exponent must be a non-negative integer");let s=S.eye(this.rows),t=this;for(let o=e;o>=1;o/=2)(o&1)!==0&&(s=s.mmul(t)),t=t.mmul(t);return s}strassen2x2(e){e=S.checkMatrix(e);let s=new S(2,2);const t=this.get(0,0),o=e.get(0,0),r=this.get(0,1),l=e.get(0,1),n=this.get(1,0),h=e.get(1,0),u=this.get(1,1),w=e.get(1,1),f=(t+u)*(o+w),a=(n+u)*o,j=t*(l-w),m=u*(h-o),p=(t+r)*w,E=(n-t)*(o+l),c=(r-u)*(h+w),M=f+m-p+c,k=j+p,I=a+m,N=f-a+j+E;return s.set(0,0,M),s.set(0,1,k),s.set(1,0,I),s.set(1,1,N),s}strassen3x3(e){e=S.checkMatrix(e);let s=new S(3,3);const t=this.get(0,0),o=this.get(0,1),r=this.get(0,2),l=this.get(1,0),n=this.get(1,1),h=this.get(1,2),u=this.get(2,0),w=this.get(2,1),f=this.get(2,2),a=e.get(0,0),j=e.get(0,1),m=e.get(0,2),p=e.get(1,0),E=e.get(1,1),c=e.get(1,2),M=e.get(2,0),k=e.get(2,1),I=e.get(2,2),N=(t+o+r-l-n-w-f)*E,V=(t-l)*(-j+E),v=n*(-a+j+p-E-c-M+I),R=(-t+l+n)*(a-j+E),F=(l+n)*(-a+j),g=t*a,d=(-t+u+w)*(a-m+c),b=(-t+u)*(m-c),y=(u+w)*(-a+m),C=(t+o+r-n-h-u-w)*c,D=w*(-a+m+p-E-c-M+k),z=(-r+w+f)*(E+M-k),P=(r-f)*(E-k),U=r*M,K=(w+f)*(-M+k),$=(-r+n+h)*(c+M-I),X=(r-h)*(c-I),A=(n+h)*(-M+I),_=o*p,O=h*k,W=l*m,Y=u*j,L=f*I,je=g+U+_,Se=N+R+F+g+z+U+K,Ee=g+d+y+C+U+$+A,be=V+v+R+g+U+$+X,ke=V+R+F+g+O,ve=U+$+X+A+W,Re=g+d+b+D+z+P+U,Ie=z+P+U+K+Y,Ne=g+d+b+y+L;return s.set(0,0,je),s.set(0,1,Se),s.set(0,2,Ee),s.set(1,0,be),s.set(1,1,ke),s.set(1,2,ve),s.set(2,0,Re),s.set(2,1,Ie),s.set(2,2,Ne),s}mmulStrassen(e){e=S.checkMatrix(e);let s=this.clone(),t=s.rows,o=s.columns,r=e.rows,l=e.columns;o!==r&&console.warn(`Multiplying ${t} x ${o} and ${r} x ${l} matrix: dimensions do not match.`);function n(f,a,j){let m=f.rows,p=f.columns;if(m===a&&p===j)return f;{let E=T.zeros(a,j);return E=E.setSubMatrix(f,0,0),E}}let h=Math.max(t,r),u=Math.max(o,l);s=n(s,h,u),e=n(e,h,u);function w(f,a,j,m){if(j<=512||m<=512)return f.mmul(a);j%2===1&&m%2===1?(f=n(f,j+1,m+1),a=n(a,j+1,m+1)):j%2===1?(f=n(f,j+1,m),a=n(a,j+1,m)):m%2===1&&(f=n(f,j,m+1),a=n(a,j,m+1));let p=parseInt(f.rows/2,10),E=parseInt(f.columns/2,10),c=f.subMatrix(0,p-1,0,E-1),M=a.subMatrix(0,p-1,0,E-1),k=f.subMatrix(0,p-1,E,f.columns-1),I=a.subMatrix(0,p-1,E,a.columns-1),N=f.subMatrix(p,f.rows-1,0,E-1),V=a.subMatrix(p,a.rows-1,0,E-1),v=f.subMatrix(p,f.rows-1,E,f.columns-1),R=a.subMatrix(p,a.rows-1,E,a.columns-1),F=w(T.add(c,v),T.add(M,R),p,E),g=w(T.add(N,v),M,p,E),d=w(c,T.sub(I,R),p,E),b=w(v,T.sub(V,M),p,E),y=w(T.add(c,k),R,p,E),C=w(T.sub(N,c),T.add(M,I),p,E),D=w(T.sub(k,v),T.add(V,R),p,E),z=T.add(F,b);z.sub(y),z.add(D);let P=T.add(d,y),U=T.add(g,b),K=T.sub(F,g);K.add(d),K.add(C);let $=T.zeros(2*z.rows,2*z.columns);return $=$.setSubMatrix(z,0,0),$=$.setSubMatrix(P,z.rows,0),$=$.setSubMatrix(U,0,z.columns),$=$.setSubMatrix(K,z.rows,z.columns),$.subMatrix(0,j-1,0,m-1)}return w(s,e,h,u)}scaleRows(e={}){if(typeof e!="object")throw new TypeError("options must be an object");const{min:s=0,max:t=1}=e;if(!Number.isFinite(s))throw new TypeError("min must be a number");if(!Number.isFinite(t))throw new TypeError("max must be a number");if(s>=t)throw new RangeError("min must be smaller than max");let o=new S(this.rows,this.columns);for(let r=0;r<this.rows;r++){const l=this.getRow(r);l.length>0&&dt(l,{min:s,max:t,output:l}),o.setRow(r,l)}return o}scaleColumns(e={}){if(typeof e!="object")throw new TypeError("options must be an object");const{min:s=0,max:t=1}=e;if(!Number.isFinite(s))throw new TypeError("min must be a number");if(!Number.isFinite(t))throw new TypeError("max must be a number");if(s>=t)throw new RangeError("min must be smaller than max");let o=new S(this.rows,this.columns);for(let r=0;r<this.columns;r++){const l=this.getColumn(r);l.length&&dt(l,{min:s,max:t,output:l}),o.setColumn(r,l)}return o}flipRows(){const e=Math.ceil(this.columns/2);for(let s=0;s<this.rows;s++)for(let t=0;t<e;t++){let o=this.get(s,t),r=this.get(s,this.columns-1-t);this.set(s,t,r),this.set(s,this.columns-1-t,o)}return this}flipColumns(){const e=Math.ceil(this.rows/2);for(let s=0;s<this.columns;s++)for(let t=0;t<e;t++){let o=this.get(t,s),r=this.get(this.rows-1-t,s);this.set(t,s,r),this.set(this.rows-1-t,s,o)}return this}kroneckerProduct(e){e=S.checkMatrix(e);let s=this.rows,t=this.columns,o=e.rows,r=e.columns,l=new S(s*o,t*r);for(let n=0;n<s;n++)for(let h=0;h<t;h++)for(let u=0;u<o;u++)for(let w=0;w<r;w++)l.set(o*n+u,r*h+w,this.get(n,h)*e.get(u,w));return l}kroneckerSum(e){if(e=S.checkMatrix(e),!this.isSquare()||!e.isSquare())throw new Error("Kronecker Sum needs two Square Matrices");let s=this.rows,t=e.rows,o=this.kroneckerProduct(S.eye(t,t)),r=S.eye(s,s).kroneckerProduct(e);return o.add(r)}transpose(){let e=new S(this.columns,this.rows);for(let s=0;s<this.rows;s++)for(let t=0;t<this.columns;t++)e.set(t,s,this.get(s,t));return e}sortRows(e=St){for(let s=0;s<this.rows;s++)this.setRow(s,this.getRow(s).sort(e));return this}sortColumns(e=St){for(let s=0;s<this.columns;s++)this.setColumn(s,this.getColumn(s).sort(e));return this}subMatrix(e,s,t,o){mt(this,e,s,t,o);let r=new S(s-e+1,o-t+1);for(let l=e;l<=s;l++)for(let n=t;n<=o;n++)r.set(l-e,n-t,this.get(l,n));return r}subMatrixRow(e,s,t){if(s===void 0&&(s=0),t===void 0&&(t=this.columns-1),s>t||s<0||s>=this.columns||t<0||t>=this.columns)throw new RangeError("Argument out of range");let o=new S(e.length,t-s+1);for(let r=0;r<e.length;r++)for(let l=s;l<=t;l++){if(e[r]<0||e[r]>=this.rows)throw new RangeError(`Row index out of range: ${e[r]}`);o.set(r,l-s,this.get(e[r],l))}return o}subMatrixColumn(e,s,t){if(s===void 0&&(s=0),t===void 0&&(t=this.rows-1),s>t||s<0||s>=this.rows||t<0||t>=this.rows)throw new RangeError("Argument out of range");let o=new S(t-s+1,e.length);for(let r=0;r<e.length;r++)for(let l=s;l<=t;l++){if(e[r]<0||e[r]>=this.columns)throw new RangeError(`Column index out of range: ${e[r]}`);o.set(l-s,r,this.get(l,e[r]))}return o}setSubMatrix(e,s,t){if(e=S.checkMatrix(e),e.isEmpty())return this;let o=s+e.rows-1,r=t+e.columns-1;mt(this,s,o,t,r);for(let l=0;l<e.rows;l++)for(let n=0;n<e.columns;n++)this.set(s+l,t+n,e.get(l,n));return this}selection(e,s){gt(this,e),wt(this,s);let t=new S(e.length,s.length);for(let o=0;o<e.length;o++){let r=e[o];for(let l=0;l<s.length;l++){let n=s[l];t.set(o,l,this.get(r,n))}}return t}trace(){let e=Math.min(this.rows,this.columns),s=0;for(let t=0;t<e;t++)s+=this.get(t,t);return s}clone(){return this.constructor.copy(this,new S(this.rows,this.columns))}static copy(e,s){for(const[t,o,r]of e.entries())s.set(t,o,r);return s}sum(e){switch(e){case"row":return Vt(this);case"column":return Ct(this);case void 0:return Pt(this);default:throw new Error(`invalid option: ${e}`)}}product(e){switch(e){case"row":return $t(this);case"column":return Lt(this);case void 0:return Bt(this);default:throw new Error(`invalid option: ${e}`)}}mean(e){const s=this.sum(e);switch(e){case"row":{for(let t=0;t<this.rows;t++)s[t]/=this.columns;return s}case"column":{for(let t=0;t<this.columns;t++)s[t]/=this.rows;return s}case void 0:return s/this.size;default:throw new Error(`invalid option: ${e}`)}}variance(e,s={}){if(typeof e=="object"&&(s=e,e=void 0),typeof s!="object")throw new TypeError("options must be an object");const{unbiased:t=!0,mean:o=this.mean(e)}=s;if(typeof t!="boolean")throw new TypeError("unbiased must be a boolean");switch(e){case"row":{if(!B(o))throw new TypeError("mean must be an array");return Ut(this,t,o)}case"column":{if(!B(o))throw new TypeError("mean must be an array");return Ot(this,t,o)}case void 0:{if(typeof o!="number")throw new TypeError("mean must be a number");return Qt(this,t,o)}default:throw new Error(`invalid option: ${e}`)}}standardDeviation(e,s){typeof e=="object"&&(s=e,e=void 0);const t=this.variance(e,s);if(e===void 0)return Math.sqrt(t);for(let o=0;o<t.length;o++)t[o]=Math.sqrt(t[o]);return t}center(e,s={}){if(typeof e=="object"&&(s=e,e=void 0),typeof s!="object")throw new TypeError("options must be an object");const{center:t=this.mean(e)}=s;switch(e){case"row":{if(!B(t))throw new TypeError("center must be an array");return Wt(this,t),this}case"column":{if(!B(t))throw new TypeError("center must be an array");return Yt(this,t),this}case void 0:{if(typeof t!="number")throw new TypeError("center must be a number");return Jt(this,t),this}default:throw new Error(`invalid option: ${e}`)}}scale(e,s={}){if(typeof e=="object"&&(s=e,e=void 0),typeof s!="object")throw new TypeError("options must be an object");let t=s.scale;switch(e){case"row":{if(t===void 0)t=Gt(this);else if(!B(t))throw new TypeError("scale must be an array");return Kt(this,t),this}case"column":{if(t===void 0)t=Xt(this);else if(!B(t))throw new TypeError("scale must be an array");return Zt(this,t),this}case void 0:{if(t===void 0)t=xt(this);else if(typeof t!="number")throw new TypeError("scale must be a number");return At(this,t),this}default:throw new Error(`invalid option: ${e}`)}}toString(e){return Mt(this,e)}[Symbol.iterator](){return this.entries()}*entries(){for(let e=0;e<this.rows;e++)for(let s=0;s<this.columns;s++)yield[e,s,this.get(e,s)]}*values(){for(let e=0;e<this.rows;e++)for(let s=0;s<this.columns;s++)yield this.get(e,s)}}T.prototype.klass="Matrix",typeof Symbol<"u"&&(T.prototype[Symbol.for("nodejs.util.inspect.custom")]=_t);function St(i,e){return i-e}function Ht(i){return i.every(e=>typeof e=="number")}T.random=T.rand,T.randomInt=T.randInt,T.diagonal=T.diag,T.prototype.diagonal=T.prototype.diag,T.identity=T.eye,T.prototype.negate=T.prototype.neg,T.prototype.tensorProduct=T.prototype.kroneckerProduct;class S extends T{data;#t(e,s){if(this.data=[],Number.isInteger(s)&&s>=0)for(let t=0;t<e;t++)this.data.push(new Float64Array(s));else throw new TypeError("nColumns must be a positive integer");this.rows=e,this.columns=s}constructor(e,s){if(super(),S.isMatrix(e))this.#t(e.rows,e.columns),S.copy(e,this);else if(Number.isInteger(e)&&e>=0)this.#t(e,s);else if(B(e)){const t=e;if(e=t.length,s=e?t[0].length:0,typeof s!="number")throw new TypeError("Data must be a 2D array with at least one element");this.data=[];for(let o=0;o<e;o++){if(t[o].length!==s)throw new RangeError("Inconsistent array dimensions");if(!Ht(t[o]))throw new TypeError("Input data contains non-numeric values");this.data.push(Float64Array.from(t[o]))}this.rows=e,this.columns=s}else throw new TypeError("First argument must be a positive number or an array")}set(e,s,t){return this.data[e][s]=t,this}get(e,s){return this.data[e][s]}removeRow(e){return J(this,e),this.data.splice(e,1),this.rows-=1,this}addRow(e,s){return s===void 0&&(s=e,e=this.rows),J(this,e,!0),s=Float64Array.from(tt(this,s)),this.data.splice(e,0,s),this.rows+=1,this}removeColumn(e){G(this,e);for(let s=0;s<this.rows;s++){const t=new Float64Array(this.columns-1);for(let o=0;o<e;o++)t[o]=this.data[s][o];for(let o=e+1;o<this.columns;o++)t[o-1]=this.data[s][o];this.data[s]=t}return this.columns-=1,this}addColumn(e,s){typeof s>"u"&&(s=e,e=this.columns),G(this,e,!0),s=et(this,s);for(let t=0;t<this.rows;t++){const o=new Float64Array(this.columns+1);let r=0;for(;r<e;r++)o[r]=this.data[t][r];for(o[r++]=s[t];r<this.columns+1;r++)o[r]=this.data[t][r-1];this.data[t]=o}return this.columns+=1,this}}zt(T,S);class H extends T{#t;get size(){return this.#t.size}get rows(){return this.#t.rows}get columns(){return this.#t.columns}get diagonalSize(){return this.rows}static isSymmetricMatrix(e){return S.isMatrix(e)&&e.klassType==="SymmetricMatrix"}static zeros(e){return new this(e)}static ones(e){return new this(e).fill(1)}constructor(e){if(super(),S.isMatrix(e)){if(!e.isSymmetric())throw new TypeError("not symmetric data");this.#t=S.copy(e,new S(e.rows,e.rows))}else if(Number.isInteger(e)&&e>=0)this.#t=new S(e,e);else if(this.#t=new S(e),!this.isSymmetric())throw new TypeError("not symmetric data")}clone(){const e=new H(this.diagonalSize);for(const[s,t,o]of this.upperRightEntries())e.set(s,t,o);return e}toMatrix(){return new S(this)}get(e,s){return this.#t.get(e,s)}set(e,s,t){return this.#t.set(e,s,t),this.#t.set(s,e,t),this}removeCross(e){return this.#t.removeRow(e),this.#t.removeColumn(e),this}addCross(e,s){s===void 0&&(s=e,e=this.diagonalSize);const t=s.slice();return t.splice(e,1),this.#t.addRow(e,t),this.#t.addColumn(e,s),this}applyMask(e){if(e.length!==this.diagonalSize)throw new RangeError("Mask size do not match with matrix size");const s=[];for(const[t,o]of e.entries())o||s.push(t);s.reverse();for(const t of s)this.removeCross(t);return this}toCompact(){const{diagonalSize:e}=this,s=new Array(e*(e+1)/2);for(let t=0,o=0,r=0;r<s.length;r++)s[r]=this.get(o,t),++t>=e&&(t=++o);return s}static fromCompact(e){const s=e.length,t=(Math.sqrt(8*s+1)-1)/2;if(!Number.isInteger(t))throw new TypeError(`This array is not a compact representation of a Symmetric Matrix, ${JSON.stringify(e)}`);const o=new H(t);for(let r=0,l=0,n=0;n<s;n++)o.set(r,l,e[n]),++r>=t&&(r=++l);return o}*upperRightEntries(){for(let e=0,s=0;e<this.diagonalSize;void 0){const t=this.get(e,s);yield[e,s,t],++s>=this.diagonalSize&&(s=++e)}}*upperRightValues(){for(let e=0,s=0;e<this.diagonalSize;void 0)yield this.get(e,s),++s>=this.diagonalSize&&(s=++e)}}H.prototype.klassType="SymmetricMatrix";class lt extends H{static isDistanceMatrix(e){return H.isSymmetricMatrix(e)&&e.klassSubType==="DistanceMatrix"}constructor(e){if(super(e),!this.isDistance())throw new TypeError("Provided arguments do no produce a distance matrix")}set(e,s,t){return e===s&&(t=0),super.set(e,s,t)}addCross(e,s){return s===void 0&&(s=e,e=this.diagonalSize),s=s.slice(),s[e]=0,super.addCross(e,s)}toSymmetricMatrix(){return new H(this)}clone(){const e=new lt(this.diagonalSize);for(const[s,t,o]of this.upperRightEntries())s!==t&&e.set(s,t,o);return e}toCompact(){const{diagonalSize:e}=this,s=(e-1)*e/2,t=new Array(s);for(let o=1,r=0,l=0;l<t.length;l++)t[l]=this.get(r,o),++o>=e&&(o=++r+1);return t}static fromCompact(e){const s=e.length;if(s===0)return new this(0);const t=(Math.sqrt(8*s+1)+1)/2;if(!Number.isInteger(t))throw new TypeError(`This array is not a compact representation of a DistanceMatrix, ${JSON.stringify(e)}`);const o=new this(t);for(let r=1,l=0,n=0;n<s;n++)o.set(r,l,e[n]),++r>=t&&(r=++l+1);return o}}lt.prototype.klassSubType="DistanceMatrix";class Z extends T{constructor(e,s,t){super(),this.matrix=e,this.rows=s,this.columns=t}}class te extends Z{constructor(e,s){G(e,s),super(e,e.rows,1),this.column=s}set(e,s,t){return this.matrix.set(e,this.column,t),this}get(e){return this.matrix.get(e,this.column)}}class ee extends Z{constructor(e,s){wt(e,s),super(e,e.rows,s.length),this.columnIndices=s}set(e,s,t){return this.matrix.set(e,this.columnIndices[s],t),this}get(e,s){return this.matrix.get(e,this.columnIndices[s])}}class se extends Z{constructor(e){super(e,e.rows,e.columns)}set(e,s,t){return this.matrix.set(e,this.columns-s-1,t),this}get(e,s){return this.matrix.get(e,this.columns-s-1)}}class oe extends Z{constructor(e){super(e,e.rows,e.columns)}set(e,s,t){return this.matrix.set(this.rows-e-1,s,t),this}get(e,s){return this.matrix.get(this.rows-e-1,s)}}class re extends Z{constructor(e,s){J(e,s),super(e,1,e.columns),this.row=s}set(e,s,t){return this.matrix.set(this.row,s,t),this}get(e,s){return this.matrix.get(this.row,s)}}class ie extends Z{constructor(e,s){gt(e,s),super(e,s.length,e.columns),this.rowIndices=s}set(e,s,t){return this.matrix.set(this.rowIndices[e],s,t),this}get(e,s){return this.matrix.get(this.rowIndices[e],s)}}class ht extends Z{constructor(e,s,t){gt(e,s),wt(e,t),super(e,s.length,t.length),this.rowIndices=s,this.columnIndices=t}set(e,s,t){return this.matrix.set(this.rowIndices[e],this.columnIndices[s],t),this}get(e,s){return this.matrix.get(this.rowIndices[e],this.columnIndices[s])}}class ne extends Z{constructor(e,s,t,o,r){mt(e,s,t,o,r),super(e,t-s+1,r-o+1),this.startRow=s,this.startColumn=o}set(e,s,t){return this.matrix.set(this.startRow+e,this.startColumn+s,t),this}get(e,s){return this.matrix.get(this.startRow+e,this.startColumn+s)}}class le extends Z{constructor(e){super(e,e.columns,e.rows)}set(e,s,t){return this.matrix.set(s,e,t),this}get(e,s){return this.matrix.get(s,e)}}class Et extends T{constructor(e,s={}){const{rows:t=1}=s;if(e.length%t!==0)throw new Error("the data length is not divisible by the number of rows");super(),this.rows=t,this.columns=e.length/t,this.data=e}set(e,s,t){let o=this._calculateIndex(e,s);return this.data[o]=t,this}get(e,s){let t=this._calculateIndex(e,s);return this.data[t]}_calculateIndex(e,s){return e*this.columns+s}}class Q extends T{constructor(e){super(),this.data=e,this.rows=e.length,this.columns=e[0].length}set(e,s,t){return this.data[e][s]=t,this}get(e,s){return this.data[e][s]}}function he(i,e){if(B(i))return i[0]&&B(i[0])?new Q(i):new Et(i,e);throw new Error("the argument is not an array")}class ut{constructor(e){e=Q.checkMatrix(e);let s=e.clone(),t=s.rows,o=s.columns,r=new Float64Array(t),l=1,n,h,u,w,f,a,j,m,p;for(n=0;n<t;n++)r[n]=n;for(m=new Float64Array(t),h=0;h<o;h++){for(n=0;n<t;n++)m[n]=s.get(n,h);for(n=0;n<t;n++){for(p=Math.min(n,h),f=0,u=0;u<p;u++)f+=s.get(n,u)*m[u];m[n]-=f,s.set(n,h,m[n])}for(w=h,n=h+1;n<t;n++)Math.abs(m[n])>Math.abs(m[w])&&(w=n);if(w!==h){for(u=0;u<o;u++)a=s.get(w,u),s.set(w,u,s.get(h,u)),s.set(h,u,a);j=r[w],r[w]=r[h],r[h]=j,l=-l}if(h<t&&s.get(h,h)!==0)for(n=h+1;n<t;n++)s.set(n,h,s.get(n,h)/s.get(h,h))}this.LU=s,this.pivotVector=r,this.pivotSign=l}isSingular(){let e=this.LU,s=e.columns;for(let t=0;t<s;t++)if(e.get(t,t)===0)return!0;return!1}solve(e){e=S.checkMatrix(e);let s=this.LU;if(s.rows!==e.rows)throw new Error("Invalid matrix dimensions");if(this.isSingular())throw new Error("LU matrix is singular");let o=e.columns,r=e.subMatrixRow(this.pivotVector,0,o-1),l=s.columns,n,h,u;for(u=0;u<l;u++)for(n=u+1;n<l;n++)for(h=0;h<o;h++)r.set(n,h,r.get(n,h)-r.get(u,h)*s.get(n,u));for(u=l-1;u>=0;u--){for(h=0;h<o;h++)r.set(u,h,r.get(u,h)/s.get(u,u));for(n=0;n<u;n++)for(h=0;h<o;h++)r.set(n,h,r.get(n,h)-r.get(u,h)*s.get(n,u))}return r}get determinant(){let e=this.LU;if(!e.isSquare())throw new Error("Matrix must be square");let s=this.pivotSign,t=e.columns;for(let o=0;o<t;o++)s*=e.get(o,o);return s}get lowerTriangularMatrix(){let e=this.LU,s=e.rows,t=e.columns,o=new S(s,t);for(let r=0;r<s;r++)for(let l=0;l<t;l++)r>l?o.set(r,l,e.get(r,l)):r===l?o.set(r,l,1):o.set(r,l,0);return o}get upperTriangularMatrix(){let e=this.LU,s=e.rows,t=e.columns,o=new S(s,t);for(let r=0;r<s;r++)for(let l=0;l<t;l++)r<=l?o.set(r,l,e.get(r,l)):o.set(r,l,0);return o}get pivotPermutationVector(){return Array.from(this.pivotVector)}}function at(i){const e=i.data,s=i.rows;for(let t=0;t<s;t++){const o=e[t];for(let r=t+1;r<s;r++){const l=o[r];o[r]=e[r][t],e[r][t]=l}}return i}function x(i,e){let s=0;return Math.abs(i)>Math.abs(e)?(s=e/i,Math.abs(i)*Math.sqrt(1+s*s)):e!==0?(s=i/e,Math.abs(e)*Math.sqrt(1+s*s)):0}class pt{constructor(e){e=Q.checkMatrix(e);let s=e.clone(),t=e.rows,o=e.columns,r=new Float64Array(o),l,n,h,u;for(h=0;h<o;h++){let w=0;for(l=h;l<t;l++)w=x(w,s.get(l,h));if(w!==0){for(s.get(h,h)<0&&(w=-w),l=h;l<t;l++)s.set(l,h,s.get(l,h)/w);for(s.set(h,h,s.get(h,h)+1),n=h+1;n<o;n++){for(u=0,l=h;l<t;l++)u+=s.get(l,h)*s.get(l,n);for(u=-u/s.get(h,h),l=h;l<t;l++)s.set(l,n,s.get(l,n)+u*s.get(l,h))}}r[h]=-w}this.QR=s,this.Rdiag=r}solve(e){e=S.checkMatrix(e);let s=this.QR,t=s.rows;if(e.rows!==t)throw new Error("Matrix row dimensions must agree");if(!this.isFullRank())throw new Error("Matrix is rank deficient");let o=e.columns,r=e.clone(),l=s.columns,n,h,u,w;for(u=0;u<l;u++)for(h=0;h<o;h++){for(w=0,n=u;n<t;n++)w+=s.get(n,u)*r.get(n,h);for(w=-w/s.get(u,u),n=u;n<t;n++)r.set(n,h,r.get(n,h)+w*s.get(n,u))}for(u=l-1;u>=0;u--){for(h=0;h<o;h++)r.set(u,h,r.get(u,h)/this.Rdiag[u]);for(n=0;n<u;n++)for(h=0;h<o;h++)r.set(n,h,r.get(n,h)-r.get(u,h)*s.get(n,u))}return r.subMatrix(0,l-1,0,o-1)}isFullRank(){let e=this.QR.columns;for(let s=0;s<e;s++)if(this.Rdiag[s]===0)return!1;return!0}get upperTriangularMatrix(){let e=this.QR,s=e.columns,t=new S(s,s),o,r;for(o=0;o<s;o++)for(r=0;r<s;r++)o<r?t.set(o,r,e.get(o,r)):o===r?t.set(o,r,this.Rdiag[o]):t.set(o,r,0);return t}get orthogonalMatrix(){let e=this.QR,s=e.rows,t=e.columns,o=new S(s,t),r,l,n,h;for(n=t-1;n>=0;n--){for(r=0;r<s;r++)o.set(r,n,0);for(o.set(n,n,1),l=n;l<t;l++)if(e.get(n,n)!==0){for(h=0,r=n;r<s;r++)h+=e.get(r,n)*o.get(r,l);for(h=-h/e.get(n,n),r=n;r<s;r++)o.set(r,l,o.get(r,l)+h*e.get(r,n))}}return o}}class ot{constructor(e,s={}){if(e=Q.checkMatrix(e),e.isEmpty())throw new Error("Matrix must be non-empty");let t=e.rows,o=e.columns;const{computeLeftSingularVectors:r=!0,computeRightSingularVectors:l=!0,autoTranspose:n=!1}=s;let h=!!r,u=!!l,w=!1,f;if(t<o)if(!n)console.warn("Computing SVD on a matrix with more columns than rows. Consider enabling autoTranspose"),f=e.transpose();else{f=e.clone(),t=e.columns,o=e.rows,w=!0;let g=h;h=u,u=g}else f=e.transpose();let a=Math.min(t,o),j=Math.min(t+1,o),m=new Float64Array(j),p=new S(a,t),E=new S(o,o),c=new Float64Array(o),M=new Float64Array(t),k=new Float64Array(j);for(let g=0;g<j;g++)k[g]=g;let I=Math.min(t-1,o),N=Math.max(0,Math.min(o-2,t)),V=Math.max(I,N);for(let g=0;g<V;g++){if(g<I){m[g]=0;for(let d=g;d<t;d++)m[g]=x(m[g],f.get(g,d));if(m[g]!==0){f.get(g,g)<0&&(m[g]=-m[g]);for(let d=g;d<t;d++)f.set(g,d,f.get(g,d)/m[g]);f.set(g,g,f.get(g,g)+1)}m[g]=-m[g]}for(let d=g+1;d<o;d++){if(g<I&&m[g]!==0){let b=0;for(let y=g;y<t;y++)b+=f.get(g,y)*f.get(d,y);b=-b/f.get(g,g);for(let y=g;y<t;y++)f.set(d,y,f.get(d,y)+b*f.get(g,y))}c[d]=f.get(d,g)}if(h&&g<I)for(let d=g;d<t;d++)p.set(g,d,f.get(g,d));if(g<N){c[g]=0;for(let d=g+1;d<o;d++)c[g]=x(c[g],c[d]);if(c[g]!==0){c[g+1]<0&&(c[g]=0-c[g]);for(let d=g+1;d<o;d++)c[d]/=c[g];c[g+1]+=1}if(c[g]=-c[g],g+1<t&&c[g]!==0){for(let d=g+1;d<t;d++)M[d]=0;for(let d=g+1;d<t;d++)for(let b=g+1;b<o;b++)M[d]+=c[b]*f.get(b,d);for(let d=g+1;d<o;d++){let b=-c[d]/c[g+1];for(let y=g+1;y<t;y++)f.set(d,y,f.get(d,y)+b*M[y])}}if(u)for(let d=g+1;d<o;d++)E.set(g,d,c[d])}}let v=Math.min(o,t+1);if(I<o&&(m[I]=f.get(I,I)),t<v&&(m[v-1]=0),N+1<v&&(c[N]=f.get(v-1,N)),c[v-1]=0,h){for(let g=I;g<a;g++){for(let d=0;d<t;d++)p.set(g,d,0);p.set(g,g,1)}for(let g=I-1;g>=0;g--)if(m[g]!==0){for(let d=g+1;d<a;d++){let b=0;for(let y=g;y<t;y++)b+=p.get(g,y)*p.get(d,y);b=-b/p.get(g,g);for(let y=g;y<t;y++)p.set(d,y,p.get(d,y)+b*p.get(g,y))}for(let d=g;d<t;d++)p.set(g,d,-p.get(g,d));p.set(g,g,1+p.get(g,g));for(let d=0;d<g-1;d++)p.set(g,d,0)}else{for(let d=0;d<t;d++)p.set(g,d,0);p.set(g,g,1)}}if(u)for(let g=o-1;g>=0;g--){if(g<N&&c[g]!==0)for(let d=g+1;d<o;d++){let b=0;for(let y=g+1;y<o;y++)b+=E.get(g,y)*E.get(d,y);b=-b/E.get(g,g+1);for(let y=g+1;y<o;y++)E.set(d,y,E.get(d,y)+b*E.get(g,y))}for(let d=0;d<o;d++)E.set(g,d,0);E.set(g,g,1)}let R=v-1,F=Number.EPSILON;for(;v>0;){let g,d;for(g=v-2;g>=-1&&g!==-1;g--){const b=Number.MIN_VALUE+F*Math.abs(m[g]+Math.abs(m[g+1]));if(Math.abs(c[g])<=b||Number.isNaN(c[g])){c[g]=0;break}}if(g===v-2)d=4;else{let b;for(b=v-1;b>=g&&b!==g;b--){let y=(b!==v?Math.abs(c[b]):0)+(b!==g+1?Math.abs(c[b-1]):0);if(Math.abs(m[b])<=F*y){m[b]=0;break}}b===g?d=3:b===v-1?d=1:(d=2,g=b)}switch(g++,d){case 1:{let b=c[v-2];c[v-2]=0;for(let y=v-2;y>=g;y--){let C=x(m[y],b),D=m[y]/C,z=b/C;if(m[y]=C,y!==g&&(b=-z*c[y-1],c[y-1]=D*c[y-1]),u)for(let P=0;P<o;P++)C=D*E.get(y,P)+z*E.get(v-1,P),E.set(v-1,P,-z*E.get(y,P)+D*E.get(v-1,P)),E.set(y,P,C)}break}case 2:{let b=c[g-1];c[g-1]=0;for(let y=g;y<v;y++){let C=x(m[y],b),D=m[y]/C,z=b/C;if(m[y]=C,b=-z*c[y],c[y]=D*c[y],h)for(let P=0;P<t;P++)C=D*p.get(y,P)+z*p.get(g-1,P),p.set(g-1,P,-z*p.get(y,P)+D*p.get(g-1,P)),p.set(y,P,C)}break}case 3:{const b=Math.max(Math.abs(m[v-1]),Math.abs(m[v-2]),Math.abs(c[v-2]),Math.abs(m[g]),Math.abs(c[g])),y=m[v-1]/b,C=m[v-2]/b,D=c[v-2]/b,z=m[g]/b,P=c[g]/b,U=((C+y)*(C-y)+D*D)/2,K=y*D*(y*D);let $=0;(U!==0||K!==0)&&(U<0?$=0-Math.sqrt(U*U+K):$=Math.sqrt(U*U+K),$=K/(U+$));let X=(z+y)*(z-y)+$,A=z*P;for(let _=g;_<v-1;_++){let O=x(X,A);O===0&&(O=Number.MIN_VALUE);let W=X/O,Y=A/O;if(_!==g&&(c[_-1]=O),X=W*m[_]+Y*c[_],c[_]=W*c[_]-Y*m[_],A=Y*m[_+1],m[_+1]=W*m[_+1],u)for(let L=0;L<o;L++)O=W*E.get(_,L)+Y*E.get(_+1,L),E.set(_+1,L,-Y*E.get(_,L)+W*E.get(_+1,L)),E.set(_,L,O);if(O=x(X,A),O===0&&(O=Number.MIN_VALUE),W=X/O,Y=A/O,m[_]=O,X=W*c[_]+Y*m[_+1],m[_+1]=-Y*c[_]+W*m[_+1],A=Y*c[_+1],c[_+1]=W*c[_+1],h&&_<t-1)for(let L=0;L<t;L++)O=W*p.get(_,L)+Y*p.get(_+1,L),p.set(_+1,L,-Y*p.get(_,L)+W*p.get(_+1,L)),p.set(_,L,O)}c[v-2]=X;break}case 4:{if(m[g]<=0&&(m[g]=m[g]<0?-m[g]:0,u))for(let b=0;b<=R;b++)E.set(g,b,-E.get(g,b));for(;g<R&&!(m[g]>=m[g+1]);){let b=m[g];if(m[g]=m[g+1],m[g+1]=b,u&&g<o-1)for(let y=0;y<o;y++)b=E.get(g+1,y),E.set(g+1,y,E.get(g,y)),E.set(g,y,b);if(h&&g<t-1)for(let y=0;y<t;y++)b=p.get(g+1,y),p.set(g+1,y,p.get(g,y)),p.set(g,y,b);g++}v--;break}}}if(p=p.isSquare()?at(p):p.transpose(),E=at(E),w){let g=E;E=p,p=g}this.m=t,this.n=o,this.s=m,this.U=p,this.V=E}solve(e){let s=e,t=this.threshold,o=this.s.length,r=S.zeros(o,o);for(let a=0;a<o;a++)Math.abs(this.s[a])<=t?r.set(a,a,0):r.set(a,a,1/this.s[a]);let l=this.U,n=this.rightSingularVectors,h=n.mmul(r),u=n.rows,w=l.rows,f=S.zeros(u,w);for(let a=0;a<u;a++)for(let j=0;j<w;j++){let m=0;for(let p=0;p<o;p++)m+=h.get(a,p)*l.get(j,p);f.set(a,j,m)}return f.mmul(s)}solveForDiagonal(e){return this.solve(S.diag(e))}inverse(){let e=this.V,s=this.threshold,t=e.rows,o=e.columns,r=new S(t,this.s.length);for(let w=0;w<t;w++)for(let f=0;f<o;f++)Math.abs(this.s[f])>s&&r.set(w,f,e.get(w,f)/this.s[f]);let l=this.U,n=l.rows,h=l.columns,u=new S(t,n);for(let w=0;w<t;w++)for(let f=0;f<n;f++){let a=0;for(let j=0;j<h;j++)a+=r.get(w,j)*l.get(f,j);u.set(w,f,a)}return u}get condition(){return this.s[0]/this.s[Math.min(this.m,this.n)-1]}get norm2(){return this.s[0]}get rank(){let e=Math.max(this.m,this.n)*this.s[0]*Number.EPSILON,s=0,t=this.s;for(let o=0,r=t.length;o<r;o++)t[o]>e&&s++;return s}get diagonal(){return Array.from(this.s)}get threshold(){return Number.EPSILON/2*Math.max(this.m,this.n)*this.s[0]}get leftSingularVectors(){return this.U}get rightSingularVectors(){return this.V}get diagonalMatrix(){return S.diag(this.s)}}function ue(i,e=!1){return i=Q.checkMatrix(i),e?new ot(i).inverse():bt(i,S.eye(i.rows))}function bt(i,e,s=!1){return i=Q.checkMatrix(i),e=Q.checkMatrix(e),s?new ot(i).solve(e):i.isSquare()?new ut(i).solve(e):new pt(i).solve(e)}function ft(i){if(i=S.checkMatrix(i),i.isSquare()){if(i.columns===0)return 1;let e,s,t,o;if(i.columns===2)return e=i.get(0,0),s=i.get(0,1),t=i.get(1,0),o=i.get(1,1),e*o-s*t;if(i.columns===3){let r,l,n;return r=new ht(i,[1,2],[1,2]),l=new ht(i,[1,2],[0,2]),n=new ht(i,[1,2],[0,1]),e=i.get(0,0),s=i.get(0,1),t=i.get(0,2),e*ft(r)-s*ft(l)+t*ft(n)}else return new ut(i).determinant}else throw Error("determinant can only be calculated for a square matrix")}function fe(i,e){let s=[];for(let t=0;t<i;t++)t!==e&&s.push(t);return s}function ce(i,e,s,t=1e-9,o=1e-9){if(i>o)return new Array(e.rows+1).fill(0);{let r=e.addRow(s,[0]);for(let l=0;l<r.rows;l++)Math.abs(r.get(l,0))<t&&r.set(l,0,0);return r.to1DArray()}}function ge(i,e={}){const{thresholdValue:s=1e-9,thresholdError:t=1e-9}=e;i=S.checkMatrix(i);let o=i.rows,r=new S(o,o);for(let l=0;l<o;l++){let n=S.columnVector(i.getRow(l)),h=i.subMatrixRow(fe(o,l)).transpose(),w=new ot(h).solve(n),f=S.sub(n,h.mmul(w)).abs().max();r.setRow(l,ce(f,w,l,s,t))}return r}function we(i,e=Number.EPSILON){if(i=S.checkMatrix(i),i.isEmpty())return i.transpose();let s=new ot(i,{autoTranspose:!0}),t=s.leftSingularVectors,o=s.rightSingularVectors,r=s.diagonal;for(let l=0;l<r.length;l++)Math.abs(r[l])>e?r[l]=1/r[l]:r[l]=0;return o.mmul(S.diag(r).mmul(t.transpose()))}function me(i,e=i,s={}){i=new S(i);let t=!1;if(typeof e=="object"&&!S.isMatrix(e)&&!B(e)?(s=e,e=i,t=!0):e=new S(e),i.rows!==e.rows)throw new TypeError("Both matrices must have the same number of rows");const{center:o=!0}=s;o&&(i=i.center("column"),t||(e=e.center("column")));const r=i.transpose().mmul(e);for(let l=0;l<r.rows;l++)for(let n=0;n<r.columns;n++)r.set(l,n,r.get(l,n)*(1/(i.rows-1)));return r}function ae(i,e=i,s={}){i=new S(i);let t=!1;if(typeof e=="object"&&!S.isMatrix(e)&&!B(e)?(s=e,e=i,t=!0):e=new S(e),i.rows!==e.rows)throw new TypeError("Both matrices must have the same number of rows");const{center:o=!0,scale:r=!0}=s;o&&(i.center("column"),t||e.center("column")),r&&(i.scale("column"),t||e.scale("column"));const l=i.standardDeviation("column",{unbiased:!0}),n=t?l:e.standardDeviation("column",{unbiased:!0}),h=i.transpose().mmul(e);for(let u=0;u<h.rows;u++)for(let w=0;w<h.columns;w++)h.set(u,w,h.get(u,w)*(1/(l[u]*n[w]))*(1/(i.rows-1)));return h}class kt{constructor(e,s={}){const{assumeSymmetric:t=!1}=s;if(e=Q.checkMatrix(e),!e.isSquare())throw new Error("Matrix is not a square matrix");if(e.isEmpty())throw new Error("Matrix must be non-empty");let o=e.columns,r=new S(o,o),l=new Float64Array(o),n=new Float64Array(o),h=e,u,w,f=!1;if(t?f=!0:f=e.isSymmetric(),f){for(u=0;u<o;u++)for(w=0;w<o;w++)r.set(w,u,h.get(u,w));pe(o,n,l,r),de(o,n,l,r),at(r)}else{let a=new S(o,o),j=new Float64Array(o);for(w=0;w<o;w++)for(u=0;u<o;u++)a.set(u,w,h.get(u,w));ye(o,a,j,r),Me(o,n,l,r,a)}this.n=o,this.e=n,this.d=l,this.V=r}get realEigenvalues(){return Array.from(this.d)}get imaginaryEigenvalues(){return Array.from(this.e)}get eigenvectorMatrix(){return this.V}get diagonalMatrix(){let e=this.n,s=this.e,t=this.d,o=new S(e,e),r,l;for(r=0;r<e;r++){for(l=0;l<e;l++)o.set(r,l,0);o.set(r,r,t[r]),s[r]>0?o.set(r,r+1,s[r]):s[r]<0&&o.set(r,r-1,s[r])}return o}}function pe(i,e,s,t){let o,r,l,n,h,u,w,f;for(h=0;h<i;h++)s[h]=t.get(h,i-1);for(n=i-1;n>0;n--){for(f=0,l=0,u=0;u<n;u++)f=f+Math.abs(s[u]);if(f===0)for(e[n]=s[n-1],h=0;h<n;h++)s[h]=t.get(h,n-1),t.set(h,n,0),t.set(n,h,0);else{for(u=0;u<n;u++)s[u]/=f,l+=s[u]*s[u];for(o=s[n-1],r=Math.sqrt(l),o>0&&(r=-r),e[n]=f*r,l=l-o*r,s[n-1]=o-r,h=0;h<n;h++)e[h]=0;for(h=0;h<n;h++){for(o=s[h],t.set(n,h,o),r=e[h]+t.get(h,h)*o,u=h+1;u<=n-1;u++)r+=t.get(h,u)*s[u],e[u]+=t.get(h,u)*o;e[h]=r}for(o=0,h=0;h<n;h++)e[h]/=l,o+=e[h]*s[h];for(w=o/(l+l),h=0;h<n;h++)e[h]-=w*s[h];for(h=0;h<n;h++){for(o=s[h],r=e[h],u=h;u<=n-1;u++)t.set(h,u,t.get(h,u)-(o*e[u]+r*s[u]));s[h]=t.get(h,n-1),t.set(h,n,0)}}s[n]=l}for(n=0;n<i-1;n++){if(t.set(n,i-1,t.get(n,n)),t.set(n,n,1),l=s[n+1],l!==0){for(u=0;u<=n;u++)s[u]=t.get(n+1,u)/l;for(h=0;h<=n;h++){for(r=0,u=0;u<=n;u++)r+=t.get(n+1,u)*t.get(h,u);for(u=0;u<=n;u++)t.set(h,u,t.get(h,u)-r*s[u])}}for(u=0;u<=n;u++)t.set(n+1,u,0)}for(h=0;h<i;h++)s[h]=t.get(h,i-1),t.set(h,i-1,0);t.set(i-1,i-1,1),e[0]=0}function de(i,e,s,t){let o,r,l,n,h,u,w,f,a,j,m,p,E,c,M,k;for(l=1;l<i;l++)e[l-1]=e[l];e[i-1]=0;let I=0,N=0,V=Number.EPSILON;for(u=0;u<i;u++){for(N=Math.max(N,Math.abs(s[u])+Math.abs(e[u])),w=u;w<i&&!(Math.abs(e[w])<=V*N);)w++;if(w>u)do{for(o=s[u],f=(s[u+1]-o)/(2*e[u]),a=x(f,1),f<0&&(a=-a),s[u]=e[u]/(f+a),s[u+1]=e[u]*(f+a),j=s[u+1],r=o-s[u],l=u+2;l<i;l++)s[l]-=r;for(I=I+r,f=s[w],m=1,p=m,E=m,c=e[u+1],M=0,k=0,l=w-1;l>=u;l--)for(E=p,p=m,k=M,o=m*e[l],r=m*f,a=x(f,e[l]),e[l+1]=M*a,M=e[l]/a,m=f/a,f=m*s[l]-M*o,s[l+1]=r+M*(m*o+M*s[l]),h=0;h<i;h++)r=t.get(l+1,h),t.set(l+1,h,M*t.get(l,h)+m*r),t.set(l,h,m*t.get(l,h)-M*r);f=-M*k*E*c*e[u]/j,e[u]=M*f,s[u]=m*f}while(Math.abs(e[u])>V*N);s[u]=s[u]+I,e[u]=0}for(l=0;l<i-1;l++){for(h=l,f=s[l],n=l+1;n<i;n++)s[n]<f&&(h=n,f=s[n]);if(h!==l)for(s[h]=s[l],s[l]=f,n=0;n<i;n++)f=t.get(l,n),t.set(l,n,t.get(h,n)),t.set(h,n,f)}}function ye(i,e,s,t){let o=0,r=i-1,l,n,h,u,w,f,a;for(f=o+1;f<=r-1;f++){for(a=0,u=f;u<=r;u++)a=a+Math.abs(e.get(u,f-1));if(a!==0){for(h=0,u=r;u>=f;u--)s[u]=e.get(u,f-1)/a,h+=s[u]*s[u];for(n=Math.sqrt(h),s[f]>0&&(n=-n),h=h-s[f]*n,s[f]=s[f]-n,w=f;w<i;w++){for(l=0,u=r;u>=f;u--)l+=s[u]*e.get(u,w);for(l=l/h,u=f;u<=r;u++)e.set(u,w,e.get(u,w)-l*s[u])}for(u=0;u<=r;u++){for(l=0,w=r;w>=f;w--)l+=s[w]*e.get(u,w);for(l=l/h,w=f;w<=r;w++)e.set(u,w,e.get(u,w)-l*s[w])}s[f]=a*s[f],e.set(f,f-1,a*n)}}for(u=0;u<i;u++)for(w=0;w<i;w++)t.set(u,w,u===w?1:0);for(f=r-1;f>=o+1;f--)if(e.get(f,f-1)!==0){for(u=f+1;u<=r;u++)s[u]=e.get(u,f-1);for(w=f;w<=r;w++){for(n=0,u=f;u<=r;u++)n+=s[u]*t.get(u,w);for(n=n/s[f]/e.get(f,f-1),u=f;u<=r;u++)t.set(u,w,t.get(u,w)+n*s[u])}}}function Me(i,e,s,t,o){let r=i-1,l=0,n=i-1,h=Number.EPSILON,u=0,w=0,f=0,a=0,j=0,m=0,p=0,E=0,c,M,k,I,N,V,v,R,F,g,d,b,y,C,D;for(c=0;c<i;c++)for((c<l||c>n)&&(s[c]=o.get(c,c),e[c]=0),M=Math.max(c-1,0);M<i;M++)w=w+Math.abs(o.get(c,M));for(;r>=l;){for(I=r;I>l&&(m=Math.abs(o.get(I-1,I-1))+Math.abs(o.get(I,I)),m===0&&(m=w),!(Math.abs(o.get(I,I-1))<h*m));)I--;if(I===r)o.set(r,r,o.get(r,r)+u),s[r]=o.get(r,r),e[r]=0,r--,E=0;else if(I===r-1){if(v=o.get(r,r-1)*o.get(r-1,r),f=(o.get(r-1,r-1)-o.get(r,r))/2,a=f*f+v,p=Math.sqrt(Math.abs(a)),o.set(r,r,o.get(r,r)+u),o.set(r-1,r-1,o.get(r-1,r-1)+u),R=o.get(r,r),a>=0){for(p=f>=0?f+p:f-p,s[r-1]=R+p,s[r]=s[r-1],p!==0&&(s[r]=R-v/p),e[r-1]=0,e[r]=0,R=o.get(r,r-1),m=Math.abs(R)+Math.abs(p),f=R/m,a=p/m,j=Math.sqrt(f*f+a*a),f=f/j,a=a/j,M=r-1;M<i;M++)p=o.get(r-1,M),o.set(r-1,M,a*p+f*o.get(r,M)),o.set(r,M,a*o.get(r,M)-f*p);for(c=0;c<=r;c++)p=o.get(c,r-1),o.set(c,r-1,a*p+f*o.get(c,r)),o.set(c,r,a*o.get(c,r)-f*p);for(c=l;c<=n;c++)p=t.get(c,r-1),t.set(c,r-1,a*p+f*t.get(c,r)),t.set(c,r,a*t.get(c,r)-f*p)}else s[r-1]=R+f,s[r]=R+f,e[r-1]=p,e[r]=-p;r=r-2,E=0}else{if(R=o.get(r,r),F=0,v=0,I<r&&(F=o.get(r-1,r-1),v=o.get(r,r-1)*o.get(r-1,r)),E===10){for(u+=R,c=l;c<=r;c++)o.set(c,c,o.get(c,c)-R);m=Math.abs(o.get(r,r-1))+Math.abs(o.get(r-1,r-2)),R=F=.75*m,v=-.4375*m*m}if(E===30&&(m=(F-R)/2,m=m*m+v,m>0)){for(m=Math.sqrt(m),F<R&&(m=-m),m=R-v/((F-R)/2+m),c=l;c<=r;c++)o.set(c,c,o.get(c,c)-m);u+=m,R=F=v=.964}for(E=E+1,N=r-2;N>=I&&(p=o.get(N,N),j=R-p,m=F-p,f=(j*m-v)/o.get(N+1,N)+o.get(N,N+1),a=o.get(N+1,N+1)-p-j-m,j=o.get(N+2,N+1),m=Math.abs(f)+Math.abs(a)+Math.abs(j),f=f/m,a=a/m,j=j/m,!(N===I||Math.abs(o.get(N,N-1))*(Math.abs(a)+Math.abs(j))<h*(Math.abs(f)*(Math.abs(o.get(N-1,N-1))+Math.abs(p)+Math.abs(o.get(N+1,N+1))))));)N--;for(c=N+2;c<=r;c++)o.set(c,c-2,0),c>N+2&&o.set(c,c-3,0);for(k=N;k<=r-1&&(C=k!==r-1,k!==N&&(f=o.get(k,k-1),a=o.get(k+1,k-1),j=C?o.get(k+2,k-1):0,R=Math.abs(f)+Math.abs(a)+Math.abs(j),R!==0&&(f=f/R,a=a/R,j=j/R)),R!==0);k++)if(m=Math.sqrt(f*f+a*a+j*j),f<0&&(m=-m),m!==0){for(k!==N?o.set(k,k-1,-m*R):I!==N&&o.set(k,k-1,-o.get(k,k-1)),f=f+m,R=f/m,F=a/m,p=j/m,a=a/f,j=j/f,M=k;M<i;M++)f=o.get(k,M)+a*o.get(k+1,M),C&&(f=f+j*o.get(k+2,M),o.set(k+2,M,o.get(k+2,M)-f*p)),o.set(k,M,o.get(k,M)-f*R),o.set(k+1,M,o.get(k+1,M)-f*F);for(c=0;c<=Math.min(r,k+3);c++)f=R*o.get(c,k)+F*o.get(c,k+1),C&&(f=f+p*o.get(c,k+2),o.set(c,k+2,o.get(c,k+2)-f*j)),o.set(c,k,o.get(c,k)-f),o.set(c,k+1,o.get(c,k+1)-f*a);for(c=l;c<=n;c++)f=R*t.get(c,k)+F*t.get(c,k+1),C&&(f=f+p*t.get(c,k+2),t.set(c,k+2,t.get(c,k+2)-f*j)),t.set(c,k,t.get(c,k)-f),t.set(c,k+1,t.get(c,k+1)-f*a)}}}if(w!==0){for(r=i-1;r>=0;r--)if(f=s[r],a=e[r],a===0)for(I=r,o.set(r,r,1),c=r-1;c>=0;c--){for(v=o.get(c,c)-f,j=0,M=I;M<=r;M++)j=j+o.get(c,M)*o.get(M,r);if(e[c]<0)p=v,m=j;else if(I=c,e[c]===0?o.set(c,r,v!==0?-j/v:-j/(h*w)):(R=o.get(c,c+1),F=o.get(c+1,c),a=(s[c]-f)*(s[c]-f)+e[c]*e[c],V=(R*m-p*j)/a,o.set(c,r,V),o.set(c+1,r,Math.abs(R)>Math.abs(p)?(-j-v*V)/R:(-m-F*V)/p)),V=Math.abs(o.get(c,r)),h*V*V>1)for(M=c;M<=r;M++)o.set(M,r,o.get(M,r)/V)}else if(a<0)for(I=r-1,Math.abs(o.get(r,r-1))>Math.abs(o.get(r-1,r))?(o.set(r-1,r-1,a/o.get(r,r-1)),o.set(r-1,r,-(o.get(r,r)-f)/o.get(r,r-1))):(D=ct(0,-o.get(r-1,r),o.get(r-1,r-1)-f,a),o.set(r-1,r-1,D[0]),o.set(r-1,r,D[1])),o.set(r,r-1,0),o.set(r,r,1),c=r-2;c>=0;c--){for(g=0,d=0,M=I;M<=r;M++)g=g+o.get(c,M)*o.get(M,r-1),d=d+o.get(c,M)*o.get(M,r);if(v=o.get(c,c)-f,e[c]<0)p=v,j=g,m=d;else if(I=c,e[c]===0?(D=ct(-g,-d,v,a),o.set(c,r-1,D[0]),o.set(c,r,D[1])):(R=o.get(c,c+1),F=o.get(c+1,c),b=(s[c]-f)*(s[c]-f)+e[c]*e[c]-a*a,y=(s[c]-f)*2*a,b===0&&y===0&&(b=h*w*(Math.abs(v)+Math.abs(a)+Math.abs(R)+Math.abs(F)+Math.abs(p))),D=ct(R*j-p*g+a*d,R*m-p*d-a*g,b,y),o.set(c,r-1,D[0]),o.set(c,r,D[1]),Math.abs(R)>Math.abs(p)+Math.abs(a)?(o.set(c+1,r-1,(-g-v*o.get(c,r-1)+a*o.get(c,r))/R),o.set(c+1,r,(-d-v*o.get(c,r)-a*o.get(c,r-1))/R)):(D=ct(-j-F*o.get(c,r-1),-m-F*o.get(c,r),p,a),o.set(c+1,r-1,D[0]),o.set(c+1,r,D[1]))),V=Math.max(Math.abs(o.get(c,r-1)),Math.abs(o.get(c,r))),h*V*V>1)for(M=c;M<=r;M++)o.set(M,r-1,o.get(M,r-1)/V),o.set(M,r,o.get(M,r)/V)}for(c=0;c<i;c++)if(c<l||c>n)for(M=c;M<i;M++)t.set(c,M,o.get(c,M));for(M=i-1;M>=l;M--)for(c=l;c<=n;c++){for(p=0,k=l;k<=Math.min(M,n);k++)p=p+t.get(c,k)*o.get(k,M);t.set(c,M,p)}}}function ct(i,e,s,t){let o,r;return Math.abs(s)>Math.abs(t)?(o=t/s,r=s+o*t,[(i+o*e)/r,(e-o*i)/r]):(o=s/t,r=t+o*s,[(o*i+e)/r,(o*e-i)/r])}class vt{constructor(e){if(e=Q.checkMatrix(e),!e.isSymmetric())throw new Error("Matrix is not symmetric");let s=e,t=s.rows,o=new S(t,t),r=!0,l,n,h;for(n=0;n<t;n++){let u=0;for(h=0;h<n;h++){let w=0;for(l=0;l<h;l++)w+=o.get(h,l)*o.get(n,l);w=(s.get(n,h)-w)/o.get(h,h),o.set(n,h,w),u=u+w*w}for(u=s.get(n,n)-u,r&&=u>0,o.set(n,n,Math.sqrt(Math.max(u,0))),h=n+1;h<t;h++)o.set(n,h,0)}this.L=o,this.positiveDefinite=r}isPositiveDefinite(){return this.positiveDefinite}solve(e){e=Q.checkMatrix(e);let s=this.L,t=s.rows;if(e.rows!==t)throw new Error("Matrix dimensions do not match");if(this.isPositiveDefinite()===!1)throw new Error("Matrix is not positive definite");let o=e.columns,r=e.clone(),l,n,h;for(h=0;h<t;h++)for(n=0;n<o;n++){for(l=0;l<h;l++)r.set(h,n,r.get(h,n)-r.get(l,n)*s.get(h,l));r.set(h,n,r.get(h,n)/s.get(h,h))}for(h=t-1;h>=0;h--)for(n=0;n<o;n++){for(l=h+1;l<t;l++)r.set(h,n,r.get(h,n)-r.get(l,n)*s.get(l,h));r.set(h,n,r.get(h,n)/s.get(h,h))}return r}get lowerTriangularMatrix(){return this.L}}class Rt{constructor(e,s={}){e=Q.checkMatrix(e);let{Y:t}=s;const{scaleScores:o=!1,maxIterations:r=1e3,terminationCriteria:l=1e-10}=s;let n;if(t){if(B(t)&&typeof t[0]=="number"?t=S.columnVector(t):t=Q.checkMatrix(t),t.rows!==e.rows)throw new Error("Y should have the same number of rows as X");n=t.getColumnVector(0)}else n=e.getColumnVector(0);let h=1,u,w,f,a;for(let j=0;j<r&&h>l;j++)f=e.transpose().mmul(n).div(n.transpose().mmul(n).get(0,0)),f=f.div(f.norm()),u=e.mmul(f).div(f.transpose().mmul(f).get(0,0)),j>0&&(h=u.clone().sub(a).pow(2).sum()),a=u.clone(),t?(w=t.transpose().mmul(u).div(u.transpose().mmul(u).get(0,0)),w=w.div(w.norm()),n=t.mmul(w).div(w.transpose().mmul(w).get(0,0))):n=u;if(t){let j=e.transpose().mmul(u).div(u.transpose().mmul(u).get(0,0));j=j.div(j.norm());let m=e.clone().sub(u.clone().mmul(j.transpose())),p=n.transpose().mmul(u).div(u.transpose().mmul(u).get(0,0)),E=t.clone().sub(u.clone().mulS(p.get(0,0)).mmul(w.transpose()));this.t=u,this.p=j.transpose(),this.w=f.transpose(),this.q=w,this.u=n,this.s=u.transpose().mmul(u),this.xResidual=m,this.yResidual=E,this.betas=p}else this.w=f.transpose(),this.s=u.transpose().mmul(u).sqrt(),o?this.t=u.clone().div(this.s.get(0,0)):this.t=u,this.xResidual=e.sub(u.mmul(f.transpose()))}}return q.__exports.AbstractMatrix=T,q.__exports.CHO=vt,q.__exports.CholeskyDecomposition=vt,q.__exports.DistanceMatrix=lt,q.__exports.EVD=kt,q.__exports.EigenvalueDecomposition=kt,q.__exports.LU=ut,q.__exports.LuDecomposition=ut,q.__exports.Matrix=S,q.__exports.MatrixColumnSelectionView=ee,q.__exports.MatrixColumnView=te,q.__exports.MatrixFlipColumnView=se,q.__exports.MatrixFlipRowView=oe,q.__exports.MatrixRowSelectionView=ie,q.__exports.MatrixRowView=re,q.__exports.MatrixSelectionView=ht,q.__exports.MatrixSubView=ne,q.__exports.MatrixTransposeView=le,q.__exports.NIPALS=Rt,q.__exports.Nipals=Rt,q.__exports.QR=pt,q.__exports.QrDecomposition=pt,q.__exports.SVD=ot,q.__exports.SingularValueDecomposition=ot,q.__exports.SymmetricMatrix=H,q.__exports.WrapperMatrix1D=Et,q.__exports.WrapperMatrix2D=Q,q.__exports.correlation=ae,q.__exports.covariance=me,q.__exports.default=S,q.__exports.determinant=ft,q.__exports.inverse=ue,q.__exports.linearDependencies=ge,q.__exports.pseudoInverse=we,q.__exports.solve=bt,q.__exports.wrap=he,q.__exports}exports.__require=Te;
+//# sourceMappingURL=index.cjs564.js.map