numpy-ts 0.4.0 → 0.6.0

package/dist/numpy-ts.node-io.cjs

@@ -0,0 +1,3 @@
+"use strict";var Gr=Object.defineProperty;var ss=Object.getOwnPropertyDescriptor;var is=Object.getOwnPropertyNames;var us=Object.prototype.hasOwnProperty;var $=(r,t)=>()=>(r&&(t=r(r=0)),t);var Hr=(r,t)=>{for(var e in t)Gr(r,e,{get:t[e],enumerable:!0})},cs=(r,t,e,n)=>{if(t&&typeof t=="object"||typeof t=="function")for(let o of is(t))!us.call(r,o)&&o!==e&&Gr(r,o,{get:()=>t[o],enumerable:!(n=ss(t,o))||n.enumerable});return r};var Yr=r=>cs(Gr({},"__esModule",{value:!0}),r);function Mt(r){if(!r.includes(":")){if(r.includes("."))throw new Error(`Invalid slice index: "${r}" (must be integer)`);let a=parseInt(r,10);if(isNaN(a))throw new Error(`Invalid slice index: "${r}"`);return{start:a,stop:null,step:1,isIndex:!0}}let t=r.split(":");if(t.length>3)throw new Error(`Invalid slice notation: "${r}" (too many colons)`);let e=t[0]===""?null:parseInt(t[0],10),n=t[1]===""||t[1]===void 0?null:parseInt(t[1],10),o=t[2]===""||t[2]===void 0?1:parseInt(t[2],10);if(e!==null&&isNaN(e))throw new Error(`Invalid start index in slice: "${r}"`);if(n!==null&&isNaN(n))throw new Error(`Invalid stop index in slice: "${r}"`);if(isNaN(o))throw new Error(`Invalid step in slice: "${r}"`);if(o===0)throw new Error("Slice step cannot be zero");return{start:e,stop:n,step:o,isIndex:!1}}function Ct(r,t){let{start:e,stop:n}=r,{step:o,isIndex:a}=r;if(a){if(e===null)throw new Error("Index cannot be null");let i=e<0?t+e:e;if(i<0||i>=t)throw new Error(`Index ${e} is out of bounds for size ${t}`);return{start:i,stop:i+1,step:1,isIndex:!0}}return o>0?(e===null&&(e=0),n===null&&(n=t)):(e===null&&(e=t-1),n===null&&(n=-t-1)),e<0&&(e=t+e),n<0&&(n=t+n),e=Math.max(0,Math.min(e,t)),n=Math.max(-1,Math.min(n,t)),{start:e,stop:n,step:o,isIndex:!1}}var kt=$(()=>{"use strict"});function _(r){switch(r){case"float64":return Float64Array;case"float32":return Float32Array;case"int64":return BigInt64Array;case"int32":return Int32Array;case"int16":return Int16Array;case"int8":return Int8Array;case"uint64":return BigUint64Array;case"uint32":return Uint32Array;case"uint16":return Uint16Array;case"uint8":return Uint8Array;case"bool":return Uint8Array;default:throw new Error(`Unknown dtype: ${r}`)}}function Vt(r){switch(r){case"float64":case"int64":case"uint64":return 8;case"float32":case"int32":case"uint32":return 4;case"int16":case"uint16":return 2;case"int8":case"uint8":case"bool":return 1;default:throw new Error(`Unknown dtype: ${r}`)}}function Pt(r){return r==="float64"||r==="float32"}function N(r){return r==="int64"||r==="uint64"}function M(r,t){if(r===t)return r;if(r==="bool")return t;if(t==="bool")return r;if(Pt(r)||Pt(t)){if(r==="float64"||t==="float64")return"float64";if(r==="float32"){let c=t;return c==="int32"||c==="int64"||c==="uint32"||c==="uint64"?"float64":"float32"}if(t==="float32"){let c=r;return c==="int32"||c==="int64"||c==="uint32"||c==="uint64"?"float64":"float32"}return"float32"}let e=r.startsWith("int"),n=t.startsWith("int"),o=r.startsWith("uint"),a=t.startsWith("uint"),i=c=>c.includes("64")?64:c.includes("32")?32:c.includes("16")?16:c.includes("8")?8:0,s=i(r),u=i(t);if(r==="int64"&&t==="uint64"||r==="uint64"&&t==="int64")return"float64";if(e&&a&&s===u){if(s===8)return"int16";if(s===16)return"int32";if(s===32)return"int64"}if(o&&n&&s===u){if(u===8)return"int16";if(u===16)return"int32";if(u===32)return"int64"}if(e&&n||o&&a){let c=Math.max(s,u);return e?c===64?"int64":c===32?"int32":c===16?"int16":"int8":c===64?"uint64":c===32?"uint32":c===16?"uint16":"uint8"}return e&&a?s>u?r:u===8?"int16":u===16?"int32":u===32?"int64":"float64":o&&n?u>s?t:s===8?"int16":s===16?"int32":s===32?"int64":"float64":"float64"}var v,V=$(()=>{"use strict";v="float64"});function Z(r){let t=new Array(r.length),e=1;for(let n=r.length-1;n>=0;n--)t[n]=e,e*=r[n];return t}var D,F=$(()=>{"use strict";V();D=class r{constructor(t,e,n,o,a){this._data=t,this._shape=e,this._strides=n,this._offset=o,this._dtype=a}get shape(){return this._shape}get ndim(){return this._shape.length}get size(){return this._shape.reduce((t,e)=>t*e,1)}get dtype(){return this._dtype}get data(){return this._data}get strides(){return this._strides}get offset(){return this._offset}get isCContiguous(){let t=this._shape,e=this._strides,n=t.length;if(n===0)return!0;if(n===1)return e[0]===1;let o=1;for(let a=n-1;a>=0;a--){if(e[a]!==o)return!1;o*=t[a]}return!0}get isFContiguous(){let t=this._shape,e=this._strides,n=t.length;if(n===0)return!0;if(n===1)return e[0]===1;let o=1;for(let a=0;a<n;a++){if(e[a]!==o)return!1;o*=t[a]}return!0}iget(t){let e=this._shape,n=this._strides,o=e.length;if(o===0)return this._data[this._offset];let a=t,i=this._offset;for(let s=0;s<o;s++){let u=1;for(let l=s+1;l<o;l++)u*=e[l];let c=Math.floor(a/u);a=a%u,i+=c*n[s]}return this._data[i]}iset(t,e){let n=this._shape,o=this._strides,a=n.length;if(a===0){this._data[this._offset]=e;return}let i=t,s=this._offset;for(let u=0;u<a;u++){let c=1;for(let y=u+1;y<a;y++)c*=n[y];let l=Math.floor(i/c);i=i%c,s+=l*o[u]}this._data[s]=e}get(...t){let e=this._strides,n=this._offset;for(let o=0;o<t.length;o++)n+=t[o]*e[o];return this._data[n]}set(t,e){let n=this._strides,o=this._offset;for(let a=0;a<t.length;a++)o+=t[a]*n[a];this._data[o]=e}copy(){let t=Array.from(this._shape),e=this._dtype,n=this.size,o=_(e);if(!o)throw new Error(`Cannot copy array with dtype ${e}`);let a=new o(n);if(this.isCContiguous&&this._offset===0)if(N(e)){let i=this._data,s=a;for(let u=0;u<n;u++)s[u]=i[u]}else a.set(this._data);else if(N(e)){let i=a;for(let s=0;s<n;s++)i[s]=this.iget(s)}else for(let i=0;i<n;i++)a[i]=this.iget(i);return new r(a,t,r._computeStrides(t),0,e)}static fromData(t,e,n,o,a){let i=o??r._computeStrides(e),s=a??0;return new r(t,e,i,s,n)}static zeros(t,e=v){let n=t.reduce((i,s)=>i*s,1),o=_(e);if(!o)throw new Error(`Cannot create array with dtype ${e}`);let a=new o(n);return new r(a,t,r._computeStrides(t),0,e)}static ones(t,e=v){let n=t.reduce((i,s)=>i*s,1),o=_(e);if(!o)throw new Error(`Cannot create array with dtype ${e}`);let a=new o(n);return N(e)?a.fill(BigInt(1)):a.fill(1),new r(a,t,r._computeStrides(t),0,e)}static _computeStrides(t){let e=new Array(t.length),n=1;for(let o=t.length-1;o>=0;o--)e[o]=n,n*=t[o];return e}}});function Ft(r,t){let e=r.length,n=t.length,o=Math.max(e,n),a=new Array(o);for(let i=0;i<o;i++){let s=i<o-e?1:r[i-(o-e)],u=i<o-n?1:t[i-(o-n)];if(s===u)a[i]=s;else if(s===1)a[i]=u;else if(u===1)a[i]=s;else throw new Error(`operands could not be broadcast together with shapes ${JSON.stringify(Array.from(r))} ${JSON.stringify(Array.from(t))}`)}return a}function fs(r,t,e){let n=r.length,o=e.length,a=new Array(o).fill(0);for(let i=0;i<n;i++){let s=o-n+i,u=r[i],c=e[s];if(u===c)a[s]=t[i];else if(u===1)a[s]=0;else throw new Error("Invalid broadcast")}return a}function Dr(r,t){let e=fs(r.shape,r.strides,t);return D.fromData(r.data,Array.from(t),r.dtype,e,r.offset)}function q(r,t,e,n){let o=Ft(r.shape,t.shape),a=Dr(r,o),i=Dr(t,o),s=M(r.dtype,t.dtype),u=D.zeros(o,s),c=u.data,l=u.size;if(N(s)){let y=c;for(let f=0;f<l;f++){let p=a.iget(f),m=i.iget(f),d=typeof p=="bigint"?p:BigInt(Math.round(p)),g=typeof m=="bigint"?m:BigInt(Math.round(m));n==="add"?y[f]=d+g:n==="subtract"?y[f]=d-g:n==="multiply"?y[f]=d*g:n==="divide"?y[f]=d/g:y[f]=BigInt(Math.round(e(Number(d),Number(g))))}}else{let y=N(r.dtype)||N(t.dtype);for(let f=0;f<l;f++){let p=a.iget(f),m=i.iget(f),d=Number(p),g=Number(m);c[f]=e(d,g)}}return u}function G(r,t,e){let n=Ft(r.shape,t.shape),o=Dr(r,n),a=Dr(t,n),i=n.reduce((c,l)=>c*l,1),s=new Uint8Array(i),u=N(r.dtype)||N(t.dtype);for(let c=0;c<i;c++){let l=o.iget(c),y=a.iget(c),f=Number(l),p=Number(y);s[c]=e(f,p)?1:0}return D.fromData(s,n,"bool")}function C(r,t,e=!0){let n=r.dtype,o=Array.from(r.shape),a=r.size,s=e?n:n!=="float32"&&n!=="float64"?"float64":n,u=D.zeros(o,s),c=u.data,l=r.data;if(N(n))if(N(s)){let y=c;for(let f=0;f<a;f++){let p=Number(l[f]);y[f]=BigInt(Math.round(t(p)))}}else for(let y=0;y<a;y++)c[y]=t(Number(l[y]));else for(let y=0;y<a;y++)c[y]=t(Number(l[y]));return u}var tr=$(()=>{"use strict";F();V()});function Jr(r,t){return r.isCContiguous&&t.isCContiguous&&r.shape.length===t.shape.length&&r.shape.every((e,n)=>e===t.shape[n])}function Rt(r,t){return typeof t=="number"?ms(r,t):Jr(r,t)?ls(r,t):q(r,t,(e,n)=>e+n,"add")}function ls(r,t){let e=M(r.dtype,t.dtype),n=D.zeros(Array.from(r.shape),e),o=r.size,a=r.data,i=t.data,s=n.data;if(N(e)){let u=s;if(!N(r.dtype)||!N(t.dtype))for(let l=0;l<o;l++){let y=typeof a[l]=="bigint"?a[l]:BigInt(Math.round(Number(a[l]))),f=typeof i[l]=="bigint"?i[l]:BigInt(Math.round(Number(i[l])));u[l]=y+f}else{let l=a,y=i;for(let f=0;f<o;f++)u[f]=l[f]+y[f]}}else if(N(r.dtype)||N(t.dtype))for(let c=0;c<o;c++){let l=typeof a[c]=="bigint"?Number(a[c]):a[c],y=typeof i[c]=="bigint"?Number(i[c]):i[c];s[c]=l+y}else for(let c=0;c<o;c++)s[c]=a[c]+i[c];return n}function qt(r,t){return typeof t=="number"?ds(r,t):Jr(r,t)?ps(r,t):q(r,t,(e,n)=>e-n,"subtract")}function ps(r,t){let e=M(r.dtype,t.dtype),n=D.zeros(Array.from(r.shape),e),o=r.size,a=r.data,i=t.data,s=n.data;if(N(e)){let u=s;if(!N(r.dtype)||!N(t.dtype))for(let l=0;l<o;l++){let y=typeof a[l]=="bigint"?a[l]:BigInt(Math.round(Number(a[l]))),f=typeof i[l]=="bigint"?i[l]:BigInt(Math.round(Number(i[l])));u[l]=y-f}else{let l=a,y=i;for(let f=0;f<o;f++)u[f]=l[f]-y[f]}}else if(N(r.dtype)||N(t.dtype))for(let c=0;c<o;c++){let l=typeof a[c]=="bigint"?Number(a[c]):a[c],y=typeof i[c]=="bigint"?Number(i[c]):i[c];s[c]=l-y}else for(let c=0;c<o;c++)s[c]=a[c]-i[c];return n}function Lt(r,t){return typeof t=="number"?gs(r,t):Jr(r,t)?ys(r,t):q(r,t,(e,n)=>e*n,"multiply")}function ys(r,t){let e=M(r.dtype,t.dtype),n=D.zeros(Array.from(r.shape),e),o=r.size,a=r.data,i=t.data,s=n.data;if(N(e)){let u=s;if(!N(r.dtype)||!N(t.dtype))for(let l=0;l<o;l++){let y=typeof a[l]=="bigint"?a[l]:BigInt(Math.round(Number(a[l]))),f=typeof i[l]=="bigint"?i[l]:BigInt(Math.round(Number(i[l])));u[l]=y*f}else{let l=a,y=i;for(let f=0;f<o;f++)u[f]=l[f]*y[f]}}else if(N(r.dtype)||N(t.dtype))for(let c=0;c<o;c++){let l=typeof a[c]=="bigint"?Number(a[c]):a[c],y=typeof i[c]=="bigint"?Number(i[c]):i[c];s[c]=l*y}else for(let c=0;c<o;c++)s[c]=a[c]*i[c];return n}function jt(r,t){if(typeof t=="number")return hs(r,t);let e=r.dtype==="float64",n=t.dtype==="float64",o=r.dtype==="float32",a=t.dtype==="float32";if(e||n){let u=e?r:er(r,"float64"),c=n?t:er(t,"float64");return q(u,c,(l,y)=>l/y,"divide")}if(o||a){let u=o?r:er(r,"float32"),c=a?t:er(t,"float32");return q(u,c,(l,y)=>l/y,"divide")}let i=er(r,"float64"),s=er(t,"float64");return q(i,s,(u,c)=>u/c,"divide")}function er(r,t){let e=D.zeros(Array.from(r.shape),t),n=r.size,o=r.data,a=e.data;for(let i=0;i<n;i++)a[i]=Number(o[i]);return e}function ms(r,t){let e=r.dtype,n=Array.from(r.shape),o=r.data,a=r.size,i=D.zeros(n,e),s=i.data;if(N(e)){let u=o,c=s,l=BigInt(Math.round(t));for(let y=0;y<a;y++)c[y]=u[y]+l}else for(let u=0;u<a;u++)s[u]=Number(o[u])+t;return i}function ds(r,t){let e=r.dtype,n=Array.from(r.shape),o=r.data,a=r.size,i=D.zeros(n,e),s=i.data;if(N(e)){let u=o,c=s,l=BigInt(Math.round(t));for(let y=0;y<a;y++)c[y]=u[y]-l}else for(let u=0;u<a;u++)s[u]=Number(o[u])-t;return i}function gs(r,t){let e=r.dtype,n=Array.from(r.shape),o=r.data,a=r.size,i=D.zeros(n,e),s=i.data;if(N(e)){let u=o,c=s,l=BigInt(Math.round(t));for(let y=0;y<a;y++)c[y]=u[y]*l}else for(let u=0;u<a;u++)s[u]=Number(o[u])*t;return i}function hs(r,t){let e=r.dtype,n=Array.from(r.shape),o=r.data,a=r.size,s=e!=="float32"&&e!=="float64"?"float64":e,u=D.zeros(n,s),c=u.data;if(N(e))for(let l=0;l<a;l++)c[l]=Number(o[l])/t;else for(let l=0;l<a;l++)c[l]=Number(o[l])/t;return u}function Zt(r){let t=r.dtype,e=Array.from(r.shape),n=r.data,o=r.size,a=D.zeros(e,t),i=a.data;if(N(t)){let s=n,u=i;for(let c=0;c<o;c++){let l=s[c];u[c]=l<0n?-l:l}}else for(let s=0;s<o;s++)i[s]=Math.abs(Number(n[s]));return a}function Wt(r){let t=r.dtype,e=Array.from(r.shape),n=r.data,o=r.size,a=D.zeros(e,t),i=a.data;if(N(t)){let s=n,u=i;for(let c=0;c<o;c++)u[c]=-s[c]}else for(let s=0;s<o;s++)i[s]=-Number(n[s]);return a}function Gt(r){let t=r.dtype,e=Array.from(r.shape),n=r.data,o=r.size,a=D.zeros(e,t),i=a.data;if(N(t)){let s=n,u=i;for(let c=0;c<o;c++){let l=s[c];u[c]=l>0n?1n:l<0n?-1n:0n}}else for(let s=0;s<o;s++){let u=Number(n[s]);i[s]=u>0?1:u<0?-1:0}return a}function Nr(r,t){return typeof t=="number"?As(r,t):q(r,t,(e,n)=>(e%n+n)%n,"mod")}function As(r,t){let e=r.dtype,n=Array.from(r.shape),o=r.data,a=r.size,i=D.zeros(n,e),s=i.data;if(N(e)){let u=o,c=s,l=BigInt(Math.round(t));for(let y=0;y<a;y++){let f=u[y];c[y]=(f%l+l)%l}}else for(let u=0;u<a;u++){let c=Number(o[u]);s[u]=(c%t+t)%t}return i}function Qr(r,t){return typeof t=="number"?bs(r,t):q(r,t,(e,n)=>Math.floor(e/n),"floor_divide")}function bs(r,t){let e=r.dtype,n=Array.from(r.shape),o=r.data,a=r.size,i=D.zeros(n,e),s=i.data;if(N(e)){let u=o,c=s,l=BigInt(Math.round(t));for(let y=0;y<a;y++)c[y]=u[y]/l}else for(let u=0;u<a;u++)s[u]=Math.floor(Number(o[u])/t);return i}function Ht(r){let t=r.dtype,e=Array.from(r.shape),n=r.data,o=r.size,a=D.zeros(e,t),i=a.data;for(let s=0;s<o;s++)i[s]=n[s];return a}function Yt(r){let t=r.dtype,e=Array.from(r.shape),n=r.data,o=r.size,i=t!=="float32"&&t!=="float64"?"float64":t,s=D.zeros(e,i),u=s.data;if(N(t))for(let c=0;c<o;c++)u[c]=1/Number(n[c]);else for(let c=0;c<o;c++)u[c]=1/Number(n[c]);return s}function Jt(r){let t=r.dtype,e=Array.from(r.shape),n=r.data,o=r.size,i=t!=="float32"&&t!=="float64"?"float64":t,s=D.zeros(e,i),u=s.data;for(let c=0;c<o;c++)u[c]=Math.cbrt(Number(n[c]));return s}function Qt(r){let t=r.dtype,e=Array.from(r.shape),n=r.data,o=r.size,a=t==="float32"?"float32":"float64",i=D.zeros(e,a),s=i.data;for(let u=0;u<o;u++)s[u]=Math.abs(Number(n[u]));return i}function Xt(r,t){let e=Qr(r,t),n=Nr(r,t);return[e,n]}function Kt(r){let t=r.dtype,e=Array.from(r.shape),n=r.data,o=r.size,a=D.zeros(e,t),i=a.data;if(N(t)){let s=n,u=i;for(let c=0;c<o;c++)u[c]=s[c]*s[c]}else for(let s=0;s<o;s++){let u=Number(n[s]);i[s]=u*u}return a}function re(r,t){return Nr(r,t)}function te(r,t){let e=r.dtype,n=Array.from(r.shape),o=r.size,a=e==="float32"?"float32":"float64",i=D.zeros(n,a),s=i.data;if(typeof t=="number")for(let u=0;u<o;u++){let c=Number(r.data[u]);c<0?s[u]=0:c===0?s[u]=t:s[u]=1}else{let u=t.data,c=t.shape;if(n.every((l,y)=>l===c[y]))for(let l=0;l<o;l++){let y=Number(r.data[l]);y<0?s[l]=0:y===0?s[l]=Number(u[l]):s[l]=1}else for(let l=0;l<o;l++){let y=Number(r.data[l]),f=l%t.size;y<0?s[l]=0:y===0?s[l]=Number(u[f]):s[l]=1}}return i}var ee=$(()=>{"use strict";F();V();tr()});function J(r){if(r.length===0)return[];if(r.length===1)return Array.from(r[0]);let t=Math.max(...r.map(n=>n.length)),e=new Array(t);for(let n=0;n<t;n++){let o=1;for(let a of r){let i=a.length-t+n,s=i<0?1:a[i];if(s!==1){if(o===1)o=s;else if(o!==s)return null}}e[n]=o}return e}function Ds(r,t,e){let n=r.length,o=e.length,a=new Array(o).fill(0);for(let i=0;i<n;i++){let s=o-n+i,u=r[i],c=e[s];if(u===c)a[s]=t[i];else if(u===1)a[s]=0;else throw new Error("Invalid broadcast")}return a}function H(r,t){let e=Ds(r.shape,r.strides,t);return D.fromData(r.data,Array.from(t),r.dtype,e,r.offset)}function Xr(...r){let t=J(r);if(t===null){let e=r.map(n=>`(${n.join(",")})`).join(" ");throw new Error(`shape mismatch: objects cannot be broadcast to a single shape. Mismatch is between ${e}`)}return t}var Kr=$(()=>{"use strict";F()});function ne(r,t){return typeof t=="number"?Ns(r,t):G(r,t,(e,n)=>e>n)}function oe(r,t){return typeof t=="number"?ws(r,t):G(r,t,(e,n)=>e>=n)}function ae(r,t){return typeof t=="number"?xs(r,t):G(r,t,(e,n)=>e<n)}function se(r,t){return typeof t=="number"?Is(r,t):G(r,t,(e,n)=>e<=n)}function ie(r,t){return typeof t=="number"?zs(r,t):G(r,t,(e,n)=>e===n)}function ue(r,t){return typeof t=="number"?Bs(r,t):G(r,t,(e,n)=>e!==n)}function rt(r,t,e=1e-5,n=1e-8){return typeof t=="number"?_s(r,t,e,n):G(r,t,(o,a)=>{let i=Math.abs(o-a),s=n+e*Math.abs(a);return i<=s})}function ce(r,t,e=1e-5,n=1e-8){let o=rt(r,t,e,n),a=o.data;for(let i=0;i<o.size;i++)if(a[i]===0)return!1;return!0}function fe(r,t){let e=[Array.from(r.shape),Array.from(t.shape)],n=J(e);if(n===null)return!1;let o=H(r,n),a=H(t,n),i=n.length,s=n.reduce((l,y)=>l*y,1),u=N(o.dtype),c=N(a.dtype);for(let l=0;l<s;l++){let y=l,f=new Array(i);for(let d=i-1;d>=0;d--)f[d]=y%n[d],y=Math.floor(y/n[d]);let p=o.get(...f),m=a.get(...f);if(u||c){let d=typeof p=="bigint"?p:BigInt(Number(p)),g=typeof m=="bigint"?m:BigInt(Number(m));if(d!==g)return!1}else if(p!==m)return!1}return!0}function Ns(r,t){let e=new Uint8Array(r.size),n=r.data;for(let o=0;o<r.size;o++)e[o]=n[o]>t?1:0;return D.fromData(e,Array.from(r.shape),"bool")}function ws(r,t){let e=new Uint8Array(r.size),n=r.data;for(let o=0;o<r.size;o++)e[o]=n[o]>=t?1:0;return D.fromData(e,Array.from(r.shape),"bool")}function xs(r,t){let e=new Uint8Array(r.size),n=r.data;for(let o=0;o<r.size;o++)e[o]=n[o]<t?1:0;return D.fromData(e,Array.from(r.shape),"bool")}function Is(r,t){let e=new Uint8Array(r.size),n=r.data;for(let o=0;o<r.size;o++)e[o]=n[o]<=t?1:0;return D.fromData(e,Array.from(r.shape),"bool")}function zs(r,t){let e=new Uint8Array(r.size),n=r.data,o=r.dtype;if(N(o)){let a=BigInt(Math.round(t)),i=n;for(let s=0;s<r.size;s++)e[s]=i[s]===a?1:0}else for(let a=0;a<r.size;a++)e[a]=n[a]===t?1:0;return D.fromData(e,Array.from(r.shape),"bool")}function Bs(r,t){let e=new Uint8Array(r.size),n=r.data;for(let o=0;o<r.size;o++)e[o]=n[o]!==t?1:0;return D.fromData(e,Array.from(r.shape),"bool")}function _s(r,t,e,n){let o=new Uint8Array(r.size),a=r.data,i=r.dtype;if(N(i)){let s=a;for(let u=0;u<r.size;u++){let c=Number(s[u]),l=Math.abs(c-t),y=n+e*Math.abs(t);o[u]=l<=y?1:0}}else for(let s=0;s<r.size;s++){let u=Number(a[s]),c=Math.abs(u-t),l=n+e*Math.abs(t);o[s]=c<=l?1:0}return D.fromData(o,Array.from(r.shape),"bool")}var le=$(()=>{"use strict";F();V();tr();Kr()});function U(r,t){let e=0,n=1;for(let o=r.length-1;o>=0;o--)e+=r[o]*n,n*=t[o];return e}function T(r,t,e,n){let o=n.length,a=new Array(o),i=Array.from(n).filter((u,c)=>c!==t),s=r;for(let u=i.length-1;u>=0;u--)a[u>=t?u+1:u]=s%i[u],s=Math.floor(s/i[u]);return a[t]=e,a}var pe=$(()=>{"use strict"});function ur(r,t,e=!1){let n=r.dtype,o=r.shape,a=o.length,i=r.size,s=r.data;if(t===void 0)if(N(n)){let m=s,d=BigInt(0);for(let g=0;g<i;g++)d+=m[g];return Number(d)}else{let m=0;for(let d=0;d<i;d++)m+=Number(s[d]);return m}let u=t;if(u<0&&(u=a+u),u<0||u>=a)throw new Error(`axis ${t} is out of bounds for array of dimension ${a}`);let c=Array.from(o).filter((m,d)=>d!==u);if(c.length===0)return ur(r);let l=D.zeros(c,n),y=l.data,f=o[u],p=c.reduce((m,d)=>m*d,1);if(N(n)){let m=s,d=y;for(let g=0;g<p;g++){let A=BigInt(0);for(let h=0;h<f;h++){let b=T(g,u,h,o),S=U(b,o);A+=m[S]}d[g]=A}}else for(let m=0;m<p;m++){let d=0;for(let g=0;g<f;g++){let A=T(m,u,g,o),h=U(A,o);d+=Number(s[h])}y[m]=d}if(e){let m=[...o];return m[u]=1,D.fromData(y,m,n)}return l}function wr(r,t,e=!1){let n=r.dtype,o=r.shape;if(t===void 0)return ur(r)/r.size;let a=t;if(a<0&&(a=o.length+a),a<0||a>=o.length)throw new Error(`axis ${t} is out of bounds for array of dimension ${o.length}`);let i=ur(r,t,e);if(typeof i=="number")return i/o[a];let s=o[a],u=n;(N(n)||n.startsWith("int")||n.startsWith("uint"))&&(u="float64");let c=D.zeros(Array.from(i.shape),u),l=c.data,y=i.data;if(N(n)){let f=y;for(let p=0;p<l.length;p++)l[p]=Number(f[p])/s}else for(let f=0;f<l.length;f++)l[f]=Number(y[f])/s;return c}function xr(r,t,e=!1){let n=r.dtype,o=r.shape,a=o.length,i=r.size,s=r.data;if(t===void 0){if(i===0)throw new Error("max of empty array");let m=s[0];for(let d=1;d<i;d++)s[d]>m&&(m=s[d]);return Number(m)}let u=t;if(u<0&&(u=a+u),u<0||u>=a)throw new Error(`axis ${t} is out of bounds for array of dimension ${a}`);let c=Array.from(o).filter((m,d)=>d!==u);if(c.length===0)return xr(r);let l=D.zeros(c,n),y=l.data,f=o[u],p=c.reduce((m,d)=>m*d,1);if(N(n)){let m=s,d=y;for(let g=0;g<p;g++){let A=T(g,u,0,o),h=U(A,o),b=m[h];for(let S=1;S<f;S++){let w=T(g,u,S,o),x=U(w,o),z=m[x];z>b&&(b=z)}d[g]=b}}else for(let m=0;m<p;m++){let d=-1/0;for(let g=0;g<f;g++){let A=T(m,u,g,o),h=U(A,o),b=Number(s[h]);b>d&&(d=b)}y[m]=d}if(e){let m=[...o];return m[u]=1,D.fromData(y,m,n)}return l}function tt(r,t,e=!1){let n=r.dtype,o=r.shape,a=o.length,i=r.size,s=r.data;if(t===void 0)if(N(n)){let m=s,d=BigInt(1);for(let g=0;g<i;g++)d*=m[g];return Number(d)}else{let m=1;for(let d=0;d<i;d++)m*=Number(s[d]);return m}let u=t;if(u<0&&(u=a+u),u<0||u>=a)throw new Error(`axis ${t} is out of bounds for array of dimension ${a}`);let c=Array.from(o).filter((m,d)=>d!==u);if(c.length===0)return tt(r);let l=D.zeros(c,n),y=l.data,f=o[u],p=c.reduce((m,d)=>m*d,1);if(N(n)){let m=s,d=y;for(let g=0;g<p;g++){let A=BigInt(1);for(let h=0;h<f;h++){let b=T(g,u,h,o),S=U(b,o);A*=m[S]}d[g]=A}}else for(let m=0;m<p;m++){let d=1;for(let g=0;g<f;g++){let A=T(m,u,g,o),h=U(A,o);d*=Number(s[h])}y[m]=d}if(e){let m=[...o];return m[u]=1,D.fromData(y,m,n)}return l}function Ir(r,t,e=!1){let n=r.dtype,o=r.shape,a=o.length,i=r.size,s=r.data;if(t===void 0){if(i===0)throw new Error("min of empty array");let m=s[0];for(let d=1;d<i;d++)s[d]<m&&(m=s[d]);return Number(m)}let u=t;if(u<0&&(u=a+u),u<0||u>=a)throw new Error(`axis ${t} is out of bounds for array of dimension ${a}`);let c=Array.from(o).filter((m,d)=>d!==u);if(c.length===0)return Ir(r);let l=D.zeros(c,n),y=l.data,f=o[u],p=c.reduce((m,d)=>m*d,1);if(N(n)){let m=s,d=y;for(let g=0;g<p;g++){let A=T(g,u,0,o),h=U(A,o),b=m[h];for(let S=1;S<f;S++){let w=T(g,u,S,o),x=U(w,o),z=m[x];z<b&&(b=z)}d[g]=b}}else for(let m=0;m<p;m++){let d=1/0;for(let g=0;g<f;g++){let A=T(m,u,g,o),h=U(A,o),b=Number(s[h]);b<d&&(d=b)}y[m]=d}if(e){let m=[...o];return m[u]=1,D.fromData(y,m,n)}return l}function et(r,t){let e=r.dtype,n=r.shape,o=n.length,a=r.size,i=r.data;if(t===void 0){if(a===0)throw new Error("argmin of empty array");let p=i[0],m=0;for(let d=1;d<a;d++)i[d]<p&&(p=i[d],m=d);return m}let s=t;if(s<0&&(s=o+s),s<0||s>=o)throw new Error(`axis ${t} is out of bounds for array of dimension ${o}`);let u=Array.from(n).filter((p,m)=>m!==s);if(u.length===0)return et(r);let c=D.zeros(u,"int32"),l=c.data,y=n[s],f=u.reduce((p,m)=>p*m,1);if(N(e)){let p=i;for(let m=0;m<f;m++){let d=T(m,s,0,n),g=U(d,n),A=p[g],h=0;for(let b=1;b<y;b++){let S=T(m,s,b,n),w=U(S,n),x=p[w];x<A&&(A=x,h=b)}l[m]=h}}else for(let p=0;p<f;p++){let m=1/0,d=0;for(let g=0;g<y;g++){let A=T(p,s,g,n),h=U(A,n),b=Number(i[h]);b<m&&(m=b,d=g)}l[p]=d}return c}function nt(r,t){let e=r.dtype,n=r.shape,o=n.length,a=r.size,i=r.data;if(t===void 0){if(a===0)throw new Error("argmax of empty array");let p=i[0],m=0;for(let d=1;d<a;d++)i[d]>p&&(p=i[d],m=d);return m}let s=t;if(s<0&&(s=o+s),s<0||s>=o)throw new Error(`axis ${t} is out of bounds for array of dimension ${o}`);let u=Array.from(n).filter((p,m)=>m!==s);if(u.length===0)return nt(r);let c=D.zeros(u,"int32"),l=c.data,y=n[s],f=u.reduce((p,m)=>p*m,1);if(N(e)){let p=i;for(let m=0;m<f;m++){let d=T(m,s,0,n),g=U(d,n),A=p[g],h=0;for(let b=1;b<y;b++){let S=T(m,s,b,n),w=U(S,n),x=p[w];x>A&&(A=x,h=b)}l[m]=h}}else for(let p=0;p<f;p++){let m=-1/0,d=0;for(let g=0;g<y;g++){let A=T(p,s,g,n),h=U(A,n),b=Number(i[h]);b>m&&(m=b,d=g)}l[p]=d}return c}function ot(r,t,e=0,n=!1){let o=r.shape,a=o.length,i=r.size,s=r.data,u=wr(r,t,n);if(t===void 0){let A=u,h=0;for(let b=0;b<i;b++){let S=Number(s[b])-A;h+=S*S}return h/(i-e)}let c=t;if(c<0&&(c=a+c),c<0||c>=a)throw new Error(`axis ${t} is out of bounds for array of dimension ${a}`);let l=o[c],y=u,f=y.data,p=n?y.shape:Array.from(o).filter((A,h)=>h!==c),m=D.zeros(Array.from(p),"float64"),d=m.data,g=p.reduce((A,h)=>A*h,1);for(let A=0;A<g;A++){let h=0,b=Number(f[A]);for(let S=0;S<l;S++){let w=T(A,c,S,o),x=U(w,o),z=Number(s[x])-b;h+=z*z}d[A]=h/(l-e)}return m}function ye(r,t,e=0,n=!1){let o=ot(r,t,e,n);if(typeof o=="number")return Math.sqrt(o);let a=D.zeros(Array.from(o.shape),"float64"),i=o.data,s=a.data;for(let u=0;u<i.length;u++)s[u]=Math.sqrt(Number(i[u]));return a}function at(r,t,e=!1){let n=r.shape,o=n.length,a=r.size,i=r.data;if(t===void 0){for(let p=0;p<a;p++)if(!i[p])return!1;return!0}let s=t;if(s<0&&(s=o+s),s<0||s>=o)throw new Error(`axis ${t} is out of bounds for array of dimension ${o}`);let u=Array.from(n).filter((p,m)=>m!==s);if(u.length===0)return at(r);let c=D.zeros(u,"bool"),l=c.data,y=n[s],f=u.reduce((p,m)=>p*m,1);for(let p=0;p<f;p++){let m=!0;for(let d=0;d<y;d++){let g=T(p,s,d,n),A=U(g,n);if(!i[A]){m=!1;break}}l[p]=m?1:0}if(e){let p=[...n];return p[s]=1,D.fromData(l,p,"bool")}return c}function st(r,t,e=!1){let n=r.shape,o=n.length,a=r.size,i=r.data;if(t===void 0){for(let p=0;p<a;p++)if(i[p])return!0;return!1}let s=t;if(s<0&&(s=o+s),s<0||s>=o)throw new Error(`axis ${t} is out of bounds for array of dimension ${o}`);let u=Array.from(n).filter((p,m)=>m!==s);if(u.length===0)return st(r);let c=D.zeros(u,"bool"),l=c.data,y=n[s],f=u.reduce((p,m)=>p*m,1);for(let p=0;p<f;p++){let m=!1;for(let d=0;d<y;d++){let g=T(p,s,d,n),A=U(g,n);if(i[A]){m=!0;break}}l[p]=m?1:0}if(e){let p=[...n];return p[s]=1,D.fromData(l,p,"bool")}return c}function it(r,t){let e=r.shape,n=e.length,o=r.data;if(t===void 0){let f=r.size,p=new Float64Array(f),m=0;for(let d=0;d<f;d++)m+=Number(o[d]),p[d]=m;return D.fromData(p,[f],"float64")}let a=t;if(a<0&&(a=n+a),a<0||a>=n)throw new Error(`axis ${t} is out of bounds for array of dimension ${n}`);let i=new Float64Array(r.size),s=e[a],u=[],c=1;for(let f=n-1;f>=0;f--)u.unshift(c),c*=e[f];let l=r.size,y=u[a];for(let f=0;f<l;f++)Math.floor(f/y)%s===0?i[f]=Number(o[f]):i[f]=i[f-y]+Number(o[f]);return D.fromData(i,[...e],"float64")}function ut(r,t){let e=r.shape,n=e.length,o=r.data;if(t===void 0){let f=r.size,p=new Float64Array(f),m=1;for(let d=0;d<f;d++)m*=Number(o[d]),p[d]=m;return D.fromData(p,[f],"float64")}let a=t;if(a<0&&(a=n+a),a<0||a>=n)throw new Error(`axis ${t} is out of bounds for array of dimension ${n}`);let i=new Float64Array(r.size),s=e[a],u=[],c=1;for(let f=n-1;f>=0;f--)u.unshift(c),c*=e[f];let l=r.size,y=u[a];for(let f=0;f<l;f++)Math.floor(f/y)%s===0?i[f]=Number(o[f]):i[f]=i[f-y]*Number(o[f]);return D.fromData(i,[...e],"float64")}function ct(r,t,e=!1){let n=xr(r,t,e),o=Ir(r,t,e);if(typeof n=="number"&&typeof o=="number")return n-o;let a=n,i=o,s=a.data,u=i.data,c=new Float64Array(a.size);for(let l=0;l<a.size;l++)c[l]=Number(s[l])-Number(u[l]);return D.fromData(c,[...a.shape],"float64")}function ft(r,t,e=!1){return nr(r,.5,t,e)}function lt(r,t,e,n=!1){return nr(r,t/100,e,n)}function nr(r,t,e,n=!1){if(t<0||t>1)throw new Error("Quantile must be between 0 and 1");let o=r.shape,a=o.length,i=r.data;if(e===void 0){let p=[];for(let b=0;b<r.size;b++)p.push(Number(i[b]));p.sort((b,S)=>b-S);let m=p.length,d=t*(m-1),g=Math.floor(d),A=Math.ceil(d);if(g===A)return p[g];let h=d-g;return p[g]*(1-h)+p[A]*h}let s=e;if(s<0&&(s=a+s),s<0||s>=a)throw new Error(`axis ${e} is out of bounds for array of dimension ${a}`);let u=Array.from(o).filter((p,m)=>m!==s);if(u.length===0)return nr(r,t);let c=u.reduce((p,m)=>p*m,1),l=o[s],y=new Float64Array(c);for(let p=0;p<c;p++){let m=[];for(let b=0;b<l;b++){let S=T(p,s,b,o),w=U(S,o);m.push(Number(i[w]))}m.sort((b,S)=>b-S);let d=m.length,g=t*(d-1),A=Math.floor(g),h=Math.ceil(g);if(A===h)y[p]=m[A];else{let b=g-A;y[p]=m[A]*(1-b)+m[h]*b}}let f=D.fromData(y,u,"float64");if(n){let p=[...o];return p[s]=1,D.fromData(y,p,"float64")}return f}function zr(r,t,e,n=!1){let o=r.shape,a=o.length,i=r.data;if(e===void 0)return wr(r,t,n);if(t===void 0){let m=0,d=0,g=e.data;for(let A=0;A<r.size;A++){let h=Number(g[A%e.size]);m+=Number(i[A])*h,d+=h}return d===0?NaN:m/d}let s=t;if(s<0&&(s=a+s),s<0||s>=a)throw new Error(`axis ${t} is out of bounds for array of dimension ${a}`);let u=Array.from(o).filter((m,d)=>d!==s);if(u.length===0)return zr(r,void 0,e);let c=u.reduce((m,d)=>m*d,1),l=o[s],y=e.data,f=new Float64Array(c);for(let m=0;m<c;m++){let d=0,g=0;for(let A=0;A<l;A++){let h=T(m,s,A,o),b=U(h,o),S=Number(y[A%e.size]);d+=Number(i[b])*S,g+=S}f[m]=g===0?NaN:d/g}let p=D.fromData(f,u,"float64");if(n){let m=[...o];return m[s]=1,D.fromData(f,m,"float64")}return p}function Br(r,t,e=!1){let n=r.shape,o=n.length,a=r.data;if(t===void 0){let f=0;for(let p=0;p<r.size;p++){let m=Number(a[p]);isNaN(m)||(f+=m)}return f}let i=t;if(i<0&&(i=o+i),i<0||i>=o)throw new Error(`axis ${t} is out of bounds for array of dimension ${o}`);let s=Array.from(n).filter((f,p)=>p!==i);if(s.length===0)return Br(r);let u=s.reduce((f,p)=>f*p,1),c=n[i],l=new Float64Array(u);for(let f=0;f<u;f++){let p=0;for(let m=0;m<c;m++){let d=T(f,i,m,n),g=U(d,n),A=Number(a[g]);isNaN(A)||(p+=A)}l[f]=p}let y=D.fromData(l,s,"float64");if(e){let f=[...n];return f[i]=1,D.fromData(l,f,"float64")}return y}function _r(r,t,e=!1){let n=r.shape,o=n.length,a=r.data;if(t===void 0){let f=1;for(let p=0;p<r.size;p++){let m=Number(a[p]);isNaN(m)||(f*=m)}return f}let i=t;if(i<0&&(i=o+i),i<0||i>=o)throw new Error(`axis ${t} is out of bounds for array of dimension ${o}`);let s=Array.from(n).filter((f,p)=>p!==i);if(s.length===0)return _r(r);let u=s.reduce((f,p)=>f*p,1),c=n[i],l=new Float64Array(u);for(let f=0;f<u;f++){let p=1;for(let m=0;m<c;m++){let d=T(f,i,m,n),g=U(d,n),A=Number(a[g]);isNaN(A)||(p*=A)}l[f]=p}let y=D.fromData(l,s,"float64");if(e){let f=[...n];return f[i]=1,D.fromData(l,f,"float64")}return y}function Er(r,t,e=!1){let n=r.shape,o=n.length,a=r.data;if(t===void 0){let f=0,p=0;for(let m=0;m<r.size;m++){let d=Number(a[m]);isNaN(d)||(f+=d,p++)}return p===0?NaN:f/p}let i=t;if(i<0&&(i=o+i),i<0||i>=o)throw new Error(`axis ${t} is out of bounds for array of dimension ${o}`);let s=Array.from(n).filter((f,p)=>p!==i);if(s.length===0)return Er(r);let u=s.reduce((f,p)=>f*p,1),c=n[i],l=new Float64Array(u);for(let f=0;f<u;f++){let p=0,m=0;for(let d=0;d<c;d++){let g=T(f,i,d,n),A=U(g,n),h=Number(a[A]);isNaN(h)||(p+=h,m++)}l[f]=m===0?NaN:p/m}let y=D.fromData(l,s,"float64");if(e){let f=[...n];return f[i]=1,D.fromData(l,f,"float64")}return y}function cr(r,t,e=0,n=!1){let o=r.shape,a=o.length,i=r.data;if(t===void 0){let p=0,m=0;for(let A=0;A<r.size;A++){let h=Number(i[A]);isNaN(h)||(p+=h,m++)}if(m-e<=0)return NaN;let d=p/m,g=0;for(let A=0;A<r.size;A++){let h=Number(i[A]);isNaN(h)||(g+=(h-d)**2)}return g/(m-e)}let s=t;if(s<0&&(s=a+s),s<0||s>=a)throw new Error(`axis ${t} is out of bounds for array of dimension ${a}`);let u=Array.from(o).filter((p,m)=>m!==s);if(u.length===0)return cr(r,void 0,e);let c=u.reduce((p,m)=>p*m,1),l=o[s],y=new Float64Array(c);for(let p=0;p<c;p++){let m=0,d=0;for(let h=0;h<l;h++){let b=T(p,s,h,o),S=U(b,o),w=Number(i[S]);isNaN(w)||(m+=w,d++)}if(d-e<=0){y[p]=NaN;continue}let g=m/d,A=0;for(let h=0;h<l;h++){let b=T(p,s,h,o),S=U(b,o),w=Number(i[S]);isNaN(w)||(A+=(w-g)**2)}y[p]=A/(d-e)}let f=D.fromData(y,u,"float64");if(n){let p=[...o];return p[s]=1,D.fromData(y,p,"float64")}return f}function pt(r,t,e=0,n=!1){let o=cr(r,t,e,n);if(typeof o=="number")return Math.sqrt(o);let a=o,i=new Float64Array(a.size);for(let s=0;s<a.size;s++)i[s]=Math.sqrt(Number(a.data[s]));return D.fromData(i,[...a.shape],"float64")}function Ur(r,t,e=!1){let n=r.shape,o=n.length,a=r.data;if(t===void 0){let f=1/0;for(let p=0;p<r.size;p++){let m=Number(a[p]);!isNaN(m)&&m<f&&(f=m)}return f===1/0?NaN:f}let i=t;if(i<0&&(i=o+i),i<0||i>=o)throw new Error(`axis ${t} is out of bounds for array of dimension ${o}`);let s=Array.from(n).filter((f,p)=>p!==i);if(s.length===0)return Ur(r);let u=s.reduce((f,p)=>f*p,1),c=n[i],l=new Float64Array(u);for(let f=0;f<u;f++){let p=1/0;for(let m=0;m<c;m++){let d=T(f,i,m,n),g=U(d,n),A=Number(a[g]);!isNaN(A)&&A<p&&(p=A)}l[f]=p===1/0?NaN:p}let y=D.fromData(l,s,"float64");if(e){let f=[...n];return f[i]=1,D.fromData(l,f,"float64")}return y}function Tr(r,t,e=!1){let n=r.shape,o=n.length,a=r.data;if(t===void 0){let f=-1/0;for(let p=0;p<r.size;p++){let m=Number(a[p]);!isNaN(m)&&m>f&&(f=m)}return f===-1/0?NaN:f}let i=t;if(i<0&&(i=o+i),i<0||i>=o)throw new Error(`axis ${t} is out of bounds for array of dimension ${o}`);let s=Array.from(n).filter((f,p)=>p!==i);if(s.length===0)return Tr(r);let u=s.reduce((f,p)=>f*p,1),c=n[i],l=new Float64Array(u);for(let f=0;f<u;f++){let p=-1/0;for(let m=0;m<c;m++){let d=T(f,i,m,n),g=U(d,n),A=Number(a[g]);!isNaN(A)&&A>p&&(p=A)}l[f]=p===-1/0?NaN:p}let y=D.fromData(l,s,"float64");if(e){let f=[...n];return f[i]=1,D.fromData(l,f,"float64")}return y}function Or(r,t){let e=r.shape,n=e.length,o=r.data;if(t===void 0){let l=1/0,y=-1;for(let f=0;f<r.size;f++){let p=Number(o[f]);!isNaN(p)&&p<l&&(l=p,y=f)}return y}let a=t;if(a<0&&(a=n+a),a<0||a>=n)throw new Error(`axis ${t} is out of bounds for array of dimension ${n}`);let i=Array.from(e).filter((l,y)=>y!==a);if(i.length===0)return Or(r);let s=i.reduce((l,y)=>l*y,1),u=e[a],c=new Int32Array(s);for(let l=0;l<s;l++){let y=1/0,f=0;for(let p=0;p<u;p++){let m=T(l,a,p,e),d=U(m,e),g=Number(o[d]);!isNaN(g)&&g<y&&(y=g,f=p)}c[l]=f}return D.fromData(c,i,"int32")}function $r(r,t){let e=r.shape,n=e.length,o=r.data;if(t===void 0){let l=-1/0,y=-1;for(let f=0;f<r.size;f++){let p=Number(o[f]);!isNaN(p)&&p>l&&(l=p,y=f)}return y}let a=t;if(a<0&&(a=n+a),a<0||a>=n)throw new Error(`axis ${t} is out of bounds for array of dimension ${n}`);let i=Array.from(e).filter((l,y)=>y!==a);if(i.length===0)return $r(r);let s=i.reduce((l,y)=>l*y,1),u=e[a],c=new Int32Array(s);for(let l=0;l<s;l++){let y=-1/0,f=0;for(let p=0;p<u;p++){let m=T(l,a,p,e),d=U(m,e),g=Number(o[d]);!isNaN(g)&&g>y&&(y=g,f=p)}c[l]=f}return D.fromData(c,i,"int32")}function yt(r,t){let e=r.shape,n=e.length,o=r.data;if(t===void 0){let f=r.size,p=new Float64Array(f),m=0;for(let d=0;d<f;d++){let g=Number(o[d]);isNaN(g)||(m+=g),p[d]=m}return D.fromData(p,[f],"float64")}let a=t;if(a<0&&(a=n+a),a<0||a>=n)throw new Error(`axis ${t} is out of bounds for array of dimension ${n}`);let i=new Float64Array(r.size),s=e[a],u=[],c=1;for(let f=n-1;f>=0;f--)u.unshift(c),c*=e[f];let l=r.size,y=u[a];for(let f=0;f<l;f++){let p=Number(o[f]);Math.floor(f/y)%s===0?i[f]=isNaN(p)?0:p:i[f]=i[f-y]+(isNaN(p)?0:p)}return D.fromData(i,[...e],"float64")}function mt(r,t){let e=r.shape,n=e.length,o=r.data;if(t===void 0){let f=r.size,p=new Float64Array(f),m=1;for(let d=0;d<f;d++){let g=Number(o[d]);isNaN(g)||(m*=g),p[d]=m}return D.fromData(p,[f],"float64")}let a=t;if(a<0&&(a=n+a),a<0||a>=n)throw new Error(`axis ${t} is out of bounds for array of dimension ${n}`);let i=new Float64Array(r.size),s=e[a],u=[],c=1;for(let f=n-1;f>=0;f--)u.unshift(c),c*=e[f];let l=r.size,y=u[a];for(let f=0;f<l;f++){let p=Number(o[f]);Math.floor(f/y)%s===0?i[f]=isNaN(p)?1:p:i[f]=i[f-y]*(isNaN(p)?1:p)}return D.fromData(i,[...e],"float64")}function vr(r,t,e=!1){let n=r.shape,o=n.length,a=r.data;if(t===void 0){let f=[];for(let d=0;d<r.size;d++){let g=Number(a[d]);isNaN(g)||f.push(g)}if(f.length===0)return NaN;f.sort((d,g)=>d-g);let p=f.length,m=Math.floor(p/2);return p%2===0?(f[m-1]+f[m])/2:f[m]}let i=t;if(i<0&&(i=o+i),i<0||i>=o)throw new Error(`axis ${t} is out of bounds for array of dimension ${o}`);let s=Array.from(n).filter((f,p)=>p!==i);if(s.length===0)return vr(r);let u=s.reduce((f,p)=>f*p,1),c=n[i],l=new Float64Array(u);for(let f=0;f<u;f++){let p=[];for(let g=0;g<c;g++){let A=T(f,i,g,n),h=U(A,n),b=Number(a[h]);isNaN(b)||p.push(b)}if(p.length===0){l[f]=NaN;continue}p.sort((g,A)=>g-A);let m=p.length,d=Math.floor(m/2);m%2===0?l[f]=(p[d-1]+p[d])/2:l[f]=p[d]}let y=D.fromData(l,s,"float64");if(e){let f=[...n];return f[i]=1,D.fromData(l,f,"float64")}return y}var me=$(()=>{"use strict";F();V();pe()});function R(r,t){let e=r.size,n=r.dtype,o=t.indexOf(-1),a;if(o!==-1){let c=t.reduce((y,f,p)=>p===o?y:y*f,1),l=e/c;if(!Number.isInteger(l))throw new Error(`cannot reshape array of size ${e} into shape ${JSON.stringify(t)}`);a=t.map((y,f)=>f===o?l:y)}else a=t;if(a.reduce((c,l)=>c*l,1)!==e)throw new Error(`cannot reshape array of size ${e} into shape ${JSON.stringify(a)}`);if(r.isCContiguous){let c=r.data;return D.fromData(c,a,n,Z(a),0)}let u=r.copy().data;return D.fromData(u,a,n,Z(a),0)}function Mr(r){let t=r.size,e=r.dtype,n=_(e);if(!n)throw new Error(`Cannot flatten array with dtype ${e}`);if(r.isCContiguous){let s=r.data.slice(r.offset,r.offset+t);return D.fromData(s,[t],e,[1],0)}let o=new n(t),a=N(e);for(let i=0;i<t;i++){let s=r.iget(i);o[i]=s}return D.fromData(o,[t],e,[1],0)}function lr(r){let t=r.size,e=r.dtype;if(r.isCContiguous){let n=r.data;return D.fromData(n,[t],e,[1],0)}return Mr(r)}function pr(r,t){let e=r.shape,n=e.length,o=r.strides,a=r.data,i=r.dtype,s;if(t===void 0)s=Array.from({length:n},(y,f)=>n-1-f);else{if(t.length!==n)throw new Error(`axes must have length ${n}, got ${t.length}`);let y=new Set;for(let f of t){let p=f<0?n+f:f;if(p<0||p>=n)throw new Error(`axis ${f} is out of bounds for array of dimension ${n}`);if(y.has(p))throw new Error("repeated axis in transpose");y.add(p)}s=t.map(f=>f<0?n+f:f)}let u=s.map(y=>e[y]),c=Array.from(o),l=s.map(y=>c[y]);return D.fromData(a,u,i,l,r.offset)}function de(r,t){let e=r.shape,n=e.length,o=r.strides,a=r.data,i=r.dtype;if(t===void 0){let s=[],u=[];for(let c=0;c<n;c++)e[c]!==1&&(s.push(e[c]),u.push(o[c]));return s.length===0&&(s.push(1),u.push(1)),D.fromData(a,s,i,u,r.offset)}else{let s=t<0?n+t:t;if(s<0||s>=n)throw new Error(`axis ${t} is out of bounds for array of dimension ${n}`);if(e[s]!==1)throw new Error(`cannot select an axis which has size not equal to one (axis ${t} has size ${e[s]})`);let u=[],c=[];for(let l=0;l<n;l++)l!==s&&(u.push(e[l]),c.push(o[l]));return D.fromData(a,u,i,c,r.offset)}}function fr(r,t){let e=r.shape,n=e.length,o=r.strides,a=r.data,i=r.dtype,s=t;if(s<0&&(s=n+t+1),s<0||s>n)throw new Error(`axis ${t} is out of bounds for array of dimension ${n+1}`);let u=[...Array.from(e)];u.splice(s,0,1);let c=[...Array.from(o)],l=s<n?o[s]*(e[s]||1):1;return c.splice(s,0,l),D.fromData(a,u,i,c,r.offset)}function ge(r,t,e){let n=r.shape,o=n.length,a=r.strides,i=r.data,s=r.dtype,u=t<0?o+t:t,c=e<0?o+e:e;if(u<0||u>=o)throw new Error(`axis1 ${t} is out of bounds for array of dimension ${o}`);if(c<0||c>=o)throw new Error(`axis2 ${e} is out of bounds for array of dimension ${o}`);if(u===c)return D.fromData(i,Array.from(n),s,Array.from(a),r.offset);let l=Array.from(n),y=Array.from(a);return[l[u],l[c]]=[l[c],l[u]],[y[u],y[c]]=[y[c],y[u]],D.fromData(i,l,s,y,r.offset)}function dt(r,t,e){let n=r.ndim,o=Array.isArray(t)?t:[t],a=Array.isArray(e)?e:[e];if(o.length!==a.length)throw new Error("source and destination must have the same number of elements");let i=o.map(c=>{let l=c<0?n+c:c;if(l<0||l>=n)throw new Error(`source axis ${c} is out of bounds for array of dimension ${n}`);return l}),s=a.map(c=>{let l=c<0?n+c:c;if(l<0||l>=n)throw new Error(`destination axis ${c} is out of bounds for array of dimension ${n}`);return l});if(new Set(i).size!==i.length)throw new Error("repeated axis in source");if(new Set(s).size!==s.length)throw new Error("repeated axis in destination");let u=[];for(let c=0;c<n;c++)i.includes(c)||u.push(c);for(let c=0;c<i.length;c++){let l=s[c];u.splice(l,0,i[c])}return pr(r,u)}function Q(r,t=0){if(r.length===0)throw new Error("need at least one array to concatenate");if(r.length===1)return r[0].copy();let e=r[0],n=e.ndim,o=e.dtype,a=t<0?n+t:t;if(a<0||a>=n)throw new Error(`axis ${t} is out of bounds for array of dimension ${n}`);for(let p=1;p<r.length;p++){let m=r[p];if(m.ndim!==n)throw new Error("all the input arrays must have same number of dimensions");for(let d=0;d<n;d++)if(d!==a&&m.shape[d]!==e.shape[d])throw new Error("all the input array dimensions except for the concatenation axis must match exactly")}let i=Array.from(e.shape),s=e.shape[a];for(let p=1;p<r.length;p++)s+=r[p].shape[a];i[a]=s;let u=i.reduce((p,m)=>p*m,1),c=_(o);if(!c)throw new Error(`Cannot concatenate arrays with dtype ${o}`);let l=new c(u),y=Z(i),f=0;for(let p of r){let m=p.shape[a];Ts(p,l,i,y,a,f,o),f+=m}return D.fromData(l,i,o)}function Ts(r,t,e,n,o,a,i){let s=r.shape,u=s.length,c=r.size,l=i==="int64"||i==="uint64";if(o===0&&r.isCContiguous&&u>0){let p=a*n[0],m=r.data,d=r.offset,g=d+c;t.set(m.subarray(d,g),p);return}if(o===1&&u===2&&r.isCContiguous){let p=s[0],m=s[1],d=e[1],g=r.data,A=r.offset;for(let h=0;h<p;h++){let b=A+h*m,S=h*d+a;t.set(g.subarray(b,b+m),S)}return}let y=new Array(u).fill(0),f=a*n[o];for(let p=0;p<c;p++){let m=r.iget(p),d=f;for(let g=0;g<u;g++)d+=y[g]*n[g];t[d]=m;for(let g=u-1;g>=0&&(y[g]++,!(y[g]<s[g]));g--)y[g]=0}}function he(r,t=0){if(r.length===0)throw new Error("need at least one array to stack");let e=r[0],n=e.shape,o=e.ndim,a=t<0?o+1+t:t;if(a<0||a>o)throw new Error(`axis ${t} is out of bounds for array of dimension ${o+1}`);for(let s=1;s<r.length;s++){let u=r[s];if(u.ndim!==o)throw new Error("all input arrays must have the same shape");for(let c=0;c<o;c++)if(u.shape[c]!==n[c])throw new Error("all input arrays must have the same shape")}let i=r.map(s=>fr(s,a));return Q(i,a)}function Ae(r){if(r.length===0)throw new Error("need at least one array to stack");let t=r.map(e=>e.ndim===1?R(e,[1,e.shape[0]]):e);return Q(t,0)}function gt(r){if(r.length===0)throw new Error("need at least one array to stack");return r.every(e=>e.ndim===1)?Q(r,0):Q(r,1)}function be(r){if(r.length===0)throw new Error("need at least one array to stack");let t=r.map(e=>e.ndim===1?R(fr(R(e,[1,e.shape[0]]),2),[1,e.shape[0],1]):e.ndim===2?fr(e,2):e);return Q(t,2)}function Se(r,t,e=0){let n=r.shape,o=n.length,a=e<0?o+e:e;if(a<0||a>=o)throw new Error(`axis ${e} is out of bounds for array of dimension ${o}`);let i=n[a],s;if(typeof t=="number"){if(i%t!==0)throw new Error("array split does not result in an equal division");let u=i/t;s=[];for(let c=1;c<t;c++)s.push(c*u)}else s=t;return De(r,s,a)}function yr(r,t,e=0){let n=r.shape,o=n.length,a=e<0?o+e:e;if(a<0||a>=o)throw new Error(`axis ${e} is out of bounds for array of dimension ${o}`);let i=n[a],s;if(typeof t=="number"){let u=t,c=Math.floor(i/u),l=i%u;s=[];let y=0;for(let f=0;f<u-1;f++)y+=c+(f<l?1:0),s.push(y)}else s=t;return De(r,s,a)}function De(r,t,e){let n=r.shape,o=n[e],a=[0,...t,o],i=[];for(let s=0;s<a.length-1;s++){let u=a[s],c=a[s+1];if(u>c)throw new Error("split indices must be in ascending order");let l=Array.from(n);l[e]=c-u;let y=r.offset+u*r.strides[e];i.push(D.fromData(r.data,l,r.dtype,Array.from(r.strides),y))}return i}function Ne(r,t){if(r.ndim<2)throw new Error("vsplit only works on arrays of 2 or more dimensions");return yr(r,t,0)}function we(r,t){if(r.ndim<1)throw new Error("hsplit only works on arrays of 1 or more dimensions");let e=r.ndim===1?0:1;return yr(r,t,e)}function xe(r,t){let e=r.shape,n=e.length,o=r.dtype,a=Array.isArray(t)?t:[t],i=Math.max(n,a.length),s=new Array(i).fill(1),u=new Array(i).fill(1);for(let h=0;h<n;h++)s[i-n+h]=e[h];for(let h=0;h<a.length;h++)u[i-a.length+h]=a[h];let c=s.map((h,b)=>h*u[b]),l=c.reduce((h,b)=>h*b,1),y=_(o);if(!y)throw new Error(`Cannot tile array with dtype ${o}`);let f=new y(l),p=Z(c),m=r;n<i&&(m=R(r,s));let d=o==="int64"||o==="uint64",g=m.strides,A=new Array(i).fill(0);for(let h=0;h<l;h++){let b=m.offset;for(let x=0;x<i;x++){let z=A[x]%s[x];b+=z*g[x]}let S=m.data[b],w=0;for(let x=0;x<i;x++)w+=A[x]*p[x];f[w]=S;for(let x=i-1;x>=0&&(A[x]++,!(A[x]<c[x]));x--)A[x]=0}return D.fromData(f,c,o)}function Ie(r,t,e){let n=r.shape,o=n.length,a=r.dtype,i=r.size;if(e===void 0){let h=i,b=Array.isArray(t)?t:new Array(h).fill(t);if(b.length!==h)throw new Error(`operands could not be broadcast together with shape (${h},) (${b.length},)`);let S=b.reduce((B,E)=>B+E,0),w=_(a);if(!w)throw new Error(`Cannot repeat array with dtype ${a}`);let x=new w(S),z=0;for(let B=0;B<h;B++){let E=r.iget(B),O=b[B];for(let Sr=0;Sr<O;Sr++)x[z++]=E}return D.fromData(x,[S],a)}let s=e<0?o+e:e;if(s<0||s>=o)throw new Error(`axis ${e} is out of bounds for array of dimension ${o}`);let u=n[s],c=Array.isArray(t)?t:new Array(u).fill(t);if(c.length!==u)throw new Error(`operands could not be broadcast together with shape (${u},) (${c.length},)`);let l=Array.from(n);l[s]=c.reduce((h,b)=>h+b,0);let y=l.reduce((h,b)=>h*b,1),f=_(a);if(!f)throw new Error(`Cannot repeat array with dtype ${a}`);let p=new f(y),m=Z(l),d=new Array(o).fill(0),g=a==="int64"||a==="uint64",A=[0];for(let h=0;h<u;h++)A.push(A[h]+c[h]);for(let h=0;h<i;h++){let b=r.iget(h),S=d[s],w=c[S],x=0;for(let E=0;E<o;E++)E!==s&&(x+=d[E]*m[E]);let z=m[s],B=A[S];for(let E=0;E<w;E++){let O=x+(B+E)*z;p[O]=b}for(let E=o-1;E>=0&&(d[E]++,!(d[E]<n[E]));E--)d[E]=0}return D.fromData(p,l,a)}function ze(r,t){let e=r.shape,n=e.length,o=r.dtype,a=r.size,i;if(t===void 0)i=new Set(Array.from({length:n},(f,p)=>p));else if(typeof t=="number"){let f=t<0?n+t:t;if(f<0||f>=n)throw new Error(`axis ${t} is out of bounds for array of dimension ${n}`);i=new Set([f])}else i=new Set(t.map(f=>{let p=f<0?n+f:f;if(p<0||p>=n)throw new Error(`axis ${f} is out of bounds for array of dimension ${n}`);return p}));let s=_(o);if(!s)throw new Error(`Cannot flip array with dtype ${o}`);let u=new s(a),c=N(o);if(n===1&&r.isCContiguous){let f=r.data,p=r.offset;for(let m=0;m<a;m++)u[m]=f[p+a-1-m];return D.fromData(u,[...e],o)}if(n===2&&r.isCContiguous){let f=e[0],p=e[1],m=r.data,d=r.offset;if(i.size===2){for(let g=0;g<a;g++)u[g]=m[d+a-1-g];return D.fromData(u,[...e],o)}if(i.size===1){if(i.has(0)){for(let g=0;g<f;g++){let A=d+(f-1-g)*p,h=g*p;for(let b=0;b<p;b++)u[h+b]=m[A+b]}return D.fromData(u,[...e],o)}else if(i.has(1)){for(let g=0;g<f;g++){let A=d+g*p,h=g*p;for(let b=0;b<p;b++)u[h+b]=m[A+p-1-b]}return D.fromData(u,[...e],o)}}}let l=new Array(n),y=new Array(n).fill(0);for(let f=0;f<a;f++){for(let d=0;d<n;d++)l[d]=i.has(d)?e[d]-1-y[d]:y[d];let p=r.offset;for(let d=0;d<n;d++)p+=l[d]*r.strides[d];let m=r.data[p];u[f]=m;for(let d=n-1;d>=0&&(y[d]++,!(y[d]<e[d]));d--)y[d]=0}return D.fromData(u,[...e],o)}function Be(r,t=1,e=[0,1]){let n=r.shape,o=n.length,a=r.dtype;if(o<2)throw new Error("Input must be at least 2-D");let i=e[0]<0?o+e[0]:e[0],s=e[1]<0?o+e[1]:e[1];if(i<0||i>=o||s<0||s>=o)throw new Error(`Axes are out of bounds for array of dimension ${o}`);if(i===s)throw new Error("Axes must be different");if(t=(t%4+4)%4,t===0)return r.copy();let u=_(a);if(!u)throw new Error(`Cannot rotate array with dtype ${a}`);let c=[...n];(t===1||t===3)&&([c[i],c[s]]=[c[s],c[i]]);let l=c.reduce((g,A)=>g*A,1),y=new u(l),f=Z(c),p=N(a),m=new Array(o).fill(0),d=new Array(o);for(let g=0;g<r.size;g++){for(let w=0;w<o;w++)d[w]=m[w];let A,h;t===1?(A=n[s]-1-m[s],h=m[i]):t===2?(A=n[i]-1-m[i],h=n[s]-1-m[s],d[i]=A,d[s]=h):(A=m[s],h=n[i]-1-m[i]),t!==2&&(d[i]=A,d[s]=h);let b=0;for(let w=0;w<o;w++)b+=d[w]*f[w];let S=r.iget(g);y[b]=S;for(let w=o-1;w>=0&&(m[w]++,!(m[w]<n[w]));w--)m[w]=0}return D.fromData(y,c,a)}function _e(r,t,e){let n=r.shape,o=n.length,a=r.dtype,i=r.size;if(e===void 0){let m=Array.isArray(t)?t.reduce((b,S)=>b+S,0):t,d=Mr(r),g=_(a);if(!g)throw new Error(`Cannot roll array with dtype ${a}`);let A=new g(i),h=N(a);for(let b=0;b<i;b++){let S=((b-m)%i+i)%i,w=d.iget(S);A[b]=w}return D.fromData(A,[...n],a)}let s=Array.isArray(t)?t:[t],u=Array.isArray(e)?e:[e];if(s.length!==u.length)throw new Error("shift and axis must have the same length");let c=u.map(m=>{let d=m<0?o+m:m;if(d<0||d>=o)throw new Error(`axis ${m} is out of bounds for array of dimension ${o}`);return d}),l=_(a);if(!l)throw new Error(`Cannot roll array with dtype ${a}`);let y=new l(i),f=N(a),p=new Array(o).fill(0);for(let m=0;m<i;m++){let d=[...p];for(let h=0;h<c.length;h++){let b=c[h],S=n[b],w=s[h];d[b]=((d[b]-w)%S+S)%S}let g=r.offset;for(let h=0;h<o;h++)g+=d[h]*r.strides[h];let A=r.data[g];y[m]=A;for(let h=o-1;h>=0&&(p[h]++,!(p[h]<n[h]));h--)p[h]=0}return D.fromData(y,[...n],a)}function Ee(r,t,e=0){let n=r.ndim,o=t<0?n+t:t;if(o<0||o>=n)throw new Error(`axis ${t} is out of bounds for array of dimension ${n}`);let a=e<0?n+e:e;if(a<0||a>n)throw new Error(`start ${e} is out of bounds`);return o<a&&a--,o===a?D.fromData(r.data,Array.from(r.shape),r.dtype,Array.from(r.strides),r.offset):dt(r,o,a)}function Ue(r,t){if(r.ndim<3)throw new Error("dsplit only works on arrays of 3 or more dimensions");return yr(r,t,2)}function Te(r){if(r.length===0)throw new Error("need at least one array to stack");let t=r.map(e=>e.ndim===1?R(e,[e.shape[0],1]):e);return gt(t)}function Oe(r,t){let e=r.dtype,n=t.reduce((u,c)=>u*c,1),o=r.size,a=_(e);if(!a)throw new Error(`Cannot resize array with dtype ${e}`);let i=new a(n),s=N(e);for(let u=0;u<n;u++){let c=u%o,l=r.iget(c);i[u]=l}return D.fromData(i,t,e)}function $e(r){return r.map(t=>t.ndim===0?R(t,[1]):t)}function ve(r){return r.map(t=>t.ndim===0?R(t,[1,1]):t.ndim===1?R(t,[1,t.shape[0]]):t)}function Me(r){return r.map(t=>t.ndim===0?R(t,[1,1,1]):t.ndim===1?R(t,[1,t.shape[0],1]):t.ndim===2?R(t,[t.shape[0],t.shape[1],1]):t)}var ht=$(()=>{"use strict";F();V()});function Os(r,t,e,n,o,a,i,s,u,c,l,y,f,p){if(y===0)for(let A=0;A<n*o;A++)f[A]=0;else if(y!==1)for(let A=0;A<n*o;A++)f[A]=(f[A]??0)*y;let m=r==="row-major",d=t==="transpose",g=e==="transpose";if(m&&!d&&!g)for(let A=0;A<n;A++)for(let h=0;h<o;h++){let b=0;for(let S=0;S<a;S++)b+=(s[A*u+S]??0)*(c[S*l+h]??0);f[A*p+h]=(f[A*p+h]??0)+i*b}else if(m&&d&&!g)for(let A=0;A<n;A++)for(let h=0;h<o;h++){let b=0;for(let S=0;S<a;S++)b+=(s[S*u+A]??0)*(c[S*l+h]??0);f[A*p+h]=(f[A*p+h]??0)+i*b}else if(m&&!d&&g)for(let A=0;A<n;A++)for(let h=0;h<o;h++){let b=0;for(let S=0;S<a;S++)b+=(s[A*u+S]??0)*(c[h*l+S]??0);f[A*p+h]=(f[A*p+h]??0)+i*b}else if(m&&d&&g)for(let A=0;A<n;A++)for(let h=0;h<o;h++){let b=0;for(let S=0;S<a;S++)b+=(s[S*u+A]??0)*(c[h*l+S]??0);f[A*p+h]=(f[A*p+h]??0)+i*b}else if(!m&&!d&&!g)for(let A=0;A<n;A++)for(let h=0;h<o;h++){let b=0;for(let S=0;S<a;S++)b+=(s[S*u+A]??0)*(c[h*l+S]??0);f[h*p+A]=(f[h*p+A]??0)+i*b}else if(!m&&d&&!g)for(let A=0;A<n;A++)for(let h=0;h<o;h++){let b=0;for(let S=0;S<a;S++)b+=(s[A*u+S]??0)*(c[h*l+S]??0);f[h*p+A]=(f[h*p+A]??0)+i*b}else if(!m&&!d&&g)for(let A=0;A<n;A++)for(let h=0;h<o;h++){let b=0;for(let S=0;S<a;S++)b+=(s[S*u+A]??0)*(c[S*l+h]??0);f[h*p+A]=(f[h*p+A]??0)+i*b}else for(let A=0;A<n;A++)for(let h=0;h<o;h++){let b=0;for(let S=0;S<a;S++)b+=(s[A*u+S]??0)*(c[S*l+h]??0);f[h*p+A]=(f[h*p+A]??0)+i*b}}function bt(r,t){let e=r.ndim,n=t.ndim;if(e===0||n===0){let o=e===0?r.get():null,a=n===0?t.get():null;if(e===0&&n===0)return typeof o=="bigint"&&typeof a=="bigint"?o*a:Number(o)*Number(a);if(e===0){let i=M(r.dtype,t.dtype),s=D.zeros([...t.shape],i);for(let u=0;u<t.size;u++){let c=t.data[u+t.offset];typeof o=="bigint"&&typeof c=="bigint"?s.data[u]=o*c:s.data[u]=Number(o)*Number(c)}return s}else{let i=M(r.dtype,t.dtype),s=D.zeros([...r.shape],i);for(let u=0;u<r.size;u++){let c=r.data[u+r.offset];typeof c=="bigint"&&typeof a=="bigint"?s.data[u]=c*a:s.data[u]=Number(c)*Number(a)}return s}}if(e===1&&n===1){if(r.shape[0]!==t.shape[0])throw new Error(`dot: incompatible shapes (${r.shape[0]},) and (${t.shape[0]},)`);let o=r.shape[0],a=0;for(let i=0;i<o;i++){let s=r.get(i),u=t.get(i);typeof s=="bigint"&&typeof u=="bigint"?a=Number(a)+Number(s*u):a+=Number(s)*Number(u)}return a}if(e===2&&n===2)return X(r,t);if(e===2&&n===1){let[o,a]=r.shape,i=t.shape[0];if(a!==i)throw new Error(`dot: incompatible shapes (${o},${a}) and (${i},)`);let s=M(r.dtype,t.dtype),u=D.zeros([o],s);for(let c=0;c<o;c++){let l=0;for(let y=0;y<a;y++){let f=r.get(c,y),p=t.get(y);typeof f=="bigint"&&typeof p=="bigint"?l=Number(l)+Number(f*p):l+=Number(f)*Number(p)}u.set([c],l)}return u}if(e===1&&n===2){let o=r.shape[0],[a,i]=t.shape;if(o!==a)throw new Error(`dot: incompatible shapes (${o},) and (${a},${i})`);let s=M(r.dtype,t.dtype),u=D.zeros([i],s);for(let c=0;c<i;c++){let l=0;for(let y=0;y<o;y++){let f=r.get(y),p=t.get(y,c);typeof f=="bigint"&&typeof p=="bigint"?l=Number(l)+Number(f*p):l+=Number(f)*Number(p)}u.set([c],l)}return u}if(e>2&&n===1){let o=r.shape[e-1],a=t.shape[0];if(o!==a)throw new Error(`dot: incompatible shapes ${JSON.stringify(r.shape)} and (${a},)`);let i=[...r.shape.slice(0,-1)],s=M(r.dtype,t.dtype),u=D.zeros(i,s),c=i.reduce((l,y)=>l*y,1);for(let l=0;l<c;l++){let y=0,f=l,p=[];for(let m=i.length-1;m>=0;m--)p[m]=f%i[m],f=Math.floor(f/i[m]);for(let m=0;m<o;m++){let d=[...p,m],g=r.get(...d),A=t.get(m);typeof g=="bigint"&&typeof A=="bigint"?y=Number(y)+Number(g*A):y+=Number(g)*Number(A)}u.set(p,y)}return u}if(e===1&&n>2){let o=r.shape[0],a=1,i=t.shape[a];if(o!==i)throw new Error(`dot: incompatible shapes (${o},) and ${JSON.stringify(t.shape)}`);let s=[...t.shape.slice(0,a),...t.shape.slice(a+1)],u=M(r.dtype,t.dtype),c=D.zeros(s,u),l=s.reduce((y,f)=>y*f,1);for(let y=0;y<l;y++){let f=y,p=[];for(let A=s.length-1;A>=0;A--)p[A]=f%s[A],f=Math.floor(f/s[A]);let m=p.slice(0,a),d=p.slice(a),g=0;for(let A=0;A<o;A++){let h=r.get(A),b=[...m,A,...d],S=t.get(...b);typeof h=="bigint"&&typeof S=="bigint"?g=Number(g)+Number(h*S):g+=Number(h)*Number(S)}c.set(p,g)}return c}if(e>=2&&n>=2&&!(e===2&&n===2)){let o=r.shape[e-1],a=t.shape[n-2];if(o!==a)throw new Error(`dot: incompatible shapes ${JSON.stringify(r.shape)} and ${JSON.stringify(t.shape)}`);let i=[...r.shape.slice(0,-1),...t.shape.slice(0,-2),t.shape[n-1]],s=M(r.dtype,t.dtype),u=D.zeros(i,s),c=r.shape.slice(0,-1).reduce((p,m)=>p*m,1),l=t.shape.slice(0,-2).reduce((p,m)=>p*m,1),y=t.shape[n-1],f=o;for(let p=0;p<c;p++)for(let m=0;m<l;m++)for(let d=0;d<y;d++){let g=0;for(let h=0;h<f;h++){let b=p*f+h,S=r.data[b+r.offset],w=m*f*y+h*y+d,x=t.data[w+t.offset];typeof S=="bigint"&&typeof x=="bigint"?g=Number(g)+Number(S*x):g+=Number(S)*Number(x)}let A=p*l*y+m*y+d;u.data[A]=g}return u}throw new Error(`dot: unexpected combination of dimensions ${e}D \xB7 ${n}D`)}function X(r,t){if(r.ndim!==2||t.ndim!==2)throw new Error("matmul requires 2D arrays");let[e=0,n=0]=r.shape,[o=0,a=0]=t.shape;if(n!==o)throw new Error(`matmul shape mismatch: (${e},${n}) @ (${o},${a})`);let i=M(r.dtype,t.dtype),s=i.startsWith("int")||i.startsWith("uint")||i==="bool"?"float64":i;if(s!=="float64")throw new Error(`matmul currently only supports float64, got ${s}`);let u=r.dtype==="float64"?r.data:Float64Array.from(Array.from(r.data).map(Number)),c=t.dtype==="float64"?t.data:Float64Array.from(Array.from(t.data).map(Number));r.offset>0&&(u=u.subarray(r.offset)),t.offset>0&&(c=c.subarray(t.offset));let[l=0,y=0]=r.strides,[f=0,p=0]=t.strides,m=y>l,d=p>f,g=m?"transpose":"no-transpose",A=d?"transpose":"no-transpose",h,b;m?h=y:h=l,d?b=p:b=f;let S=D.zeros([e,a],"float64");return Os("row-major",g,A,e,a,n,1,u,h,c,b,0,S.data,a),S}function Pe(r){if(r.ndim!==2)throw new Error(`trace requires 2D array, got ${r.ndim}D`);let[t=0,e=0]=r.shape,n=Math.min(t,e),o=0;for(let a=0;a<n;a++){let i=r.get(a,a);typeof i=="bigint"?o=(typeof o=="bigint"?o:BigInt(o))+i:o=(typeof o=="bigint"?Number(o):o)+i}return o}function ke(r,t){return pr(r,t)}function Ve(r,t){let e=r.ndim,n=t.ndim,o=r.shape[e-1],a=t.shape[n-1];if(o!==a)throw new Error(`inner: incompatible shapes - last dimensions ${o} and ${a} don't match`);if(e===1&&n===1)return bt(r,t);let i=[...r.shape.slice(0,-1),...t.shape.slice(0,-1)],s=M(r.dtype,t.dtype),u=D.zeros(i,s),c=e===1?1:r.shape.slice(0,-1).reduce((f,p)=>f*p,1),l=n===1?1:t.shape.slice(0,-1).reduce((f,p)=>f*p,1),y=o;for(let f=0;f<c;f++)for(let p=0;p<l;p++){let m=0;for(let g=0;g<y;g++){let A=e===1?g:f*y+g,h=n===1?g:p*y+g,b=r.data[A+r.offset],S=t.data[h+t.offset];typeof b=="bigint"&&typeof S=="bigint"?m=Number(m)+Number(b*S):m+=Number(b)*Number(S)}if(i.length===0)return m;let d=c===1?p:f*l+p;u.data[d]=m}return u}function St(r,t){let e=r.ndim===1?r:lr(r),n=t.ndim===1?t:lr(t),o=e.size,a=n.size,i=M(r.dtype,t.dtype),s=D.zeros([o,a],i);for(let u=0;u<o;u++)for(let c=0;c<a;c++){let l=e.get(u),y=n.get(c),f;typeof l=="bigint"&&typeof y=="bigint"?f=l*y:f=Number(l)*Number(y),s.set([u,c],f)}return s}function Fe(r,t,e){let n,o;if(typeof e=="number"){let f=e;if(f<0)throw new Error("tensordot: axes must be non-negative");if(f>r.ndim||f>t.ndim)throw new Error("tensordot: axes exceeds array dimensions");n=Array.from({length:f},(p,m)=>r.ndim-f+m),o=Array.from({length:f},(p,m)=>m)}else if([n,o]=e,n.length!==o.length)throw new Error("tensordot: axes lists must have same length");for(let f=0;f<n.length;f++){let p=n[f],m=o[f];if(p<0||p>=r.ndim||m<0||m>=t.ndim)throw new Error("tensordot: axis out of bounds");if(r.shape[p]!==t.shape[m])throw new Error(`tensordot: shape mismatch on axes ${p} and ${m}: ${r.shape[p]} != ${t.shape[m]}`)}let a=[],i=[];for(let f=0;f<r.ndim;f++)n.includes(f)||a.push(f);for(let f=0;f<t.ndim;f++)o.includes(f)||i.push(f);let s=[...a.map(f=>r.shape[f]),...i.map(f=>t.shape[f])];if(s.length===0){let f=0,p=n.map(m=>r.shape[m]).reduce((m,d)=>m*d,1);for(let m=0;m<p;m++){let d=m,g=new Array(n.length);for(let w=n.length-1;w>=0;w--){let x=n[w];g[w]=d%r.shape[x],d=Math.floor(d/r.shape[x])}let A=new Array(r.ndim),h=new Array(t.ndim);for(let w=0;w<n.length;w++)A[n[w]]=g[w];for(let w=0;w<o.length;w++)h[o[w]]=g[w];let b=r.get(...A),S=t.get(...h);typeof b=="bigint"&&typeof S=="bigint"?f=Number(f)+Number(b*S):f+=Number(b)*Number(S)}return f}let u=M(r.dtype,t.dtype),c=D.zeros(s,u),l=s.reduce((f,p)=>f*p,1),y=n.map(f=>r.shape[f]).reduce((f,p)=>f*p,1);for(let f=0;f<l;f++){let p=f,m=[];for(let h=s.length-1;h>=0;h--)m[h]=p%s[h],p=Math.floor(p/s[h]);let d=m.slice(0,a.length),g=m.slice(a.length),A=0;for(let h=0;h<y;h++){p=h;let b=[];for(let B=n.length-1;B>=0;B--){let E=n[B];b[B]=p%r.shape[E],p=Math.floor(p/r.shape[E])}let S=new Array(r.ndim),w=new Array(t.ndim);for(let B=0;B<a.length;B++)S[a[B]]=d[B];for(let B=0;B<i.length;B++)w[i[B]]=g[B];for(let B=0;B<n.length;B++)S[n[B]]=b[B],w[o[B]]=b[B];let x=r.get(...S),z=t.get(...w);typeof x=="bigint"&&typeof z=="bigint"?A=Number(A)+Number(x*z):A+=Number(x)*Number(z)}c.set(m,A)}return c}function Re(r,t=0,e=0,n=1){let o=r.shape,a=o.length;if(a<2)throw new Error("diagonal requires an array of at least two dimensions");let i=e<0?a+e:e,s=n<0?a+n:n;if(i<0||i>=a||s<0||s>=a)throw new Error("axis out of bounds");if(i===s)throw new Error("axis1 and axis2 cannot be the same");let u=o[i],c=o[s],l;t>=0?l=Math.max(0,Math.min(u,c-t)):l=Math.max(0,Math.min(u+t,c));let y=[];for(let d=0;d<a;d++)d!==i&&d!==s&&y.push(o[d]);y.push(l);let f=D.zeros(y,r.dtype),p=o.filter((d,g)=>g!==i&&g!==s),m=p.reduce((d,g)=>d*g,1);for(let d=0;d<m;d++){let g=d,A=[];for(let h=p.length-1;h>=0;h--)A.unshift(g%p[h]),g=Math.floor(g/p[h]);for(let h=0;h<l;h++){let b=new Array(a),S=0;for(let z=0;z<a;z++)z===i?b[z]=t>=0?h:h-t:z===s?b[z]=t>=0?h+t:h:b[z]=A[S++];let w=[...A,h],x=r.get(...b);f.set(w,x)}}return f}function qe(r,...t){let e=r.indexOf("->"),n,o;e===-1?(n=r,o=$s(n)):(n=r.slice(0,e),o=r.slice(e+2));let a=n.split(",").map(d=>d.trim());if(a.length!==t.length)throw new Error(`einsum: expected ${a.length} operands, got ${t.length}`);let i=new Map;for(let d=0;d<t.length;d++){let g=a[d],A=t[d];if(g.length!==A.ndim)throw new Error(`einsum: operand ${d} has ${A.ndim} dimensions but subscript '${g}' has ${g.length} indices`);for(let h=0;h<g.length;h++){let b=g[h],S=A.shape[h];if(i.has(b)){if(i.get(b)!==S)throw new Error(`einsum: size mismatch for index '${b}': ${i.get(b)} vs ${S}`)}else i.set(b,S)}}for(let d of o)if(!i.has(d))throw new Error(`einsum: output subscript contains unknown index '${d}'`);let s=new Set(o),u=new Set;for(let d of a)for(let g of d)u.add(g);let c=[];for(let d of u)s.has(d)||c.push(d);if(t.length===2&&a.length===2){let[d,g]=a,[A,h]=t;if(d.length===2&&g.length===2&&o.length===2&&A.ndim===2&&h.ndim===2){let[b,S]=[d[0],d[1]],[w,x]=[g[0],g[1]],[z,B]=[o[0],o[1]];if(b===z&&x===B&&S===w&&c.length===1&&c[0]===S||b===z&&x===B&&S===w&&c.length===1&&c[0]===S)return X(A,h);if(S===z&&x===B&&b===w&&c.length===1&&c[0]===b){let E=ke(A);return X(E,h)}if(b===z&&w===B&&S===x&&c.length===1&&c[0]===S){let E=ke(h);return X(A,E)}}if(d.length===1&&g.length===1&&d===g&&o.length===0&&A.ndim===1&&h.ndim===1)return At(t,a,c,i);if(d&&g&&d.length===1&&g.length===1&&o.length===2&&o===d+g&&c.length===0&&A.ndim===1&&h.ndim===1)return St(A,h)}if(t.length===1&&a[0].length===2&&o.length===0){let d=a[0];if(d[0]===d[1]&&t[0].ndim===2)return At(t,a,c,i)}let l=Array.from(o).map(d=>i.get(d));if(l.length===0)return At(t,a,c,i);let y=t[0].dtype;for(let d=1;d<t.length;d++)y=M(y,t[d].dtype);let f=D.zeros(l,y),p=l.reduce((d,g)=>d*g,1),m=1;for(let d of c)m*=i.get(d);for(let d=0;d<p;d++){let g=vs(d,l),A=new Map;for(let b=0;b<o.length;b++)A.set(o[b],g[b]);let h=0;for(let b=0;b<m;b++){let S=b;for(let x=c.length-1;x>=0;x--){let z=c[x],B=i.get(z);A.set(z,S%B),S=Math.floor(S/B)}let w=1;for(let x=0;x<t.length;x++){let z=t[x],B=a[x],E=[];for(let Sr of B)E.push(A.get(Sr));let O=z.get(...E);w*=Number(O)}h+=w}f.set(g,h)}return f}function $s(r){let t=new Map,e=r.split(",");for(let o of e)for(let a of o.trim())t.set(a,(t.get(a)||0)+1);let n=[];for(let[o,a]of t)a===1&&n.push(o);return n.sort().join("")}function At(r,t,e,n){let o=1;for(let i of e)o*=n.get(i);let a=0;for(let i=0;i<o;i++){let s=new Map,u=i;for(let l=e.length-1;l>=0;l--){let y=e[l],f=n.get(y);s.set(y,u%f),u=Math.floor(u/f)}let c=1;for(let l=0;l<r.length;l++){let y=r[l],f=t[l],p=[];for(let d of f)p.push(s.get(d));let m=y.get(...p);c*=Number(m)}a+=c}return a}function vs(r,t){let e=new Array(t.length),n=r;for(let o=t.length-1;o>=0;o--)e[o]=n%t[o],n=Math.floor(n/t[o]);return e}function Le(r,t){let e=r.shape,n=t.shape,o=e.length,a=n.length,i=M(r.dtype,t.dtype),s=Math.max(o,a),u=new Array(s),c=new Array(s).fill(1),l=new Array(s).fill(1);for(let m=0;m<o;m++)c[s-o+m]=e[m];for(let m=0;m<a;m++)l[s-a+m]=n[m];for(let m=0;m<s;m++)u[m]=c[m]*l[m];let y=D.zeros(u,i),f=e.reduce((m,d)=>m*d,1),p=n.reduce((m,d)=>m*d,1);for(let m=0;m<f;m++){let d=m,g=new Array(o);for(let b=o-1;b>=0;b--)g[b]=d%e[b],d=Math.floor(d/e[b]);let A=new Array(s).fill(0);for(let b=0;b<o;b++)A[s-o+b]=g[b];let h=r.get(...g);for(let b=0;b<p;b++){let S=b,w=new Array(a);for(let O=a-1;O>=0;O--)w[O]=S%n[O],S=Math.floor(S/n[O]);let x=new Array(s).fill(0);for(let O=0;O<a;O++)x[s-a+O]=w[O];let z=t.get(...w),B=new Array(s);for(let O=0;O<s;O++)B[O]=A[O]*l[O]+x[O];let E=typeof h=="bigint"||typeof z=="bigint"?BigInt(Number(h))*BigInt(Number(z)):Number(h)*Number(z);y.set(B,E)}}return y}var je=$(()=>{"use strict";F();V();ht()});function Ze(r){return C(r,Math.sqrt,!1)}function We(r,t){return typeof t=="number"?Cs(r,t):q(r,t,Math.pow,"power")}function Cs(r,t){let e=r.dtype,n=Array.from(r.shape),o=r.data,a=r.size,u=e!=="float32"&&e!=="float64"&&(t<0||!Number.isInteger(t))?"float64":e,c=D.zeros(n,u),l=c.data;if(N(e))if(N(u)&&Number.isInteger(t)&&t>=0){let y=o,f=l;for(let p=0;p<a;p++)f[p]=y[p]**BigInt(t)}else for(let y=0;y<a;y++)l[y]=Math.pow(Number(o[y]),t);else for(let y=0;y<a;y++)l[y]=Math.pow(Number(o[y]),t);return c}var Ge=$(()=>{"use strict";F();tr();V()});function He(r){return C(r,Math.sin,!1)}function Ye(r){return C(r,Math.cos,!1)}function Je(r){return C(r,Math.tan,!1)}function Qe(r){return C(r,Math.asin,!1)}function Xe(r){return C(r,Math.acos,!1)}function Ke(r){return C(r,Math.atan,!1)}function rn(r,t){return typeof t=="number"?Vs(r,t):Ps(r,t)}function Ps(r,t){let e=Array.from(r.shape),n=r.size,o=r.dtype,a=t.dtype,i=o==="float32"&&a==="float32"?"float32":"float64",s=D.zeros(e,i),u=s.data;for(let c=0;c<n;c++){let l=(N(o),Number(r.data[c])),y=(N(a),Number(t.data[c]));u[c]=Math.atan2(l,y)}return s}function Vs(r,t){let e=r.dtype,n=Array.from(r.shape),o=r.data,a=r.size,i=e==="float32"?"float32":"float64",s=D.zeros(n,i),u=s.data;if(N(e))for(let c=0;c<a;c++)u[c]=Math.atan2(Number(o[c]),t);else for(let c=0;c<a;c++)u[c]=Math.atan2(Number(o[c]),t);return s}function tn(r,t){return typeof t=="number"?Rs(r,t):Fs(r,t)}function Fs(r,t){let e=Array.from(r.shape),n=r.size,o=r.dtype,a=t.dtype,i=o==="float32"&&a==="float32"?"float32":"float64",s=D.zeros(e,i),u=s.data;for(let c=0;c<n;c++){let l=(N(o),Number(r.data[c])),y=(N(a),Number(t.data[c]));u[c]=Math.hypot(l,y)}return s}function Rs(r,t){let e=r.dtype,n=Array.from(r.shape),o=r.data,a=r.size,i=e==="float32"?"float32":"float64",s=D.zeros(n,i),u=s.data;if(N(e))for(let c=0;c<a;c++)u[c]=Math.hypot(Number(o[c]),t);else for(let c=0;c<a;c++)u[c]=Math.hypot(Number(o[c]),t);return s}function en(r){let t=180/Math.PI;return C(r,e=>e*t,!1)}function nn(r){let t=Math.PI/180;return C(r,e=>e*t,!1)}var on=$(()=>{"use strict";F();tr();V()});function an(r){return C(r,Math.sinh,!1)}function sn(r){return C(r,Math.cosh,!1)}function un(r){return C(r,Math.tanh,!1)}function cn(r){return C(r,Math.asinh,!1)}function fn(r){return C(r,Math.acosh,!1)}function ln(r){return C(r,Math.atanh,!1)}var pn=$(()=>{"use strict";tr()});function Dt(r,t){let e=r.shape,n=e.length,o=t.length;if(o<n)throw new Error("input operand has more dimensions than allowed by the axis remapping");let a=J([Array.from(e),t]);if(a===null)throw new Error(`operands could not be broadcast together with shape (${e.join(",")}) (${t.join(",")})`);for(let i=0;i<o;i++)if(a[i]!==t[i])throw new Error(`operands could not be broadcast together with shape (${e.join(",")}) (${t.join(",")})`);return H(r,t)}function yn(r){if(r.length===0)return[];if(r.length===1)return[r[0]];let t=r.map(n=>Array.from(n.shape)),e=J(t);if(e===null)throw new Error(`operands could not be broadcast together with shapes ${t.map(n=>`(${n.join(",")})`).join(" ")}`);return r.map(n=>H(n,e))}function mn(r,t,e){let n=r.shape,o=n.length,a=r.dtype;if(e===void 0){let m=r.size;for(let h of t){let b=h<0?m+h:h;if(b<0||b>=m)throw new Error(`index ${h} is out of bounds for axis 0 with size ${m}`)}let d=t.length,g=_(a);if(!g)throw new Error(`Cannot take from array with dtype ${a}`);let A=new g(d);for(let h=0;h<d;h++){let b=t[h];b<0&&(b=m+b);let S=r.iget(b);N(a),A[h]=S}return D.fromData(A,[d],a)}let i=e<0?o+e:e;if(i<0||i>=o)throw new Error(`axis ${e} is out of bounds for array of dimension ${o}`);let s=n[i];for(let m of t){let d=m<0?s+m:m;if(d<0||d>=s)throw new Error(`index ${m} is out of bounds for axis ${i} with size ${s}`)}let u=Array.from(n);u[i]=t.length;let c=u.reduce((m,d)=>m*d,1),l=_(a);if(!l)throw new Error(`Cannot take from array with dtype ${a}`);let y=new l(c),f=Z(u),p=new Array(o).fill(0);for(let m=0;m<c;m++){let d=[...p],g=p[i],A=t[g];A<0&&(A=s+A),d[i]=A;let h=r.get(...d),b=0;for(let S=0;S<o;S++)b+=p[S]*f[S];N(a),y[b]=h;for(let S=o-1;S>=0&&(p[S]++,!(p[S]<u[S]));S--)p[S]=0}return D.fromData(y,u,a)}function dn(r,t,e){let n=r.size,o=r.dtype,a;if(typeof e=="number"||typeof e=="bigint")a=new Array(t.length).fill(e);else{a=[];for(let i=0;i<e.size;i++)a.push(e.iget(i));if(a.length===1)a=new Array(t.length).fill(a[0]);else if(a.length!==t.length){let i=[...a];a=[];for(let s=0;s<t.length;s++)a.push(i[s%i.length])}}for(let i=0;i<t.length;i++){let s=t[i];if(s<0&&(s=n+s),s<0||s>=n)throw new Error(`index ${t[i]} is out of bounds for axis 0 with size ${n}`);let u=a[i];N(o)?typeof u!="bigint"&&(u=BigInt(Math.round(Number(u)))):typeof u=="bigint"&&(u=Number(u)),r.iset(s,u)}}function gn(r,t){if(t.length===0)throw new Error("choices cannot be empty");let e=r.shape,n=t.length,o=t[0].dtype,a=t.map(f=>Array.from(f.shape));a.unshift(Array.from(e));let i=J(a);if(i===null)throw new Error("operands could not be broadcast together");let s=H(r,i),u=t.map(f=>H(f,i)),c=i.reduce((f,p)=>f*p,1),l=_(o);if(!l)throw new Error(`Cannot choose with dtype ${o}`);let y=new l(c);for(let f=0;f<c;f++){let p=Number(s.iget(f));if(p<0||p>=n)throw new Error(`index ${p} is out of bounds for axis 0 with size ${n}`);let m=u[p].iget(f);N(o),y[f]=m}return D.fromData(y,i,o)}function hn(r,t,e=!1){if(r.ndim!==t.ndim)return!1;for(let o=0;o<r.ndim;o++)if(r.shape[o]!==t.shape[o])return!1;let n=r.size;for(let o=0;o<n;o++){let a=r.iget(o),i=t.iget(o);if(e){let s=typeof a=="number"&&Number.isNaN(a),u=typeof i=="number"&&Number.isNaN(i);if(s&&u)continue}if(a!==i)return!1}return!0}var An=$(()=>{"use strict";F();V();Kr()});function L(r,t=v){let e=D.zeros(r,t);return new I(e)}function Nt(r,t=v){let e=D.ones(r,t);return new I(e)}function Zs(r){let t=[],e=r;for(;Array.isArray(e);)t.push(e.length),e=e[0];return t}function bn(r){return typeof r=="bigint"?!0:Array.isArray(r)?r.some(t=>bn(t)):!1}function Ws(r){let t=[];function e(n){Array.isArray(n)?n.forEach(o=>e(o)):t.push(n)}return e(r),t}function k(r,t){if(r instanceof I)return!t||r.dtype===t?r.copy():r.astype(t);let e=bn(r),n=Zs(r),o=n.reduce((l,y)=>l*y,1),a=t;a||(e?a="int64":a=v);let i=_(a);if(!i)throw new Error(`Cannot create array with dtype ${a}`);let s=new i(o),u=Ws(r);if(N(a)){let l=s;for(let y=0;y<o;y++){let f=u[y];l[y]=typeof f=="bigint"?f:BigInt(Math.round(Number(f)))}}else if(a==="bool"){let l=s;for(let y=0;y<o;y++)l[y]=u[y]?1:0}else{let l=s;for(let y=0;y<o;y++){let f=u[y];l[y]=Number(f)}}let c=D.fromData(s,n,a);return new I(c)}function Sn(r,t,e=1,n=v){let o=r,a=t;if(t===void 0&&(o=0,a=r),a===void 0)throw new Error("stop is required");let i=Math.max(0,Math.ceil((a-o)/e)),s=_(n);if(!s)throw new Error(`Cannot create arange array with dtype ${n}`);let u=new s(i);if(N(n))for(let l=0;l<i;l++)u[l]=BigInt(Math.round(o+l*e));else if(n==="bool")for(let l=0;l<i;l++)u[l]=o+l*e!==0?1:0;else for(let l=0;l<i;l++)u[l]=o+l*e;let c=D.fromData(u,[i],n);return new I(c)}function Dn(r,t,e=50,n=v){if(e<0)throw new Error("num must be non-negative");if(e===0)return k([],n);if(e===1)return k([r],n);let o=_(n);if(!o)throw new Error(`Cannot create linspace array with dtype ${n}`);let a=new o(e),i=(t-r)/(e-1);if(N(n))for(let u=0;u<e;u++)a[u]=BigInt(Math.round(r+u*i));else if(n==="bool")for(let u=0;u<e;u++)a[u]=r+u*i!==0?1:0;else for(let u=0;u<e;u++)a[u]=r+u*i;let s=D.fromData(a,[e],n);return new I(s)}function Nn(r,t,e=50,n=10,o=v){if(e<0)throw new Error("num must be non-negative");if(e===0)return k([],o);if(e===1)return k([Math.pow(n,r)],o);let a=_(o);if(!a)throw new Error(`Cannot create logspace array with dtype ${o}`);let i=new a(e),s=(t-r)/(e-1);if(N(o))for(let c=0;c<e;c++){let l=r+c*s;i[c]=BigInt(Math.round(Math.pow(n,l)))}else if(o==="bool")for(let c=0;c<e;c++){let l=r+c*s;i[c]=Math.pow(n,l)!==0?1:0}else for(let c=0;c<e;c++){let l=r+c*s;i[c]=Math.pow(n,l)}let u=D.fromData(i,[e],o);return new I(u)}function wn(r,t,e=50,n=v){if(e<0)throw new Error("num must be non-negative");if(r===0||t===0)throw new Error("Geometric sequence cannot include zero");if(e===0)return k([],n);if(e===1)return k([r],n);let o=Math.sign(r),a=Math.sign(t);if(o!==a)throw new Error("Geometric sequence cannot contain both positive and negative values");let i=_(n);if(!i)throw new Error(`Cannot create geomspace array with dtype ${n}`);let s=new i(e),u=Math.log(Math.abs(r)),l=(Math.log(Math.abs(t))-u)/(e-1);if(N(n))for(let f=0;f<e;f++){let p=o*Math.exp(u+f*l);s[f]=BigInt(Math.round(p))}else if(n==="bool")for(let f=0;f<e;f++){let p=o*Math.exp(u+f*l);s[f]=p!==0?1:0}else for(let f=0;f<e;f++){let p=o*Math.exp(u+f*l);s[f]=p}let y=D.fromData(s,[e],n);return new I(y)}function wt(r,t,e=0,n=v){let o=t??r,a=L([r,o],n),i=a.data;if(N(n)){let s=i;for(let u=0;u<r;u++){let c=u+e;c>=0&&c<o&&(s[u*o+c]=BigInt(1))}}else{let s=i;for(let u=0;u<r;u++){let c=u+e;c>=0&&c<o&&(s[u*o+c]=1)}}return a}function xt(r,t=v){return L(r,t)}function It(r,t,e){let n=e;n||(typeof t=="bigint"?n="int64":typeof t=="boolean"?n="bool":Number.isInteger(t)?n="int32":n=v);let o=_(n);if(!o)throw new Error(`Cannot create full array with dtype ${n}`);let a=r.reduce((u,c)=>u*c,1),i=new o(a);if(N(n)){let u=typeof t=="bigint"?t:BigInt(Math.round(Number(t)));i.fill(u)}else n==="bool"?i.fill(t?1:0):i.fill(Number(t));let s=D.fromData(i,r,n);return new I(s)}function xn(r,t=v){return wt(r,r,0,t)}function mr(r,t){return r instanceof I?!t||r.dtype===t?r:r.astype(t):k(r,t)}function In(r){return r.copy()}function zn(r,t){return L(Array.from(r.shape),t??r.dtype)}function Bn(r,t){return Nt(Array.from(r.shape),t??r.dtype)}function _n(r,t){return xt(Array.from(r.shape),t??r.dtype)}function En(r,t,e){return It(Array.from(r.shape),t,e??r.dtype)}function Un(r,t){return mr(r,t)}function Tn(r,t){let e=mr(r,t);return e.flags.C_CONTIGUOUS?e:e.copy()}function On(r,t){return mr(r,t).copy()}function zt(r,t=0){if(r.ndim===1){let e=r.size,n=e+Math.abs(t),o=L([n,n],r.dtype);for(let a=0;a<e;a++){let i=t>=0?a:a-t,s=t>=0?a+t:a;o.set([i,s],r.get([a]))}return o}else if(r.ndim===2){let[e,n]=r.shape,o,a,i;if(t>=0?(o=0,a=t,i=Math.min(e,n-t)):(o=-t,a=0,i=Math.min(e+t,n)),i<=0)return L([0],r.dtype);let s=_(r.dtype),u=new s(i);for(let l=0;l<i;l++){let y=r.get([o+l,a+l]);N(r.dtype)?u[l]=typeof y=="bigint"?y:BigInt(y):u[l]=y}let c=D.fromData(u,[i],r.dtype);return new I(c)}else throw new Error("Input must be 1-D or 2-D")}function $n(r,t=0){let e=r.flatten();return zt(e,t)}function vn(r,t,e=v){let n=t.reduce((c,l)=>c*l,1),o=_(e);if(!o)throw new Error(`Cannot create array with dtype ${e}`);let a=new o(n),i=t.length,s=new Array(i).fill(0);for(let c=0;c<n;c++){let l=r(...s);N(e)?a[c]=typeof l=="bigint"?l:BigInt(Number(l)):e==="bool"?a[c]=l?1:0:a[c]=Number(l);for(let y=i-1;y>=0&&(s[y]++,!(s[y]<t[y]));y--)s[y]=0}let u=D.fromData(a,t,e);return new I(u)}function Mn(...r){let t=[],e="xy";for(let s of r)s instanceof I?t.push(s):typeof s=="object"&&"indexing"in s&&(e=s.indexing||"xy");if(t.length===0)return[];if(t.length===1)return[t[0].copy()];let n=t.map(s=>s.size);e==="xy"&&t.length>=2&&(t=[t[1],t[0],...t.slice(2)],[n[0],n[1]]=[n[1],n[0]]);let o=n,a=o.length,i=[];for(let s=0;s<t.length;s++){let u=t[s],c=u.size,l=new Array(a).fill(1);l[s]=c;let y=u.reshape(...l),f=Dt(y.storage,o),p=I._fromStorage(f.copy());i.push(p)}return e==="xy"&&i.length>=2&&([i[0],i[1]]=[i[1],i[0]]),i}function Cn(r,t,e=0,n=v){let o=t??r,a=L([r,o],n);for(let i=0;i<r;i++)for(let s=0;s<=i+e&&s<o;s++)s>=0&&a.set([i,s],1);return a}function kn(r,t=0){if(r.ndim<2)throw new Error("Input must have at least 2 dimensions");let e=r.copy(),n=e.shape,o=n[n.length-2],a=n[n.length-1],i=n.slice(0,-2).reduce((s,u)=>s*u,1);for(let s=0;s<i;s++)for(let u=0;u<o;u++)for(let c=0;c<a;c++)if(c>u+t){let l=[],y=s;for(let f=n.length-3;f>=0;f--)l.unshift(y%n[f]),y=Math.floor(y/n[f]);l.push(u,c),e.set(l,0)}return e}function Pn(r,t=0){if(r.ndim<2)throw new Error("Input must have at least 2 dimensions");let e=r.copy(),n=e.shape,o=n[n.length-2],a=n[n.length-1],i=n.slice(0,-2).reduce((s,u)=>s*u,1);for(let s=0;s<i;s++)for(let u=0;u<o;u++)for(let c=0;c<a;c++)if(c<u+t){let l=[],y=s;for(let f=n.length-3;f>=0;f--)l.unshift(y%n[f]),y=Math.floor(y/n[f]);l.push(u,c),e.set(l,0)}return e}function Vn(r,t,e=!1){if(r.ndim!==1)throw new Error("Input must be 1-D");let n=r.size,o=t??n;if(o<0)throw new Error("N must be non-negative");let a=L([n,o],r.dtype);for(let i=0;i<n;i++){let s=r.get([i]);for(let u=0;u<o;u++){let c=e?u:o-1-u;a.set([i,u],Math.pow(s,c))}}return a}function Fn(r,t=v,e=-1,n=0){let o,a=n;r instanceof ArrayBuffer?o=r:(o=r.buffer,a+=r.byteOffset);let i=Gs(t),s=o.byteLength-a,u=Math.floor(s/i),c=e<0?u:Math.min(e,u);if(c<=0)return k([],t);let l=_(t);if(!l)throw new Error(`Unsupported dtype: ${t}`);let y=new l(o,a,c),f=D.fromData(y,[c],t);return new I(f)}function Rn(r,t=v,e=-1){let n=[],o=0;for(let a of r){if(e>=0&&o>=e)break;n.push(a),o++}return k(n,t)}function qn(r,t=v,e=-1){let n=[],o=0;for(let a of r){if(e>=0&&o>=e)break;n.push(a),o++}return k(n,t)}function Ln(r,t=v,e=-1,n=""){let o;n===""?o=r.trim().split(/\s+/):o=r.split(n);let a=[],i=0;for(let s of o){if(e>=0&&i>=e)break;let u=s.trim();u!==""&&(N(t)?a.push(BigInt(u)):a.push(parseFloat(u)),i++)}return k(a,t)}function Gs(r){switch(r){case"int8":case"uint8":case"bool":return 1;case"int16":case"uint16":return 2;case"int32":case"uint32":case"float32":return 4;case"int64":case"uint64":case"float64":return 8;default:return 8}}function jn(r){return r.sqrt()}function Zn(r,t){return r.power(t)}function Wn(r){return r.absolute()}function Gn(r){return r.negative()}function Hn(r){return r.sign()}function Yn(r,t){return r.mod(t)}function Jn(r,t){return r.floor_divide(t)}function Qn(r){return r.positive()}function Xn(r){return r.reciprocal()}function Kn(r,t){return r.dot(t)}function ro(r){return r.trace()}function to(r,t=0,e=0,n=1){let o=Re(r.storage,t,e,n);return I._fromStorage(o)}function eo(r,t){let e=Le(r.storage,t.storage);return I._fromStorage(e)}function no(r,t){return r.transpose(t)}function oo(r,t){return r.inner(t)}function ao(r,t){return r.outer(t)}function so(r,t,e=2){return r.tensordot(t,e)}function io(r){return r.sin()}function uo(r){return r.cos()}function co(r){return r.tan()}function fo(r){return r.arcsin()}function lo(r){return r.arccos()}function po(r){return r.arctan()}function yo(r,t){return r.arctan2(t)}function mo(r,t){return r.hypot(t)}function go(r){return r.degrees()}function ho(r){return r.radians()}function Ao(r){return r.radians()}function bo(r){return r.degrees()}function So(r){return r.sinh()}function Do(r){return r.cosh()}function No(r){return r.tanh()}function wo(r){return r.arcsinh()}function xo(r){return r.arccosh()}function Io(r){return r.arctanh()}function zo(r,t,e){return r.swapaxes(t,e)}function Bo(r,t,e){return r.moveaxis(t,e)}function K(r,t=0){if(r.length===0)throw new Error("need at least one array to concatenate");let e=r.map(o=>o.storage),n=Q(e,t);return I._fromStorage(n)}function _o(r,t=0){if(r.length===0)throw new Error("need at least one array to stack");let e=r.map(o=>o.storage),n=he(e,t);return I._fromStorage(n)}function Bt(r){if(r.length===0)throw new Error("need at least one array to stack");let t=r.map(n=>n.storage),e=Ae(t);return I._fromStorage(e)}function Eo(r){if(r.length===0)throw new Error("need at least one array to stack");let t=r.map(n=>n.storage),e=gt(t);return I._fromStorage(e)}function Uo(r){if(r.length===0)throw new Error("need at least one array to stack");let t=r.map(n=>n.storage),e=be(t);return I._fromStorage(e)}function To(r,t,e=0){return Se(r.storage,t,e).map(o=>I._fromStorage(o,r.base??r))}function Oo(r,t,e=0){return yr(r.storage,t,e).map(o=>I._fromStorage(o,r.base??r))}function $o(r,t){return Ne(r.storage,t).map(n=>I._fromStorage(n,r.base??r))}function vo(r,t){return we(r.storage,t).map(n=>I._fromStorage(n,r.base??r))}function Mo(r,t){let e=xe(r.storage,t);return I._fromStorage(e)}function Co(r,t,e){return r.repeat(t,e)}function ko(r){return r.ravel()}function Po(r,t){return r.reshape(...t)}function Vo(r,t){return r.squeeze(t)}function Fo(r,t){return r.expand_dims(t)}function Cr(r,t){let e=ze(r.storage,t);return I._fromStorage(e)}function Ro(r){if(r.ndim<2)throw new Error("Input must be at least 2-D");return Cr(r,1)}function qo(r){if(r.ndim<2)throw new Error("Input must be at least 2-D");return Cr(r,0)}function Lo(r,t=1,e=[0,1]){let n=Be(r.storage,t,e);return I._fromStorage(n)}function jo(r,t,e){let n=_e(r.storage,t,e);return I._fromStorage(n)}function Zo(r,t,e=0){let n=Ee(r.storage,t,e);return I._fromStorage(n,r.base??r)}function Wo(...r){let t=r.map(o=>o.storage),n=$e(t).map((o,a)=>o===t[a]?r[a]:I._fromStorage(o));return n.length===1?n[0]:n}function Go(...r){let t=r.map(o=>o.storage),n=ve(t).map((o,a)=>o===t[a]?r[a]:I._fromStorage(o));return n.length===1?n[0]:n}function Ho(...r){let t=r.map(o=>o.storage),n=Me(t).map((o,a)=>o===t[a]?r[a]:I._fromStorage(o));return n.length===1?n[0]:n}function Yo(r,t){return Ue(r.storage,t).map(n=>I._fromStorage(n,r.base??r))}function Jo(r){if(r.length===0)throw new Error("need at least one array to stack");let t=r.map(n=>n.storage),e=Te(t);return I._fromStorage(e)}function Qo(r){return Bt(r)}function Xo(r,t){let e=Oe(r.storage,t);return I._fromStorage(e)}function Ko(r,t,e){let n=t instanceof I?t:k(t,r.dtype);if(e===void 0){let o=r.flatten(),a=n.flatten();return K([o,a])}return K([r,n],e)}function ra(r,t,e){let n=r.dtype;if(e===void 0){let p=r.flatten(),d=(Array.isArray(t)?t:[t]).map(S=>S<0?p.size+S:S),g=[];for(let S=0;S<p.size;S++)d.includes(S)||g.push(S);let A=_(n),h=new A(g.length);for(let S=0;S<g.length;S++){let w=p.get([g[S]]);N(n)?h[S]=typeof w=="bigint"?w:BigInt(w):h[S]=w}let b=D.fromData(h,[g.length],n);return new I(b)}let o=r.shape,a=o.length,i=e<0?a+e:e;if(i<0||i>=a)throw new Error(`axis ${e} is out of bounds for array of dimension ${a}`);let s=o[i],u=Array.isArray(t)?t:[t],c=new Set(u.map(p=>p<0?s+p:p)),l=[],y=0;for(let p=0;p<=s;p++)(c.has(p)||p===s)&&(p>y&&l.push([y,p]),y=p+1);if(l.length===0){let p=[...o];return p[i]=0,L(p,n)}let f=[];for(let[p,m]of l){let d=o.map(()=>":");d[i]=`${p}:${m}`,f.push(r.slice(...d))}return K(f,i)}function ta(r,t,e,n){let o=e instanceof I?e:k(e,r.dtype);if(n===void 0){let y=r.flatten(),f=o.flatten(),p=t<0?y.size+t:t;if(p<0||p>y.size)throw new Error(`index ${t} is out of bounds for array of size ${y.size}`);let m=p>0?y.slice(`0:${p}`):null,d=p<y.size?y.slice(`${p}:`):null,g=[];return m&&g.push(m),g.push(f),d&&g.push(d),K(g)}let a=r.shape,i=a.length,s=n<0?i+n:n;if(s<0||s>=i)throw new Error(`axis ${n} is out of bounds for array of dimension ${i}`);let u=a[s],c=t<0?u+t:t;if(c<0||c>u)throw new Error(`index ${t} is out of bounds for axis ${n} with size ${u}`);let l=[];if(c>0){let y=a.map(()=>":");y[s]=`0:${c}`,l.push(r.slice(...y))}if(l.push(o),c<u){let y=a.map(()=>":");y[s]=`${c}:`,l.push(r.slice(...y))}return K(l,s)}function ea(r,t,e="constant",n=0){let o=r.shape,a=o.length,i=r.dtype,s;if(typeof t=="number"?s=o.map(()=>[t,t]):Array.isArray(t)&&typeof t[0]=="number"?s=o.map(()=>t):s=t,s.length!==a)throw new Error(`pad_width must have ${a} elements`);let u=o.map((d,g)=>d+s[g][0]+s[g][1]),c=u.reduce((d,g)=>d*g,1),l=_(i),y=new l(c),f=N(i);e==="constant"&&(f?y.fill(BigInt(n)):y.fill(n));let p=new Array(a).fill(0);for(let d=0;d<c;d++){let g=!0,A=[];for(let b=0;b<a;b++){let[S]=s[b],w=p[b]-S;if(w<0||w>=o[b]){g=!1;break}A.push(w)}let h;if(g)h=r.get(A);else if(e==="constant"){for(let b=a-1;b>=0&&(p[b]++,!(p[b]<u[b]));b--)p[b]=0;continue}else{let b=[];for(let S=0;S<a;S++){let[w]=s[S],x=p[S]-w,z=o[S];x<0?e==="edge"?x=0:e==="reflect"?(x=-x,x>=z&&(x=z-1)):e==="symmetric"?(x=-x-1,x>=z&&(x=z-1),x<0&&(x=0)):e==="wrap"&&(x=(x%z+z)%z):x>=z&&(e==="edge"?x=z-1:e==="reflect"?(x=2*z-x-2,x<0&&(x=0)):e==="symmetric"?(x=2*z-x-1,x<0&&(x=0)):e==="wrap"&&(x=x%z)),b.push(Math.max(0,Math.min(z-1,x)))}h=r.get(b)}f?y[d]=typeof h=="bigint"?h:BigInt(Number(h)):y[d]=Number(h);for(let b=a-1;b>=0&&(p[b]++,!(p[b]<u[b]));b--)p[b]=0}let m=D.fromData(y,u,i);return new I(m)}function na(r,t){let e=Dt(r.storage,t);return I._fromStorage(e,r.base??r)}function oa(...r){let t=r.map(n=>n.storage);return yn(t).map((n,o)=>I._fromStorage(n,r[o].base??r[o]))}function aa(...r){return Xr(...r)}function sa(r,t,e){return r.take(t,e)}function ia(r,t,e){r.put(t,e)}function ua(r,t){let e=t.map(o=>o.storage),n=gn(r.storage,e);return I._fromStorage(n)}function ca(r,t,e=!1){return hn(r.storage,t.storage,e)}function fa(r,t){return fe(r.storage,t.storage)}function la(r,t){return I._fromStorage(it(r.storage,t))}function pa(r,t){return I._fromStorage(ut(r.storage,t))}function ya(r,t,e=!1){let n=ct(r.storage,t,e);return typeof n=="number"?n:I._fromStorage(n)}function ma(r,t,e=!1){let n=ft(r.storage,t,e);return typeof n=="number"?n:I._fromStorage(n)}function da(r,t,e,n=!1){let o=lt(r.storage,t,e,n);return typeof o=="number"?o:I._fromStorage(o)}function ga(r,t,e,n=!1){let o=nr(r.storage,t,e,n);return typeof o=="number"?o:I._fromStorage(o)}function ha(r,t,e,n=!1){let o=e?e.storage:void 0,a=zr(r.storage,t,o,n);return typeof a=="number"?a:I._fromStorage(a)}function Aa(r,t,e=!1){let n=Br(r.storage,t,e);return typeof n=="number"?n:I._fromStorage(n)}function ba(r,t,e=!1){let n=_r(r.storage,t,e);return typeof n=="number"?n:I._fromStorage(n)}function Sa(r,t,e=!1){let n=Er(r.storage,t,e);return typeof n=="number"?n:I._fromStorage(n)}function Da(r,t,e=0,n=!1){let o=cr(r.storage,t,e,n);return typeof o=="number"?o:I._fromStorage(o)}function Na(r,t,e=0,n=!1){let o=pt(r.storage,t,e,n);return typeof o=="number"?o:I._fromStorage(o)}function wa(r,t,e=!1){let n=Ur(r.storage,t,e);return typeof n=="number"?n:I._fromStorage(n)}function xa(r,t,e=!1){let n=Tr(r.storage,t,e);return typeof n=="number"?n:I._fromStorage(n)}function Ia(r,t){let e=Or(r.storage,t);return typeof e=="number"?e:I._fromStorage(e)}function za(r,t){let e=$r(r.storage,t);return typeof e=="number"?e:I._fromStorage(e)}function Ba(r,t){return I._fromStorage(yt(r.storage,t))}function _a(r,t){return I._fromStorage(mt(r.storage,t))}function Ea(r,t,e=!1){let n=vr(r.storage,t,e);return typeof n=="number"?n:I._fromStorage(n)}function Ua(r){return r.cbrt()}function Ta(r){return r.fabs()}function Oa(r,t){return r.divmod(t)}function $a(r){return r.square()}function va(r,t){return r.remainder(t)}function Ma(r,t){return r.heaviside(t)}function Ca(r,...t){let e=t.map(o=>o.storage),n=qe(r,...e);return typeof n=="number"||typeof n=="bigint"?n:I._fromStorage(n)}var I,_t=$(()=>{"use strict";kt();V();F();ee();le();me();ht();je();Ge();on();pn();An();I=class r{constructor(t,e){this._storage=t,this._base=e}get storage(){return this._storage}static _fromStorage(t,e){return new r(t,e)}get shape(){return this._storage.shape}get ndim(){return this._storage.ndim}get size(){return this._storage.size}get dtype(){return this._storage.dtype}get data(){return this._storage.data}get strides(){return this._storage.strides}get flags(){return{C_CONTIGUOUS:this._storage.isCContiguous,F_CONTIGUOUS:this._storage.isFContiguous,OWNDATA:this._base===void 0}}get base(){return this._base??null}get(t){if(t.length!==this.ndim)throw new Error(`Index has ${t.length} dimensions, but array has ${this.ndim} dimensions`);let e=t.map((n,o)=>{let a=n;if(a<0&&(a=this.shape[o]+a),a<0||a>=this.shape[o])throw new Error(`Index ${n} is out of bounds for axis ${o} with size ${this.shape[o]}`);return a});return this._storage.get(...e)}set(t,e){if(t.length!==this.ndim)throw new Error(`Index has ${t.length} dimensions, but array has ${this.ndim} dimensions`);let n=t.map((i,s)=>{let u=i;if(u<0&&(u=this.shape[s]+u),u<0||u>=this.shape[s])throw new Error(`Index ${i} is out of bounds for axis ${s} with size ${this.shape[s]}`);return u}),o=this.dtype,a;N(o)?a=typeof e=="bigint"?e:BigInt(Math.round(e)):o==="bool"?a=e?1:0:a=Number(e),this._storage.set(n,a)}copy(){return new r(this._storage.copy())}astype(t,e=!0){let n=this.dtype;if(n===t&&!e)return this;if(n===t&&e)return this.copy();let o=Array.from(this.shape),a=this.size,i=_(t);if(!i)throw new Error(`Cannot convert to dtype ${t}`);let s=new i(a),u=this.data;if(N(n)&&!N(t)){let l=u;if(t==="bool")for(let y=0;y<a;y++)s[y]=l[y]!==BigInt(0)?1:0;else for(let y=0;y<a;y++)s[y]=Number(l[y])}else if(!N(n)&&N(t)){let l=u;for(let y=0;y<a;y++)s[y]=BigInt(Math.round(Number(l[y])))}else if(t==="bool"){let l=u;for(let y=0;y<a;y++)s[y]=l[y]!==0?1:0}else if(n==="bool"&&!N(t)){let l=u;for(let y=0;y<a;y++)s[y]=l[y]}else if(!N(n)&&!N(t)){let l=u;for(let y=0;y<a;y++)s[y]=l[y]}else{let l=u;for(let y=0;y<a;y++)s[y]=l[y]}let c=D.fromData(s,o,t);return new r(c)}add(t){let e=typeof t=="number"?t:t._storage,n=Rt(this._storage,e);return r._fromStorage(n)}subtract(t){let e=typeof t=="number"?t:t._storage,n=qt(this._storage,e);return r._fromStorage(n)}multiply(t){let e=typeof t=="number"?t:t._storage,n=Lt(this._storage,e);return r._fromStorage(n)}divide(t){let e=typeof t=="number"?t:t._storage,n=jt(this._storage,e);return r._fromStorage(n)}mod(t){let e=typeof t=="number"?t:t._storage,n=Nr(this._storage,e);return r._fromStorage(n)}floor_divide(t){let e=typeof t=="number"?t:t._storage,n=Qr(this._storage,e);return r._fromStorage(n)}positive(){let t=Ht(this._storage);return r._fromStorage(t)}reciprocal(){let t=Yt(this._storage);return r._fromStorage(t)}sqrt(){let t=Ze(this._storage);return r._fromStorage(t)}power(t){let e=typeof t=="number"?t:t._storage,n=We(this._storage,e);return r._fromStorage(n)}absolute(){let t=Zt(this._storage);return r._fromStorage(t)}negative(){let t=Wt(this._storage);return r._fromStorage(t)}sign(){let t=Gt(this._storage);return r._fromStorage(t)}sin(){let t=He(this._storage);return r._fromStorage(t)}cos(){let t=Ye(this._storage);return r._fromStorage(t)}tan(){let t=Je(this._storage);return r._fromStorage(t)}arcsin(){let t=Qe(this._storage);return r._fromStorage(t)}arccos(){let t=Xe(this._storage);return r._fromStorage(t)}arctan(){let t=Ke(this._storage);return r._fromStorage(t)}arctan2(t){let e=typeof t=="number"?t:t._storage,n=rn(this._storage,e);return r._fromStorage(n)}hypot(t){let e=typeof t=="number"?t:t._storage,n=tn(this._storage,e);return r._fromStorage(n)}degrees(){let t=en(this._storage);return r._fromStorage(t)}radians(){let t=nn(this._storage);return r._fromStorage(t)}sinh(){let t=an(this._storage);return r._fromStorage(t)}cosh(){let t=sn(this._storage);return r._fromStorage(t)}tanh(){let t=un(this._storage);return r._fromStorage(t)}arcsinh(){let t=cn(this._storage);return r._fromStorage(t)}arccosh(){let t=fn(this._storage);return r._fromStorage(t)}arctanh(){let t=ln(this._storage);return r._fromStorage(t)}greater(t){let e=typeof t=="number"?t:t._storage,n=ne(this._storage,e);return r._fromStorage(n)}greater_equal(t){let e=typeof t=="number"?t:t._storage,n=oe(this._storage,e);return r._fromStorage(n)}less(t){let e=typeof t=="number"?t:t._storage,n=ae(this._storage,e);return r._fromStorage(n)}less_equal(t){let e=typeof t=="number"?t:t._storage,n=se(this._storage,e);return r._fromStorage(n)}equal(t){let e=typeof t=="number"?t:t._storage,n=ie(this._storage,e);return r._fromStorage(n)}not_equal(t){let e=typeof t=="number"?t:t._storage,n=ue(this._storage,e);return r._fromStorage(n)}isclose(t,e=1e-5,n=1e-8){let o=typeof t=="number"?t:t._storage,a=rt(this._storage,o,e,n);return r._fromStorage(a)}allclose(t,e=1e-5,n=1e-8){let o=typeof t=="number"?t:t._storage;return ce(this._storage,o,e,n)}sum(t,e=!1){let n=ur(this._storage,t,e);return typeof n=="number"?n:r._fromStorage(n)}mean(t,e=!1){let n=wr(this._storage,t,e);return typeof n=="number"?n:r._fromStorage(n)}max(t,e=!1){let n=xr(this._storage,t,e);return typeof n=="number"?n:r._fromStorage(n)}min(t,e=!1){let n=Ir(this._storage,t,e);return typeof n=="number"?n:r._fromStorage(n)}prod(t,e=!1){let n=tt(this._storage,t,e);return typeof n=="number"?n:r._fromStorage(n)}argmin(t){let e=et(this._storage,t);return typeof e=="number"?e:r._fromStorage(e)}argmax(t){let e=nt(this._storage,t);return typeof e=="number"?e:r._fromStorage(e)}var(t,e=0,n=!1){let o=ot(this._storage,t,e,n);return typeof o=="number"?o:r._fromStorage(o)}std(t,e=0,n=!1){let o=ye(this._storage,t,e,n);return typeof o=="number"?o:r._fromStorage(o)}all(t,e=!1){let n=at(this._storage,t,e);return typeof n=="boolean"?n:r._fromStorage(n)}any(t,e=!1){let n=st(this._storage,t,e);return typeof n=="boolean"?n:r._fromStorage(n)}cumsum(t){return r._fromStorage(it(this._storage,t))}cumprod(t){return r._fromStorage(ut(this._storage,t))}ptp(t,e=!1){let n=ct(this._storage,t,e);return typeof n=="number"?n:r._fromStorage(n)}median(t,e=!1){let n=ft(this._storage,t,e);return typeof n=="number"?n:r._fromStorage(n)}percentile(t,e,n=!1){let o=lt(this._storage,t,e,n);return typeof o=="number"?o:r._fromStorage(o)}quantile(t,e,n=!1){let o=nr(this._storage,t,e,n);return typeof o=="number"?o:r._fromStorage(o)}average(t,e){let n=zr(this._storage,e,t?.storage);return typeof n=="number"?n:r._fromStorage(n)}nansum(t,e=!1){let n=Br(this._storage,t,e);return typeof n=="number"?n:r._fromStorage(n)}nanprod(t,e=!1){let n=_r(this._storage,t,e);return typeof n=="number"?n:r._fromStorage(n)}nanmean(t,e=!1){let n=Er(this._storage,t,e);return typeof n=="number"?n:r._fromStorage(n)}nanvar(t,e=0,n=!1){let o=cr(this._storage,t,e,n);return typeof o=="number"?o:r._fromStorage(o)}nanstd(t,e=0,n=!1){let o=pt(this._storage,t,e,n);return typeof o=="number"?o:r._fromStorage(o)}nanmin(t,e=!1){let n=Ur(this._storage,t,e);return typeof n=="number"?n:r._fromStorage(n)}nanmax(t,e=!1){let n=Tr(this._storage,t,e);return typeof n=="number"?n:r._fromStorage(n)}nanargmin(t){let e=Or(this._storage,t);return typeof e=="number"?e:r._fromStorage(e)}nanargmax(t){let e=$r(this._storage,t);return typeof e=="number"?e:r._fromStorage(e)}nancumsum(t){return r._fromStorage(yt(this._storage,t))}nancumprod(t){return r._fromStorage(mt(this._storage,t))}nanmedian(t,e=!1){let n=vr(this._storage,t,e);return typeof n=="number"?n:r._fromStorage(n)}reshape(...t){let e=t.length===1&&Array.isArray(t[0])?t[0]:t,n=R(this._storage,e),a=n.data===this.data?this._base??this:void 0;return r._fromStorage(n,a)}flatten(){let t=Mr(this._storage);return r._fromStorage(t)}ravel(){let t=lr(this._storage),n=t.data===this.data?this._base??this:void 0;return r._fromStorage(t,n)}transpose(t){let e=pr(this._storage,t),n=this._base??this;return r._fromStorage(e,n)}squeeze(t){let e=de(this._storage,t),n=this._base??this;return r._fromStorage(e,n)}expand_dims(t){let e=fr(this._storage,t),n=this._base??this;return r._fromStorage(e,n)}swapaxes(t,e){let n=ge(this._storage,t,e),o=this._base??this;return r._fromStorage(n,o)}moveaxis(t,e){let n=dt(this._storage,t,e),o=this._base??this;return r._fromStorage(n,o)}repeat(t,e){let n=Ie(this._storage,t,e);return r._fromStorage(n)}take(t,e){let n=mn(this._storage,t,e);return r._fromStorage(n)}put(t,e){let n=e instanceof r?e._storage:e;dn(this._storage,t,n)}matmul(t){let e=X(this._storage,t._storage);return r._fromStorage(e)}dot(t){let e=bt(this._storage,t._storage);return typeof e=="number"||typeof e=="bigint"?e:r._fromStorage(e)}trace(){return Pe(this._storage)}inner(t){let e=Ve(this._storage,t._storage);return typeof e=="number"||typeof e=="bigint"?e:r._fromStorage(e)}outer(t){let e=St(this._storage,t._storage);return r._fromStorage(e)}tensordot(t,e=2){let n=Fe(this._storage,t._storage,e);return typeof n=="number"||typeof n=="bigint"?n:r._fromStorage(n)}cbrt(){let t=Jt(this._storage);return r._fromStorage(t)}fabs(){let t=Qt(this._storage);return r._fromStorage(t)}divmod(t){let e=typeof t=="number"?t:t._storage,[n,o]=Xt(this._storage,e);return[r._fromStorage(n),r._fromStorage(o)]}square(){let t=Kt(this._storage);return r._fromStorage(t)}remainder(t){let e=typeof t=="number"?t:t._storage,n=re(this._storage,e);return r._fromStorage(n)}heaviside(t){let e=typeof t=="number"?t:t._storage,n=te(this._storage,e);return r._fromStorage(n)}slice(...t){if(t.length===0)return this;if(t.length>this.ndim)throw new Error(`Too many indices for array: array is ${this.ndim}-dimensional, but ${t.length} were indexed`);let e=t.map((u,c)=>{let l=Mt(u);return Ct(l,this.shape[c])});for(;e.length<this.ndim;)e.push({start:0,stop:this.shape[e.length],step:1,isIndex:!1});let n=[],o=[],a=this._storage.offset;for(let u=0;u<e.length;u++){let c=e[u],l=this._storage.strides[u];if(a+=c.start*l,!c.isIndex){let y;c.step>0?y=Math.max(0,Math.ceil((c.stop-c.start)/c.step)):y=Math.max(0,Math.ceil((c.start-c.stop)/Math.abs(c.step))),n.push(y),o.push(l*c.step)}}let i=D.fromData(this._storage.data,n,this._storage.dtype,o,a),s=this._base??this;return new r(i,s)}row(t){if(this.ndim<2)throw new Error("row() requires at least 2 dimensions");return this.slice(String(t),":")}col(t){if(this.ndim<2)throw new Error("col() requires at least 2 dimensions");return this.slice(":",String(t))}rows(t,e){if(this.ndim<2)throw new Error("rows() requires at least 2 dimensions");return this.slice(`${t}:${e}`,":")}cols(t,e){if(this.ndim<2)throw new Error("cols() requires at least 2 dimensions");return this.slice(":",`${t}:${e}`)}toString(){return`NDArray(shape=${JSON.stringify(this.shape)}, dtype=${this.dtype})`}toArray(){if(this.ndim===0)return this._storage.iget(0);let t=this.shape,e=t.length,n=(o,a)=>{if(a===e)return this._storage.get(...o);let i=[];for(let s=0;s<t[a];s++)o[a]=s,i.push(n(o,a+1));return i};return n(new Array(e),0)}}});function Et(){let r=new ArrayBuffer(2);return new DataView(r).setInt16(0,256,!0),new Int16Array(r)[0]===256}function Pa(r){if(r.startsWith("[")||r.startsWith("("))throw new j(`Structured/compound dtypes are not supported: ${r}`);let t="",e=r;(r[0]==="<"||r[0]===">"||r[0]==="="||r[0]==="|")&&(t=r[0],e=r.slice(1));let n=e[0];if(n&&n in ka)throw new j(`Unsupported dtype: ${ka[n]} (${r}). Use the 'force' parameter to skip arrays with unsupported dtypes.`);let o=Hs[e];if(!o)throw new j(`Unknown or unsupported dtype descriptor: ${r}. Supported types: ${kr.join(", ")}. Use the 'force' parameter to skip arrays with unsupported dtypes.`);let a=Et(),i=t==="<"||t==="|"||t==="="&&a,s=t===">"||t==="="&&!a,u=parseInt(e.slice(1),10),c=u>1&&(s&&a||i&&!a);return{dtype:o,needsByteSwap:c,itemsize:u}}var dr,kr,Hs,gr,ka,j,P,hr=$(()=>{"use strict";dr=new Uint8Array([147,78,85,77,80,89]),kr=["float64","float32","int64","int32","int16","int8","uint64","uint32","uint16","uint8","bool"];Hs={f8:"float64",f4:"float32",i8:"int64",i4:"int32",i2:"int16",i1:"int8",u8:"uint64",u4:"uint32",u2:"uint16",u1:"uint8",b1:"bool"},gr={float64:"<f8",float32:"<f4",int64:"<i8",int32:"<i4",int16:"<i2",int8:"|i1",uint64:"<u8",uint32:"<u4",uint16:"<u2",uint8:"|u1",bool:"|b1"},ka={c:"complex numbers",S:"byte strings",U:"Unicode strings",O:"Python objects",V:"structured arrays (void)",M:"datetime64",m:"timedelta64"};j=class extends Error{constructor(t){super(t),this.name="UnsupportedDTypeError"}},P=class extends Error{constructor(t){super(t),this.name="InvalidNpyError"}}});function Y(r){let t=r instanceof ArrayBuffer?new Uint8Array(r):r,e=Pr(t);return Vr(t,e)}function Pr(r){if(r.length<10)throw new P("File too small to be a valid NPY file");for(let c=0;c<dr.length;c++)if(r[c]!==dr[c])throw new P("Invalid NPY magic number");let t=r[6],e=r[7];if(t!==1&&t!==2&&t!==3)throw new P(`Unsupported NPY version: ${t}.${e}`);let n,o;t===1?(n=r[8]|r[9]<<8,o=10):(n=r[8]|r[9]<<8|r[10]<<16|r[11]<<24,o=12);let a=o+n;if(r.length<a)throw new P("File truncated: header extends beyond file");let i=r.slice(o,a),s=new TextDecoder("utf-8").decode(i).trim(),u=Ys(s);return{version:{major:t,minor:e},header:u,dataOffset:a}}function Vr(r,t){let{header:e,dataOffset:n}=t,{dtype:o,needsByteSwap:a,itemsize:i}=Pa(e.descr),s=e.shape.reduce((d,g)=>d*g,1),u=s*i,c=r.length-n;if(c<u)throw new P(`File truncated: expected ${u} bytes of data, got ${c}`);let l=new ArrayBuffer(u);new Uint8Array(l).set(r.subarray(n,n+u));let f=Js(l,o,s,a,i),p=e.shape,m;if(e.fortran_order&&p.length>1){let d=[...p].reverse(),g=D.fromData(f,d,o);m=Qs(g,d),p=e.shape}else m=D.fromData(f,[...p],o);return new I(m)}function Ys(r){let t=r.match(/'descr'\s*:\s*'([^']+)'/),e=r.match(/'fortran_order'\s*:\s*(True|False)/),n=r.match(/'shape'\s*:\s*\(([^)]*)\)/);if(!t||!e||!n)throw new P(`Failed to parse NPY header: ${r}`);let o=t[1],a=e[1]==="True",i=n[1].trim(),s;return i===""?s=[]:s=i.split(",").map(u=>u.trim()).filter(u=>u!=="").map(u=>{let c=parseInt(u,10);if(isNaN(c))throw new P(`Invalid shape value: ${u}`);return c}),{descr:o,fortran_order:a,shape:s}}function Js(r,t,e,n,o){let a=_(t);if(!a)throw new P(`Cannot create array for dtype: ${t}`);if(!n)return new a(r,0,e);let i=new Uint8Array(r),s=new Uint8Array(r.byteLength);for(let u=0;u<e;u++){let c=u*o;for(let l=0;l<o;l++)s[c+l]=i[c+o-1-l]}return new a(s.buffer,0,e)}function Qs(r,t){let e=t.length,n=r.size,o=r.dtype,a=_(o);if(!a)throw new P(`Cannot create array for dtype: ${o}`);let i=new a(n),s=[...t].reverse(),u=Va(t),c=Va(s),l=new Array(e).fill(0);for(let y=0;y<n;y++){let f=y;for(let m=0;m<e;m++){let d=u[m];l[m]=Math.floor(f/d),f=f%d}let p=0;for(let m=0;m<e;m++)p+=l[e-1-m]*c[m];N(o),i[p]=r.iget(y)}return D.fromData(i,s,o)}function Va(r){let t=new Array(r.length),e=1;for(let n=r.length-1;n>=0;n--)t[n]=e,e*=r[n];return t}var Fr=$(()=>{"use strict";_t();F();V();hr()});function W(r){let t=r.shape,e=r.dtype,n=gr[e],o=t.length===0?"()":t.length===1?`(${t[0]},)`:`(${t.join(", ")})`,a=`{'descr': '${n}', 'fortran_order': False, 'shape': ${o}, }`,i=12,u=(64-(i+a.length+1)%64)%64;a=a+" ".repeat(u)+`
+`;let c=new TextEncoder().encode(a),l=c.length,y=r.size,f=Vt(e),p=y*f,m=i+l+p,d=new Uint8Array(m);d.set(dr,0),d[6]=3,d[7]=0,d[8]=l&255,d[9]=l>>8&255,d[10]=l>>16&255,d[11]=l>>24&255,d.set(c,i);let g=i+l;return Xs(r,d.subarray(g),f),d}function Xs(r,t,e){let n=r.dtype,o=r.size,a=Et(),i=N(n),s=r._storage;if(s.isCContiguous&&s.offset===0&&a){let c=s.data,l=new Uint8Array(c.buffer,c.byteOffset,o*e);t.set(l)}else{let c=new DataView(t.buffer,t.byteOffset,t.byteLength);for(let l=0;l<o;l++){let y=s.iget(l),f=l*e;i?Ks(c,f,y,n==="uint64"):ri(c,f,y,n)}}}function Ks(r,t,e,n){n?r.setBigUint64(t,e,!0):r.setBigInt64(t,e,!0)}function ri(r,t,e,n){switch(n){case"float64":r.setFloat64(t,e,!0);break;case"float32":r.setFloat32(t,e,!0);break;case"int32":r.setInt32(t,e,!0);break;case"int16":r.setInt16(t,e,!0);break;case"int8":r.setInt8(t,e);break;case"uint32":r.setUint32(t,e,!0);break;case"uint16":r.setUint16(t,e,!0);break;case"uint8":case"bool":r.setUint8(t,e);break;default:throw new Error(`Unsupported dtype for serialization: ${n}`)}}var Rr=$(()=>{"use strict";V();hr()});function Ut(r){let t=4294967295;for(let e=0;e<r.length;e++)t=ti[(t^r[e])&255]^t>>>8;return(t^4294967295)>>>0}var ti,Tt=$(()=>{"use strict";ti=(()=>{let r=new Uint32Array(256);for(let t=0;t<256;t++){let e=t;for(let n=0;n<8;n++)e=e&1?3988292384^e>>>1:e>>>1;r[t]=e}return r})()});async function Fa(r){let t=qa(r),e=new Map;for(let n of t){let o=await ei(n);e.set(n.name,o)}return e}function Ra(r){let t=qa(r),e=new Map;for(let n of t){if(n.compressionMethod!==0)throw new Error(`Cannot read compressed entry synchronously: ${n.name}. Use readZip() (async) for DEFLATE-compressed files.`);e.set(n.name,n.compressedData)}return e}function qa(r){let t=r instanceof ArrayBuffer?new Uint8Array(r):r,e=new DataView(t.buffer,t.byteOffset,t.byteLength),n=[],o=-1;for(let c=t.length-22;c>=0;c--)if(e.getUint32(c,!0)===101010256){o=c;break}if(o===-1)throw new Error("Invalid ZIP file: end of central directory not found");let a=e.getUint32(o+16,!0),i=e.getUint16(o+10,!0),s=[],u=a;for(let c=0;c<i&&e.getUint32(u,!0)===33639248;c++){let y=e.getUint16(u+10,!0),f=e.getUint32(u+16,!0),p=e.getUint32(u+20,!0),m=e.getUint32(u+24,!0),d=e.getUint16(u+28,!0),g=e.getUint16(u+30,!0),A=e.getUint16(u+32,!0),h=e.getUint32(u+42,!0),b=t.slice(u+46,u+46+d),S=new TextDecoder("utf-8").decode(b);s.push({name:S,compressionMethod:y,crc32:f,compressedSize:p,uncompressedSize:m,localHeaderOffset:h}),u=u+46+d+g+A}for(let c of s){let l=c.localHeaderOffset;if(e.getUint32(l,!0)!==67324752)throw new Error(`Invalid local file header at offset ${l}`);let f=e.getUint16(l+26,!0),p=e.getUint16(l+28,!0),m=l+30+f+p,d=t.slice(m,m+c.compressedSize);n.push({name:c.name,compressedData:d,compressionMethod:c.compressionMethod,crc32:c.crc32,compressedSize:c.compressedSize,uncompressedSize:c.uncompressedSize})}return n}async function ei(r){if(r.compressionMethod===0)return r.compressedData;if(r.compressionMethod===8)return await ni(r.compressedData);throw new Error(`Unsupported compression method: ${r.compressionMethod}`)}async function ni(r){if(typeof DecompressionStream>"u")throw new Error("DecompressionStream is not available. This environment does not support the Compression Streams API. Please use a modern browser or Node.js 18+.");let t=new DecompressionStream("deflate-raw"),e=new Uint8Array(r.length);e.set(r);let n=t.writable.getWriter();n.write(e),n.close();let o=t.readable.getReader(),a=[];for(;;){let{done:c,value:l}=await o.read();if(c)break;a.push(l)}let i=a.reduce((c,l)=>c+l.length,0),s=new Uint8Array(i),u=0;for(let c of a)s.set(c,u),u+=c.length;return s}var La=$(()=>{"use strict";Tt()});var Za={};Hr(Za,{loadNpz:()=>qr,loadNpzSync:()=>Lr,parseNpz:()=>rr,parseNpzSync:()=>br});async function rr(r,t={}){let e=t.force??!1,n=await Fa(r);return ja(n,e)}function br(r,t={}){let e=t.force??!1,n=Ra(r);return ja(n,e)}function ja(r,t){let e=new Map,n=[],o=new Map;for(let[a,i]of r){if(!a.endsWith(".npy"))continue;let s=a.slice(0,-4);try{let u=Y(i);e.set(s,u)}catch(u){if(u instanceof j&&t)n.push(s),o.set(s,u.message);else throw u}}return{arrays:e,skipped:n,errors:o}}async function qr(r,t={}){let e=await rr(r,t);return Object.fromEntries(e.arrays)}function Lr(r,t={}){let e=br(r,t);return Object.fromEntries(e.arrays)}var jr=$(()=>{"use strict";Fr();hr();La()});async function Wa(r,t={}){let e=t.compress??!1,n=[];for(let[f,p]of r){let m=Ut(p),d,g;e?(d=await ai(p),d.length<p.length?g=8:(d=p,g=0)):(d=p,g=0),n.push({name:f,data:p,compressedData:d,crc:m,compressionMethod:g,offset:0})}let o=0;for(let f of n){let p=new TextEncoder().encode(f.name);o+=30+p.length+f.compressedData.length}let a=0;for(let f of n){let p=new TextEncoder().encode(f.name);a+=46+p.length}let s=o+a+22,u=new Uint8Array(s),c=new DataView(u.buffer),l=0;for(let f of n)f.offset=l,l=Ha(u,c,l,f);let y=l;for(let f of n)l=Ya(u,c,l,f);return Ja(c,l,n.length,a,y),u}function Ga(r){let t=[];for(let[l,y]of r){let f=Ut(y);t.push({name:l,data:y,compressedData:y,crc:f,compressionMethod:0,offset:0})}let e=0;for(let l of t){let y=new TextEncoder().encode(l.name);e+=30+y.length+l.compressedData.length}let n=0;for(let l of t){let y=new TextEncoder().encode(l.name);n+=46+y.length}let a=e+n+22,i=new Uint8Array(a),s=new DataView(i.buffer),u=0;for(let l of t)l.offset=u,u=Ha(i,s,u,l);let c=u;for(let l of t)u=Ya(i,s,u,l);return Ja(s,u,t.length,n,c),i}function Ha(r,t,e,n){let o=new TextEncoder().encode(n.name);return t.setUint32(e,67324752,!0),e+=4,t.setUint16(e,n.compressionMethod===8?20:10,!0),e+=2,t.setUint16(e,0,!0),e+=2,t.setUint16(e,n.compressionMethod,!0),e+=2,t.setUint16(e,0,!0),e+=2,t.setUint16(e,33,!0),e+=2,t.setUint32(e,n.crc,!0),e+=4,t.setUint32(e,n.compressedData.length,!0),e+=4,t.setUint32(e,n.data.length,!0),e+=4,t.setUint16(e,o.length,!0),e+=2,t.setUint16(e,0,!0),e+=2,r.set(o,e),e+=o.length,r.set(n.compressedData,e),e+=n.compressedData.length,e}function Ya(r,t,e,n){let o=new TextEncoder().encode(n.name);return t.setUint32(e,33639248,!0),e+=4,t.setUint16(e,20,!0),e+=2,t.setUint16(e,n.compressionMethod===8?20:10,!0),e+=2,t.setUint16(e,0,!0),e+=2,t.setUint16(e,n.compressionMethod,!0),e+=2,t.setUint16(e,0,!0),e+=2,t.setUint16(e,33,!0),e+=2,t.setUint32(e,n.crc,!0),e+=4,t.setUint32(e,n.compressedData.length,!0),e+=4,t.setUint32(e,n.data.length,!0),e+=4,t.setUint16(e,o.length,!0),e+=2,t.setUint16(e,0,!0),e+=2,t.setUint16(e,0,!0),e+=2,t.setUint16(e,0,!0),e+=2,t.setUint16(e,0,!0),e+=2,t.setUint32(e,0,!0),e+=4,t.setUint32(e,n.offset,!0),e+=4,r.set(o,e),e+=o.length,e}function Ja(r,t,e,n,o){r.setUint32(t,101010256,!0),t+=4,r.setUint16(t,0,!0),t+=2,r.setUint16(t,0,!0),t+=2,r.setUint16(t,e,!0),t+=2,r.setUint16(t,e,!0),t+=2,r.setUint32(t,n,!0),t+=4,r.setUint32(t,o,!0),t+=4,r.setUint16(t,0,!0)}async function ai(r){if(typeof CompressionStream>"u")throw new Error("CompressionStream is not available. This environment does not support the Compression Streams API. Please use a modern browser or Node.js 18+.");let t=new CompressionStream("deflate-raw"),e=new Uint8Array(r.length);e.set(r);let n=t.writable.getWriter();n.write(e),n.close();let o=t.readable.getReader(),a=[];for(;;){let{done:c,value:l}=await o.read();if(c)break;a.push(l)}let i=a.reduce((c,l)=>c+l.length,0),s=new Uint8Array(i),u=0;for(let c of a)s.set(c,u),u+=c.length;return s}var Qa=$(()=>{"use strict";Tt()});var Ka={};Hr(Ka,{serializeNpz:()=>ar,serializeNpzSync:()=>Zr});async function ar(r,t={}){let e=Xa(r);return Wa(e,{compress:t.compress??!1})}function Zr(r){let t=Xa(r);return Ga(t)}function Xa(r){let t=new Map;if(Array.isArray(r)){for(let n=0;n<r.length;n++){let o=r[n],a=W(o);t.set(`arr_${n}.npy`,a)}return t}let e=r instanceof Map?r.entries():Object.entries(r);for(let[n,o]of e){if(typeof n!="string"||n.length===0)throw new Error("Array names must be non-empty strings");let a=W(o),i=n.endsWith(".npy")?n:`${n}.npy`;t.set(i,a)}return t}var Wr=$(()=>{"use strict";Rr();Qa()});var mi={};Hr(mi,{DTYPE_TO_DESCR:()=>gr,InvalidNpyError:()=>P,NDArray:()=>I,SUPPORTED_DTYPES:()=>kr,UnsupportedDTypeError:()=>j,__version__:()=>si,absolute:()=>Wn,append:()=>Ko,arange:()=>Sn,arccos:()=>lo,arccosh:()=>xo,arcsin:()=>fo,arcsinh:()=>wo,arctan:()=>po,arctan2:()=>yo,arctanh:()=>Io,array:()=>k,array_equal:()=>ca,array_equiv:()=>fa,array_split:()=>Oo,asanyarray:()=>Un,asarray:()=>mr,ascontiguousarray:()=>Tn,asfortranarray:()=>On,atleast_1d:()=>Wo,atleast_2d:()=>Go,atleast_3d:()=>Ho,average:()=>ha,broadcast_arrays:()=>oa,broadcast_shapes:()=>aa,broadcast_to:()=>na,cbrt:()=>Ua,choose:()=>ua,column_stack:()=>Jo,concatenate:()=>K,copy:()=>In,cos:()=>uo,cosh:()=>Do,cumprod:()=>pa,cumsum:()=>la,deg2rad:()=>Ao,degrees:()=>go,delete:()=>ra,diag:()=>zt,diagflat:()=>$n,diagonal:()=>to,divmod:()=>Oa,dot:()=>Kn,dsplit:()=>Yo,dstack:()=>Uo,einsum:()=>Ca,empty:()=>xt,empty_like:()=>_n,expand_dims:()=>Fo,eye:()=>wt,fabs:()=>Ta,flip:()=>Cr,fliplr:()=>Ro,flipud:()=>qo,floor_divide:()=>Jn,frombuffer:()=>Fn,fromfile:()=>Rn,fromfunction:()=>vn,fromiter:()=>qn,fromstring:()=>Ln,full:()=>It,full_like:()=>En,geomspace:()=>wn,heaviside:()=>Ma,hsplit:()=>vo,hstack:()=>Eo,hypot:()=>mo,identity:()=>xn,inner:()=>oo,insert:()=>ta,kron:()=>eo,linspace:()=>Dn,load:()=>ui,loadNpy:()=>rs,loadNpySync:()=>ts,loadNpz:()=>qr,loadNpzFile:()=>os,loadNpzFileSync:()=>as,loadNpzSync:()=>Lr,loadSync:()=>ci,logspace:()=>Nn,median:()=>ma,meshgrid:()=>Mn,mod:()=>Yn,moveaxis:()=>Bo,nanargmax:()=>za,nanargmin:()=>Ia,nancumprod:()=>_a,nancumsum:()=>Ba,nanmax:()=>xa,nanmean:()=>Sa,nanmedian:()=>Ea,nanmin:()=>wa,nanprod:()=>ba,nanstd:()=>Na,nansum:()=>Aa,nanvar:()=>Da,negative:()=>Gn,ones:()=>Nt,ones_like:()=>Bn,outer:()=>ao,pad:()=>ea,parseNpy:()=>Y,parseNpyData:()=>Vr,parseNpyHeader:()=>Pr,parseNpz:()=>rr,parseNpzSync:()=>br,percentile:()=>da,positive:()=>Qn,power:()=>Zn,ptp:()=>ya,put:()=>ia,quantile:()=>ga,rad2deg:()=>bo,radians:()=>ho,ravel:()=>ko,reciprocal:()=>Xn,remainder:()=>va,repeat:()=>Co,reshape:()=>Po,resize:()=>Xo,roll:()=>jo,rollaxis:()=>Zo,rot90:()=>Lo,row_stack:()=>Qo,save:()=>fi,saveNpy:()=>es,saveNpySync:()=>ns,saveNpz:()=>vt,saveNpzSync:()=>ii,saveSync:()=>li,savez:()=>pi,savez_compressed:()=>yi,serializeNpy:()=>W,serializeNpz:()=>ar,serializeNpzSync:()=>Zr,sign:()=>Hn,sin:()=>io,sinh:()=>So,split:()=>To,sqrt:()=>jn,square:()=>$a,squeeze:()=>Vo,stack:()=>_o,swapaxes:()=>zo,take:()=>sa,tan:()=>co,tanh:()=>No,tensordot:()=>so,tile:()=>Mo,trace:()=>ro,transpose:()=>no,tri:()=>Cn,tril:()=>kn,triu:()=>Pn,vander:()=>Vn,vsplit:()=>$o,vstack:()=>Bt,zeros:()=>L,zeros_like:()=>zn});module.exports=Yr(mi);var sr=require("node:fs"),ir=require("node:fs/promises");Fr();Rr();jr();Wr();_t();Fr();Rr();hr();jr();Wr();var si="0.6.0";async function rs(r){let t=await(0,ir.readFile)(r);return Y(t)}function ts(r){let t=(0,sr.readFileSync)(r);return Y(t)}async function es(r,t){let e=W(t);await(0,ir.writeFile)(r,e)}function ns(r,t){let e=W(t);(0,sr.writeFileSync)(r,e)}async function os(r,t={}){let e=await(0,ir.readFile)(r);return rr(e,t)}function as(r,t={}){let e=(0,sr.readFileSync)(r),{parseNpzSync:n}=(jr(),Yr(Za));return n(e,t)}async function vt(r,t,e={}){let n=await ar(t,e);await(0,ir.writeFile)(r,n)}function ii(r,t){let{serializeNpzSync:e}=(Wr(),Yr(Ka)),n=e(t);(0,sr.writeFileSync)(r,n)}async function ui(r,t={}){if(r.endsWith(".npy")){if(t.allowNpy===!1)throw new Error("Loading .npy files is disabled (allowNpy: false)");return rs(r)}else{if(r.endsWith(".npz"))return os(r,t);throw new Error(`Unknown file extension. Expected .npy or .npz, got: ${r}`)}}function ci(r,t={}){if(r.endsWith(".npy")){if(t.allowNpy===!1)throw new Error("Loading .npy files is disabled (allowNpy: false)");return ts(r)}else{if(r.endsWith(".npz"))return as(r,t);throw new Error(`Unknown file extension. Expected .npy or .npz, got: ${r}`)}}async function fi(r,t){if(!r.endsWith(".npy"))throw new Error(`save() is for .npy files. Use saveNpz() for .npz files. Got: ${r}`);return es(r,t)}function li(r,t){if(!r.endsWith(".npy"))throw new Error(`saveSync() is for .npy files. Use saveNpzSync() for .npz files. Got: ${r}`);return ns(r,t)}async function pi(r,t){return r.endsWith(".npz")||(r=r+".npz"),vt(r,t,{compress:!1})}async function yi(r,t){return r.endsWith(".npz")||(r=r+".npz"),vt(r,t,{compress:!0})}0&&(module.exports={DTYPE_TO_DESCR,InvalidNpyError,NDArray,SUPPORTED_DTYPES,UnsupportedDTypeError,__version__,absolute,append,arange,arccos,arccosh,arcsin,arcsinh,arctan,arctan2,arctanh,array,array_equal,array_equiv,array_split,asanyarray,asarray,ascontiguousarray,asfortranarray,atleast_1d,atleast_2d,atleast_3d,average,broadcast_arrays,broadcast_shapes,broadcast_to,cbrt,choose,column_stack,concatenate,copy,cos,cosh,cumprod,cumsum,deg2rad,degrees,delete:null,diag,diagflat,diagonal,divmod,dot,dsplit,dstack,einsum,empty,empty_like,expand_dims,eye,fabs,flip,fliplr,flipud,floor_divide,frombuffer,fromfile,fromfunction,fromiter,fromstring,full,full_like,geomspace,heaviside,hsplit,hstack,hypot,identity,inner,insert,kron,linspace,load,loadNpy,loadNpySync,loadNpz,loadNpzFile,loadNpzFileSync,loadNpzSync,loadSync,logspace,median,meshgrid,mod,moveaxis,nanargmax,nanargmin,nancumprod,nancumsum,nanmax,nanmean,nanmedian,nanmin,nanprod,nanstd,nansum,nanvar,negative,ones,ones_like,outer,pad,parseNpy,parseNpyData,parseNpyHeader,parseNpz,parseNpzSync,percentile,positive,power,ptp,put,quantile,rad2deg,radians,ravel,reciprocal,remainder,repeat,reshape,resize,roll,rollaxis,rot90,row_stack,save,saveNpy,saveNpySync,saveNpz,saveNpzSync,saveSync,savez,savez_compressed,serializeNpy,serializeNpz,serializeNpzSync,sign,sin,sinh,split,sqrt,square,squeeze,stack,swapaxes,take,tan,tanh,tensordot,tile,trace,transpose,tri,tril,triu,vander,vsplit,vstack,zeros,zeros_like});
+//# sourceMappingURL=numpy-ts.node-io.cjs.map