numpy-ts 0.13.1 → 1.1.0

package/dist/esm/common/wasm/reduce_argmax.js

@@ -0,0 +1 @@
+import{reduce_argmax_f32 as m,reduce_argmax_i64 as c,reduce_argmax_i32 as s,reduce_argmax_i16 as d,reduce_argmax_i8 as p,reduce_argmax_u64 as l,reduce_argmax_u32 as _,reduce_argmax_u16 as A,reduce_argmax_u8 as f}from"./bins/reduce_argmax.wasm";import{ensureMemory as g,resetAllocator as T,copyIn as E}from"./runtime";import{wasmConfig as x}from"./config";const b=64,R={float32:m,int64:c,uint64:l,int32:s,uint32:_,int16:d,uint16:A,int8:p,uint8:f},C={float64:Float64Array,float32:Float32Array,int64:BigInt64Array,uint64:BigUint64Array,int32:Int32Array,uint32:Uint32Array,int16:Int16Array,uint16:Uint16Array,int8:Int8Array,uint8:Uint8Array};function M(r){if(!r.isCContiguous)return null;const t=r.size;if(t<b*x.thresholdMultiplier)return null;const e=r.dtype,n=R[e],a=C[e];if(!n||!a)return null;const o=a.BYTES_PER_ELEMENT;g(t*o),T();const i=r.offset,u=r.data.subarray(i,i+t),y=E(u);return Number(n(y,t))}export{M as wasmReduceArgmax};

package/dist/esm/common/wasm/reduce_argmin.js

@@ -0,0 +1 @@
+import{reduce_argmin_f32 as m,reduce_argmin_i64 as c,reduce_argmin_i32 as s,reduce_argmin_i16 as d,reduce_argmin_i8 as p,reduce_argmin_u64 as l,reduce_argmin_u32 as _,reduce_argmin_u16 as A,reduce_argmin_u8 as f}from"./bins/reduce_argmin.wasm";import{ensureMemory as g,resetAllocator as T,copyIn as E}from"./runtime";import{wasmConfig as b}from"./config";const R=64,C={float32:m,int64:c,uint64:l,int32:s,uint32:_,int16:d,uint16:A,int8:p,uint8:f},P={float64:Float64Array,float32:Float32Array,int64:BigInt64Array,uint64:BigUint64Array,int32:Int32Array,uint32:Uint32Array,int16:Int16Array,uint16:Uint16Array,int8:Int8Array,uint8:Uint8Array};function w(r){if(!r.isCContiguous)return null;const n=r.size;if(n<R*b.thresholdMultiplier)return null;const t=r.dtype,e=C[t],i=P[t];if(!e||!i)return null;const o=i.BYTES_PER_ELEMENT;g(n*o),T();const a=r.offset,u=r.data.subarray(a,a+n),y=E(u);return Number(e(y,n))}export{w as wasmReduceArgmin};

package/dist/esm/common/wasm/reduce_count_nz.js

@@ -0,0 +1 @@
+import{reduce_count_nz_f64 as _,reduce_count_nz_f32 as d,reduce_count_nz_i64 as u,reduce_count_nz_i32 as a,reduce_count_nz_i16 as y,reduce_count_nz_i8 as c}from"./bins/reduce_count_nz.wasm";import{ensureMemory as A,resetAllocator as f,copyIn as m}from"./runtime";import{wasmConfig as T}from"./config";const E=64,b={float64:_,float32:d,int64:u,uint64:u,int32:a,uint32:a,int16:y,uint16:y,int8:c,uint8:c},z={float64:Float64Array,float32:Float32Array,int64:BigInt64Array,uint64:BigUint64Array,int32:Int32Array,uint32:Uint32Array,int16:Int16Array,uint16:Uint16Array,int8:Int8Array,uint8:Uint8Array};function B(r){if(!r.isCContiguous)return null;const t=r.size;if(t<E*T.thresholdMultiplier)return null;const n=r.dtype,e=b[n],o=z[n];if(!e||!o)return null;const s=o.BYTES_PER_ELEMENT;A(t*s),f();const i=r.offset,l=r.data.subarray(i,i+t),p=m(l);return Number(e(p,t))}export{B as wasmReduceCountNz};

package/dist/esm/common/wasm/reduce_max.js

@@ -0,0 +1 @@
+import{reduce_max_f32 as B,reduce_max_i64 as k,reduce_max_i32 as P,reduce_max_i16 as g,reduce_max_i8 as R,reduce_max_u64 as C,reduce_max_u32 as M,reduce_max_u16 as D,reduce_max_u8 as I,reduce_max_strided_f32 as S,reduce_max_strided_i64 as U,reduce_max_strided_i32 as F,reduce_max_strided_i16 as N,reduce_max_strided_i8 as h,reduce_max_strided_u64 as L,reduce_max_strided_u32 as O,reduce_max_strided_u16 as Y,reduce_max_strided_u8 as H}from"./bins/reduce_max.wasm";import{ensureMemory as l,resetAllocator as _,copyIn as A,alloc as v,copyOut as K}from"./runtime";import{ArrayStorage as j}from"../storage";import{wasmConfig as b}from"./config";const p=64,q={float32:B,int64:k,uint64:C,int32:P,uint32:M,int16:g,uint16:D,int8:R,uint8:I},T={float64:Float64Array,float32:Float32Array,int64:BigInt64Array,uint64:BigUint64Array,int32:Int32Array,uint32:Uint32Array,int16:Int16Array,uint16:Uint16Array,int8:Int8Array,uint8:Uint8Array};function Z(r){if(!r.isCContiguous)return null;const n=r.size;if(n<p*b.thresholdMultiplier)return null;const a=r.dtype,e=q[a],t=T[a];if(!e||!t)return null;const u=t.BYTES_PER_ELEMENT;l(n*u),_();const o=r.offset,i=r.data.subarray(o,o+n),y=A(i);return Number(e(y,n))}const G={float32:S,int64:U,uint64:L,int32:F,uint32:O,int16:N,uint16:Y,int8:h,uint8:H},J={float32:Float32Array,int64:BigInt64Array,uint64:BigUint64Array,int32:Int32Array,uint32:Uint32Array,int16:Int16Array,uint16:Uint16Array,int8:Int8Array,uint8:Uint8Array};function $(r,n,a,e){if(!r.isCContiguous)return null;const t=n*a*e;if(t<p*b.thresholdMultiplier)return null;const u=r.dtype,o=G[u],i=T[u],y=J[u];if(!o||!i||!y)return null;const s=i.BYTES_PER_ELEMENT,d=s,f=n*e;l(t*s+f*d),_();const m=r.offset,w=r.data.subarray(m,m+t),x=A(w),c=v(f*d);o(x,c,n,a,e);const E=K(c,f,y);return j.fromData(E,[f],u)}export{Z as wasmReduceMax,$ as wasmReduceMaxStrided};

package/dist/esm/common/wasm/reduce_mean.js

@@ -0,0 +1 @@
+import{reduce_mean_f64 as R,reduce_mean_f32 as g,reduce_mean_i64 as B,reduce_mean_i32 as M,reduce_mean_i16 as S,reduce_mean_i8 as C,reduce_mean_u64 as D,reduce_mean_u32 as w,reduce_mean_u16 as F,reduce_mean_u8 as N,reduce_mean_strided_f64 as k,reduce_mean_strided_f32 as I,reduce_mean_strided_i64 as L,reduce_mean_strided_i32 as h,reduce_mean_strided_i16 as Y,reduce_mean_strided_i8 as O,reduce_mean_strided_u64 as U,reduce_mean_strided_u32 as x,reduce_mean_strided_u16 as H,reduce_mean_strided_u8 as v}from"./bins/reduce_mean.wasm";import{ensureMemory as y,resetAllocator as l,copyIn as f,alloc as K,copyOut as j}from"./runtime";import{ArrayStorage as q}from"../storage";import{wasmConfig as p}from"./config";const A=64,G={float64:R,float32:g,int64:B,uint64:D,int32:M,uint32:w,int16:S,uint16:F,int8:C,uint8:N},E={float64:Float64Array,float32:Float32Array,int64:BigInt64Array,uint64:BigUint64Array,int32:Int32Array,uint32:Uint32Array,int16:Int16Array,uint16:Uint16Array,int8:Int8Array,uint8:Uint8Array};function Z(r){if(!r.isCContiguous)return null;const e=r.size;if(e<A*p.thresholdMultiplier)return null;const u=r.dtype,t=G[u],n=E[u];if(!t||!n)return null;const o=n.BYTES_PER_ELEMENT;y(e*o),l();const a=r.offset,i=r.data.subarray(a,a+e),s=f(i);return Number(t(s,e))}const J={float64:k,float32:I,int64:L,uint64:U,int32:h,uint32:x,int16:Y,uint16:H,int8:O,uint8:v};function $(r,e,u,t){if(!r.isCContiguous)return null;const n=e*u*t;if(n<A*p.thresholdMultiplier)return null;const o=r.dtype,a=J[o],i=E[o];if(!a||!i)return null;const s=i.BYTES_PER_ELEMENT,_=Float64Array.BYTES_PER_ELEMENT,d=e*t;y(n*s+d*_),l();const c=r.offset,b=r.data.subarray(c,c+n),T=f(b),m=K(d*_);a(T,m,e,u,t);const P=j(m,d,Float64Array);return q.fromData(P,[d],"float64")}export{Z as wasmReduceMean,$ as wasmReduceMeanStrided};

package/dist/esm/common/wasm/reduce_min.js

@@ -0,0 +1 @@
+import{reduce_min_f32 as k,reduce_min_i64 as P,reduce_min_i32 as g,reduce_min_i16 as R,reduce_min_i8 as C,reduce_min_u64 as M,reduce_min_u32 as D,reduce_min_u16 as I,reduce_min_u8 as S,reduce_min_strided_f32 as U,reduce_min_strided_i64 as F,reduce_min_strided_i32 as N,reduce_min_strided_i16 as h,reduce_min_strided_i8 as L,reduce_min_strided_u64 as O,reduce_min_strided_u32 as Y,reduce_min_strided_u16 as x,reduce_min_strided_u8 as H}from"./bins/reduce_min.wasm";import{ensureMemory as l,resetAllocator as _,copyIn as A,alloc as v,copyOut as K}from"./runtime";import{ArrayStorage as j}from"../storage";import{wasmConfig as b}from"./config";const p=64,q={float32:k,int64:P,uint64:M,int32:g,uint32:D,int16:R,uint16:I,int8:C,uint8:S},T={float64:Float64Array,float32:Float32Array,int64:BigInt64Array,uint64:BigUint64Array,int32:Int32Array,uint32:Uint32Array,int16:Int16Array,uint16:Uint16Array,int8:Int8Array,uint8:Uint8Array};function Z(r){if(!r.isCContiguous)return null;const n=r.size;if(n<p*b.thresholdMultiplier)return null;const a=r.dtype,e=q[a],t=T[a];if(!e||!t)return null;const u=t.BYTES_PER_ELEMENT;l(n*u),_();const i=r.offset,o=r.data.subarray(i,i+n),y=A(o);return Number(e(y,n))}const G={float32:U,int64:F,uint64:O,int32:N,uint32:Y,int16:h,uint16:x,int8:L,uint8:H},J={float32:Float32Array,int64:BigInt64Array,uint64:BigUint64Array,int32:Int32Array,uint32:Uint32Array,int16:Int16Array,uint16:Uint16Array,int8:Int8Array,uint8:Uint8Array};function $(r,n,a,e){if(!r.isCContiguous)return null;const t=n*a*e;if(t<p*b.thresholdMultiplier)return null;const u=r.dtype,i=G[u],o=T[u],y=J[u];if(!i||!o||!y)return null;const s=o.BYTES_PER_ELEMENT,d=s,f=n*e;l(t*s+f*d),_();const m=r.offset,w=r.data.subarray(m,m+t),E=A(w),c=v(f*d);i(E,c,n,a,e);const B=K(c,f,y);return j.fromData(B,[f],u)}export{Z as wasmReduceMin,$ as wasmReduceMinStrided};

package/dist/esm/common/wasm/reduce_nanmax.js

@@ -0,0 +1 @@
+import{reduce_nanmax_f64 as l,reduce_nanmax_f32 as p}from"./bins/reduce_nanmax.wasm";import{ensureMemory as m,resetAllocator as c,copyIn as f}from"./runtime";import{wasmConfig as i}from"./config";const d=64,A={float64:l,float32:p},T={float64:Float64Array,float32:Float32Array};function C(r){if(!r.isCContiguous)return null;const e=r.size;if(e<d*i.thresholdMultiplier)return null;const t=r.dtype,n=A[t],o=T[t];if(!n||!o)return null;const y=o.BYTES_PER_ELEMENT;m(e*y),c();const a=r.offset,u=r.data.subarray(a,a+e),s=f(u);return Number(n(s,e))}export{C as wasmReduceNanmax};

package/dist/esm/common/wasm/reduce_nanmin.js

@@ -0,0 +1 @@
+import{reduce_nanmin_f64 as l,reduce_nanmin_f32 as p}from"./bins/reduce_nanmin.wasm";import{ensureMemory as i,resetAllocator as m,copyIn as c}from"./runtime";import{wasmConfig as f}from"./config";const d=64,A={float64:l,float32:p},T={float64:Float64Array,float32:Float32Array};function C(r){if(!r.isCContiguous)return null;const e=r.size;if(e<d*f.thresholdMultiplier)return null;const t=r.dtype,n=A[t],o=T[t];if(!n||!o)return null;const y=o.BYTES_PER_ELEMENT;i(e*y),m();const a=r.offset,u=r.data.subarray(a,a+e),s=c(u);return Number(n(s,e))}export{C as wasmReduceNanmin};

package/dist/esm/common/wasm/reduce_nansum.js

@@ -0,0 +1 @@
+import{reduce_nansum_f64 as l,reduce_nansum_f32 as p}from"./bins/reduce_nansum.wasm";import{ensureMemory as m,resetAllocator as c,copyIn as f}from"./runtime";import{wasmConfig as i}from"./config";const d=64,A={float64:l,float32:p},T={float64:Float64Array,float32:Float32Array};function C(r){if(!r.isCContiguous)return null;const e=r.size;if(e<d*i.thresholdMultiplier)return null;const t=r.dtype,n=A[t],o=T[t];if(!n||!o)return null;const u=o.BYTES_PER_ELEMENT;m(e*u),c();const a=r.offset,s=r.data.subarray(a,a+e),y=f(s);return Number(n(y,e))}export{C as wasmReduceNansum};

package/dist/esm/common/wasm/reduce_prod.js

@@ -0,0 +1 @@
+import{reduce_prod_f64 as D,reduce_prod_f32 as R,reduce_prod_i64 as c,reduce_prod_i32 as A,reduce_prod_i16 as C,reduce_prod_i8 as I,reduce_prod_u16 as S,reduce_prod_u8 as F,reduce_prod_strided_f64 as M,reduce_prod_strided_f32 as O,reduce_prod_strided_i64 as U,reduce_prod_strided_i32 as N,reduce_prod_strided_i16 as h,reduce_prod_strided_i8 as L,reduce_prod_strided_u64 as Y,reduce_prod_strided_u32 as x,reduce_prod_strided_u16 as H,reduce_prod_strided_u8 as v}from"./bins/reduce_prod.wasm";import{ensureMemory as _,resetAllocator as b,copyIn as m,alloc as K,copyOut as j}from"./runtime";import{ArrayStorage as q}from"../storage";import{wasmConfig as T}from"./config";const w=64,G={float64:D,float32:R,int64:c,uint64:c,int32:A,uint32:A,int16:C,uint16:S,int8:I,uint8:F},B={float64:Float64Array,float32:Float32Array,int64:BigInt64Array,uint64:BigUint64Array,int32:Int32Array,uint32:Uint32Array,int16:Int16Array,uint16:Uint16Array,int8:Int8Array,uint8:Uint8Array};function rr(r){if(!r.isCContiguous)return null;const n=r.size;if(n<w*T.thresholdMultiplier)return null;const o=r.dtype,e=G[o],u=B[o];if(!e||!u)return null;const t=u.BYTES_PER_ELEMENT;_(n*t),b();const a=r.offset,i=r.data.subarray(a,a+n),d=m(i);return Number(e(d,n))}const J={float64:M,float32:O,int64:U,uint64:Y,int32:N,uint32:x,int16:h,uint16:H,int8:L,uint8:v},Q={float64:"float64",float32:"float32",int64:"int64",uint64:"uint64",int32:"int32",uint32:"uint32",int16:"int64",uint16:"uint64",int8:"int64",uint8:"uint64"},V={float64:Float64Array,float32:Float32Array,int64:BigInt64Array,uint64:BigUint64Array,int32:Int32Array,uint32:Uint32Array,int16:BigInt64Array,uint16:BigUint64Array,int8:BigInt64Array,uint8:BigUint64Array},W={float64:8,float32:4,int64:8,uint64:8,int32:4,uint32:4,int16:8,uint16:8,int8:8,uint8:8};function nr(r,n,o,e){if(!r.isCContiguous)return null;const u=n*o*e;if(u<w*T.thresholdMultiplier)return null;const t=r.dtype,a=J[t],i=B[t],d=V[t],s=Q[t],y=W[t];if(!a||!i||!d||!s||!y)return null;const E=i.BYTES_PER_ELEMENT,f=n*e;_(u*E+f*y),b();const p=r.offset,P=r.data.subarray(p,p+u),g=m(P),l=K(f*y);a(g,l,n,o,e);const k=j(l,f,d);return q.fromData(k,[f],s)}export{rr as wasmReduceProd,nr as wasmReduceProdStrided};

package/dist/esm/common/wasm/reduce_quantile.js

@@ -0,0 +1 @@
+import{reduce_quantile_f64 as a}from"./bins/reduce_quantile.wasm";import{ensureMemory as f,resetAllocator as y,copyIn as p}from"./runtime";import{isBigIntDType as m,isComplexDType as c}from"../dtype";import{wasmConfig as d}from"./config";const A=64;function S(r,u){if(!r.isCContiguous||c(r.dtype))return null;const e=r.size;if(e<A*d.thresholdMultiplier)return null;const o=new Float64Array(e),n=r.offset,s=r.data;if(r.dtype==="float64")o.set(s.subarray(n,n+e));else if(m(r.dtype)){const t=s;for(let i=0;i<e;i++)o[i]=Number(t[n+i])}else for(let t=0;t<e;t++)o[t]=Number(s[n+t]);f(e*8),y();const l=p(o);return a(l,e,u)}export{S as wasmReduceQuantile};

package/dist/esm/common/wasm/reduce_std.js

@@ -0,0 +1 @@
+import{reduce_std_f64 as s,reduce_std_f32 as y,reduce_std_i64 as c,reduce_std_i32 as _,reduce_std_i16 as l,reduce_std_i8 as p,reduce_std_u64 as A,reduce_std_u32 as f,reduce_std_u16 as m,reduce_std_u8 as T}from"./bins/reduce_std.wasm";import{ensureMemory as E,resetAllocator as b,copyIn as g}from"./runtime";import{wasmConfig as R}from"./config";const S=64,C={float64:s,float32:y,int64:c,uint64:A,int32:_,uint32:f,int16:l,uint16:m,int8:p,uint8:T},P={float64:Float64Array,float32:Float32Array,int64:BigInt64Array,uint64:BigUint64Array,int32:Int32Array,uint32:Uint32Array,int16:Int16Array,uint16:Uint16Array,int8:Int8Array,uint8:Uint8Array};function F(r){if(!r.isCContiguous)return null;const t=r.size;if(t<S*R.thresholdMultiplier)return null;const e=r.dtype,n=C[e],o=P[e];if(!n||!o)return null;const i=o.BYTES_PER_ELEMENT;E(t*i),b();const u=r.offset,a=r.data.subarray(u,u+t),d=g(a);return Number(n(d,t))}export{F as wasmReduceStd};

package/dist/esm/common/wasm/reduce_sum.js

@@ -0,0 +1 @@
+import{reduce_sum_f64 as S,reduce_sum_f32 as g,reduce_sum_i64 as m,reduce_sum_i32 as _,reduce_sum_strided_f64 as B,reduce_sum_strided_f32 as C,reduce_sum_strided_i64 as D,reduce_sum_strided_i32 as w,reduce_sum_strided_i16 as M,reduce_sum_strided_i8 as F,reduce_sum_strided_u64 as k,reduce_sum_strided_u32 as I,reduce_sum_strided_u16 as N,reduce_sum_strided_u8 as h}from"./bins/reduce_sum.wasm";import{ensureMemory as l,resetAllocator as f,copyIn as p,alloc as L,copyOut as O}from"./runtime";import{ArrayStorage as U}from"../storage";import{wasmConfig as A}from"./config";const b=64,Y={float64:S,float32:g,int64:m,uint64:m,int32:_,uint32:_},E={float64:Float64Array,float32:Float32Array,int64:BigInt64Array,uint64:BigUint64Array,int32:Int32Array,uint32:Uint32Array,int16:Int16Array,uint16:Uint16Array,int8:Int8Array,uint8:Uint8Array};function q(r){if(!r.isCContiguous)return null;const t=r.size;if(t<b*A.thresholdMultiplier)return null;const u=r.dtype,e=Y[u],n=E[u];if(!e||!n)return null;const s=n.BYTES_PER_ELEMENT;l(t*s),f();const o=r.offset,a=r.data.subarray(o,o+t),d=p(a);return Number(e(d,t))}const x={float64:B,float32:C,int64:D,uint64:k,int32:w,uint32:I,int16:M,uint16:N,int8:F,uint8:h};function G(r,t,u,e){if(!r.isCContiguous)return null;const n=t*u*e;if(n<b*A.thresholdMultiplier)return null;const s=r.dtype,o=x[s],a=E[s];if(!o||!a)return null;const d=a.BYTES_PER_ELEMENT,i=t*e;l(n*d+i*8),f();const c=r.offset,T=r.data.subarray(c,c+n),P=p(T),y=L(i*8);o(P,y,t,u,e);const R=O(y,i,Float64Array);return U.fromData(R,[i],"float64")}export{q as wasmReduceSum,G as wasmReduceSumStrided};

package/dist/esm/common/wasm/reduce_var.js

@@ -0,0 +1 @@
+import{reduce_var_f64 as c,reduce_var_f32 as s,reduce_var_i64 as d,reduce_var_i32 as _,reduce_var_i16 as l,reduce_var_i8 as p,reduce_var_u64 as A,reduce_var_u32 as f,reduce_var_u16 as m,reduce_var_u8 as T}from"./bins/reduce_var.wasm";import{ensureMemory as E,resetAllocator as b,copyIn as v}from"./runtime";import{wasmConfig as g}from"./config";const R=64,C={float64:c,float32:s,int64:d,uint64:A,int32:_,uint32:f,int16:l,uint16:m,int8:p,uint8:T},P={float64:Float64Array,float32:Float32Array,int64:BigInt64Array,uint64:BigUint64Array,int32:Int32Array,uint32:Uint32Array,int16:Int16Array,uint16:Uint16Array,int8:Int8Array,uint8:Uint8Array};function w(r){if(!r.isCContiguous)return null;const t=r.size;if(t<R*g.thresholdMultiplier)return null;const e=r.dtype,n=C[e],a=P[e];if(!n||!a)return null;const u=a.BYTES_PER_ELEMENT;E(t*u),b();const o=r.offset,i=r.data.subarray(o,o+t),y=v(i);return Number(n(y,t))}export{w as wasmReduceVar};

package/dist/esm/common/wasm/repeat.js

@@ -0,0 +1 @@
+import{repeat_f64 as E,repeat_f32 as _,repeat_i64 as l,repeat_i32 as f,repeat_i16 as c,repeat_i8 as A}from"./bins/repeat.wasm";import{ensureMemory as g,resetAllocator as B,copyIn as P,alloc as R,copyOut as S}from"./runtime";import{ArrayStorage as D}from"../storage";import{wasmConfig as w}from"./config";const C=64,I={float64:E,float32:_,int64:l,uint64:l,int32:f,uint32:f,int16:c,uint16:c,int8:A,uint8:A},M={float64:Float64Array,float32:Float32Array,int64:BigInt64Array,uint64:BigUint64Array,int32:Int32Array,uint32:Uint32Array,int16:Int16Array,uint16:Uint16Array,int8:Int8Array,uint8:Uint8Array};function z(t,a){if(!t.isCContiguous)return null;const r=t.size;if(r<C*w.thresholdMultiplier)return null;const n=t.dtype,i=I[n],e=M[n];if(!i||!e)return null;const o=r*a,y=e.BYTES_PER_ELEMENT,m=r*y,u=o*y;g(m+u),B();const p=t.offset,d=t.data.subarray(p,p+r),T=P(d),s=R(u);i(T,s,r,a);const b=S(s,o,e);return D.fromData(b,[o],n)}export{z as wasmRepeat};

package/dist/esm/common/wasm/right_shift.js

@@ -0,0 +1 @@
+import{right_shift_i64 as C,right_shift_i32 as w,right_shift_i16 as R,right_shift_i8 as M,right_shift_u64 as O,right_shift_u32 as k,right_shift_u16 as N,right_shift_u8 as I,right_shift_scalar_i64 as L,right_shift_scalar_i32 as z,right_shift_scalar_i16 as F,right_shift_scalar_i8 as U,right_shift_scalar_u64 as Y,right_shift_scalar_u32 as v,right_shift_scalar_u16 as x,right_shift_scalar_u8 as H}from"./bins/right_shift.wasm";import{ensureMemory as m,resetAllocator as g,copyIn as h,alloc as A,copyOut as b}from"./runtime";import{ArrayStorage as d}from"../storage";import{promoteDTypes as K}from"../dtype";import{wasmConfig as T}from"./config";const E=64,j={int64:C,uint64:O,int32:w,uint32:k,int16:R,uint16:N,int8:M,uint8:I},q={int64:L,uint64:Y,int32:z,uint32:v,int16:F,uint16:x,int8:U,uint8:H},P={int64:BigInt64Array,uint64:BigUint64Array,int32:Int32Array,uint32:Uint32Array,int16:Int16Array,uint16:Uint16Array,int8:Int8Array,uint8:Uint8Array};function X(r,i){if(!r.isCContiguous||!i.isCContiguous)return null;const t=r.size;if(t<E*T.thresholdMultiplier)return null;const n=K(r.dtype,i.dtype),s=j[n],a=P[n];if(!s||!a)return null;const e=a.BYTES_PER_ELEMENT,l=t*e,u=t*e,o=t*e;m(l+u+o),g();const y=r.offset,f=i.offset,c=r.data.subarray(y,y+t),_=i.data.subarray(f,f+t),S=h(c),B=h(_),p=A(o);s(S,B,p,t);const D=b(p,t,a);return d.fromData(D,Array.from(r.shape),n)}function Z(r,i){if(!r.isCContiguous)return null;const t=r.size;if(t<E*T.thresholdMultiplier)return null;const n=r.dtype,s=q[n],a=P[n];if(!s||!a)return null;const e=a.BYTES_PER_ELEMENT,l=t*e,u=t*e;m(l+u),g();const o=r.offset,y=r.data.subarray(o,o+t),f=h(y),c=A(u);s(f,c,t,i);const _=b(c,t,a);return d.fromData(_,Array.from(r.shape),n)}export{X as wasmRightShift,Z as wasmRightShiftScalar};

package/dist/esm/common/wasm/roll.js

@@ -0,0 +1 @@
+import{roll_f64 as b,roll_f32 as E,roll_i64 as u,roll_i32 as s,roll_i16 as f,roll_i8 as p}from"./bins/roll.wasm";import{ensureMemory as _,resetAllocator as g,copyIn as B,alloc as P,copyOut as R}from"./runtime";import{ArrayStorage as D}from"../storage";import{wasmConfig as S}from"./config";const w=64,C={float64:b,float32:E,int64:u,uint64:u,int32:s,uint32:s,int16:f,uint16:f,int8:p,uint8:p},h={float64:Float64Array,float32:Float32Array,int64:BigInt64Array,uint64:BigUint64Array,int32:Int32Array,uint32:Uint32Array,int16:Int16Array,uint16:Uint16Array,int8:Int8Array,uint8:Uint8Array};function U(r,A){if(!r.isCContiguous)return null;const t=r.size;if(t<w*S.thresholdMultiplier)return null;const n=r.dtype,e=C[n],o=h[n];if(!e||!o)return null;const a=o.BYTES_PER_ELEMENT,m=t*a,i=t*a;_(m+i),g();const l=r.offset,c=r.data.subarray(l,l+t),d=B(c),y=P(i);e(d,y,t,A);const T=R(y,t,o);return D.fromData(T,Array.from(r.shape),n)}export{U as wasmRoll};

package/dist/esm/common/wasm/rot90.js

@@ -0,0 +1 @@
+import{rot90_f64 as b,rot90_f32 as _,rot90_i64 as p,rot90_i32 as c,rot90_i16 as f,rot90_i8 as A}from"./bins/rot90.wasm";import{ensureMemory as g,resetAllocator as w,copyIn as B,alloc as P,copyOut as R}from"./runtime";import{ArrayStorage as h}from"../storage";import{wasmConfig as D}from"./config";const S=64,C={float64:b,float32:_,int64:p,uint64:p,int32:c,uint32:c,int16:f,uint16:f,int8:A,uint8:A},I={float64:Float64Array,float32:Float32Array,int64:BigInt64Array,uint64:BigUint64Array,int32:Int32Array,uint32:Uint32Array,int16:Int16Array,uint16:Uint16Array,int8:Int8Array,uint8:Uint8Array};function L(t){if(!t.isCContiguous||t.shape.length!==2)return null;const r=t.size;if(r<S*D.thresholdMultiplier)return null;const n=t.dtype,e=C[n],o=I[n];if(!e||!o)return null;const a=t.shape[0],i=t.shape[1],s=o.BYTES_PER_ELEMENT,m=r*s,y=r*s;g(m+y),w();const u=t.offset,d=t.data.subarray(u,u+r),T=B(d),l=P(y);e(T,l,a,i);const E=R(l,r,o);return h.fromData(E,[i,a],n)}export{L as wasmRot90};

package/dist/esm/common/wasm/runtime.js

@@ -0,0 +1 @@
+import{wasmConfig as b}from"./config";let f=null,m=0,u=0;function y(){return f||(f=new WebAssembly.Memory({initial:17})),f}function a(e){const n=y(),t=u+e,o=n.buffer.byteLength;if(t>o){const r=Math.ceil((t-o)/65536);n.grow(r)}}function l(e=u){m=e,b.wasmCallCount++}function p(e){e>u&&(u=e)}function s(e){const n=m+7&-8;return m=n+e,n}function c(e){const n=s(e.byteLength),t=y();return new Uint8Array(t.buffer,n,e.byteLength).set(new Uint8Array(e.buffer,e.byteOffset,e.byteLength)),n}function d(e,n,t){const o=y(),r=new t(new ArrayBuffer(n*t.BYTES_PER_ELEMENT),0,n);return new Uint8Array(r.buffer,0,r.byteLength).set(new Uint8Array(o.buffer,e,r.byteLength)),r}export{s as alloc,c as copyIn,d as copyOut,a as ensureMemory,y as getSharedMemory,l as resetAllocator,p as setHeapBase};

package/dist/esm/common/wasm/searchsorted.js

@@ -0,0 +1 @@
+import{searchsorted_left_f64 as b,searchsorted_right_f64 as E,searchsorted_left_f32 as S,searchsorted_right_f32 as P,searchsorted_left_i64 as B,searchsorted_right_i64 as D,searchsorted_left_u64 as C,searchsorted_right_u64 as M,searchsorted_left_i32 as I,searchsorted_right_i32 as w,searchsorted_left_u32 as R,searchsorted_right_u32 as k,searchsorted_left_i16 as F,searchsorted_right_i16 as O,searchsorted_left_u16 as U,searchsorted_right_u16 as L,searchsorted_left_i8 as N,searchsorted_right_i8 as z,searchsorted_left_u8 as H,searchsorted_right_u8 as K}from"./bins/searchsorted.wasm";import{ensureMemory as Y,resetAllocator as v,copyIn as l,alloc as x,copyOut as j}from"./runtime";import{ArrayStorage as q}from"../storage";import{wasmConfig as G}from"./config";const J=64,Q={float64:b,float32:S,int64:B,uint64:C,int32:I,uint32:R,int16:F,uint16:U,int8:N,uint8:H},V={float64:E,float32:P,int64:D,uint64:M,int32:w,uint32:k,int16:O,uint16:L,int8:z,uint8:K},W={float64:Float64Array,float32:Float32Array,int64:BigInt64Array,uint64:BigUint64Array,int32:Int32Array,uint32:Uint32Array,int16:Int16Array,uint16:Uint16Array,int8:Int8Array,uint8:Uint8Array};function er(e,r,_){if(!e.isCContiguous||!r.isCContiguous)return null;const n=e.size,t=r.size;if(t<J*G.thresholdMultiplier)return null;const o=e.dtype,a=(_==="left"?Q:V)[o],s=W[o];if(!a||!s||r.dtype!==o)return null;const i=s.BYTES_PER_ELEMENT,h=n*i,d=t*i,c=t*4;Y(h+d+c),v();const u=e.offset,p=e.data.subarray(u,u+n),A=l(p),y=r.offset,m=r.data.subarray(y,y+t),g=l(m),f=x(c);a(A,n,g,f,t);const T=j(f,t,Int32Array);return q.fromData(T,Array.from(r.shape),"int32")}export{er as wasmSearchsorted};

package/dist/esm/common/wasm/sign.js

@@ -0,0 +1 @@
+import{sign_f64 as A,sign_f32 as c,sign_i64 as m,sign_i32 as g,sign_i16 as d,sign_i8 as T}from"./bins/sign.wasm";import{ensureMemory as E,resetAllocator as b,copyIn as _,alloc as B,copyOut as P}from"./runtime";import{ArrayStorage as S}from"../storage";import{wasmConfig as D}from"./config";const w=64,C={float64:A,float32:c,int64:m,int32:g,int16:d,int8:T},U={float64:Float64Array,float32:Float32Array,int64:BigInt64Array,uint64:BigUint64Array,int32:Int32Array,uint32:Uint32Array,int16:Int16Array,uint16:Uint16Array,int8:Int8Array,uint8:Uint8Array};function k(r){if(!r.isCContiguous)return null;const t=r.size;if(t<w*D.thresholdMultiplier)return null;const n=r.dtype,e=C[n],o=U[n];if(!e||!o)return null;const a=o.BYTES_PER_ELEMENT,u=t*a,i=t*a;E(u+i),b();const y=r.offset,f=r.data.subarray(y,y+t),l=_(f),s=B(i);e(l,s,t);const p=P(s,t,o);return S.fromData(p,Array.from(r.shape),n)}export{k as wasmSign};

package/dist/esm/common/wasm/sinh.js

@@ -0,0 +1 @@
+import{sinh_f64 as w,sinh_f32 as S,sinh_i64 as k,sinh_u64 as B}from"./bins/sinh.wasm";import{ensureMemory as y,resetAllocator as i,copyIn as l,alloc as c,copyOut as m}from"./runtime";import{ArrayStorage as p}from"../storage";import{isComplexDType as C,isBigIntDType as _}from"../dtype";import{wasmConfig as v}from"./config";const x=64,z={float64:w,float32:S};function R(n){if(!n.isCContiguous)return null;const r=n.size;if(r<x*v.thresholdMultiplier)return null;const o=n.dtype;if(C(o))return null;const A=z[o];if(A){const t=o==="float32",a=t?4:8,f=t?Float32Array:Float64Array;y(r*a*2),i();const e=n.offset,u=n.data.subarray(e,e+r),F=l(u),T=c(r*a);A(F,T,r);const O=m(T,r,f);return p.fromData(O,Array.from(n.shape),o)}if(o==="int64"||o==="uint64"){y(r*16),i();const t=n.offset,a=n.data.subarray(t,t+r),f=l(a),e=c(r*8);(o==="int64"?k:B)(f,e,r);const u=m(e,r,Float64Array);return p.fromData(u,Array.from(n.shape),"float64")}const b=8;y(r*b*2),i();const d=n.offset,h=n.data,s=new Float64Array(r);if(_(o))for(let t=0;t<r;t++)s[t]=Number(h[d+t]);else for(let t=0;t<r;t++)s[t]=h[d+t];const g=l(s),D=c(r*b);w(g,D,r);const P=m(D,r,Float64Array);return p.fromData(P,Array.from(n.shape),"float64")}export{R as wasmSinh};

package/dist/esm/common/wasm/sort.js

@@ -0,0 +1 @@
+import{sort_f64 as P,sort_f32 as h,sort_i64 as x,sort_u64 as L,sort_i32 as C,sort_u32 as M,sort_i16 as R,sort_u16 as F,sort_i8 as U,sort_u8 as k,sort_slices_f64 as N,sort_slices_f32 as I,sort_slices_i64 as Y,sort_slices_u64 as O,sort_slices_i32 as v,sort_slices_u32 as H,sort_slices_i16 as K,sort_slices_u16 as j,sort_slices_i8 as q,sort_slices_u8 as G,sort_c128 as J,sort_c64 as Q,sort_slices_c128 as V,sort_slices_c64 as W}from"./bins/sort.wasm";import{ensureMemory as f,resetAllocator as _,copyIn as p,copyOut as X,getSharedMemory as T}from"./runtime";import{ArrayStorage as Z}from"../storage";import{wasmConfig as $}from"./config";const z=64,S={float64:P,float32:h,int64:x,uint64:L,int32:C,uint32:M,int16:R,uint16:F,int8:U,uint8:k,complex128:J,complex64:Q},D={float64:N,float32:I,int64:Y,uint64:O,int32:v,uint32:H,int16:K,uint16:j,int8:q,uint8:G,complex128:V,complex64:W},A={float64:Float64Array,float32:Float32Array,int64:BigInt64Array,uint64:BigUint64Array,int32:Int32Array,uint32:Uint32Array,int16:Int16Array,uint16:Uint16Array,int8:Int8Array,uint8:Uint8Array,complex128:Float64Array,complex64:Float32Array};function or(r,e,t,o,n){if(t<2)return!0;const y=D[n];if(y&&e[0]===0&&o>1&&e[1]===t){const i=A[n];if(!i)return!1;const d=i.BYTES_PER_ELEMENT,w=r.length*d;f(w),_();const E=p(r);y(E,t,o);const B=T();return new Uint8Array(r.buffer,r.byteOffset,r.byteLength).set(new Uint8Array(B.buffer,E,r.byteLength)),!0}const l=S[n],c=A[n];if(!l||!c)return!1;const b=n==="complex128"||n==="complex64",s=c.BYTES_PER_ELEMENT,m=b?s*2:s,a=r.length*s;f(a),_();const u=p(r);for(let i=0;i<o;i++)l(u+e[i]*m,t);const g=T();return new Uint8Array(r.buffer,r.byteOffset,r.byteLength).set(new Uint8Array(g.buffer,u,r.byteLength)),!0}function sr(r){if(!r.isCContiguous)return null;const e=r.size;if(e<z*$.thresholdMultiplier)return null;const t=r.dtype,o=S[t],n=A[t];if(!o||!n)return null;const y=n.BYTES_PER_ELEMENT,l=e*y;f(l),_();const c=r.offset,s=t==="complex128"||t==="complex64"?e*2:e,m=r.data.subarray(c,c+s),a=p(m);o(a,e);const u=X(a,s,n);return Z.fromData(u,Array.from(r.shape),t)}export{sr as wasmSort,or as wasmSortSlices};

package/dist/esm/common/wasm/sqrt.js

@@ -0,0 +1 @@
+import{sqrt_f64 as B,sqrt_f32 as g,sqrt_i64 as p,sqrt_i32 as l,sqrt_i16 as A,sqrt_i8 as c}from"./bins/sqrt.wasm";import{ensureMemory as m,resetAllocator as E,copyIn as T,alloc as b,copyOut as d}from"./runtime";import{ArrayStorage as _}from"../storage";import{isComplexDType as C}from"../dtype";import{wasmConfig as S}from"./config";const w=64,R={float64:B,float32:g},h={int64:p,uint64:p,int32:l,uint32:l,int16:A,uint16:A,int8:c,uint8:c},P={float64:Float64Array,float32:Float32Array,int64:BigInt64Array,uint64:BigUint64Array,int32:Int32Array,uint32:Uint32Array,int16:Int16Array,uint16:Uint16Array,int8:Int8Array,uint8:Uint8Array};function L(t){if(!t.isCContiguous)return null;const r=t.size;if(r<w*S.thresholdMultiplier)return null;const n=t.dtype;if(C(n))return null;const i=R[n];if(i){const y=P[n],o=y.BYTES_PER_ELEMENT;m(r*o*2),E();const e=t.data.subarray(t.offset,t.offset+r),u=T(e),a=b(r*o);i(u,a,r);const s=d(a,r,y);return _.fromData(s,Array.from(t.shape),n)}const f=h[n];if(f){const o=P[n].BYTES_PER_ELEMENT,e=8;m(r*o+r*e),E();const u=t.data.subarray(t.offset,t.offset+r),a=T(u),s=b(r*e);f(a,s,r);const D=d(s,r,Float64Array);return _.fromData(D,Array.from(t.shape),"float64")}return null}export{L as wasmSqrt};

package/dist/esm/common/wasm/square.js

@@ -0,0 +1 @@
+import{square_f64 as A,square_f32 as T,square_i64 as s,square_i32 as l,square_i16 as p,square_i8 as c,square_c128 as d,square_c64 as _}from"./bins/square.wasm";import{ensureMemory as E,resetAllocator as b,copyIn as q,alloc as g,copyOut as P}from"./runtime";import{ArrayStorage as S}from"../storage";import{isComplexDType as x}from"../dtype";import{wasmConfig as C}from"./config";const D=64,F={float64:A,float32:T,complex128:d,complex64:_,int64:s,uint64:s,int32:l,uint32:l,int16:p,uint16:p,int8:c,uint8:c},w={float64:Float64Array,float32:Float32Array,complex128:Float64Array,complex64:Float32Array,int64:BigInt64Array,uint64:BigUint64Array,int32:Int32Array,uint32:Uint32Array,int16:Int16Array,uint16:Uint16Array,int8:Int8Array,uint8:Uint8Array};function h(r){if(!r.isCContiguous)return null;const t=r.size;if(t<D*C.thresholdMultiplier)return null;const n=r.dtype,y=F[n],e=w[n];if(!y||!e)return null;const o=x(n)?2:1,i=e.BYTES_PER_ELEMENT,a=t*o;E(a*i*2),b();const m=q(r.data.subarray(r.offset*o,(r.offset+t)*o)),u=g(a*i);y(m,u,t);const f=P(u,a,e);return S.fromData(f,Array.from(r.shape),n)}export{h as wasmSquare};

package/dist/esm/common/wasm/sub.js

@@ -0,0 +1 @@
+import{sub_f64 as z,sub_f32 as U,sub_i64 as d,sub_i32 as E,sub_i16 as T,sub_i8 as P,sub_c128 as Y,sub_c64 as v,sub_scalar_f64 as H,sub_scalar_f32 as K,sub_scalar_i64 as S,sub_scalar_i32 as B,sub_scalar_i16 as g,sub_scalar_i8 as D,sub_scalar_c128 as j,sub_scalar_c64 as q}from"./bins/sub.wasm";import{ensureMemory as x,resetAllocator as C,copyIn as A,alloc as w,copyOut as h}from"./runtime";import{ArrayStorage as F}from"../storage";import{promoteDTypes as G}from"../dtype";import{wasmConfig as R}from"./config";const M=64,J={float64:z,float32:U,int64:d,uint64:d,int32:E,uint32:E,int16:T,uint16:T,int8:P,uint8:P,complex128:Y,complex64:v},Q={float64:H,float32:K,int64:S,uint64:S,int32:B,uint32:B,int16:g,uint16:g,int8:D,uint8:D,complex128:j,complex64:q},O={float64:Float64Array,float32:Float32Array,complex128:Float64Array,complex64:Float32Array,int64:BigInt64Array,uint64:BigUint64Array,int32:Int32Array,uint32:Uint32Array,int16:Int16Array,uint16:Uint16Array,int8:Int8Array,uint8:Uint8Array},k={complex128:2,complex64:2};function tt(t,e){if(!t.isCContiguous||!e.isCContiguous)return null;const a=t.size;if(a<M*R.thresholdMultiplier)return null;const n=G(t.dtype,e.dtype),c=J[n],o=O[n];if(!c||!o)return null;const s=k[n]??1,u=o.BYTES_PER_ELEMENT,r=a*s,b=r*u,y=r*u,i=r*u;x(b+y+i),C();const l=t.offset*s,f=e.offset*s,p=t.data.subarray(l,l+r),m=e.data.subarray(f,f+r),N=A(p),I=A(m),_=w(i);c(N,I,_,a);const L=h(_,r,o);return F.fromData(L,Array.from(t.shape),n)}function rt(t,e){if(!t.isCContiguous)return null;const a=t.size;if(a<M*R.thresholdMultiplier)return null;const n=t.dtype,c=Q[n],o=O[n];if(!c||!o)return null;const s=k[n]??1,u=o.BYTES_PER_ELEMENT,r=a*s,b=r*u,y=r*u;x(b+y),C();const i=t.offset*s,l=t.data.subarray(i,i+r),f=A(l),p=w(y);c(f,p,a,e);const m=h(p,r,o);return F.fromData(m,Array.from(t.shape),n)}export{tt as wasmSub,rt as wasmSubScalar};

package/dist/esm/common/wasm/svd.js

@@ -0,0 +1 @@
+import{svd_f64 as E}from"./bins/svd.wasm";import{ensureMemory as H,resetAllocator as O,copyIn as _,alloc as n,copyOut as i}from"./runtime";import{ArrayStorage as m}from"../storage";import{wasmConfig as v}from"./config";const D=4;function L(r){if(r.ndim!==2)return null;const o=r.shape[0],t=r.shape[1];if(o<D*v.thresholdMultiplier||t<D*v.thresholdMultiplier)return null;const u=Math.min(o,t),f=o*t,e=o*o,c=u,l=t*t,S=o*t+t*t,h=(f+e+c+l+S)*8;H(h),O();const y=new Float64Array(f);for(let a=0;a<o;a++)for(let s=0;s<t;s++)y[a*t+s]=Number(r.get(a,s));const F=_(y),p=n(e*8),A=n(c*8),g=n(l*8),d=n(S*8);E(F,p,A,g,d,o,t);const w=i(p,e,Float64Array),z=i(A,c,Float64Array),P=i(g,l,Float64Array),M=m.fromData(w,[o,o],"float64"),k=m.fromData(z,[u],"float64"),B=m.fromData(P,[t,t],"float64");return{u:M,s:k,vt:B}}export{L as wasmSvd};

package/dist/esm/common/wasm/tan.js

@@ -0,0 +1 @@
+import{tan_f64 as i,tan_f32 as C}from"./bins/tan.wasm";import{ensureMemory as m,resetAllocator as A,copyIn as b,alloc as d,copyOut as T}from"./runtime";import{ArrayStorage as D}from"../storage";import{isComplexDType as S,isBigIntDType as k}from"../dtype";import{wasmConfig as B}from"./config";const v=64,_={float64:i,float32:C};function M(o){if(!o.isCContiguous)return null;const r=o.size;if(r<v*B.thresholdMultiplier)return null;const e=o.dtype;if(S(e))return null;const a=_[e];if(a){const t=e==="float32",l=t?4:8,F=t?Float32Array:Float64Array;m(r*l*2),A();const c=o.offset,P=o.data.subarray(c,c+r),h=b(P),p=d(r*l);a(h,p,r);const O=T(p,r,F);return D.fromData(O,Array.from(o.shape),e)}const s=8;m(r*s*2),A();const f=o.offset,y=o.data,n=new Float64Array(r);if(k(e))for(let t=0;t<r;t++)n[t]=Number(y[f+t]);else for(let t=0;t<r;t++)n[t]=y[f+t];const g=b(n),u=d(r*s);i(g,u,r);const w=T(u,r,Float64Array);return D.fromData(w,Array.from(o.shape),"float64")}export{M as wasmTan};

package/dist/esm/common/wasm/tanh.js

@@ -0,0 +1 @@
+import{tanh_f64 as w,tanh_f32 as k,tanh_i64 as B,tanh_u64 as C}from"./bins/tanh.wasm";import{ensureMemory as y,resetAllocator as l,copyIn as c,alloc as i,copyOut as m}from"./runtime";import{ArrayStorage as p}from"../storage";import{isComplexDType as S,isBigIntDType as _}from"../dtype";import{wasmConfig as v}from"./config";const x=64,z={float64:w,float32:k};function R(n){if(!n.isCContiguous)return null;const r=n.size;if(r<x*v.thresholdMultiplier)return null;const o=n.dtype;if(S(o))return null;const A=z[o];if(A){const t=o==="float32",e=t?4:8,f=t?Float32Array:Float64Array;y(r*e*2),l();const a=n.offset,u=n.data.subarray(a,a+r),F=c(u),T=i(r*e);A(F,T,r);const O=m(T,r,f);return p.fromData(O,Array.from(n.shape),o)}if(o==="int64"||o==="uint64"){y(r*16),l();const t=n.offset,e=n.data.subarray(t,t+r),f=c(e),a=i(r*8);(o==="int64"?B:C)(f,a,r);const u=m(a,r,Float64Array);return p.fromData(u,Array.from(n.shape),"float64")}const b=8;y(r*b*2),l();const d=n.offset,h=n.data,s=new Float64Array(r);if(_(o))for(let t=0;t<r;t++)s[t]=Number(h[d+t]);else for(let t=0;t<r;t++)s[t]=h[d+t];const g=c(s),D=i(r*b);w(g,D,r);const P=m(D,r,Float64Array);return p.fromData(P,Array.from(n.shape),"float64")}export{R as wasmTanh};

package/dist/esm/common/wasm/tile.js

@@ -0,0 +1 @@
+import{tile_2d_f64 as g,tile_2d_f32 as B,tile_2d_i64 as m,tile_2d_i32 as A,tile_2d_i16 as d,tile_2d_i8 as T}from"./bins/tile.wasm";import{ensureMemory as P,resetAllocator as S,copyIn as h,alloc as w,copyOut as C}from"./runtime";import{ArrayStorage as I}from"../storage";import{wasmConfig as M}from"./config";const k=64,F={float64:g,float32:B,int64:m,uint64:m,int32:A,uint32:A,int16:d,uint16:d,int8:T,uint8:T},R={float64:Float64Array,float32:Float32Array,int64:BigInt64Array,uint64:BigUint64Array,int32:Int32Array,uint32:Uint32Array,int16:Int16Array,uint16:Uint16Array,int8:Int8Array,uint8:Uint8Array};function H(t,r,n){if(!t.isCContiguous||t.shape.length!==2)return null;const e=t.size;if(e<k*M.thresholdMultiplier)return null;const o=t.dtype,y=F[o],i=R[o];if(!y||!i)return null;const a=t.shape[0],u=t.shape[1],s=a*r*u*n,l=i.BYTES_PER_ELEMENT,_=e*l,c=s*l;P(_+c),S();const f=t.offset,b=t.data.subarray(f,f+e),E=h(b),p=w(c);y(E,p,a,u,r,n);const D=C(p,s,i);return I.fromData(D,[a*r,u*n],o)}export{H as wasmTile2D};

package/dist/esm/common/wasm/vdot.js

@@ -0,0 +1 @@
+import{vdot_c128 as b,vdot_c64 as A}from"./bins/vdot.wasm";import{ensureMemory as E,resetAllocator as C,copyIn as p,alloc as D,copyOut as g}from"./runtime";import{promoteDTypes as w}from"../dtype";import{Complex as x}from"../complex";import{wasmConfig as P}from"./config";const B=64,S={complex128:b,complex64:A},_={complex128:Float64Array,complex64:Float32Array};function K(r,t){if(r.ndim!==1||t.ndim!==1||!r.isCContiguous||!t.isCContiguous)return null;const e=r.shape[0];if(e!==t.shape[0]||e<B*P.thresholdMultiplier)return null;const s=w(r.dtype,t.dtype),a=S[s],o=_[s];if(!a||!o)return null;const n=o.BYTES_PER_ELEMENT,y=e*2*n,c=e*2*n,u=2*n;E(y+c+u),C();const f=r.data.subarray(r.offset*2,r.offset*2+e*2),i=t.data.subarray(t.offset*2,t.offset*2+e*2),d=p(f),T=p(i),m=D(u);a(d,T,m,e);const l=g(m,2,o);return new x(Number(l[0]),Number(l[1]))}export{K as wasmVdotComplex};

package/dist/esm/common/wasm/vecdot.js

@@ -0,0 +1 @@
+import{vecdot_f64 as D,vecdot_f32 as v,vecdot_c128 as P,vecdot_c64 as h,vecdot_i64 as l,vecdot_i32 as p,vecdot_i16 as m,vecdot_i8 as f}from"./bins/vecdot.wasm";import{ensureMemory as x,resetAllocator as C,copyIn as d,alloc as S,copyOut as w}from"./runtime";import{ArrayStorage as F}from"../storage";import{promoteDTypes as R}from"../dtype";import{wasmConfig as I}from"./config";const M=256,U={float64:D,float32:v,complex128:P,complex64:h,int64:l,uint64:l,int32:p,uint32:p,int16:m,uint16:m,int8:f,uint8:f},k={float64:Float64Array,float32:Float32Array,complex128:Float64Array,complex64:Float32Array,int64:BigInt64Array,uint64:BigUint64Array,int32:Int32Array,uint32:Uint32Array,int16:Int16Array,uint16:Uint16Array,int8:Int8Array,uint8:Uint8Array},K={complex128:2,complex64:2};function W(e,o){if(e.ndim!==2||o.ndim!==2||!e.isCContiguous||!o.isCContiguous)return null;const t=e.shape[0],n=e.shape[1];if(t!==o.shape[0]||n!==o.shape[1]||t*n<M*I.thresholdMultiplier)return null;const a=R(e.dtype,o.dtype),c=U[a],s=k[a];if(!c||!s)return null;const r=K[a]??1,y=s.BYTES_PER_ELEMENT,A=t*n*r*y,T=t*n*r*y,i=t*r*y;x(A+T+i),C();const E=e.data.subarray(e.offset*r,e.offset*r+t*n*r),_=o.data.subarray(o.offset*r,o.offset*r+t*n*r),b=d(E),g=d(_),u=S(i);c(b,g,u,t,n);const B=w(u,t*r,s);return F.fromData(B,[t],a)}export{W as wasmVecdot};

package/dist/esm/common/wasm/vecmat.js

@@ -0,0 +1 @@
+import{vecmat_f64 as v,vecmat_f32 as B,vecmat_c128 as P,vecmat_c64 as C,vecmat_i64 as u,vecmat_i32 as l,vecmat_i16 as p,vecmat_i8 as f}from"./bins/vecmat.wasm";import{ensureMemory as S,resetAllocator as h,copyIn as A,alloc as w,copyOut as F}from"./runtime";import{ArrayStorage as x}from"../storage";import{promoteDTypes as R}from"../dtype";import{wasmConfig as I}from"./config";const M=128,N={float64:v,float32:B,complex128:P,complex64:C,int64:u,uint64:u,int32:l,uint32:l,int16:p,uint16:p,int8:f,uint8:f},U={float64:Float64Array,float32:Float32Array,complex128:Float64Array,complex64:Float32Array,int64:BigInt64Array,uint64:BigUint64Array,int32:Int32Array,uint32:Uint32Array,int16:Int16Array,uint16:Uint16Array,int8:Int8Array,uint8:Uint8Array},k={complex128:2,complex64:2};function W(e,r){if(e.ndim!==1||r.ndim!==2||!e.isCContiguous||!r.isCContiguous)return null;const n=r.shape[0],a=r.shape[1];if(n!==e.shape[0]||n*a<M*I.thresholdMultiplier)return null;const o=R(e.dtype,r.dtype),c=N[o],y=U[o];if(!c||!y)return null;const t=k[o]??1,s=y.BYTES_PER_ELEMENT,d=n*t*s,T=n*a*t*s,i=a*t*s;S(d+T+i),h();const b=e.data.subarray(e.offset*t,e.offset*t+n*t),E=r.data.subarray(r.offset*t,r.offset*t+n*a*t),_=A(b),g=A(E),m=w(i);c(_,g,m,n,a);const D=F(m,a*t,y);return x.fromData(D,[a],o)}export{W as wasmVecmat};

package/dist/esm/core/creation.js

@@ -1 +1 @@
-import{ArrayStorage as g}from"../common/storage";import{NDArrayCore as b}from"../common/ndarray-core";import{Complex as M,isComplexLike as U}from"../common/complex";import{getTypedArrayConstructor as w,isBigIntDType as C,isComplexDType as S,DEFAULT_DTYPE as m}from"../common/dtype";function B(r){return new b(r)}function z(r){return!!(r instanceof b||r&&typeof r=="object"&&"storage"in r&&r.storage instanceof g)}function N(r,e=m){const t=g.zeros(r,e);return new b(t)}function F(r,e=m){const t=g.ones(r,e);return new b(t)}function k(r,e=m){return N(r,e)}function v(r,e,t){let n=t;n||(typeof e=="bigint"?n="int64":typeof e=="boolean"?n="bool":Number.isInteger(e)?n="int32":n=m);const o=w(n);if(!o)throw new Error(`Cannot create full array with dtype ${n}`);const u=r.reduce((s,a)=>s*a,1),i=new o(u);if(C(n)){const s=typeof e=="bigint"?e:BigInt(Math.round(Number(e)));i.fill(s)}else n==="bool"?i.fill(e?1:0):i.fill(Number(e));const f=g.fromData(i,r,n);return new b(f)}function j(r){const e=[];let t=r;for(;Array.isArray(t);)e.push(t.length),t=t[0];return e}function E(r){return typeof r=="bigint"?!0:Array.isArray(r)?r.some(e=>E(e)):!1}function I(r){return U(r)?!0:Array.isArray(r)?r.some(e=>I(e)):!1}function $(r){const e=[];function t(n){Array.isArray(n)?n.forEach(o=>t(o)):e.push(n)}return t(r),e}function h(r,e){if(r instanceof b)return!e||r.dtype===e?r.copy():r.astype(e);const t=E(r),n=I(r),o=j(r),u=o.reduce((A,p)=>A*p,1);let i=e;i||(n?i="complex128":t?i="int64":i=m);const f=S(i),s=w(i);if(!s)throw new Error(`Cannot create array with dtype ${i}`);const a=f?u*2:u,c=new s(a),y=$(r);if(C(i)){const A=c;for(let p=0;p<u;p++){const D=y[p];A[p]=typeof D=="bigint"?D:BigInt(Math.round(Number(D)))}}else if(i==="bool"){const A=c;for(let p=0;p<u;p++)A[p]=y[p]?1:0}else if(f){const A=c;for(let p=0;p<u;p++){const D=y[p];let d,x;D instanceof M?(d=D.re,x=D.im):typeof D=="object"&&D!==null&&"re"in D?(d=D.re,x=D.im??0):(d=Number(D),x=0),A[p*2]=d,A[p*2+1]=x}}else{const A=c;for(let p=0;p<u;p++){const D=y[p];A[p]=Number(D)}}const l=g.fromData(c,o,i);return new b(l)}function K(r,e,t=1,n=m){let o=r,u=e;if(e===void 0&&(o=0,u=r),u===void 0)throw new Error("stop is required");const i=Math.max(0,Math.ceil((u-o)/t)),f=w(n);if(!f)throw new Error(`Cannot create arange array with dtype ${n}`);const s=new f(i);if(C(n))for(let c=0;c<i;c++)s[c]=BigInt(Math.round(o+c*t));else if(n==="bool")for(let c=0;c<i;c++)s[c]=o+c*t!==0?1:0;else for(let c=0;c<i;c++)s[c]=o+c*t;const a=g.fromData(s,[i],n);return new b(a)}function O(r,e,t=50,n=m){if(t<0)throw new Error("num must be non-negative");if(t===0)return h([],n);if(t===1)return h([r],n);const o=w(n);if(!o)throw new Error(`Cannot create linspace array with dtype ${n}`);const u=new o(t),i=(e-r)/(t-1);if(C(n))for(let s=0;s<t;s++)u[s]=BigInt(Math.round(r+s*i));else if(n==="bool")for(let s=0;s<t;s++)u[s]=r+s*i!==0?1:0;else for(let s=0;s<t;s++)u[s]=r+s*i;const f=g.fromData(u,[t],n);return new b(f)}function R(r,e,t=50,n=10,o=m){if(t<0)throw new Error("num must be non-negative");if(t===0)return h([],o);if(t===1)return h([Math.pow(n,r)],o);const u=w(o);if(!u)throw new Error(`Cannot create logspace array with dtype ${o}`);const i=new u(t),f=(e-r)/(t-1);if(C(o))for(let a=0;a<t;a++){const c=r+a*f;i[a]=BigInt(Math.round(Math.pow(n,c)))}else if(o==="bool")for(let a=0;a<t;a++){const c=r+a*f;i[a]=Math.pow(n,c)!==0?1:0}else for(let a=0;a<t;a++){const c=r+a*f;i[a]=Math.pow(n,c)}const s=g.fromData(i,[t],o);return new b(s)}function H(r,e,t=50,n=m){if(t<0)throw new Error("num must be non-negative");if(r===0||e===0)throw new Error("Geometric sequence cannot include zero");if(t===0)return h([],n);if(t===1)return h([r],n);const o=Math.sign(r),u=Math.sign(e);if(o!==u)throw new Error("Geometric sequence cannot contain both positive and negative values");const i=w(n);if(!i)throw new Error(`Cannot create geomspace array with dtype ${n}`);const f=new i(t),s=Math.log(Math.abs(r)),c=(Math.log(Math.abs(e))-s)/(t-1);if(C(n))for(let l=0;l<t;l++){const A=o*Math.exp(s+l*c);f[l]=BigInt(Math.round(A))}else if(n==="bool")for(let l=0;l<t;l++){const A=o*Math.exp(s+l*c);f[l]=A!==0?1:0}else for(let l=0;l<t;l++){const A=o*Math.exp(s+l*c);f[l]=A}const y=g.fromData(f,[t],n);return new b(y)}function L(r,e,t=0,n=m){const o=e??r,u=N([r,o],n),i=u.data;if(C(n)){const f=i;for(let s=0;s<r;s++){const a=s+t;a>=0&&a<o&&(f[s*o+a]=BigInt(1))}}else{const f=i;for(let s=0;s<r;s++){const a=s+t;a>=0&&a<o&&(f[s*o+a]=1)}}return u}function J(r,e=m){return L(r,r,0,e)}function T(r,e){return z(r)?!e||r.dtype===e?r:r.astype(e):h(r,e)}function Q(r,e){return N(Array.from(r.shape),e??r.dtype)}function W(r,e){return F(Array.from(r.shape),e??r.dtype)}function X(r,e){return k(Array.from(r.shape),e??r.dtype)}function Z(r,e,t){return v(Array.from(r.shape),e,t??r.dtype)}function rr(r){return r.copy()}function er(r,e){return T(r,e)}function tr(r,e){return T(r,e).copy()}function nr(r,e){return T(r,e).copy()}function or(r,e){const t=T(r,e),n=t.data;for(let o=0;o<n.length;o++){const u=n[o];if(!Number.isFinite(u))throw new Error("array must not contain infs or NaNs")}return t}function ar(r,e,t){let n=r;return e&&e!==r.dtype&&(n=n.astype(e)),n}function q(r){const e=r.data,t=g.fromData(e.slice(),[e.length],r.dtype);return B(t)}function _(r,e=0){const t=r.shape,n=r.data,o=r.dtype;if(t.length===1){const u=t[0],i=u+Math.abs(e),f=N([i,i],o),s=f.data;for(let a=0;a<u;a++){const c=e>=0?a:a-e,y=e>=0?a+e:a;c>=0&&c<i&&y>=0&&y<i&&(s[c*i+y]=n[a])}return f}else if(t.length===2){const u=t[0],i=t[1],f=Math.min(e>=0?Math.min(u,i-e):Math.min(u+e,i),Math.max(0,e>=0?i-e:u+e));if(f<=0)return h([],o);const s=[];for(let a=0;a<f;a++){const c=e>=0?a:a-e,y=e>=0?a+e:a;c>=0&&c<u&&y>=0&&y<i&&s.push(n[c*i+y])}return h(s,o)}throw new Error("Input must be 1-D or 2-D")}function sr(r,e=0){const t=q(r);return _(t,e)}function ir(r,e,t=0,n=m){const o=e??r,u=N([r,o],n),i=u.data;for(let f=0;f<r;f++)for(let s=0;s<=Math.min(f+t,o-1);s++)s>=0&&(i[f*o+s]=1);return u}function cr(r,e=0){const t=r.shape;if(t.length<2)throw new Error("Input must be at least 2-D");const n=r.copy(),o=n.data,u=t[t.length-2],i=t[t.length-1],f=t.slice(0,-2).reduce((a,c)=>a*c,1),s=u*i;for(let a=0;a<f;a++){const c=a*s;for(let y=0;y<u;y++)for(let l=y+e+1;l<i;l++)o[c+y*i+l]=0}return n}function ur(r,e=0){const t=r.shape;if(t.length<2)throw new Error("Input must be at least 2-D");const n=r.copy(),o=n.data,u=t[t.length-2],i=t[t.length-1],f=t.slice(0,-2).reduce((a,c)=>a*c,1),s=u*i;for(let a=0;a<f;a++){const c=a*s;for(let y=0;y<u;y++)for(let l=0;l<Math.min(y+e,i);l++)o[c+y*i+l]=0}return n}function fr(r,e,t=!1){const n=r.size,o=e??n,u=r.data,i=N([n,o],r.dtype),f=i.data;for(let s=0;s<n;s++){const a=u[s];for(let c=0;c<o;c++){const y=t?c:o-1-c;f[s*o+c]=Math.pow(a,y)}}return i}function lr(r,e=m,t=-1,n=0){const o=w(e);if(!o)throw new Error(`Unsupported dtype: ${e}`);let u;if(r instanceof ArrayBuffer){const f=n,s=t<0?(r.byteLength-f)/o.BYTES_PER_ELEMENT:t;u=new o(r,f,s)}else{const f=n,s=t<0?r.length:n+t,a=Array.from(r.slice(f,s));u=new o(a.length);for(let c=0;c<a.length;c++)u[c]=a[c]}const i=g.fromData(u,[u.length],e);return B(i)}function yr(r,e,t=m){const n=e.reduce((a,c)=>a*c,1),o=w(t);if(!o)throw new Error(`Unsupported dtype: ${t}`);const u=new o(n),i=[];let f=1;for(let a=e.length-1;a>=0;a--)i.unshift(f),f*=e[a];for(let a=0;a<n;a++){const c=[];let y=a;for(let l=0;l<e.length;l++)c.push(Math.floor(y/i[l])),y=y%i[l];u[a]=r(...c)}const s=g.fromData(u,e,t);return B(s)}function pr(r,e=m,t=-1){const n=[];let o=0;for(const u of r){if(t>=0&&o>=t)break;n.push(u),o++}return h(n,e)}function Ar(r,e=m,t=-1,n){const o=n??/\s+/,i=r.split(o).filter(s=>s.trim()!=="").map(s=>parseFloat(s.trim())),f=t>=0?i.slice(0,t):i;return h(f,e)}function Dr(r,e=m,t=-1,n=""){throw new Error("fromfile requires Node.js file system access")}function mr(...r){if(r.length===0)return[];if(r.length===1)return[r[0].copy()];const e=r.map(o=>o.size),t=[...e],n=[];for(let o=0;o<r.length;o++){const u=r[o],i=u.data,f=N(t,u.dtype),s=f.data,a=[];let c=1;for(let l=t.length-1;l>=0;l--)a.unshift(c),c*=t[l];const y=t.reduce((l,A)=>l*A,1);for(let l=0;l<y;l++){const A=Math.floor(l/a[o])%e[o];s[l]=i[A]}n.push(f)}return n}export{K as arange,h as array,er as asanyarray,T as asarray,or as asarray_chkfinite,tr as ascontiguousarray,nr as asfortranarray,rr as copy,_ as diag,sr as diagflat,k as empty,X as empty_like,L as eye,lr as frombuffer,Dr as fromfile,yr as fromfunction,pr as fromiter,Ar as fromstring,v as full,Z as full_like,H as geomspace,J as identity,O as linspace,R as logspace,mr as meshgrid,F as ones,W as ones_like,ar as require,ir as tri,cr as tril,ur as triu,fr as vander,N as zeros,Q as zeros_like};
+import{ArrayStorage as D}from"../common/storage";import{NDArrayCore as h}from"../common/ndarray-core";import{Complex as E,isComplexLike as U}from"../common/complex";import{getTypedArrayConstructor as w,isBigIntDType as C,isComplexDType as S,DEFAULT_DTYPE as m}from"../common/dtype";function B(r){return new h(r)}function z(r){return!!(r instanceof h||r&&typeof r=="object"&&"storage"in r&&r.storage instanceof D)}function N(r,e=m){const t=D.zeros(r,e);return new h(t)}function v(r,e=m){const t=D.ones(r,e);return new h(t)}function k(r,e=m){return N(r,e)}function F(r,e,t){let n=t;n||(typeof e=="bigint"?n="int64":typeof e=="boolean"?n="bool":Number.isInteger(e)?n="int32":n=m);const o=w(n);if(!o)throw new Error(`Cannot create full array with dtype ${n}`);const u=r.reduce((s,c)=>s*c,1),a=new o(u);if(C(n)){const s=typeof e=="bigint"?e:BigInt(Math.round(Number(e)));a.fill(s)}else n==="bool"?a.fill(e?1:0):a.fill(Number(e));const f=D.fromData(a,r,n);return new h(f)}function $(r){const e=[];let t=r;for(;Array.isArray(t);)e.push(t.length),t=t[0];return e}function M(r){return typeof r=="bigint"?!0:Array.isArray(r)?r.some(e=>M(e)):!1}function I(r){return U(r)?!0:Array.isArray(r)?r.some(e=>I(e)):!1}function j(r){const e=[];function t(n){Array.isArray(n)?n.forEach(o=>t(o)):e.push(n)}return t(r),e}function b(r,e){if(r instanceof h)return!e||r.dtype===e?r.copy():r.astype(e);const t=M(r),n=I(r),o=$(r),u=o.reduce((A,p)=>A*p,1);let a=e;a||(n?a="complex128":t?a="int64":a=m);const f=S(a),s=w(a);if(!s)throw new Error(`Cannot create array with dtype ${a}`);const c=f?u*2:u,i=new s(c),l=j(r);if(C(a)){const A=i;for(let p=0;p<u;p++){const g=l[p];A[p]=typeof g=="bigint"?g:BigInt(Math.round(Number(g)))}}else if(a==="bool"){const A=i;for(let p=0;p<u;p++)A[p]=l[p]?1:0}else if(f){const A=i;for(let p=0;p<u;p++){const g=l[p];let d,x;g instanceof E?(d=g.re,x=g.im):typeof g=="object"&&g!==null&&"re"in g?(d=g.re,x=g.im??0):(d=Number(g),x=0),A[p*2]=d,A[p*2+1]=x}}else{const A=i;for(let p=0;p<u;p++){const g=l[p];A[p]=Number(g)}}const y=D.fromData(i,o,a);return new h(y)}function Y(r,e,t=1,n=m){let o=r,u=e;if(e===void 0&&(o=0,u=r),u===void 0)throw new Error("stop is required");const a=Math.max(0,Math.ceil((u-o)/t)),f=w(n);if(!f)throw new Error(`Cannot create arange array with dtype ${n}`);const s=new f(a);if(C(n))for(let i=0;i<a;i++)s[i]=BigInt(Math.round(o+i*t));else if(n==="bool")for(let i=0;i<a;i++)s[i]=o+i*t!==0?1:0;else for(let i=0;i<a;i++)s[i]=o+i*t;const c=D.fromData(s,[a],n);return new h(c)}function K(r,e,t=50,n=m){if(t<0)throw new Error("num must be non-negative");if(t===0)return b([],n);if(t===1)return b([r],n);const o=w(n);if(!o)throw new Error(`Cannot create linspace array with dtype ${n}`);const u=new o(t),a=(e-r)/(t-1);if(C(n))for(let s=0;s<t;s++)u[s]=BigInt(Math.round(r+s*a));else if(n==="bool")for(let s=0;s<t;s++)u[s]=r+s*a!==0?1:0;else for(let s=0;s<t;s++)u[s]=r+s*a;const f=D.fromData(u,[t],n);return new h(f)}function R(r,e,t=50,n=10,o=m){if(t<0)throw new Error("num must be non-negative");if(t===0)return b([],o);if(t===1)return b([Math.pow(n,r)],o);const u=w(o);if(!u)throw new Error(`Cannot create logspace array with dtype ${o}`);const a=new u(t),f=(e-r)/(t-1);if(C(o))for(let c=0;c<t;c++){const i=r+c*f;a[c]=BigInt(Math.round(Math.pow(n,i)))}else if(o==="bool")for(let c=0;c<t;c++){const i=r+c*f;a[c]=Math.pow(n,i)!==0?1:0}else for(let c=0;c<t;c++){const i=r+c*f;a[c]=Math.pow(n,i)}const s=D.fromData(a,[t],o);return new h(s)}function H(r,e,t=50,n=m){if(t<0)throw new Error("num must be non-negative");if(r===0||e===0)throw new Error("Geometric sequence cannot include zero");if(t===0)return b([],n);if(t===1)return b([r],n);const o=Math.sign(r),u=Math.sign(e);if(o!==u)throw new Error("Geometric sequence cannot contain both positive and negative values");const a=w(n);if(!a)throw new Error(`Cannot create geomspace array with dtype ${n}`);const f=new a(t),s=Math.log(Math.abs(r)),i=(Math.log(Math.abs(e))-s)/(t-1);if(C(n))for(let y=0;y<t;y++){const A=o*Math.exp(s+y*i);f[y]=BigInt(Math.round(A))}else if(n==="bool")for(let y=0;y<t;y++){const A=o*Math.exp(s+y*i);f[y]=A!==0?1:0}else for(let y=0;y<t;y++){const A=o*Math.exp(s+y*i);f[y]=A}const l=D.fromData(f,[t],n);return new h(l)}function L(r,e,t=0,n=m){const o=e??r,u=N([r,o],n),a=u.data;if(C(n)){const f=a;for(let s=0;s<r;s++){const c=s+t;c>=0&&c<o&&(f[s*o+c]=BigInt(1))}}else{const f=a;for(let s=0;s<r;s++){const c=s+t;c>=0&&c<o&&(f[s*o+c]=1)}}return u}function J(r,e=m){return L(r,r,0,e)}function T(r,e){return z(r)?!e||r.dtype===e?r:r.astype(e):b(r,e)}function Q(r,e){return N(Array.from(r.shape),e??r.dtype)}function W(r,e){return v(Array.from(r.shape),e??r.dtype)}function X(r,e){return k(Array.from(r.shape),e??r.dtype)}function Z(r,e,t){return F(Array.from(r.shape),e,t??r.dtype)}function rr(r){return r.copy()}function er(r,e){return T(r,e)}function tr(r,e){return T(r,e).copy()}function nr(r,e){return T(r,e).copy()}function or(r,e){const t=T(r,e),n=t.data;if(n instanceof Float64Array){const o=new Uint32Array(n.buffer,n.byteOffset,n.length*2),u=2146435072;for(let a=1;a<o.length;a+=2)if((o[a]&u)===u)throw new Error("array must not contain infs or NaNs")}else if(n instanceof Float32Array){const o=new Uint32Array(n.buffer,n.byteOffset,n.length),u=2139095040;for(let a=0;a<o.length;a++)if((o[a]&u)===u)throw new Error("array must not contain infs or NaNs")}return t}function ar(r,e,t){let n=r;return e&&e!==r.dtype&&(n=n.astype(e)),n}function q(r){const e=r.data,t=D.fromData(e.slice(),[e.length],r.dtype);return B(t)}function _(r,e=0){const t=r.shape,n=r.data,o=r.dtype;if(t.length===1){const u=t[0],a=u+Math.abs(e),f=N([a,a],o),s=f.data;for(let c=0;c<u;c++){const i=e>=0?c:c-e,l=e>=0?c+e:c;i>=0&&i<a&&l>=0&&l<a&&(s[i*a+l]=n[c])}return f}else if(t.length===2){const u=t[0],a=t[1],f=Math.min(e>=0?Math.min(u,a-e):Math.min(u+e,a),Math.max(0,e>=0?a-e:u+e));if(f<=0)return b([],o);const s=[];for(let c=0;c<f;c++){const i=e>=0?c:c-e,l=e>=0?c+e:c;i>=0&&i<u&&l>=0&&l<a&&s.push(n[i*a+l])}return b(s,o)}throw new Error("Input must be 1-D or 2-D")}function sr(r,e=0){const t=q(r);return _(t,e)}function ir(r,e,t=0,n=m){const o=e??r,u=N([r,o],n),a=u.data;for(let f=0;f<r;f++)for(let s=0;s<=Math.min(f+t,o-1);s++)s>=0&&(a[f*o+s]=1);return u}function cr(r,e=0){const t=r.shape;if(t.length<2)throw new Error("Input must be at least 2-D");const n=r.copy(),o=n.data,u=t[t.length-2],a=t[t.length-1],f=t.slice(0,-2).reduce((i,l)=>i*l,1),s=u*a,c=o instanceof BigInt64Array||o instanceof BigUint64Array;for(let i=0;i<f;i++){const l=i*s;for(let y=0;y<u;y++){const A=l+y*a+Math.max(0,Math.min(y+e+1,a)),p=l+y*a+a;A<p&&o.fill(c?0n:0,A,p)}}return n}function ur(r,e=0){const t=r.shape;if(t.length<2)throw new Error("Input must be at least 2-D");const n=r.copy(),o=n.data,u=t[t.length-2],a=t[t.length-1],f=t.slice(0,-2).reduce((i,l)=>i*l,1),s=u*a,c=o instanceof BigInt64Array||o instanceof BigUint64Array;for(let i=0;i<f;i++){const l=i*s;for(let y=0;y<u;y++){const A=l+y*a+Math.max(0,Math.min(y+e,a)),p=l+y*a;p<A&&o.fill(c?0n:0,p,A)}}return n}function fr(r,e,t=!1){const n=r.size,o=e??n,u=r.data,a=N([n,o],r.dtype),f=a.data;for(let s=0;s<n;s++){const c=u[s];for(let i=0;i<o;i++){const l=t?i:o-1-i;f[s*o+i]=Math.pow(c,l)}}return a}function yr(r,e=m,t=-1,n=0){const o=w(e);if(!o)throw new Error(`Unsupported dtype: ${e}`);let u;if(r instanceof ArrayBuffer){const f=n,s=t<0?(r.byteLength-f)/o.BYTES_PER_ELEMENT:t;u=new o(r,f,s)}else{const f=n,s=t<0?r.length:n+t,c=Array.from(r.slice(f,s));u=new o(c.length);for(let i=0;i<c.length;i++)u[i]=c[i]}const a=D.fromData(u,[u.length],e);return B(a)}function lr(r,e,t=m){const n=e.reduce((c,i)=>c*i,1),o=w(t);if(!o)throw new Error(`Unsupported dtype: ${t}`);const u=new o(n),a=[];let f=1;for(let c=e.length-1;c>=0;c--)a.unshift(f),f*=e[c];for(let c=0;c<n;c++){const i=[];let l=c;for(let y=0;y<e.length;y++)i.push(Math.floor(l/a[y])),l=l%a[y];u[c]=r(...i)}const s=D.fromData(u,e,t);return B(s)}function pr(r,e=m,t=-1){const n=[];let o=0;for(const u of r){if(t>=0&&o>=t)break;n.push(u),o++}return b(n,e)}function Ar(r,e=m,t=-1,n){const o=n??/\s+/,a=r.split(o).filter(s=>s.trim()!=="").map(s=>parseFloat(s.trim())),f=t>=0?a.slice(0,t):a;return b(f,e)}function gr(r,e=m,t=-1,n=""){throw new Error("fromfile requires Node.js file system access")}function mr(...r){if(r.length===0)return[];if(r.length===1)return[r[0].copy()];const e=r.map(o=>o.size),t=[...e],n=[];for(let o=0;o<r.length;o++){const u=r[o],a=u.data,f=N(t,u.dtype),s=f.data,c=[];let i=1;for(let y=t.length-1;y>=0;y--)c.unshift(i),i*=t[y];const l=t.reduce((y,A)=>y*A,1);for(let y=0;y<l;y++){const A=Math.floor(y/c[o])%e[o];s[y]=a[A]}n.push(f)}return n}export{Y as arange,b as array,er as asanyarray,T as asarray,or as asarray_chkfinite,tr as ascontiguousarray,nr as asfortranarray,rr as copy,_ as diag,sr as diagflat,k as empty,X as empty_like,L as eye,yr as frombuffer,gr as fromfile,lr as fromfunction,pr as fromiter,Ar as fromstring,F as full,Z as full_like,H as geomspace,J as identity,K as linspace,R as logspace,mr as meshgrid,v as ones,W as ones_like,ar as require,ir as tri,cr as tril,ur as triu,fr as vander,N as zeros,Q as zeros_like};

package/dist/esm/core/linalg.js

@@ -1 +1 @@
-import*as u from"../common/ops/linalg";import{toStorage as o,fromStorage as a,fromStorageView as D,Complex as C}from"./types";function N(r,e){const n=u.dot(o(r),o(e));return typeof n=="number"||typeof n=="bigint"||n instanceof C?n:a(n)}function i(r){return u.trace(o(r))}function s(r,e=0,n=0,t=1){return a(u.diagonal(o(r),e,n,t))}function w(r,e){return a(u.kron(o(r),o(e)))}function b(r,e){return D(u.transpose(o(r),e),r)}function l(r,e){const n=u.inner(o(r),o(e));return typeof n=="number"||typeof n=="bigint"||n instanceof C?n:a(n)}function p(r,e){return a(u.outer(o(r),o(e)))}function c(r,e,n=2){const t=u.tensordot(o(r),o(e),n);return typeof t=="number"||typeof t=="bigint"||t instanceof C?t:a(t)}function k(r,...e){const n=e.map(y=>o(y)),t=u.einsum(r,...n);return typeof t=="number"||typeof t=="bigint"||t instanceof C?t:a(t)}function q(r,...e){const n=e.map(t=>o(t));return u.einsum_path(r,...n)}function S(r,e){return u.vdot(o(r),o(e))}function f(r,e,n=-1){const t=u.vecdot(o(r),o(e),n);return typeof t=="number"||typeof t=="bigint"||t instanceof C?t:a(t)}function x(r){return D(u.matrix_transpose(o(r)),r)}function g(r,e){return D(u.permute_dims(o(r),e),r)}function O(r,e){return a(u.matvec(o(r),o(e)))}function L(r,e){return a(u.vecmat(o(r),o(e)))}function v(r,e,n=-1,t=-1,y=-1,A){const m=u.cross(o(r),o(e),n,t,y,A);return typeof m=="number"||m instanceof C?m:a(m)}function U(r,e){return a(u.matmul(o(r),o(e)))}const V={matmul:(r,e)=>a(u.matmul(o(r),o(e))),dot:(r,e)=>N(r,e),det:r=>u.det(o(r)),inv:r=>a(u.inv(o(r))),solve:(r,e)=>a(u.solve(o(r),o(e))),lstsq:(r,e,n)=>{const t=u.lstsq(o(r),o(e),n);return{x:a(t.x),residuals:a(t.residuals),rank:t.rank,s:a(t.s)}},norm:(r,e,n,t)=>{const y=u.norm(o(r),e,n,t);return typeof y=="number"?y:a(y)},cond:(r,e)=>u.cond(o(r),e),matrix_rank:(r,e)=>u.matrix_rank(o(r),e),matrix_power:(r,e)=>a(u.matrix_power(o(r),e)),pinv:(r,e)=>a(u.pinv(o(r),e)),qr:(r,e)=>{const n=u.qr(o(r),e);return"h"in n&&"tau"in n?{h:a(n.h),tau:a(n.tau)}:"q"in n&&"r"in n?{q:a(n.q),r:a(n.r)}:a(n)},cholesky:(r,e)=>a(u.cholesky(o(r),e)),svd:(r,e,n)=>{const t=u.svd(o(r),e,n);return"u"in t&&"s"in t&&"vt"in t?{u:a(t.u),s:a(t.s),vt:a(t.vt)}:a(t)},eig:r=>{const e=u.eig(o(r));return{w:a(e.w),v:a(e.v)}},eigh:(r,e)=>{const n=u.eigh(o(r),e);return{w:a(n.w),v:a(n.v)}},eigvals:r=>a(u.eigvals(o(r))),eigvalsh:(r,e)=>a(u.eigvalsh(o(r),e)),slogdet:r=>u.slogdet(o(r)),svdvals:r=>a(u.svdvals(o(r))),multi_dot:r=>a(u.multi_dot(r.map(e=>o(e)))),tensorinv:(r,e)=>a(u.tensorinv(o(r),e)),tensorsolve:(r,e,n)=>a(u.tensorsolve(o(r),o(e),n)),vector_norm:(r,e,n,t)=>{const y=u.vector_norm(o(r),e,n,t);return typeof y=="number"?y:a(y)},matrix_norm:(r,e,n)=>{const t=u.matrix_norm(o(r),e,n);return typeof t=="number"?t:a(t)},cross:(r,e,n,t,y,A)=>v(r,e,n,t,y,A),matrix_transpose:r=>x(r),permute_dims:(r,e)=>g(r,e),trace:r=>i(r),diagonal:(r,e,n,t)=>s(r,e,n,t),outer:(r,e)=>p(r,e),inner:(r,e)=>l(r,e),tensordot:(r,e,n)=>c(r,e,n),vecdot:(r,e,n)=>f(r,e,n),transpose:(r,e)=>b(r,e)};export{v as cross,s as diagonal,N as dot,k as einsum,q as einsum_path,l as inner,w as kron,V as linalg,U as matmul,x as matrix_transpose,O as matvec,p as outer,g as permute_dims,c as tensordot,i as trace,b as transpose,S as vdot,f as vecdot,L as vecmat};
+import*as u from"../common/ops/linalg";import{toStorage as n,fromStorage as a,fromStorageView as A,Complex as C,ArrayStorage as s}from"./types";function D(r,e){const o=u.dot(n(r),n(e));return typeof o=="number"||typeof o=="bigint"||o instanceof C?o:a(o)}function N(r,e=0,o=0,t=1){const y=u.trace(n(r),e,o,t);return y instanceof s?a(y):y}function b(r,e=0,o=0,t=1){return a(u.diagonal(n(r),e,o,t))}function w(r,e){return a(u.kron(n(r),n(e)))}function l(r,e){return A(u.transpose(n(r),e),r)}function c(r,e){const o=u.inner(n(r),n(e));return typeof o=="number"||typeof o=="bigint"||o instanceof C?o:a(o)}function p(r,e){return a(u.outer(n(r),n(e)))}function f(r,e,o=2){const t=u.tensordot(n(r),n(e),o);return typeof t=="number"||typeof t=="bigint"||t instanceof C?t:a(t)}function k(r,...e){const o=e.map(y=>n(y)),t=u.einsum(r,...o);return typeof t=="number"||typeof t=="bigint"||t instanceof C?t:a(t)}function q(r,...e){const o=e.map(t=>n(t));return u.einsum_path(r,...o)}function S(r,e){return u.vdot(n(r),n(e))}function g(r,e,o=-1){const t=u.vecdot(n(r),n(e),o);return typeof t=="number"||typeof t=="bigint"||t instanceof C?t:a(t)}function x(r){return A(u.matrix_transpose(n(r)),r)}function d(r,e){return A(u.permute_dims(n(r),e),r)}function O(r,e){return a(u.matvec(n(r),n(e)))}function L(r,e){return a(u.vecmat(n(r),n(e)))}function v(r,e,o=-1,t=-1,y=-1,i){const m=u.cross(n(r),n(e),o,t,y,i);return typeof m=="number"||typeof m=="bigint"||m instanceof C?m:a(m)}function U(r,e){return a(u.matmul(n(r),n(e)))}const V={matmul:(r,e)=>a(u.matmul(n(r),n(e))),dot:(r,e)=>D(r,e),det:r=>{const e=u.det(n(r));return typeof e=="number"?e:a(e)},inv:r=>a(u.inv(n(r))),solve:(r,e)=>a(u.solve(n(r),n(e))),lstsq:(r,e,o)=>{const t=u.lstsq(n(r),n(e),o);return{x:a(t.x),residuals:a(t.residuals),rank:t.rank,s:a(t.s)}},norm:(r,e,o,t)=>{const y=u.norm(n(r),e,o,t);return typeof y=="number"?y:a(y)},cond:(r,e)=>u.cond(n(r),e),matrix_rank:(r,e)=>u.matrix_rank(n(r),e),matrix_power:(r,e)=>a(u.matrix_power(n(r),e)),pinv:(r,e)=>a(u.pinv(n(r),e)),qr:(r,e)=>{const o=u.qr(n(r),e);return"h"in o&&"tau"in o?{h:a(o.h),tau:a(o.tau)}:"q"in o&&"r"in o?{q:a(o.q),r:a(o.r)}:a(o)},cholesky:(r,e)=>a(u.cholesky(n(r),e)),svd:(r,e,o)=>{const t=u.svd(n(r),e,o);return"u"in t&&"s"in t&&"vt"in t?{u:a(t.u),s:a(t.s),vt:a(t.vt)}:a(t)},eig:r=>{const e=u.eig(n(r));return{w:a(e.w),v:a(e.v)}},eigh:(r,e)=>{const o=u.eigh(n(r),e);return{w:a(o.w),v:a(o.v)}},eigvals:r=>a(u.eigvals(n(r))),eigvalsh:(r,e)=>a(u.eigvalsh(n(r),e)),slogdet:r=>{const e=u.slogdet(n(r)),o=e.sign instanceof s?a(e.sign):e.sign,t=e.logabsdet instanceof s?a(e.logabsdet):e.logabsdet;return{sign:o,logabsdet:t}},svdvals:r=>a(u.svdvals(n(r))),multi_dot:r=>a(u.multi_dot(r.map(e=>n(e)))),tensorinv:(r,e)=>a(u.tensorinv(n(r),e)),tensorsolve:(r,e,o)=>a(u.tensorsolve(n(r),n(e),o)),vector_norm:(r,e,o,t)=>{const y=u.vector_norm(n(r),e,o,t);return typeof y=="number"?y:a(y)},matrix_norm:(r,e,o)=>{const t=u.matrix_norm(n(r),e,o);return typeof t=="number"?t:a(t)},cross:(r,e,o,t,y,i)=>v(r,e,o,t,y,i),matrix_transpose:r=>x(r),permute_dims:(r,e)=>d(r,e),trace:(r,e,o,t)=>N(r,e,o,t),diagonal:(r,e,o,t)=>b(r,e,o,t),outer:(r,e)=>p(r,e),inner:(r,e)=>c(r,e),tensordot:(r,e,o)=>f(r,e,o),vecdot:(r,e,o)=>g(r,e,o),transpose:(r,e)=>l(r,e)};export{v as cross,b as diagonal,D as dot,k as einsum,q as einsum_path,c as inner,w as kron,V as linalg,U as matmul,x as matrix_transpose,O as matvec,p as outer,d as permute_dims,f as tensordot,N as trace,l as transpose,S as vdot,g as vecdot,L as vecmat};

package/dist/esm/core/polynomial.js

@@ -1 +1 @@
-import{NDArrayCore as z}from"../common/ndarray-core";import{ArrayStorage as j}from"../common/storage";import{array as A}from"./creation";function D(n){return n instanceof z?n:A(n)}function L(n){const m=D(n),e=m.data,a=m.size;if(a===0)return A([1]);let r=[1];for(let i=0;i<a;i++){const l=e[i],o=new Array(r.length+1).fill(0);for(let c=0;c<r.length;c++)o[c]+=r[c],o[c+1]-=r[c]*l;r=o}return A(r)}function S(n,m){const e=D(n),a=D(m),r=e.data,i=a.data,l=Math.max(e.size,a.size),o=new Array(l).fill(0);for(let s=0;s<e.size;s++)o[l-e.size+s]+=r[s];for(let s=0;s<a.size;s++)o[l-a.size+s]+=i[s];let c=0;for(;c<o.length-1&&o[c]===0;)c++;return A(o.slice(c))}function T(n,m=1){let e=D(n);for(let a=0;a<m;a++){const r=e.data,i=e.size;if(i<=1)return A([0]);const l=[];for(let o=0;o<i-1;o++){const c=i-1-o;l.push(r[o]*c)}e=A(l)}return e}function B(n,m){const e=[...D(n).data],a=[...D(m).data];if(a.length===0||a.length===1&&a[0]===0)throw new Error("Division by zero polynomial");for(;e.length>1&&e[0]===0;)e.shift();for(;a.length>1&&a[0]===0;)a.shift();if(e.length<a.length)return[A([0]),A(e)];const r=[],i=[...e];for(;i.length>=a.length;){const l=i[0]/a[0];r.push(l);for(let o=0;o<a.length;o++)i[o]-=l*a[o];i.shift()}for(;i.length>1&&Math.abs(i[0])<1e-15;)i.shift();return[A(r.length>0?r:[0]),A(i.length>0?i:[0])]}function F(n,m,e){const a=n.data,r=m.data,i=n.size;if(e>=i)throw new Error("polyfit: degree must be less than number of points");const l=[];for(let u=0;u<i;u++){const p=[];for(let t=e;t>=0;t--)p.push(Math.pow(a[u],t));l.push(p)}const o=[],c=[];for(let u=0;u<=e;u++){o.push([]);for(let t=0;t<=e;t++){let h=0;for(let y=0;y<i;y++)h+=l[y][u]*l[y][t];o[u].push(h)}let p=0;for(let t=0;t<i;t++)p+=l[t][u]*r[t];c.push(p)}const s=e+1,b=o.map((u,p)=>[...u,c[p]]);for(let u=0;u<s;u++){let p=u;for(let t=u+1;t<s;t++)Math.abs(b[t][u])>Math.abs(b[p][u])&&(p=t);[b[u],b[p]]=[b[p],b[u]];for(let t=u+1;t<s;t++){const h=b[t][u]/b[u][u];for(let y=u;y<=s;y++)b[t][y]-=h*b[u][y]}}const f=new Array(s).fill(0);for(let u=s-1;u>=0;u--){let p=b[u][s];for(let t=u+1;t<s;t++)p-=b[u][t]*f[t];f[u]=p/b[u][u]}return A(f)}function P(n,m=1,e=0){let a=D(n);const r=Array.isArray(e)?e:[e];for(let i=0;i<m;i++){const l=a.data,o=a.size,c=[];for(let b=0;b<o;b++){const f=o-b;c.push(l[b]/f)}const s=i<r.length?r[i]:0;c.push(s),a=A(c)}return a}function Q(n,m){const e=D(n),a=D(m),r=e.data,i=a.data,l=e.size+a.size-1,o=new Array(l).fill(0);for(let c=0;c<e.size;c++)for(let s=0;s<a.size;s++)o[c+s]+=r[c]*i[s];return A(o)}function V(n,m){const e=D(n),a=D(m),r=e.data,i=a.data,l=Math.max(e.size,a.size),o=new Array(l).fill(0);for(let s=0;s<e.size;s++)o[l-e.size+s]+=r[s];for(let s=0;s<a.size;s++)o[l-a.size+s]-=i[s];let c=0;for(;c<o.length-1&&o[c]===0;)c++;return A(o.slice(c))}function Z(n,m){const e=D(n),a=e.data;if(typeof m=="number"){let o=a[0];for(let c=1;c<e.size;c++)o=o*m+a[c];return o}const r=m instanceof z?m:A(m),i=r.data,l=[];for(let o=0;o<r.size;o++){const c=i[o];let s=a[0];for(let b=1;b<e.size;b++)s=s*c+a[b];l.push(s)}return A(l)}function G(n){const e=[...D(n).data];for(;e.length>1&&e[0]===0;)e.shift();let a=0;for(;e.length>1&&e[e.length-1]===0;)e.pop(),a++;const r=e.length-1;if(r+a===0)return q([],[]);let l=[],o=[];if(r===1)l.push(-e[1]/e[0]),o.push(0);else if(r===2){const f=e[0],u=e[1],p=e[2],t=u*u-4*f*p;if(t>=0){const h=Math.sqrt(t);l.push((-u+h)/(2*f),(-u-h)/(2*f)),o.push(0,0)}else{const h=Math.sqrt(-t);l.push(-u/(2*f),-u/(2*f)),o.push(h/(2*f),-h/(2*f))}}else if(r>=3){const f=k(e,r);for(const u of f)l.push(u.re),o.push(u.im)}for(let f=0;f<a;f++)l.push(0),o.push(0);const c=l.map((f,u)=>u);c.sort((f,u)=>{const p=Math.sqrt(l[f]**2+o[f]**2),t=Math.sqrt(l[u]**2+o[u]**2);return Math.abs(p-t)>1e-10?t-p:Math.abs(l[f]-l[u])>1e-10?l[u]-l[f]:o[u]-o[f]});const s=c.map(f=>l[f]),b=c.map(f=>o[f]);return q(s,b)}function q(n,m){const e=n.length,a=new Float64Array(2*e);for(let i=0;i<e;i++)a[2*i]=n[i],a[2*i+1]=m[i];const r=j.fromData(a,[e],"complex128");return new z(r)}function k(n,m){const e=Array.from({length:m},()=>new Array(m).fill(0)),a=n[0];for(let r=0;r<m;r++)e[0][r]=-n[r+1]/a;for(let r=1;r<m;r++)e[r][r-1]=1;return R(e,m)}function R(n,m){const e=[];let r=m,i=0;const l=100*m;let o=m,c=0;for(;r>0&&i<l;){if(i++,r===o?c++:(c=0,o=r),r===1){e.push({re:n[0][0],im:0}),r=0;break}if(r===2){e.push(...v(n[0][0],n[0][1],n[1][0],n[1][1])),r=0;break}let s=r-1;for(;s>0;){const t=Math.abs(n[s-1][s-1])+Math.abs(n[s][s]),h=222e-18*(t===0?1:t);if(Math.abs(n[s][s-1])<=h){n[s][s-1]=0;break}s--}const b=r-s;if(b===1){e.push({re:n[r-1][r-1],im:0}),r--;continue}if(b===2){e.push(...v(n[r-2][r-2],n[r-2][r-1],n[r-1][r-2],n[r-1][r-1])),r-=2;continue}let f;if(c>0&&c%10===0)f=Math.abs(n[r-1][r-2])+Math.abs(n[r-2][r-3]);else{const t=n[r-2][r-2],h=n[r-2][r-1],y=n[r-1][r-2],C=n[r-1][r-1],d=t+C,g=t*C-h*y,w=d*d-4*g;if(w>=0){const N=Math.sqrt(w),M=(d+N)/2,x=(d-N)/2;f=Math.abs(M-C)<Math.abs(x-C)?M:x}else f=C}for(let t=s;t<r;t++)n[t][t]=n[t][t]-f;const u=[],p=[];for(let t=s;t<r-1;t++){const h=Math.hypot(n[t][t],n[t+1][t]),y=h===0?1:n[t][t]/h,C=h===0?0:n[t+1][t]/h;u.push(y),p.push(C);for(let d=t;d<r;d++){const g=n[t][d],w=n[t+1][d];n[t][d]=y*g+C*w,n[t+1][d]=-C*g+y*w}}for(let t=0;t<u.length;t++){const h=s+t,y=u[t],C=p[t],d=Math.min(h+2,r-1);for(let g=s;g<=d;g++){const w=n[g][h],N=n[g][h+1];n[g][h]=y*w+C*N,n[g][h+1]=-C*w+y*N}}for(let t=s;t<r;t++)n[t][t]=n[t][t]+f}if(r>0)for(let s=0;s<r;s++)e.push({re:n[s][s],im:0});return e}function v(n,m,e,a){const r=n+a,i=n*a-m*e,l=r*r-4*i;if(l>=0){const o=Math.sqrt(l);return[{re:(r+o)/2,im:0},{re:(r-o)/2,im:0}]}else{const o=Math.sqrt(-l);return[{re:r/2,im:o/2},{re:r/2,im:-o/2}]}}export{L as poly,S as polyadd,T as polyder,B as polydiv,F as polyfit,P as polyint,Q as polymul,V as polysub,Z as polyval,G as roots};
+import{NDArrayCore as z}from"../common/ndarray-core";import{ArrayStorage as q}from"../common/storage";import{array as A}from"./creation";function D(n){return n instanceof z?n:A(n)}function L(n){const f=D(n),e=f.data,a=f.size;if(a===0)return A([1]);let r=[1];for(let i=0;i<a;i++){const l=e[i],u=new Array(r.length+1).fill(0);for(let c=0;c<r.length;c++)u[c]+=r[c],u[c+1]-=r[c]*l;r=u}return A(r)}function S(n,f){const e=D(n),a=D(f),r=e.data,i=a.data,l=Math.max(e.size,a.size),u=new Array(l).fill(0);for(let o=0;o<e.size;o++)u[l-e.size+o]+=r[o];for(let o=0;o<a.size;o++)u[l-a.size+o]+=i[o];let c=0;for(;c<u.length-1&&u[c]===0;)c++;return A(u.slice(c))}function F(n,f=1){let e=D(n);for(let a=0;a<f;a++){const r=e.data,i=e.size;if(i<=1)return A([0]);const l=[];for(let u=0;u<i-1;u++){const c=i-1-u;l.push(r[u]*c)}e=A(l)}return e}function T(n,f){const e=[...D(n).data],a=[...D(f).data];if(a.length===0||a.length===1&&a[0]===0)throw new Error("Division by zero polynomial");for(;e.length>1&&e[0]===0;)e.shift();for(;a.length>1&&a[0]===0;)a.shift();if(e.length<a.length)return[A([0]),A(e)];const r=[],i=[...e];for(;i.length>=a.length;){const l=i[0]/a[0];r.push(l);for(let u=0;u<a.length;u++)i[u]-=l*a[u];i.shift()}for(;i.length>1&&Math.abs(i[0])<1e-15;)i.shift();return[A(r.length>0?r:[0]),A(i.length>0?i:[0])]}function B(n,f,e){const a=n.data,r=f.data,i=n.size;if(e>=i)throw new Error("polyfit: degree must be less than number of points");const l=[];for(let s=0;s<i;s++){const p=[];for(let t=e;t>=0;t--)p.push(Math.pow(a[s],t));l.push(p)}const u=[],c=[];for(let s=0;s<=e;s++){u.push([]);for(let t=0;t<=e;t++){let b=0;for(let y=0;y<i;y++)b+=l[y][s]*l[y][t];u[s].push(b)}let p=0;for(let t=0;t<i;t++)p+=l[t][s]*r[t];c.push(p)}const o=e+1,h=u.map((s,p)=>[...s,c[p]]);for(let s=0;s<o;s++){let p=s;for(let t=s+1;t<o;t++)Math.abs(h[t][s])>Math.abs(h[p][s])&&(p=t);[h[s],h[p]]=[h[p],h[s]];for(let t=s+1;t<o;t++){const b=h[t][s]/h[s][s];for(let y=s;y<=o;y++)h[t][y]-=b*h[s][y]}}const m=new Array(o).fill(0);for(let s=o-1;s>=0;s--){let p=h[s][o];for(let t=s+1;t<o;t++)p-=h[s][t]*m[t];m[s]=p/h[s][s]}return A(m)}function P(n,f=1,e=0){let a=D(n);const r=Array.isArray(e)?e:[e];for(let i=0;i<f;i++){const l=a.data,u=a.size,c=[];for(let h=0;h<u;h++){const m=u-h;c.push(l[h]/m)}const o=i<r.length?r[i]:0;c.push(o),a=A(c)}return a}function Q(n,f){const e=D(n),a=D(f),r=e.data,i=a.data,l=e.size+a.size-1,u=new Array(l).fill(0);for(let c=0;c<e.size;c++)for(let o=0;o<a.size;o++)u[c+o]+=r[c]*i[o];return A(u)}function V(n,f){const e=D(n),a=D(f),r=e.data,i=a.data,l=Math.max(e.size,a.size),u=new Array(l).fill(0);for(let o=0;o<e.size;o++)u[l-e.size+o]+=r[o];for(let o=0;o<a.size;o++)u[l-a.size+o]-=i[o];let c=0;for(;c<u.length-1&&u[c]===0;)c++;return A(u.slice(c))}function Z(n,f){const e=D(n),a=e.data;if(typeof f=="number"){let o=a[0];for(let h=1;h<e.size;h++)o=o*f+a[h];return o}const r=f instanceof z?f:A(f),i=r.data,l=r.size,u=e.size,c=new Float64Array(l);for(let o=0;o<l;o++){const h=Number(i[o]);let m=Number(a[0]);for(let s=1;s<u;s++)m=m*h+Number(a[s]);c[o]=m}return new z(q.fromData(c,Array.from(r.shape),"float64"))}function G(n){const e=[...D(n).data];for(;e.length>1&&e[0]===0;)e.shift();let a=0;for(;e.length>1&&e[e.length-1]===0;)e.pop(),a++;const r=e.length-1;if(r+a===0)return v([],[]);let l=[],u=[];if(r===1)l.push(-e[1]/e[0]),u.push(0);else if(r===2){const m=e[0],s=e[1],p=e[2],t=s*s-4*m*p;if(t>=0){const b=Math.sqrt(t);l.push((-s+b)/(2*m),(-s-b)/(2*m)),u.push(0,0)}else{const b=Math.sqrt(-t);l.push(-s/(2*m),-s/(2*m)),u.push(b/(2*m),-b/(2*m))}}else if(r>=3){const m=k(e,r);for(const s of m)l.push(s.re),u.push(s.im)}for(let m=0;m<a;m++)l.push(0),u.push(0);const c=l.map((m,s)=>s);c.sort((m,s)=>{const p=Math.sqrt(l[m]**2+u[m]**2),t=Math.sqrt(l[s]**2+u[s]**2);return Math.abs(p-t)>1e-10?t-p:Math.abs(l[m]-l[s])>1e-10?l[s]-l[m]:u[s]-u[m]});const o=c.map(m=>l[m]),h=c.map(m=>u[m]);return v(o,h)}function v(n,f){const e=n.length,a=new Float64Array(2*e);for(let i=0;i<e;i++)a[2*i]=n[i],a[2*i+1]=f[i];const r=q.fromData(a,[e],"complex128");return new z(r)}function k(n,f){const e=Array.from({length:f},()=>new Array(f).fill(0)),a=n[0];for(let r=0;r<f;r++)e[0][r]=-n[r+1]/a;for(let r=1;r<f;r++)e[r][r-1]=1;return R(e,f)}function R(n,f){const e=[];let r=f,i=0;const l=100*f;let u=f,c=0;for(;r>0&&i<l;){if(i++,r===u?c++:(c=0,u=r),r===1){e.push({re:n[0][0],im:0}),r=0;break}if(r===2){e.push(...j(n[0][0],n[0][1],n[1][0],n[1][1])),r=0;break}let o=r-1;for(;o>0;){const t=Math.abs(n[o-1][o-1])+Math.abs(n[o][o]),b=222e-18*(t===0?1:t);if(Math.abs(n[o][o-1])<=b){n[o][o-1]=0;break}o--}const h=r-o;if(h===1){e.push({re:n[r-1][r-1],im:0}),r--;continue}if(h===2){e.push(...j(n[r-2][r-2],n[r-2][r-1],n[r-1][r-2],n[r-1][r-1])),r-=2;continue}let m;if(c>0&&c%10===0)m=Math.abs(n[r-1][r-2])+Math.abs(n[r-2][r-3]);else{const t=n[r-2][r-2],b=n[r-2][r-1],y=n[r-1][r-2],w=n[r-1][r-1],d=t+w,g=t*w-b*y,C=d*d-4*g;if(C>=0){const N=Math.sqrt(C),M=(d+N)/2,x=(d-N)/2;m=Math.abs(M-w)<Math.abs(x-w)?M:x}else m=w}for(let t=o;t<r;t++)n[t][t]=n[t][t]-m;const s=[],p=[];for(let t=o;t<r-1;t++){const b=Math.hypot(n[t][t],n[t+1][t]),y=b===0?1:n[t][t]/b,w=b===0?0:n[t+1][t]/b;s.push(y),p.push(w);for(let d=t;d<r;d++){const g=n[t][d],C=n[t+1][d];n[t][d]=y*g+w*C,n[t+1][d]=-w*g+y*C}}for(let t=0;t<s.length;t++){const b=o+t,y=s[t],w=p[t],d=Math.min(b+2,r-1);for(let g=o;g<=d;g++){const C=n[g][b],N=n[g][b+1];n[g][b]=y*C+w*N,n[g][b+1]=-w*C+y*N}}for(let t=o;t<r;t++)n[t][t]=n[t][t]+m}if(r>0)for(let o=0;o<r;o++)e.push({re:n[o][o],im:0});return e}function j(n,f,e,a){const r=n+a,i=n*a-f*e,l=r*r-4*i;if(l>=0){const u=Math.sqrt(l);return[{re:(r+u)/2,im:0},{re:(r-u)/2,im:0}]}else{const u=Math.sqrt(-l);return[{re:r/2,im:u/2},{re:r/2,im:-u/2}]}}export{L as poly,S as polyadd,F as polyder,T as polydiv,B as polyfit,P as polyint,Q as polymul,V as polysub,Z as polyval,G as roots};

package/dist/esm/core/reduction.js

@@ -1 +1 @@
-import*as a from"../common/ops/reduction";import{toStorage as u,fromStorage as m,Complex as b}from"./types";function A(n,o,e){const r=a.sum(u(n),o,e);return typeof r=="number"||typeof r=="bigint"||r instanceof b?r:m(r)}function D(n,o,e){const r=a.mean(u(n),o,e);return typeof r=="number"||r instanceof b?r:m(r)}function N(n,o,e){const r=a.prod(u(n),o,e);return typeof r=="number"||typeof r=="bigint"||r instanceof b?r:m(r)}function f(n,o,e){const r=a.max(u(n),o,e);return typeof r=="number"||r instanceof b?r:m(r)}const x=f;function p(n,o,e){const r=a.min(u(n),o,e);return typeof r=="number"||r instanceof b?r:m(r)}const d=p;function g(n,o,e){const r=a.ptp(u(n),o,e);return typeof r=="number"||r instanceof b?r:m(r)}function v(n,o){const e=a.argmin(u(n),o);return typeof e=="number"?e:m(e)}function q(n,o){const e=a.argmax(u(n),o);return typeof e=="number"?e:m(e)}function y(n,o,e,r){const t=a.variance(u(n),o,e,r);return typeof t=="number"?t:m(t)}const S=y;function _(n,o,e,r){const t=a.std(u(n),o,e,r);return typeof t=="number"?t:m(t)}function O(n,o,e){const r=a.median(u(n),o,e);return typeof r=="number"?r:m(r)}function h(n,o,e,r){const t=a.percentile(u(n),o,e,r);return typeof t=="number"?t:m(t)}function j(n,o,e,r){const t=a.quantile(u(n),o,e,r);return typeof t=="number"?t:m(t)}function k(n,o,e,r){const t=e?u(e):void 0,i=a.average(u(n),o,t,r);return typeof i=="number"||i instanceof b?i:m(i)}function w(n,o,e){const r=a.all(u(n),o,e);return typeof r=="boolean"?r:m(r)}function z(n,o,e){const r=a.any(u(n),o,e);return typeof r=="boolean"?r:m(r)}function c(n,o){return m(a.cumsum(u(n),o))}const B=c;function s(n,o){return m(a.cumprod(u(n),o))}const E=s;function F(n,o,e){const r=a.nansum(u(n),o,e);return typeof r=="number"||r instanceof b?r:m(r)}function G(n,o,e){const r=a.nanprod(u(n),o,e);return typeof r=="number"||r instanceof b?r:m(r)}function H(n,o,e){const r=a.nanmean(u(n),o,e);return typeof r=="number"||r instanceof b?r:m(r)}function I(n,o,e,r){const t=a.nanvar(u(n),o,e,r);return typeof t=="number"?t:m(t)}function J(n,o,e,r){const t=a.nanstd(u(n),o,e,r);return typeof t=="number"?t:m(t)}function K(n,o,e){const r=a.nanmin(u(n),o,e);return typeof r=="number"||r instanceof b?r:m(r)}function L(n,o,e){const r=a.nanmax(u(n),o,e);return typeof r=="number"||r instanceof b?r:m(r)}function M(n,o){const e=a.nanargmin(u(n),o);return typeof e=="number"?e:m(e)}function P(n,o){const e=a.nanargmax(u(n),o);return typeof e=="number"?e:m(e)}function Q(n,o){return m(a.nancumsum(u(n),o))}function R(n,o){return m(a.nancumprod(u(n),o))}function T(n,o,e){const r=a.nanmedian(u(n),o,e);return typeof r=="number"?r:m(r)}function U(n,o,e,r){const t=a.nanquantile(u(n),o,e,r);return typeof t=="number"?t:m(t)}function V(n,o,e,r){const t=a.nanpercentile(u(n),o,e,r);return typeof t=="number"?t:m(t)}export{w as all,x as amax,d as amin,z as any,q as argmax,v as argmin,k as average,s as cumprod,c as cumsum,E as cumulative_prod,B as cumulative_sum,f as max,D as mean,O as median,p as min,P as nanargmax,M as nanargmin,R as nancumprod,Q as nancumsum,L as nanmax,H as nanmean,T as nanmedian,K as nanmin,V as nanpercentile,G as nanprod,U as nanquantile,J as nanstd,F as nansum,I as nanvar,h as percentile,N as prod,g as ptp,j as quantile,_ as std,A as sum,S as var_,y as variance};
+import*as u from"../common/ops/reduction";import*as N from"../common/ops/shape";import{toStorage as b,fromStorage as y,Complex as A,ArrayStorage as p}from"./types";function f(o,r,n,e){const t=b(o).ndim,a=r.map(s=>s<0?t+s:s),m=[...a].sort((s,i)=>i-s),l=[...a].sort((s,i)=>s-i);let c=b(o);for(const s of m){const i=e(c,s,!1);if(!(i instanceof p)){if(n){const D=Array(t).fill(1),C=p.zeros(D,c.dtype);return C.iset(0,i),y(C)}return i}c=i}if(n)for(const s of l)c=N.expandDims(c,s);return y(c)}function _(o,r,n){if(Array.isArray(r))return f(o,r,n??!1,(t,a,m)=>u.sum(t,a,m));const e=u.sum(b(o),r,n);return typeof e=="number"||typeof e=="bigint"||e instanceof A?e:y(e)}function h(o,r,n){if(Array.isArray(r))return f(o,r,n??!1,(t,a,m)=>u.mean(t,a,m));const e=u.mean(b(o),r,n);return typeof e=="number"||e instanceof A?e:y(e)}function O(o,r,n){if(Array.isArray(r))return f(o,r,n??!1,(t,a,m)=>u.prod(t,a,m));const e=u.prod(b(o),r,n);return typeof e=="number"||typeof e=="bigint"||e instanceof A?e:y(e)}function x(o,r,n){if(Array.isArray(r))return f(o,r,n??!1,(t,a,m)=>u.max(t,a,m));const e=u.max(b(o),r,n);return typeof e=="number"||e instanceof A?e:y(e)}const z=x;function g(o,r,n){if(Array.isArray(r))return f(o,r,n??!1,(t,a,m)=>u.min(t,a,m));const e=u.min(b(o),r,n);return typeof e=="number"||e instanceof A?e:y(e)}const M=g;function j(o,r,n){const e=u.ptp(b(o),r,n);return typeof e=="number"||e instanceof A?e:y(e)}function w(o,r){const n=u.argmin(b(o),r);return typeof n=="number"?n:y(n)}function B(o,r){const n=u.argmax(b(o),r);return typeof n=="number"?n:y(n)}function v(o,r,n,e){if(Array.isArray(r))return f(o,r,e??!1,(a,m,l)=>u.variance(a,m,n,l));const t=u.variance(b(o),r,n,e);return typeof t=="number"?t:y(t)}const E=v;function F(o,r,n,e){if(Array.isArray(r))return f(o,r,e??!1,(a,m,l)=>u.std(a,m,n,l));const t=u.std(b(o),r,n,e);return typeof t=="number"?t:y(t)}function G(o,r,n){if(Array.isArray(r))return f(o,r,n??!1,(t,a,m)=>u.median(t,a,m));const e=u.median(b(o),r,n);return typeof e=="number"?e:y(e)}function H(o,r,n,e){const t=u.percentile(b(o),r,n,e);return typeof t=="number"?t:y(t)}function I(o,r,n,e){const t=u.quantile(b(o),r,n,e);return typeof t=="number"?t:y(t)}function J(o,r,n,e){const t=n?b(n):void 0,a=u.average(b(o),r,t,e);return typeof a=="number"||a instanceof A?a:y(a)}function K(o,r,n){if(Array.isArray(r))return f(o,r,n??!1,(t,a,m)=>u.all(t,a,m));const e=u.all(b(o),r,n);return typeof e=="boolean"?e:y(e)}function L(o,r,n){if(Array.isArray(r))return f(o,r,n??!1,(t,a,m)=>u.any(t,a,m));const e=u.any(b(o),r,n);return typeof e=="boolean"?e:y(e)}function S(o,r){return y(u.cumsum(b(o),r))}const P=S;function d(o,r){return y(u.cumprod(b(o),r))}const Q=d;function R(o,r,n){if(Array.isArray(r))return f(o,r,n??!1,(t,a,m)=>u.nansum(t,a,m));const e=u.nansum(b(o),r,n);return typeof e=="number"||e instanceof A?e:y(e)}function U(o,r,n){if(Array.isArray(r))return f(o,r,n??!1,(t,a,m)=>u.nanprod(t,a,m));const e=u.nanprod(b(o),r,n);return typeof e=="number"||e instanceof A?e:y(e)}function V(o,r,n){if(Array.isArray(r))return f(o,r,n??!1,(t,a,m)=>u.nanmean(t,a,m));const e=u.nanmean(b(o),r,n);return typeof e=="number"||e instanceof A?e:y(e)}function W(o,r,n,e){if(Array.isArray(r))return f(o,r,e??!1,(a,m,l)=>u.nanvar(a,m,n,l));const t=u.nanvar(b(o),r,n,e);return typeof t=="number"?t:y(t)}function X(o,r,n,e){if(Array.isArray(r))return f(o,r,e??!1,(a,m,l)=>u.nanstd(a,m,n,l));const t=u.nanstd(b(o),r,n,e);return typeof t=="number"?t:y(t)}function Y(o,r,n){if(Array.isArray(r))return f(o,r,n??!1,(t,a,m)=>u.nanmin(t,a,m));const e=u.nanmin(b(o),r,n);return typeof e=="number"||e instanceof A?e:y(e)}function Z(o,r,n){if(Array.isArray(r))return f(o,r,n??!1,(t,a,m)=>u.nanmax(t,a,m));const e=u.nanmax(b(o),r,n);return typeof e=="number"||e instanceof A?e:y(e)}function $(o,r){const n=u.nanargmin(b(o),r);return typeof n=="number"?n:y(n)}function k(o,r){const n=u.nanargmax(b(o),r);return typeof n=="number"?n:y(n)}function rr(o,r){return y(u.nancumsum(b(o),r))}function er(o,r){return y(u.nancumprod(b(o),r))}function nr(o,r,n){const e=u.nanmedian(b(o),r,n);return typeof e=="number"?e:y(e)}function or(o,r,n,e){const t=u.nanquantile(b(o),r,n,e);return typeof t=="number"?t:y(t)}function tr(o,r,n,e){const t=u.nanpercentile(b(o),r,n,e);return typeof t=="number"?t:y(t)}export{K as all,z as amax,M as amin,L as any,B as argmax,w as argmin,J as average,d as cumprod,S as cumsum,Q as cumulative_prod,P as cumulative_sum,x as max,h as mean,G as median,g as min,k as nanargmax,$ as nanargmin,er as nancumprod,rr as nancumsum,Z as nanmax,V as nanmean,nr as nanmedian,Y as nanmin,tr as nanpercentile,U as nanprod,or as nanquantile,X as nanstd,R as nansum,W as nanvar,H as percentile,O as prod,j as ptp,I as quantile,F as std,_ as sum,E as var_,v as variance};

package/dist/esm/core/sets.js

@@ -1 +1 @@
-import*as i from"../common/ops/sets";import{toStorage as o,fromStorage as n}from"./types";function A(e,r=!1,s=!1,t=!1){const u=i.unique(o(e),r,s,t);return"values"in u?{values:n(u.values),indices:u.indices?n(u.indices):void 0,inverse:u.inverse?n(u.inverse):void 0,counts:u.counts?n(u.counts):void 0}:n(u)}function C(e,r){return n(i.in1d(o(e),o(r)))}function D(e,r){return n(i.intersect1d(o(e),o(r)))}function N(e,r){return n(i.isin(o(e),o(r)))}function d(e,r){return n(i.setdiff1d(o(e),o(r)))}function l(e,r){return n(i.setxor1d(o(e),o(r)))}function f(e,r){return n(i.union1d(o(e),o(r)))}function v(e,r){return n(i.trim_zeros(o(e),r))}function _(e){const r=i.unique_all(o(e));return{values:n(r.values),indices:n(r.indices),inverse_indices:n(r.inverse_indices),counts:n(r.counts)}}function p(e){const r=i.unique_counts(o(e));return{values:n(r.values),counts:n(r.counts)}}function x(e){const r=i.unique_inverse(o(e));return{values:n(r.values),inverse_indices:n(r.inverse_indices)}}function q(e){return n(i.unique_values(o(e)))}export{C as in1d,D as intersect1d,N as isin,d as setdiff1d,l as setxor1d,v as trim_zeros,f as union1d,A as unique,_ as unique_all,p as unique_counts,x as unique_inverse,q as unique_values};
+import*as u from"../common/ops/sets";import{toStorage as o,fromStorage as n}from"./types";function C(e,r=!1,s=!1,t=!1,a){const i=u.unique(o(e),r,s,t,a);return"values"in i?{values:n(i.values),indices:i.indices?n(i.indices):void 0,inverse:i.inverse?n(i.inverse):void 0,counts:i.counts?n(i.counts):void 0}:n(i)}function D(e,r){return n(u.in1d(o(e),o(r)))}function N(e,r){return n(u.intersect1d(o(e),o(r)))}function d(e,r){return n(u.isin(o(e),o(r)))}function l(e,r){return n(u.setdiff1d(o(e),o(r)))}function f(e,r){return n(u.setxor1d(o(e),o(r)))}function v(e,r){return n(u.union1d(o(e),o(r)))}function _(e,r){return n(u.trim_zeros(o(e),r))}function p(e){const r=u.unique_all(o(e));return{values:n(r.values),indices:n(r.indices),inverse_indices:n(r.inverse_indices),counts:n(r.counts)}}function x(e){const r=u.unique_counts(o(e));return{values:n(r.values),counts:n(r.counts)}}function q(e){const r=u.unique_inverse(o(e));return{values:n(r.values),inverse_indices:n(r.inverse_indices)}}function m(e){return n(u.unique_values(o(e)))}export{D as in1d,N as intersect1d,d as isin,l as setdiff1d,f as setxor1d,_ as trim_zeros,v as union1d,C as unique,p as unique_all,x as unique_counts,q as unique_inverse,m as unique_values};

package/dist/esm/core/shape-extra.js

@@ -1 +1 @@
-import{NDArrayCore as z}from"../common/ndarray-core";import{array as C,zeros as I}from"./creation";import{concatenate as T,flatten as S}from"./shape";function F(l,n,a){const b=n instanceof z?n:C(Array.isArray(n)?n:[n]);if(a===void 0){const r=S(l),x=S(b);return T([r,x],0)}return T([l,b],a)}function k(l,n,a){const b=Array.isArray(n)?n:[n],r=[...l.shape],x=l.storage;if(a===void 0){const e=S(l),f=e.storage,o=[],w=new Set(b.map(d=>d<0?e.size+d:d));for(let d=0;d<e.size;d++)w.has(d)||o.push(f.iget(d));return C(o,l.dtype)}const i=a<0?r.length+a:a,p=r[i],A=new Set(b.map(e=>e<0?p+e:e)),D=p-A.size,c=[...r];c[i]=D;const u=I(c,l.dtype),s=u.storage,m=[];let g=1;for(let e=r.length-1;e>=0;e--)m.unshift(g),g*=r[e];const t=[];g=1;for(let e=c.length-1;e>=0;e--)t.unshift(g),g*=c[e];const h=r.reduce((e,f)=>e*f,1);for(let e=0;e<h;e++){const f=[];let o=e;for(let y=0;y<r.length;y++)f.push(Math.floor(o/m[y])),o=o%m[y];if(A.has(f[i]))continue;let w=0;for(const y of A)y<f[i]&&w++;const d=[...f];d[i]-=w;let N=0;for(let y=0;y<c.length;y++)N+=d[y]*t[y];s.iset(N,x.iget(e))}return u}function V(l,n,a,b){const r=Array.isArray(n)?n:[n],x=a instanceof z?a:C(Array.isArray(a)?a:[a]);if(b===void 0){const i=S(l),p=S(x),A=i.data,D=p.data,c=Array.from(A);if(r.length===1){const u=r[0],s=u<0?i.size+u:u,m=Array.from(D);c.splice(s,0,...m)}else{const u=r.map((s,m)=>({idx:s<0?i.size+s:s,valIdx:m})).sort((s,m)=>s.idx-m.idx);for(let s=0;s<u.length;s++){const{idx:m,valIdx:g}=u[s],t=D[g%D.length];c.splice(m+s,0,t)}}return C(c,l.dtype)}throw new Error("insert along axis not fully implemented in standalone")}function B(l,n,a="constant",b=0){const r=[...l.shape],x=r.length;let i;typeof n=="number"?i=Array(x).fill([n,n]):Array.isArray(n)&&typeof n[0]=="number"?i=Array(x).fill(n):i=n;const p=r.map((t,h)=>t+i[h][0]+i[h][1]);if(a!=="constant")throw new Error(`pad mode '${a}' not fully implemented in standalone`);const A=I(p,l.dtype),D=A.storage;if(b!==0)for(let t=0;t<A.size;t++)D.iset(t,b);const c=l.storage,u=[];let s=1;for(let t=r.length-1;t>=0;t--)u.unshift(s),s*=r[t];const m=[];s=1;for(let t=p.length-1;t>=0;t--)m.unshift(s),s*=p[t];const g=r.reduce((t,h)=>t*h,1);for(let t=0;t<g;t++){const h=[];let e=t;for(let o=0;o<r.length;o++)h.push(Math.floor(e/u[o])),e=e%u[o];let f=0;for(let o=0;o<p.length;o++)f+=(h[o]+i[o][0])*m[o];D.iset(f,c.iget(t))}return A}export{F as append,k as delete,k as delete_,V as insert,B as pad};
+import{NDArrayCore as N}from"../common/ndarray-core";import{array as C,zeros as I}from"./creation";import{concatenate as T,flatten as S}from"./shape";import{wasmPad2D as k}from"../common/wasm/pad";function V(i,n,o){const c=n instanceof N?n:C(Array.isArray(n)?n:[n]);if(o===void 0){const r=S(i),p=S(c);return T([r,p],0)}return T([i,c],o)}function L(i,n,o){const c=Array.isArray(n)?n:[n],r=[...i.shape],p=i.storage;if(o===void 0){const t=S(i),y=t.storage,a=[],w=new Set(c.map(d=>d<0?t.size+d:d));for(let d=0;d<t.size;d++)w.has(d)||a.push(y.iget(d));return C(a,i.dtype)}const l=o<0?r.length+o:o,A=r[l],g=new Set(c.map(t=>t<0?A+t:t)),D=A-g.size,f=[...r];f[l]=D;const u=I(f,i.dtype),s=u.storage,m=[];let h=1;for(let t=r.length-1;t>=0;t--)m.unshift(h),h*=r[t];const e=[];h=1;for(let t=f.length-1;t>=0;t--)e.unshift(h),h*=f[t];const x=r.reduce((t,y)=>t*y,1);for(let t=0;t<x;t++){const y=[];let a=t;for(let b=0;b<r.length;b++)y.push(Math.floor(a/m[b])),a=a%m[b];if(g.has(y[l]))continue;let w=0;for(const b of g)b<y[l]&&w++;const d=[...y];d[l]-=w;let z=0;for(let b=0;b<f.length;b++)z+=d[b]*e[b];s.iset(z,p.iget(t))}return u}function B(i,n,o,c){const r=Array.isArray(n)?n:[n],p=o instanceof N?o:C(Array.isArray(o)?o:[o]);if(c===void 0){const l=S(i),A=S(p),g=l.data,D=A.data,f=Array.from(g);if(r.length===1){const u=r[0],s=u<0?l.size+u:u,m=Array.from(D);f.splice(s,0,...m)}else{const u=r.map((s,m)=>({idx:s<0?l.size+s:s,valIdx:m})).sort((s,m)=>s.idx-m.idx);for(let s=0;s<u.length;s++){const{idx:m,valIdx:h}=u[s],e=D[h%D.length];f.splice(m+s,0,e)}}return C(f,i.dtype)}throw new Error("insert along axis not fully implemented in standalone")}function W(i,n,o="constant",c=0){const r=[...i.shape],p=r.length;let l;typeof n=="number"?l=Array(p).fill([n,n]):Array.isArray(n)&&typeof n[0]=="number"?l=Array(p).fill(n):l=n;const A=r.map((e,x)=>e+l[x][0]+l[x][1]);if(o==="constant"&&c===0&&p===2&&typeof n=="number"){const e=k(i.storage,n);if(e)return new N(e)}if(o!=="constant")throw new Error(`pad mode '${o}' not fully implemented in standalone`);const g=I(A,i.dtype),D=g.storage;if(c!==0)for(let e=0;e<g.size;e++)D.iset(e,c);const f=i.storage,u=[];let s=1;for(let e=r.length-1;e>=0;e--)u.unshift(s),s*=r[e];const m=[];s=1;for(let e=A.length-1;e>=0;e--)m.unshift(s),s*=A[e];const h=r.reduce((e,x)=>e*x,1);for(let e=0;e<h;e++){const x=[];let t=e;for(let a=0;a<r.length;a++)x.push(Math.floor(t/u[a])),t=t%u[a];let y=0;for(let a=0;a<A.length;a++)y+=(x[a]+l[a][0])*m[a];D.iset(y,f.iget(e))}return g}export{V as append,L as delete,L as delete_,B as insert,W as pad};

package/dist/esm/core/sorting.js

@@ -1 +1 @@
-import*as t from"../common/ops/sorting";import{toStorage as o,fromStorage as a,fromStorageArray as y}from"./types";function D(r,e=-1){return a(t.sort(o(r),e))}function N(r,e=-1){return a(t.argsort(o(r),e))}function i(r){return a(t.lexsort(r.map(e=>o(e))))}function f(r,e,n=-1){return a(t.partition(o(r),e,n))}function p(r,e,n=-1){return a(t.argpartition(o(r),e,n))}function c(r){return a(t.sort_complex(o(r)))}function s(r){return y(t.nonzero(o(r)))}function m(r){return a(t.argwhere(o(r)))}function x(r){return a(t.flatnonzero(o(r)))}function g(r,e,n){const u=t.where(o(r),e?o(e):void 0,n?o(n):void 0);return Array.isArray(u)?y(u):a(u)}function l(r,e,n="left"){return a(t.searchsorted(o(r),o(e),n))}function b(r,e){return a(t.extract(o(r),o(e)))}function h(r,e){const n=t.count_nonzero(o(r),e);return typeof n=="number"?n:a(n)}export{p as argpartition,N as argsort,m as argwhere,h as count_nonzero,b as extract,x as flatnonzero,i as lexsort,s as nonzero,f as partition,l as searchsorted,D as sort,c as sort_complex,g as where};
+import*as t from"../common/ops/sorting";import{toStorage as o,fromStorage as n,fromStorageArray as D}from"./types";function p(r,e=-1){return n(t.sort(o(r),e))}function m(r,e=-1){return n(t.argsort(o(r),e))}function x(r){return n(t.lexsort(r.map(e=>o(e))))}function g(r,e,a=-1){return n(t.partition(o(r),e,a))}function b(r,e,a=-1){return n(t.argpartition(o(r),e,a))}function l(r){return n(t.sort_complex(o(r)))}function d(r){return D(t.nonzero(o(r)))}function h(r){return n(t.argwhere(o(r)))}function z(r){return n(t.flatnonzero(o(r)))}function S(r,e,a){const A=t.where(o(r),e?o(e):void 0,a?o(a):void 0);return Array.isArray(A)?D(A):n(A)}function _(r,e,a="left"){return n(t.searchsorted(o(r),o(e),a))}function w(r,e){return n(t.extract(o(r),o(e)))}function O(r,e){if(Array.isArray(e)){const A=o(r).ndim,C=[...e.map(u=>u<0?A+u:u)].sort((u,y)=>y-u);let i=o(r);for(const u of C){const y=t.count_nonzero(i,u);if(typeof y=="number")return y;i=y}return n(i)}const a=t.count_nonzero(o(r),e);return typeof a=="number"?a:n(a)}export{b as argpartition,m as argsort,h as argwhere,O as count_nonzero,w as extract,z as flatnonzero,x as lexsort,d as nonzero,g as partition,_ as searchsorted,p as sort,l as sort_complex,S as where};

package/dist/esm/full/index.js

@@ -1 +1 @@
-import*as o from"../core";import{NDArray as N}from"./ndarray";import{NDArrayCore as i}from"../common/ndarray-core";import{ArrayStorage as s}from"../common/storage";import{Complex as l}from"../common/complex";import{DEFAULT_DTYPE as y,getTypedArrayConstructor as b,promoteDTypes as d,isComplexDType as _}from"../common/dtype";const n=r=>{if(r instanceof N)return r;const e=r.base?n(r.base):void 0;return N.fromStorage(r.storage,e)};import{NDArray as Bt,meshgrid as Lt}from"./ndarray";import{NDArrayCore as Yt}from"../common/ndarray-core";import{Complex as Ht}from"../common/complex";function T(r,e){return n(o.broadcast_to(r,e))}function q(r,e,t){return n(o.take(r,e,t))}function z(r,e,t){return n(o.take_along_axis(r,e,t))}function j(r,e){return n(o.choose(r,e))}function S(r,e,t){return n(o.compress(r,e,t))}function F(r,e,t=0){return n(o.iindex(r,e,t))}function I(r,e,t){return n(o.bindex(r,e,t))}function P(r,e,t=0){return n(o.select(r,e,t))}function E(r,e="int32"){return n(o.indices(r,e))}function M(r,e,t="raise"){return n(o.ravel_multi_index(r,e,t))}function B(r,e){return o.unravel_index(r,e).map(n)}function L(r,e=2){return o.diag_indices(r,e).map(n)}function U(r){return o.diag_indices_from(r).map(n)}function Y(r,e=0,t){return o.tril_indices(r,e,t).map(n)}function G(r,e=0){return o.tril_indices_from(r,e).map(n)}function H(r,e=0,t){return o.triu_indices(r,e,t).map(n)}function J(r,e=0){return o.triu_indices_from(r,e).map(n)}function K(r,e,t=0){return o.mask_indices(r,e,t).map(n)}function O(r,e,t){const a=u=>r(n(u));return n(o.apply_along_axis(a,e,t))}function Q(r,e,t){const a=(u,m)=>r(n(u),m);return n(o.apply_over_axes(a,e,t))}function R(r,e){return n(o.add(r,e))}function V(r,e){return n(o.subtract(r,e))}function W(r,e){return n(o.multiply(r,e))}function X(r){return n(o.sqrt(r))}function Z(r,e){return n(o.power(r,e))}function $(r){return n(o.exp(r))}function rr(r){return n(o.exp2(r))}function er(r){return n(o.expm1(r))}function nr(r){return n(o.log(r))}function or(r){return n(o.log2(r))}function tr(r){return n(o.log10(r))}function ar(r){return n(o.log1p(r))}function ur(r,e){return n(o.logaddexp(r,e))}function ir(r,e){return n(o.logaddexp2(r,e))}function yr(r){return n(o.absolute(r))}function Dr(r){return n(o.negative(r))}function Ar(r){return n(o.sign(r))}function Nr(r,e){return n(o.mod(r,e))}function cr(r,e){return n(o.divide(r,e))}function pr(r,e){return n(o.floor_divide(r,e))}function mr(r){return n(o.positive(r))}function fr(r){return n(o.reciprocal(r))}function sr(r){return n(o.cbrt(r))}function br(r){return n(o.fabs(r))}function Cr(r,e){const t=o.divmod(r,e);return[n(t[0]),n(t[1])]}function xr(r){return n(o.square(r))}function lr(r,e){return n(o.remainder(r,e))}function dr(r,e){return n(o.heaviside(r,e))}function _r(r,e){return n(o.float_power(r,e))}function gr(r,e){return n(o.fmod(r,e))}function vr(r){const e=o.frexp(r);return[n(e[0]),n(e[1])]}function hr(r,e){return n(o.gcd(r,e))}function kr(r,e){return n(o.lcm(r,e))}function wr(r,e){return n(o.ldexp(r,e))}function Tr(r){const e=o.modf(r);return[n(e[0]),n(e[1])]}function qr(r,e,t){return n(o.clip(r,e,t))}function zr(r,e){return n(o.maximum(r,e))}function jr(r,e){return n(o.minimum(r,e))}function Sr(r,e){return n(o.fmax(r,e))}function Fr(r,e){return n(o.fmin(r,e))}function Ir(r,e=0,t,a){return n(o.nan_to_num(r,e,t,a))}function Pr(r,e,t,a,u){return n(o.interp(r,e,t,a,u))}function Er(r,e=Math.PI,t=-1,a=2*Math.PI){return n(o.unwrap(r,e,t,a))}function Mr(r){return n(o.sinc(r))}function Br(r){return n(o.i0(r))}function Lr(r,e){return n(o.bitwise_and(r,e))}function Ur(r,e){return n(o.bitwise_or(r,e))}function Yr(r,e){return n(o.bitwise_xor(r,e))}function Gr(r){return n(o.bitwise_not(r))}function Hr(r){return n(o.invert(r))}function Jr(r,e){return n(o.left_shift(r,e))}function Kr(r,e){return n(o.right_shift(r,e))}function Or(r,e,t){return n(o.packbits(r,e,t))}function Qr(r,e,t,a){return n(o.unpackbits(r,e,t,a))}function Rr(r){return n(o.bitwise_count(r))}function Vr(r){return n(o.bitwise_invert(r))}function Wr(r,e){return n(o.bitwise_left_shift(r,e))}function Xr(r,e){return n(o.bitwise_right_shift(r,e))}function Zr(r){return n(o.real(r))}function $r(r){return n(o.imag(r))}function re(r){return n(o.conj(r))}function ee(r,e){return n(o.angle(r,e))}function ne(r,e=y){return n(o.zeros(r,e))}function oe(r,e=y){return n(o.ones(r,e))}function te(r,e=y){return n(o.empty(r,e))}function ae(r,e,t){return n(o.full(r,e,t))}function ue(r,e){return n(o.array(r,e))}function ie(r,e,t=1,a=y){return n(o.arange(r,e,t,a))}function ye(r,e,t=50,a=y){return n(o.linspace(r,e,t,a))}function De(r,e,t=50,a=10,u=y){return n(o.logspace(r,e,t,a,u))}function Ae(r,e,t=50,a=y){return n(o.geomspace(r,e,t,a))}function Ne(r,e,t=0,a=y){return n(o.eye(r,e,t,a))}function ce(r,e=y){return n(o.identity(r,e))}function pe(r,e){return n(o.asarray(r,e))}function me(r,e){return n(o.zeros_like(r,e))}function fe(r,e){return n(o.ones_like(r,e))}function se(r,e){return n(o.empty_like(r,e))}function be(r,e,t){return n(o.full_like(r,e,t))}function Ce(r){return n(o.copy(r))}function xe(r,e){return n(o.asanyarray(r,e))}function le(r,e){return n(o.ascontiguousarray(r,e))}function de(r,e){return n(o.asfortranarray(r,e))}function _e(r,e){return n(o.asarray_chkfinite(r,e))}function ge(r,e,t){return n(o.require(r,e,t))}function ve(r,e=0){return n(o.diag(r,e))}function he(r,e=0){return n(o.diagflat(r,e))}function ke(r,e,t=0,a=y){return n(o.tri(r,e,t,a))}function we(r,e=0){return n(o.tril(r,e))}function Te(r,e=0){return n(o.triu(r,e))}function qe(r,e,t=!1){return n(o.vander(r,e,t))}function ze(r,e=y,t=-1,a=0){return n(o.frombuffer(r,e,t,a))}function je(r,e,t=y){return n(o.fromfunction(r,e,t))}function Se(r,e=y,t=-1){return n(o.fromiter(r,e,t))}function Fe(r,e=y,t=-1,a){return n(o.fromstring(r,e,t,a))}function Ie(r,e=y,t=-1,a=""){return n(o.fromfile(r,e,t,a))}function Pe(r,e,t){return n(o.diff(r,e,t))}function Ee(r,e,t){return n(o.ediff1d(r,e,t))}function Me(r,e){const t=o.dot(r,e);return t instanceof i?n(t):t}function Be(r,e=0,t=0,a=1){return n(o.diagonal(r,e,t,a))}function Le(r,e){return n(o.kron(r,e))}function Ue(r,e){return n(o.transpose(r,e))}function Ye(r,e){const t=o.inner(r,e);return t instanceof i?n(t):t}function Ge(r,e){return n(o.outer(r,e))}function He(r,e,t=2){const a=o.tensordot(r,e,t);return a instanceof i?n(a):a}function Je(r,e,t=-1){const a=o.vecdot(r,e,t);return a instanceof i?n(a):a}function Ke(r){return n(o.matrix_transpose(r))}function Oe(r,e){return n(o.permute_dims(r,e))}function Qe(r,e){return n(o.matvec(r,e))}function Re(r,e){return n(o.vecmat(r,e))}function Ve(r,e,t=-1,a=-1,u=-1,m){const D=o.cross(r,e,t,a,u,m),A=d(r.dtype,e.dtype);if(D instanceof l){const c=b(A==="complex64"?"float32":"float64"),p=new c(2);return p[0]=D.re,p[1]=D.im,N.fromStorage(s.fromData(p,[],A))}if(typeof D=="number"){if(_(A)){const x=b(A==="complex64"?"float32":"float64"),f=new x(2);return f[0]=D,f[1]=0,N.fromStorage(s.fromData(f,[],A))}const C=b(A),c=new C(1);return c[0]=D,N.fromStorage(s.fromData(c,[],A))}return n(D)}function We(r,e){return n(o.matmul(r,e))}function Xe(r,e){return n(o.logical_and(r,e))}function Ze(r,e){return n(o.logical_or(r,e))}function $e(r){return n(o.logical_not(r))}function rn(r,e){return n(o.logical_xor(r,e))}function en(r){return n(o.isfinite(r))}function nn(r){return n(o.isinf(r))}function on(r){return n(o.isnan(r))}function tn(r){return n(o.isnat(r))}function an(r){return n(o.isneginf(r))}function un(r){return n(o.isposinf(r))}function yn(r){return n(o.iscomplex(r))}function Dn(r){return n(o.isreal(r))}function An(r,e){return n(o.real_if_close(r,e))}function Nn(r,e){return n(o.copysign(r,e))}function cn(r){return n(o.signbit(r))}function pn(r,e){return n(o.nextafter(r,e))}function mn(r){return n(o.spacing(r))}function fn(r,e){return n(o.greater(r,e))}function sn(r,e){return n(o.greater_equal(r,e))}function bn(r,e){return n(o.less(r,e))}function Cn(r,e){return n(o.less_equal(r,e))}function xn(r,e){return n(o.equal(r,e))}function ln(r,e){return n(o.not_equal(r,e))}function dn(r,e,t=1e-5,a=1e-8){return n(o.isclose(r,e,t,a))}function _n(r){return n(o.poly(r))}function gn(r,e){return n(o.polyadd(r,e))}function vn(r,e=1){return n(o.polyder(r,e))}function hn(r,e){const t=o.polydiv(r,e);return[n(t[0]),n(t[1])]}function kn(r,e,t){return n(o.polyfit(r,e,t))}function wn(r,e=1,t=0){return n(o.polyint(r,e,t))}function Tn(r,e){return n(o.polymul(r,e))}function qn(r,e){return n(o.polysub(r,e))}function zn(r,e){const t=o.polyval(r,e);return t instanceof i?n(t):t}function jn(r){return n(o.roots(r))}function Sn(r,e,t){const a=o.sum(r,e,t);return a instanceof i?n(a):a}function Fn(r,e,t){const a=o.mean(r,e,t);return a instanceof i?n(a):a}function In(r,e,t){const a=o.prod(r,e,t);return a instanceof i?n(a):a}function Pn(r,e,t){const a=o.max(r,e,t);return a instanceof i?n(a):a}function En(r,e,t){const a=o.min(r,e,t);return a instanceof i?n(a):a}function Mn(r,e,t){const a=o.ptp(r,e,t);return a instanceof i?n(a):a}function Bn(r,e){const t=o.argmin(r,e);return t instanceof i?n(t):t}function Ln(r,e){const t=o.argmax(r,e);return t instanceof i?n(t):t}function Un(r,e,t,a){const u=o.variance(r,e,t,a);return u instanceof i?n(u):u}function Yn(r,e,t,a){const u=o.std(r,e,t,a);return u instanceof i?n(u):u}function Gn(r,e,t){const a=o.median(r,e,t);return a instanceof i?n(a):a}function Hn(r,e,t,a){const u=o.percentile(r,e,t,a);return u instanceof i?n(u):u}function Jn(r,e,t,a){const u=o.quantile(r,e,t,a);return u instanceof i?n(u):u}function Kn(r,e,t,a){const u=o.average(r,e,t,a);return u instanceof i?n(u):u}function On(r,e,t){const a=o.all(r,e,t);return a instanceof i?n(a):a}function Qn(r,e,t){const a=o.any(r,e,t);return a instanceof i?n(a):a}function Rn(r,e){return n(o.cumsum(r,e))}function Vn(r,e){return n(o.cumprod(r,e))}function Wn(r,e,t){const a=o.nansum(r,e,t);return a instanceof i?n(a):a}function Xn(r,e,t){const a=o.nanprod(r,e,t);return a instanceof i?n(a):a}function Zn(r,e,t){const a=o.nanmean(r,e,t);return a instanceof i?n(a):a}function $n(r,e,t,a){const u=o.nanvar(r,e,t,a);return u instanceof i?n(u):u}function ro(r,e,t,a){const u=o.nanstd(r,e,t,a);return u instanceof i?n(u):u}function eo(r,e,t){const a=o.nanmin(r,e,t);return a instanceof i?n(a):a}function no(r,e,t){const a=o.nanmax(r,e,t);return a instanceof i?n(a):a}function oo(r,e){const t=o.nanargmin(r,e);return t instanceof i?n(t):t}function to(r,e){const t=o.nanargmax(r,e);return t instanceof i?n(t):t}function ao(r,e){return n(o.nancumsum(r,e))}function uo(r,e){return n(o.nancumprod(r,e))}function io(r,e,t){const a=o.nanmedian(r,e,t);return a instanceof i?n(a):a}function yo(r,e,t,a){const u=o.nanquantile(r,e,t,a);return u instanceof i?n(u):u}function Do(r,e,t,a){const u=o.nanpercentile(r,e,t,a);return u instanceof i?n(u):u}function Ao(r,e=0){return n(o.around(r,e))}function No(r,e=0){return n(o.round(r,e))}function co(r){return n(o.ceil(r))}function po(r){return n(o.fix(r))}function mo(r){return n(o.floor(r))}function fo(r){return n(o.rint(r))}function so(r){return n(o.trunc(r))}function bo(r,e){return n(o.in1d(r,e))}function Co(r,e){return n(o.intersect1d(r,e))}function xo(r,e){return n(o.isin(r,e))}function lo(r,e){return n(o.setdiff1d(r,e))}function _o(r,e){return n(o.setxor1d(r,e))}function go(r,e){return n(o.union1d(r,e))}function vo(r,e){return n(o.trim_zeros(r,e))}function ho(r){return n(o.unique_values(r))}function ko(r,e,t){return n(o.append(r,e,t))}function wo(r,e,t){return n(o.delete(r,e,t))}function To(r,e,t,a){return n(o.insert(r,e,t,a))}function qo(r,e,t="constant",a=0){return n(o.pad(r,e,t,a))}function zo(r,e){return n(o.reshape(r,e))}function jo(r){return n(o.flatten(r))}function So(r){return n(o.ravel(r))}function Fo(r,e){return n(o.squeeze(r,e))}function Io(r,e){return n(o.expand_dims(r,e))}function Po(r,e,t){return n(o.swapaxes(r,e,t))}function Eo(r,e,t){return n(o.moveaxis(r,e,t))}function Mo(r,e,t=0){return n(o.rollaxis(r,e,t))}function Bo(r,e=0){return n(o.concatenate(r,e))}function Lo(r,e=0){return n(o.stack(r,e))}function Uo(r){return n(o.vstack(r))}function Yo(r){return n(o.hstack(r))}function Go(r){return n(o.dstack(r))}function Ho(r,e=0){return n(o.concat(r,e))}function Jo(r){return n(o.column_stack(r))}function Ko(r){return n(o.block(r))}function Oo(r,e=0){return o.unstack(r,e).map(n)}function Qo(r,e){return n(o.tile(r,e))}function Ro(r,e,t){return n(o.repeat(r,e,t))}function Vo(r,e){return n(o.flip(r,e))}function Wo(r){return n(o.fliplr(r))}function Xo(r){return n(o.flipud(r))}function Zo(r,e=1,t=[0,1]){return n(o.rot90(r,e,t))}function $o(r,e,t){return n(o.roll(r,e,t))}function rt(r,e){return n(o.resize(r,e))}function et(r,e=-1){return n(o.sort(r,e))}function nt(r,e=-1){return n(o.argsort(r,e))}function ot(r){return n(o.lexsort(r))}function tt(r,e,t=-1){return n(o.partition(r,e,t))}function at(r,e,t=-1){return n(o.argpartition(r,e,t))}function ut(r){return n(o.sort_complex(r))}function it(r){return o.nonzero(r).map(n)}function yt(r){return n(o.argwhere(r))}function Dt(r){return n(o.flatnonzero(r))}function At(r,e){return n(o.extract(r,e))}function Nt(r,e,t){return n(o.bincount(r,e,t))}function ct(r,e,t){return n(o.digitize(r,e,t))}function pt(r,e,t){return n(o.correlate(r,e,t))}function mt(r,e,t){return n(o.convolve(r,e,t))}function ft(r,e,t,a,u){return n(o.cov(r,e,t,a,u))}function st(r,e,t){return n(o.corrcoef(r,e,t))}function bt(r,e,t,a){const u=o.trapezoid(r,e,t,a);return u instanceof i?n(u):u}function Ct(r){return n(o.sin(r))}function xt(r){return n(o.cos(r))}function lt(r){return n(o.tan(r))}function dt(r){return n(o.arcsin(r))}function _t(r){return n(o.arccos(r))}function gt(r){return n(o.arctan(r))}function vt(r,e){return n(o.arctan2(r,e))}function ht(r,e){return n(o.hypot(r,e))}function kt(r){return n(o.degrees(r))}function wt(r){return n(o.radians(r))}function Tt(r){return n(o.deg2rad(r))}function qt(r){return n(o.rad2deg(r))}function zt(r){return n(o.sinh(r))}function jt(r){return n(o.cosh(r))}function St(r){return n(o.tanh(r))}function Ft(r){return n(o.arcsinh(r))}function It(r){return n(o.arccosh(r))}function Pt(r){return n(o.arctanh(r))}import{broadcast_arrays as Kt,broadcast_shapes as Ot,put as Qt,put_along_axis as Rt,place as Vt,putmask as Wt,copyto as Xt,ix_ as Zt,fill_diagonal as $t,array_equal as ra,array_equiv as ea,may_share_memory as na,shares_memory as oa,geterr as ta,seterr as aa,pow as ua,abs as ia,true_divide as ya,conjugate as Da,array2string as Aa,array_repr as Na,array_str as ca,set_printoptions as pa,get_printoptions as ma,printoptions as fa,format_float_positional as sa,format_float_scientific as ba,base_repr as Ca,binary_repr as xa,gradient as la,trace as da,einsum as _a,einsum_path as ga,vdot as va,linalg as ha,iscomplexobj as ka,isrealobj as wa,isfortran as Ta,isscalar as qa,iterable as za,isdtype as ja,promote_types as Sa,allclose as Fa,amax as Ia,amin as Pa,var_ as Ea,cumulative_sum as Ma,cumulative_prod as Ba,round_ as La,unique as Ua,unique_all as Ya,unique_counts as Ga,unique_inverse as Ha,split as Ja,array_split as Ka,vsplit as Oa,hsplit as Qa,dsplit as Ra,atleast_1d as Va,atleast_2d as Wa,atleast_3d as Xa,row_stack as Za,where as $a,searchsorted as ru,count_nonzero as eu,histogram as nu,histogram2d as ou,histogramdd as tu,histogram_bin_edges as au,asin as uu,acos as iu,atan as yu,atan2 as Du,asinh as Au,acosh as Nu,atanh as cu,can_cast as pu,common_type as mu,result_type as fu,min_scalar_type as su,issubdtype as bu,typename as Cu,mintypecode as xu,ndim as lu,shape as du,size as _u,item as gu,tolist as vu,tobytes as hu,byteswap as ku,view as wu,tofile as Tu,fill as qu}from"../core";export{Ht as Complex,Bt as NDArray,Yt as NDArrayCore,ia as abs,yr as absolute,iu as acos,Nu as acosh,R as add,On as all,Fa as allclose,Ia as amax,Pa as amin,ee as angle,Qn as any,ko as append,O as apply_along_axis,Q as apply_over_axes,ie as arange,_t as arccos,It as arccosh,dt as arcsin,Ft as arcsinh,gt as arctan,vt as arctan2,Pt as arctanh,Ln as argmax,Bn as argmin,at as argpartition,nt as argsort,yt as argwhere,Ao as around,ue as array,Aa as array2string,ra as array_equal,ea as array_equiv,Na as array_repr,Ka as array_split,ca as array_str,xe as asanyarray,pe as asarray,_e as asarray_chkfinite,le as ascontiguousarray,de as asfortranarray,uu as asin,Au as asinh,yu as atan,Du as atan2,cu as atanh,Va as atleast_1d,Wa as atleast_2d,Xa as atleast_3d,Kn as average,Ca as base_repr,xa as binary_repr,Nt as bincount,I as bindex,Lr as bitwise_and,Rr as bitwise_count,Vr as bitwise_invert,Wr as bitwise_left_shift,Gr as bitwise_not,Ur as bitwise_or,Xr as bitwise_right_shift,Yr as bitwise_xor,Ko as block,Kt as broadcast_arrays,Ot as broadcast_shapes,T as broadcast_to,ku as byteswap,pu as can_cast,sr as cbrt,co as ceil,j as choose,qr as clip,Jo as column_stack,mu as common_type,S as compress,Ho as concat,Bo as concatenate,re as conj,Da as conjugate,mt as convolve,Ce as copy,Nn as copysign,Xt as copyto,st as corrcoef,pt as correlate,xt as cos,jt as cosh,eu as count_nonzero,ft as cov,Ve as cross,Vn as cumprod,Rn as cumsum,Ba as cumulative_prod,Ma as cumulative_sum,Tt as deg2rad,kt as degrees,wo as delete_,ve as diag,L as diag_indices,U as diag_indices_from,he as diagflat,Be as diagonal,Pe as diff,ct as digitize,cr as divide,Cr as divmod,Me as dot,Ra as dsplit,Go as dstack,Ee as ediff1d,_a as einsum,ga as einsum_path,te as empty,se as empty_like,xn as equal,$ as exp,rr as exp2,Io as expand_dims,er as expm1,At as extract,Ne as eye,br as fabs,qu as fill,$t as fill_diagonal,po as fix,Dt as flatnonzero,jo as flatten,Vo as flip,Wo as fliplr,Xo as flipud,_r as float_power,mo as floor,pr as floor_divide,Sr as fmax,Fr as fmin,gr as fmod,sa as format_float_positional,ba as format_float_scientific,vr as frexp,ze as frombuffer,Ie as fromfile,je as fromfunction,Se as fromiter,Fe as fromstring,ae as full,be as full_like,hr as gcd,Ae as geomspace,ma as get_printoptions,ta as geterr,la as gradient,fn as greater,sn as greater_equal,dr as heaviside,nu as histogram,ou as histogram2d,au as histogram_bin_edges,tu as histogramdd,Qa as hsplit,Yo as hstack,ht as hypot,Br as i0,ce as identity,F as iindex,$r as imag,bo as in1d,E as indices,Ye as inner,To as insert,Pr as interp,Co as intersect1d,Hr as invert,dn as isclose,yn as iscomplex,ka as iscomplexobj,ja as isdtype,en as isfinite,Ta as isfortran,xo as isin,nn as isinf,on as isnan,tn as isnat,an as isneginf,un as isposinf,Dn as isreal,wa as isrealobj,qa as isscalar,bu as issubdtype,gu as item,za as iterable,Zt as ix_,Le as kron,kr as lcm,wr as ldexp,Jr as left_shift,bn as less,Cn as less_equal,ot as lexsort,ha as linalg,ye as linspace,nr as log,tr as log10,ar as log1p,or as log2,ur as logaddexp,ir as logaddexp2,Xe as logical_and,$e as logical_not,Ze as logical_or,rn as logical_xor,De as logspace,K as mask_indices,We as matmul,Ke as matrix_transpose,Qe as matvec,Pn as max,zr as maximum,na as may_share_memory,Fn as mean,Gn as median,Lt as meshgrid,En as min,su as min_scalar_type,jr as minimum,xu as mintypecode,Nr as mod,Tr as modf,Eo as moveaxis,W as multiply,Ir as nan_to_num,to as nanargmax,oo as nanargmin,uo as nancumprod,ao as nancumsum,no as nanmax,Zn as nanmean,io as nanmedian,eo as nanmin,Do as nanpercentile,Xn as nanprod,yo as nanquantile,ro as nanstd,Wn as nansum,$n as nanvar,lu as ndim,Dr as negative,pn as nextafter,it as nonzero,ln as not_equal,oe as ones,fe as ones_like,Ge as outer,Or as packbits,qo as pad,tt as partition,Hn as percentile,Oe as permute_dims,Vt as place,_n as poly,gn as polyadd,vn as polyder,hn as polydiv,kn as polyfit,wn as polyint,Tn as polymul,qn as polysub,zn as polyval,mr as positive,ua as pow,Z as power,fa as printoptions,In as prod,Sa as promote_types,Mn as ptp,Qt as put,Rt as put_along_axis,Wt as putmask,Jn as quantile,qt as rad2deg,wt as radians,So as ravel,M as ravel_multi_index,Zr as real,An as real_if_close,fr as reciprocal,lr as remainder,Ro as repeat,ge as require,zo as reshape,rt as resize,fu as result_type,Kr as right_shift,fo as rint,$o as roll,Mo as rollaxis,jn as roots,Zo as rot90,No as round,La as round_,Za as row_stack,ru as searchsorted,P as select,pa as set_printoptions,lo as setdiff1d,aa as seterr,_o as setxor1d,du as shape,oa as shares_memory,Ar as sign,cn as signbit,Ct as sin,Mr as sinc,zt as sinh,_u as size,et as sort,ut as sort_complex,mn as spacing,Ja as split,X as sqrt,xr as square,Fo as squeeze,Lo as stack,Yn as std,V as subtract,Sn as sum,Po as swapaxes,q as take,z as take_along_axis,lt as tan,St as tanh,He as tensordot,Qo as tile,hu as tobytes,Tu as tofile,vu as tolist,da as trace,Ue as transpose,bt as trapezoid,ke as tri,we as tril,Y as tril_indices,G as tril_indices_from,vo as trim_zeros,Te as triu,H as triu_indices,J as triu_indices_from,ya as true_divide,so as trunc,Cu as typename,go as union1d,Ua as unique,Ya as unique_all,Ga as unique_counts,Ha as unique_inverse,ho as unique_values,Qr as unpackbits,B as unravel_index,Oo as unstack,Er as unwrap,qe as vander,Ea as var_,Un as variance,va as vdot,Je as vecdot,Re as vecmat,wu as view,Oa as vsplit,Uo as vstack,$a as where,ne as zeros,me as zeros_like};
+import*as o from"../core";import{NDArray as N}from"./ndarray";import{NDArrayCore as i}from"../common/ndarray-core";import{ArrayStorage as s}from"../common/storage";import{Complex as l}from"../common/complex";import{DEFAULT_DTYPE as y,getTypedArrayConstructor as b,promoteDTypes as d,isComplexDType as _}from"../common/dtype";const n=r=>{if(r instanceof N)return r;const e=r.base?n(r.base):void 0;return N.fromStorage(r.storage,e)};import{NDArray as Lt,meshgrid as Ut}from"./ndarray";import{NDArrayCore as Gt}from"../common/ndarray-core";import{Complex as Jt}from"../common/complex";function T(r,e){return n(o.broadcast_to(r,e))}function q(r,e,t){return n(o.take(r,e,t))}function z(r,e,t){return n(o.take_along_axis(r,e,t))}function j(r,e){return n(o.choose(r,e))}function S(r,e,t){return n(o.compress(r,e,t))}function F(r,e,t=0){return n(o.iindex(r,e,t))}function I(r,e,t){return n(o.bindex(r,e,t))}function P(r,e,t=0){return n(o.select(r,e,t))}function E(r,e="int32"){return n(o.indices(r,e))}function M(r,e,t="raise"){return n(o.ravel_multi_index(r,e,t))}function B(r,e){return o.unravel_index(r,e).map(n)}function L(r,e=2){return o.diag_indices(r,e).map(n)}function U(r){return o.diag_indices_from(r).map(n)}function Y(r,e=0,t){return o.tril_indices(r,e,t).map(n)}function G(r,e=0){return o.tril_indices_from(r,e).map(n)}function H(r,e=0,t){return o.triu_indices(r,e,t).map(n)}function J(r,e=0){return o.triu_indices_from(r,e).map(n)}function K(r,e,t=0){return o.mask_indices(r,e,t).map(n)}function O(r,e,t){const a=u=>r(n(u));return n(o.apply_along_axis(a,e,t))}function Q(r,e,t){const a=(u,p)=>r(n(u),p);return n(o.apply_over_axes(a,e,t))}function R(r,e){return n(o.add(r,e))}function V(r,e){return n(o.subtract(r,e))}function W(r,e){return n(o.multiply(r,e))}function X(r){return n(o.sqrt(r))}function Z(r,e){return n(o.power(r,e))}function $(r){return n(o.exp(r))}function rr(r){return n(o.exp2(r))}function er(r){return n(o.expm1(r))}function nr(r){return n(o.log(r))}function or(r){return n(o.log2(r))}function tr(r){return n(o.log10(r))}function ar(r){return n(o.log1p(r))}function ur(r,e){return n(o.logaddexp(r,e))}function ir(r,e){return n(o.logaddexp2(r,e))}function yr(r){return n(o.absolute(r))}function Dr(r){return n(o.negative(r))}function Ar(r){return n(o.sign(r))}function Nr(r,e){return n(o.mod(r,e))}function cr(r,e){return n(o.divide(r,e))}function mr(r,e){return n(o.floor_divide(r,e))}function pr(r){return n(o.positive(r))}function fr(r){return n(o.reciprocal(r))}function sr(r){return n(o.cbrt(r))}function br(r){return n(o.fabs(r))}function Cr(r,e){const t=o.divmod(r,e);return[n(t[0]),n(t[1])]}function xr(r){return n(o.square(r))}function lr(r,e){return n(o.remainder(r,e))}function dr(r,e){return n(o.heaviside(r,e))}function _r(r,e){return n(o.float_power(r,e))}function gr(r,e){return n(o.fmod(r,e))}function vr(r){const e=o.frexp(r);return[n(e[0]),n(e[1])]}function hr(r,e){return n(o.gcd(r,e))}function kr(r,e){return n(o.lcm(r,e))}function wr(r,e){return n(o.ldexp(r,e))}function Tr(r){const e=o.modf(r);return[n(e[0]),n(e[1])]}function qr(r,e,t){return n(o.clip(r,e,t))}function zr(r,e){return n(o.maximum(r,e))}function jr(r,e){return n(o.minimum(r,e))}function Sr(r,e){return n(o.fmax(r,e))}function Fr(r,e){return n(o.fmin(r,e))}function Ir(r,e=0,t,a){return n(o.nan_to_num(r,e,t,a))}function Pr(r,e,t,a,u){return n(o.interp(r,e,t,a,u))}function Er(r,e=Math.PI,t=-1,a=2*Math.PI){return n(o.unwrap(r,e,t,a))}function Mr(r){return n(o.sinc(r))}function Br(r){return n(o.i0(r))}function Lr(r,e){return n(o.bitwise_and(r,e))}function Ur(r,e){return n(o.bitwise_or(r,e))}function Yr(r,e){return n(o.bitwise_xor(r,e))}function Gr(r){return n(o.bitwise_not(r))}function Hr(r){return n(o.invert(r))}function Jr(r,e){return n(o.left_shift(r,e))}function Kr(r,e){return n(o.right_shift(r,e))}function Or(r,e,t){return n(o.packbits(r,e,t))}function Qr(r,e,t,a){return n(o.unpackbits(r,e,t,a))}function Rr(r){return n(o.bitwise_count(r))}function Vr(r){return n(o.bitwise_invert(r))}function Wr(r,e){return n(o.bitwise_left_shift(r,e))}function Xr(r,e){return n(o.bitwise_right_shift(r,e))}function Zr(r){return n(o.real(r))}function $r(r){return n(o.imag(r))}function re(r){return n(o.conj(r))}function ee(r,e){return n(o.angle(r,e))}function ne(r,e=y){return n(o.zeros(r,e))}function oe(r,e=y){return n(o.ones(r,e))}function te(r,e=y){return n(o.empty(r,e))}function ae(r,e,t){return n(o.full(r,e,t))}function ue(r,e){return n(o.array(r,e))}function ie(r,e,t=1,a=y){return n(o.arange(r,e,t,a))}function ye(r,e,t=50,a=y){return n(o.linspace(r,e,t,a))}function De(r,e,t=50,a=10,u=y){return n(o.logspace(r,e,t,a,u))}function Ae(r,e,t=50,a=y){return n(o.geomspace(r,e,t,a))}function Ne(r,e,t=0,a=y){return n(o.eye(r,e,t,a))}function ce(r,e=y){return n(o.identity(r,e))}function me(r,e){return n(o.asarray(r,e))}function pe(r,e){return n(o.zeros_like(r,e))}function fe(r,e){return n(o.ones_like(r,e))}function se(r,e){return n(o.empty_like(r,e))}function be(r,e,t){return n(o.full_like(r,e,t))}function Ce(r){return n(o.copy(r))}function xe(r,e){return n(o.asanyarray(r,e))}function le(r,e){return n(o.ascontiguousarray(r,e))}function de(r,e){return n(o.asfortranarray(r,e))}function _e(r,e){return n(o.asarray_chkfinite(r,e))}function ge(r,e,t){return n(o.require(r,e,t))}function ve(r,e=0){return n(o.diag(r,e))}function he(r,e=0){return n(o.diagflat(r,e))}function ke(r,e,t=0,a=y){return n(o.tri(r,e,t,a))}function we(r,e=0){return n(o.tril(r,e))}function Te(r,e=0){return n(o.triu(r,e))}function qe(r,e,t=!1){return n(o.vander(r,e,t))}function ze(r,e=y,t=-1,a=0){return n(o.frombuffer(r,e,t,a))}function je(r,e,t=y){return n(o.fromfunction(r,e,t))}function Se(r,e=y,t=-1){return n(o.fromiter(r,e,t))}function Fe(r,e=y,t=-1,a){return n(o.fromstring(r,e,t,a))}function Ie(r,e=y,t=-1,a=""){return n(o.fromfile(r,e,t,a))}function Pe(r,e,t){return n(o.diff(r,e,t))}function Ee(r,e,t){return n(o.ediff1d(r,e,t))}function Me(r,e){const t=o.dot(r,e);return t instanceof i?n(t):t}function Be(r,e=0,t=0,a=1){const u=o.trace(r,e,t,a);return u instanceof i?n(u):u}function Le(r,e=0,t=0,a=1){return n(o.diagonal(r,e,t,a))}function Ue(r,e){return n(o.kron(r,e))}function Ye(r,e){return n(o.transpose(r,e))}function Ge(r,e){const t=o.inner(r,e);return t instanceof i?n(t):t}function He(r,e){return n(o.outer(r,e))}function Je(r,e,t=2){const a=o.tensordot(r,e,t);return a instanceof i?n(a):a}function Ke(r,e,t=-1){const a=o.vecdot(r,e,t);return a instanceof i?n(a):a}function Oe(r){return n(o.matrix_transpose(r))}function Qe(r,e){return n(o.permute_dims(r,e))}function Re(r,e){return n(o.matvec(r,e))}function Ve(r,e){return n(o.vecmat(r,e))}function We(r,e,t=-1,a=-1,u=-1,p){const D=o.cross(r,e,t,a,u,p),A=d(r.dtype,e.dtype);if(D instanceof l){const c=b(A==="complex64"?"float32":"float64"),m=new c(2);return m[0]=D.re,m[1]=D.im,N.fromStorage(s.fromData(m,[],A))}if(typeof D=="number"||typeof D=="bigint"){if(_(A)){const x=b(A==="complex64"?"float32":"float64"),f=new x(2);return f[0]=Number(D),f[1]=0,N.fromStorage(s.fromData(f,[],A))}const C=b(A),c=new C(1);return c[0]=D,N.fromStorage(s.fromData(c,[],A))}return n(D)}function Xe(r,e){return n(o.matmul(r,e))}function Ze(r,e){return n(o.logical_and(r,e))}function $e(r,e){return n(o.logical_or(r,e))}function rn(r){return n(o.logical_not(r))}function en(r,e){return n(o.logical_xor(r,e))}function nn(r){return n(o.isfinite(r))}function on(r){return n(o.isinf(r))}function tn(r){return n(o.isnan(r))}function an(r){return n(o.isnat(r))}function un(r){return n(o.isneginf(r))}function yn(r){return n(o.isposinf(r))}function Dn(r){return n(o.iscomplex(r))}function An(r){return n(o.isreal(r))}function Nn(r,e){return n(o.real_if_close(r,e))}function cn(r,e){return n(o.copysign(r,e))}function mn(r){return n(o.signbit(r))}function pn(r,e){return n(o.nextafter(r,e))}function fn(r){return n(o.spacing(r))}function sn(r,e){return n(o.greater(r,e))}function bn(r,e){return n(o.greater_equal(r,e))}function Cn(r,e){return n(o.less(r,e))}function xn(r,e){return n(o.less_equal(r,e))}function ln(r,e){return n(o.equal(r,e))}function dn(r,e){return n(o.not_equal(r,e))}function _n(r,e,t=1e-5,a=1e-8){return n(o.isclose(r,e,t,a))}function gn(r){return n(o.poly(r))}function vn(r,e){return n(o.polyadd(r,e))}function hn(r,e=1){return n(o.polyder(r,e))}function kn(r,e){const t=o.polydiv(r,e);return[n(t[0]),n(t[1])]}function wn(r,e,t){return n(o.polyfit(r,e,t))}function Tn(r,e=1,t=0){return n(o.polyint(r,e,t))}function qn(r,e){return n(o.polymul(r,e))}function zn(r,e){return n(o.polysub(r,e))}function jn(r,e){const t=o.polyval(r,e);return t instanceof i?n(t):t}function Sn(r){return n(o.roots(r))}function Fn(r,e,t){const a=o.sum(r,e,t);return a instanceof i?n(a):a}function In(r,e,t){const a=o.mean(r,e,t);return a instanceof i?n(a):a}function Pn(r,e,t){const a=o.prod(r,e,t);return a instanceof i?n(a):a}function En(r,e,t){const a=o.max(r,e,t);return a instanceof i?n(a):a}function Mn(r,e,t){const a=o.min(r,e,t);return a instanceof i?n(a):a}function Bn(r,e,t){const a=o.ptp(r,e,t);return a instanceof i?n(a):a}function Ln(r,e){const t=o.argmin(r,e);return t instanceof i?n(t):t}function Un(r,e){const t=o.argmax(r,e);return t instanceof i?n(t):t}function Yn(r,e,t,a){const u=o.variance(r,e,t,a);return u instanceof i?n(u):u}function Gn(r,e,t,a){const u=o.std(r,e,t,a);return u instanceof i?n(u):u}function Hn(r,e,t){const a=o.median(r,e,t);return a instanceof i?n(a):a}function Jn(r,e,t,a){const u=o.percentile(r,e,t,a);return u instanceof i?n(u):u}function Kn(r,e,t,a){const u=o.quantile(r,e,t,a);return u instanceof i?n(u):u}function On(r,e,t,a){const u=o.average(r,e,t,a);return u instanceof i?n(u):u}function Qn(r,e,t){const a=o.all(r,e,t);return a instanceof i?n(a):a}function Rn(r,e,t){const a=o.any(r,e,t);return a instanceof i?n(a):a}function Vn(r,e){return n(o.cumsum(r,e))}function Wn(r,e){return n(o.cumprod(r,e))}function Xn(r,e,t){const a=o.nansum(r,e,t);return a instanceof i?n(a):a}function Zn(r,e,t){const a=o.nanprod(r,e,t);return a instanceof i?n(a):a}function $n(r,e,t){const a=o.nanmean(r,e,t);return a instanceof i?n(a):a}function ro(r,e,t,a){const u=o.nanvar(r,e,t,a);return u instanceof i?n(u):u}function eo(r,e,t,a){const u=o.nanstd(r,e,t,a);return u instanceof i?n(u):u}function no(r,e,t){const a=o.nanmin(r,e,t);return a instanceof i?n(a):a}function oo(r,e,t){const a=o.nanmax(r,e,t);return a instanceof i?n(a):a}function to(r,e){const t=o.nanargmin(r,e);return t instanceof i?n(t):t}function ao(r,e){const t=o.nanargmax(r,e);return t instanceof i?n(t):t}function uo(r,e){return n(o.nancumsum(r,e))}function io(r,e){return n(o.nancumprod(r,e))}function yo(r,e,t){const a=o.nanmedian(r,e,t);return a instanceof i?n(a):a}function Do(r,e,t,a){const u=o.nanquantile(r,e,t,a);return u instanceof i?n(u):u}function Ao(r,e,t,a){const u=o.nanpercentile(r,e,t,a);return u instanceof i?n(u):u}function No(r,e=0){return n(o.around(r,e))}function co(r,e=0){return n(o.round(r,e))}function mo(r){return n(o.ceil(r))}function po(r){return n(o.fix(r))}function fo(r){return n(o.floor(r))}function so(r){return n(o.rint(r))}function bo(r){return n(o.trunc(r))}function Co(r,e){return n(o.in1d(r,e))}function xo(r,e){return n(o.intersect1d(r,e))}function lo(r,e){return n(o.isin(r,e))}function _o(r,e){return n(o.setdiff1d(r,e))}function go(r,e){return n(o.setxor1d(r,e))}function vo(r,e){return n(o.union1d(r,e))}function ho(r,e){return n(o.trim_zeros(r,e))}function ko(r){return n(o.unique_values(r))}function wo(r,e,t){return n(o.append(r,e,t))}function To(r,e,t){return n(o.delete(r,e,t))}function qo(r,e,t,a){return n(o.insert(r,e,t,a))}function zo(r,e,t="constant",a=0){return n(o.pad(r,e,t,a))}function jo(r,e){return n(o.reshape(r,e))}function So(r){return n(o.flatten(r))}function Fo(r){return n(o.ravel(r))}function Io(r,e){return n(o.squeeze(r,e))}function Po(r,e){return n(o.expand_dims(r,e))}function Eo(r,e,t){return n(o.swapaxes(r,e,t))}function Mo(r,e,t){return n(o.moveaxis(r,e,t))}function Bo(r,e,t=0){return n(o.rollaxis(r,e,t))}function Lo(r,e=0){return n(o.concatenate(r,e))}function Uo(r,e=0){return n(o.stack(r,e))}function Yo(r){return n(o.vstack(r))}function Go(r){return n(o.hstack(r))}function Ho(r){return n(o.dstack(r))}function Jo(r,e=0){return n(o.concat(r,e))}function Ko(r){return n(o.column_stack(r))}function Oo(r){return n(o.block(r))}function Qo(r,e=0){return o.unstack(r,e).map(n)}function Ro(r,e){return n(o.tile(r,e))}function Vo(r,e,t){return n(o.repeat(r,e,t))}function Wo(r,e){return n(o.flip(r,e))}function Xo(r){return n(o.fliplr(r))}function Zo(r){return n(o.flipud(r))}function $o(r,e=1,t=[0,1]){return n(o.rot90(r,e,t))}function rt(r,e,t){return n(o.roll(r,e,t))}function et(r,e){return n(o.resize(r,e))}function nt(r,e=-1){return n(o.sort(r,e))}function ot(r,e=-1){return n(o.argsort(r,e))}function tt(r){return n(o.lexsort(r))}function at(r,e,t=-1){return n(o.partition(r,e,t))}function ut(r,e,t=-1){return n(o.argpartition(r,e,t))}function it(r){return n(o.sort_complex(r))}function yt(r){return o.nonzero(r).map(n)}function Dt(r){return n(o.argwhere(r))}function At(r){return n(o.flatnonzero(r))}function Nt(r,e){return n(o.extract(r,e))}function ct(r,e,t){return n(o.bincount(r,e,t))}function mt(r,e,t){return n(o.digitize(r,e,t))}function pt(r,e,t){return n(o.correlate(r,e,t))}function ft(r,e,t){return n(o.convolve(r,e,t))}function st(r,e,t,a,u){return n(o.cov(r,e,t,a,u))}function bt(r,e,t){return n(o.corrcoef(r,e,t))}function Ct(r,e,t,a){const u=o.trapezoid(r,e,t,a);return u instanceof i?n(u):u}function xt(r){return n(o.sin(r))}function lt(r){return n(o.cos(r))}function dt(r){return n(o.tan(r))}function _t(r){return n(o.arcsin(r))}function gt(r){return n(o.arccos(r))}function vt(r){return n(o.arctan(r))}function ht(r,e){return n(o.arctan2(r,e))}function kt(r,e){return n(o.hypot(r,e))}function wt(r){return n(o.degrees(r))}function Tt(r){return n(o.radians(r))}function qt(r){return n(o.deg2rad(r))}function zt(r){return n(o.rad2deg(r))}function jt(r){return n(o.sinh(r))}function St(r){return n(o.cosh(r))}function Ft(r){return n(o.tanh(r))}function It(r){return n(o.arcsinh(r))}function Pt(r){return n(o.arccosh(r))}function Et(r){return n(o.arctanh(r))}import{broadcast_arrays as Ot,broadcast_shapes as Qt,put as Rt,put_along_axis as Vt,place as Wt,putmask as Xt,copyto as Zt,ix_ as $t,fill_diagonal as ra,array_equal as ea,array_equiv as na,may_share_memory as oa,shares_memory as ta,geterr as aa,seterr as ua,pow as ia,abs as ya,true_divide as Da,conjugate as Aa,array2string as Na,array_repr as ca,array_str as ma,set_printoptions as pa,get_printoptions as fa,printoptions as sa,format_float_positional as ba,format_float_scientific as Ca,base_repr as xa,binary_repr as la,gradient as da,einsum as _a,einsum_path as ga,vdot as va,linalg as ha,iscomplexobj as ka,isrealobj as wa,isfortran as Ta,isscalar as qa,iterable as za,isdtype as ja,promote_types as Sa,allclose as Fa,amax as Ia,amin as Pa,var_ as Ea,cumulative_sum as Ma,cumulative_prod as Ba,round_ as La,unique as Ua,unique_all as Ya,unique_counts as Ga,unique_inverse as Ha,split as Ja,array_split as Ka,vsplit as Oa,hsplit as Qa,dsplit as Ra,atleast_1d as Va,atleast_2d as Wa,atleast_3d as Xa,row_stack as Za,where as $a,searchsorted as ru,count_nonzero as eu,histogram as nu,histogram2d as ou,histogramdd as tu,histogram_bin_edges as au,asin as uu,acos as iu,atan as yu,atan2 as Du,asinh as Au,acosh as Nu,atanh as cu,can_cast as mu,common_type as pu,result_type as fu,min_scalar_type as su,issubdtype as bu,typename as Cu,mintypecode as xu,ndim as lu,shape as du,size as _u,item as gu,tolist as vu,tobytes as hu,byteswap as ku,view as wu,tofile as Tu,fill as qu}from"../core";export{Jt as Complex,Lt as NDArray,Gt as NDArrayCore,ya as abs,yr as absolute,iu as acos,Nu as acosh,R as add,Qn as all,Fa as allclose,Ia as amax,Pa as amin,ee as angle,Rn as any,wo as append,O as apply_along_axis,Q as apply_over_axes,ie as arange,gt as arccos,Pt as arccosh,_t as arcsin,It as arcsinh,vt as arctan,ht as arctan2,Et as arctanh,Un as argmax,Ln as argmin,ut as argpartition,ot as argsort,Dt as argwhere,No as around,ue as array,Na as array2string,ea as array_equal,na as array_equiv,ca as array_repr,Ka as array_split,ma as array_str,xe as asanyarray,me as asarray,_e as asarray_chkfinite,le as ascontiguousarray,de as asfortranarray,uu as asin,Au as asinh,yu as atan,Du as atan2,cu as atanh,Va as atleast_1d,Wa as atleast_2d,Xa as atleast_3d,On as average,xa as base_repr,la as binary_repr,ct as bincount,I as bindex,Lr as bitwise_and,Rr as bitwise_count,Vr as bitwise_invert,Wr as bitwise_left_shift,Gr as bitwise_not,Ur as bitwise_or,Xr as bitwise_right_shift,Yr as bitwise_xor,Oo as block,Ot as broadcast_arrays,Qt as broadcast_shapes,T as broadcast_to,ku as byteswap,mu as can_cast,sr as cbrt,mo as ceil,j as choose,qr as clip,Ko as column_stack,pu as common_type,S as compress,Jo as concat,Lo as concatenate,re as conj,Aa as conjugate,ft as convolve,Ce as copy,cn as copysign,Zt as copyto,bt as corrcoef,pt as correlate,lt as cos,St as cosh,eu as count_nonzero,st as cov,We as cross,Wn as cumprod,Vn as cumsum,Ba as cumulative_prod,Ma as cumulative_sum,qt as deg2rad,wt as degrees,To as delete_,ve as diag,L as diag_indices,U as diag_indices_from,he as diagflat,Le as diagonal,Pe as diff,mt as digitize,cr as divide,Cr as divmod,Me as dot,Ra as dsplit,Ho as dstack,Ee as ediff1d,_a as einsum,ga as einsum_path,te as empty,se as empty_like,ln as equal,$ as exp,rr as exp2,Po as expand_dims,er as expm1,Nt as extract,Ne as eye,br as fabs,qu as fill,ra as fill_diagonal,po as fix,At as flatnonzero,So as flatten,Wo as flip,Xo as fliplr,Zo as flipud,_r as float_power,fo as floor,mr as floor_divide,Sr as fmax,Fr as fmin,gr as fmod,ba as format_float_positional,Ca as format_float_scientific,vr as frexp,ze as frombuffer,Ie as fromfile,je as fromfunction,Se as fromiter,Fe as fromstring,ae as full,be as full_like,hr as gcd,Ae as geomspace,fa as get_printoptions,aa as geterr,da as gradient,sn as greater,bn as greater_equal,dr as heaviside,nu as histogram,ou as histogram2d,au as histogram_bin_edges,tu as histogramdd,Qa as hsplit,Go as hstack,kt as hypot,Br as i0,ce as identity,F as iindex,$r as imag,Co as in1d,E as indices,Ge as inner,qo as insert,Pr as interp,xo as intersect1d,Hr as invert,_n as isclose,Dn as iscomplex,ka as iscomplexobj,ja as isdtype,nn as isfinite,Ta as isfortran,lo as isin,on as isinf,tn as isnan,an as isnat,un as isneginf,yn as isposinf,An as isreal,wa as isrealobj,qa as isscalar,bu as issubdtype,gu as item,za as iterable,$t as ix_,Ue as kron,kr as lcm,wr as ldexp,Jr as left_shift,Cn as less,xn as less_equal,tt as lexsort,ha as linalg,ye as linspace,nr as log,tr as log10,ar as log1p,or as log2,ur as logaddexp,ir as logaddexp2,Ze as logical_and,rn as logical_not,$e as logical_or,en as logical_xor,De as logspace,K as mask_indices,Xe as matmul,Oe as matrix_transpose,Re as matvec,En as max,zr as maximum,oa as may_share_memory,In as mean,Hn as median,Ut as meshgrid,Mn as min,su as min_scalar_type,jr as minimum,xu as mintypecode,Nr as mod,Tr as modf,Mo as moveaxis,W as multiply,Ir as nan_to_num,ao as nanargmax,to as nanargmin,io as nancumprod,uo as nancumsum,oo as nanmax,$n as nanmean,yo as nanmedian,no as nanmin,Ao as nanpercentile,Zn as nanprod,Do as nanquantile,eo as nanstd,Xn as nansum,ro as nanvar,lu as ndim,Dr as negative,pn as nextafter,yt as nonzero,dn as not_equal,oe as ones,fe as ones_like,He as outer,Or as packbits,zo as pad,at as partition,Jn as percentile,Qe as permute_dims,Wt as place,gn as poly,vn as polyadd,hn as polyder,kn as polydiv,wn as polyfit,Tn as polyint,qn as polymul,zn as polysub,jn as polyval,pr as positive,ia as pow,Z as power,sa as printoptions,Pn as prod,Sa as promote_types,Bn as ptp,Rt as put,Vt as put_along_axis,Xt as putmask,Kn as quantile,zt as rad2deg,Tt as radians,Fo as ravel,M as ravel_multi_index,Zr as real,Nn as real_if_close,fr as reciprocal,lr as remainder,Vo as repeat,ge as require,jo as reshape,et as resize,fu as result_type,Kr as right_shift,so as rint,rt as roll,Bo as rollaxis,Sn as roots,$o as rot90,co as round,La as round_,Za as row_stack,ru as searchsorted,P as select,pa as set_printoptions,_o as setdiff1d,ua as seterr,go as setxor1d,du as shape,ta as shares_memory,Ar as sign,mn as signbit,xt as sin,Mr as sinc,jt as sinh,_u as size,nt as sort,it as sort_complex,fn as spacing,Ja as split,X as sqrt,xr as square,Io as squeeze,Uo as stack,Gn as std,V as subtract,Fn as sum,Eo as swapaxes,q as take,z as take_along_axis,dt as tan,Ft as tanh,Je as tensordot,Ro as tile,hu as tobytes,Tu as tofile,vu as tolist,Be as trace,Ye as transpose,Ct as trapezoid,ke as tri,we as tril,Y as tril_indices,G as tril_indices_from,ho as trim_zeros,Te as triu,H as triu_indices,J as triu_indices_from,Da as true_divide,bo as trunc,Cu as typename,vo as union1d,Ua as unique,Ya as unique_all,Ga as unique_counts,Ha as unique_inverse,ko as unique_values,Qr as unpackbits,B as unravel_index,Qo as unstack,Er as unwrap,qe as vander,Ea as var_,Yn as variance,va as vdot,Ke as vecdot,Ve as vecmat,wu as view,Oa as vsplit,Yo as vstack,$a as where,ne as zeros,pe as zeros_like};