numpy-ts 1.1.0 → 1.3.0

package/dist/esm/core/index.js

@@ -1 +1 @@
-import{NDArrayCore as a,Complex as t,ArrayStorage as i}from"./types";import{zeros as s,ones as n,empty as p,full as l,array as m,arange as d,linspace as c,logspace as _,geomspace as f,eye as u,identity as y,asarray as x,zeros_like as g,ones_like as v,empty_like as h,full_like as b,copy as z,asanyarray as k,ascontiguousarray as N,asfortranarray as w,asarray_chkfinite as q,require as T,diag as S,diagflat as D,tri as P,tril as E,triu as O,vander as j,frombuffer as A,fromfunction as C,fromiter as R,fromstring as H,fromfile as I,meshgrid as U}from"./creation";import{add as M,subtract as V,multiply as B,sqrt as F,power as G,pow as J,exp as K,exp2 as L,expm1 as Q,log as W,log2 as X,log10 as Z,log1p as $,logaddexp as ee,logaddexp2 as re,absolute as ae,abs as te,negative as ie,sign as oe,mod as se,divide as ne,true_divide as pe,floor_divide as le,positive as me,reciprocal as de,cbrt as ce,fabs as _e,divmod as fe,square as ue,remainder as ye,heaviside as xe,float_power as ge,fmod as ve,frexp as he,gcd as be,lcm as ze,ldexp as ke,modf as Ne,clip as we,maximum as qe,minimum as Te,fmax as Se,fmin as De,nan_to_num as Pe,interp as Ee,unwrap as Oe,sinc as je,i0 as Ae}from"./arithmetic";import{sin as Re,cos as He,tan as Ie,arcsin as Ue,asin as Ye,arccos as Me,acos as Ve,arctan as Be,atan as Fe,arctan2 as Ge,atan2 as Je,hypot as Ke,degrees as Le,radians as Qe,deg2rad as We,rad2deg as Xe,sinh as Ze,cosh as $e,tanh as er,arcsinh as rr,asinh as ar,arccosh as tr,acosh as ir,arctanh as or,atanh as sr}from"./trig";import{dot as pr,matmul as lr,trace as mr,diagonal as dr,kron as cr,transpose as _r,inner as fr,outer as ur,tensordot as yr,einsum as xr,einsum_path as gr,vdot as vr,vecdot as hr,matrix_transpose as br,permute_dims as zr,matvec as kr,vecmat as Nr,cross as wr,linalg as qr}from"./linalg";import{reshape as Sr,flatten as Dr,ravel as Pr,squeeze as Er,expand_dims as Or,swapaxes as jr,moveaxis as Ar,rollaxis as Cr,concatenate as Rr,stack as Hr,vstack as Ir,hstack as Ur,dstack as Yr,concat as Mr,column_stack as Vr,row_stack as Br,block as Fr,split as Gr,array_split as Jr,vsplit as Kr,hsplit as Lr,dsplit as Qr,unstack as Wr,tile as Xr,repeat as Zr,flip as $r,fliplr as ea,flipud as ra,rot90 as aa,roll as ta,resize as ia,atleast_1d as oa,atleast_2d as sa,atleast_3d as na}from"./shape";import{sum as la,mean as ma,prod as da,max as ca,amax as _a,min as fa,amin as ua,ptp as ya,argmin as xa,argmax as ga,variance as va,var_ as ha,variance as ba,std as za,median as ka,percentile as Na,quantile as wa,average as qa,all as Ta,any as Sa,cumsum as Da,cumulative_sum as Pa,cumprod as Ea,cumulative_prod as Oa,nansum as ja,nanprod as Aa,nanmean as Ca,nanvar as Ra,nanstd as Ha,nanmin as Ia,nanmax as Ua,nanargmin as Ya,nanargmax as Ma,nancumsum as Va,nancumprod as Ba,nanmedian as Fa,nanquantile as Ga,nanpercentile as Ja}from"./reduction";import{logical_and as La,logical_or as Qa,logical_not as Wa,logical_xor as Xa,isfinite as Za,isinf as $a,isnan as et,isnat as rt,isneginf as at,isposinf as tt,iscomplex as it,iscomplexobj as ot,isreal as st,isrealobj as nt,real_if_close as pt,isfortran as lt,isscalar as mt,iterable as dt,isdtype as ct,promote_types as _t,copysign as ft,signbit as ut,nextafter as yt,spacing as xt,greater as gt,greater_equal as vt,less as ht,less_equal as bt,equal as zt,not_equal as kt,isclose as Nt,allclose as wt}from"./logic";import{sort as Tt,argsort as St,lexsort as Dt,partition as Pt,argpartition as Et,sort_complex as Ot,nonzero as jt,argwhere as At,flatnonzero as Ct,where as Rt,searchsorted as Ht,extract as It,count_nonzero as Ut}from"./sorting";import{bitwise_and as Mt,bitwise_or as Vt,bitwise_xor as Bt,bitwise_not as Ft,invert as Gt,left_shift as Jt,right_shift as Kt,packbits as Lt,unpackbits as Qt,bitwise_count as Wt,bitwise_invert as Xt,bitwise_left_shift as Zt,bitwise_right_shift as $t}from"./bitwise";import{around as ri,round_ as ai,round as ti,ceil as ii,fix as oi,floor as si,rint as ni,trunc as pi}from"./rounding";import{unique as mi,in1d as di,intersect1d as ci,isin as _i,setdiff1d as fi,setxor1d as ui,union1d as yi,trim_zeros as xi,unique_all as gi,unique_counts as vi,unique_inverse as hi,unique_values as bi}from"./sets";import{bincount as ki,digitize as Ni,histogram as wi,histogram2d as qi,histogramdd as Ti,correlate as Si,convolve as Di,cov as Pi,corrcoef as Ei,histogram_bin_edges as Oi,trapezoid as ji}from"./statistics";import{diff as Ci,ediff1d as Ri,gradient as Hi}from"./gradient";import{real as Ui,imag as Yi,conj as Mi,conjugate as Vi,angle as Bi}from"./complex";import{broadcast_to as Gi,broadcast_arrays as Ji,broadcast_shapes as Ki,take as Li,put as Qi,take_along_axis as Wi,put_along_axis as Xi,choose as Zi,compress as $i,select as eo,place as ro,putmask as ao,copyto as to,indices as io,ix_ as oo,ravel_multi_index as so,unravel_index as no,diag_indices as po,diag_indices_from as lo,fill_diagonal as mo,tril_indices as co,tril_indices_from as _o,triu_indices as fo,triu_indices_from as uo,mask_indices as yo,array_equal as xo,array_equiv as go,apply_along_axis as vo,apply_over_axes as ho,may_share_memory as bo,shares_memory as zo,geterr as ko,seterr as No,iindex as wo,bindex as qo}from"./advanced";import{set_printoptions as So,get_printoptions as Do,printoptions as Po,format_float_positional as Eo,format_float_scientific as Oo,base_repr as jo,binary_repr as Ao,array2string as Co,array_repr as Ro,array_str as Ho}from"./formatting";import{ndim as Uo,shape as Yo,size as Mo,item as Vo,tolist as Bo,tobytes as Fo,byteswap as Go,view as Jo,tofile as Ko,fill as Lo}from"./utility";import{append as Wo,delete_ as Xo,delete_ as Zo,insert as $o,pad as es}from"./shape-extra";import{can_cast as as,common_type as ts,result_type as is,min_scalar_type as os,issubdtype as ss,typename as ns,mintypecode as ps}from"./typechecking";import{poly as ms,polyadd as ds,polyder as cs,polydiv as _s,polyfit as fs,polyint as us,polymul as ys,polysub as xs,polyval as gs,roots as vs}from"./polynomial";import{parseNpy as bs,parseNpyHeader as zs,parseNpyData as ks}from"../io/npy/parser";import{serializeNpy as ws}from"../io/npy/serializer";import{UnsupportedDTypeError as Ts,InvalidNpyError as Ss,SUPPORTED_DTYPES as Ds,DTYPE_TO_DESCR as Ps}from"../io/npy/format";import{parseNpz as Os,parseNpzSync as js,loadNpz as As,loadNpzSync as Cs}from"../io/npz/parser";import{serializeNpz as Hs,serializeNpzSync as Is}from"../io/npz/serializer";import{parseTxt as Ys,genfromtxt as Ms,fromregex as Vs}from"../io/txt/parser";import{serializeTxt as Fs}from"../io/txt/serializer";export{i as ArrayStorage,t as Complex,Ps as DTYPE_TO_DESCR,Ss as InvalidNpyError,a as NDArrayCore,Ds as SUPPORTED_DTYPES,Ts as UnsupportedDTypeError,te as abs,ae as absolute,Ve as acos,ir as acosh,M as add,Ta as all,wt as allclose,_a as amax,ua as amin,Bi as angle,Sa as any,Wo as append,vo as apply_along_axis,ho as apply_over_axes,d as arange,Me as arccos,tr as arccosh,Ue as arcsin,rr as arcsinh,Be as arctan,Ge as arctan2,or as arctanh,ga as argmax,xa as argmin,Et as argpartition,St as argsort,At as argwhere,ri as around,m as array,Co as array2string,xo as array_equal,go as array_equiv,Ro as array_repr,Jr as array_split,Ho as array_str,k as asanyarray,x as asarray,q as asarray_chkfinite,N as ascontiguousarray,w as asfortranarray,Ye as asin,ar as asinh,Fe as atan,Je as atan2,sr as atanh,oa as atleast_1d,sa as atleast_2d,na as atleast_3d,qa as average,jo as base_repr,Ao as binary_repr,ki as bincount,qo as bindex,Mt as bitwise_and,Wt as bitwise_count,Xt as bitwise_invert,Zt as bitwise_left_shift,Ft as bitwise_not,Vt as bitwise_or,$t as bitwise_right_shift,Bt as bitwise_xor,Fr as block,Ji as broadcast_arrays,Ki as broadcast_shapes,Gi as broadcast_to,Go as byteswap,as as can_cast,ce as cbrt,ii as ceil,Zi as choose,we as clip,Vr as column_stack,ts as common_type,$i as compress,Mr as concat,Rr as concatenate,Mi as conj,Vi as conjugate,Di as convolve,z as copy,ft as copysign,to as copyto,Ei as corrcoef,Si as correlate,He as cos,$e as cosh,Ut as count_nonzero,Pi as cov,wr as cross,Ea as cumprod,Da as cumsum,Oa as cumulative_prod,Pa as cumulative_sum,We as deg2rad,Le as degrees,Zo as delete,Xo as delete_,S as diag,po as diag_indices,lo as diag_indices_from,D as diagflat,dr as diagonal,Ci as diff,Ni as digitize,ne as divide,fe as divmod,pr as dot,Qr as dsplit,Yr as dstack,Ri as ediff1d,xr as einsum,gr as einsum_path,p as empty,h as empty_like,zt as equal,K as exp,L as exp2,Or as expand_dims,Q as expm1,It as extract,u as eye,_e as fabs,Lo as fill,mo as fill_diagonal,oi as fix,Ct as flatnonzero,Dr as flatten,$r as flip,ea as fliplr,ra as flipud,ge as float_power,si as floor,le as floor_divide,Se as fmax,De as fmin,ve as fmod,Eo as format_float_positional,Oo as format_float_scientific,he as frexp,A as frombuffer,I as fromfile,C as fromfunction,R as fromiter,Vs as fromregex,H as fromstring,l as full,b as full_like,be as gcd,Ms as genfromtxt,f as geomspace,Do as get_printoptions,ko as geterr,Hi as gradient,gt as greater,vt as greater_equal,xe as heaviside,wi as histogram,qi as histogram2d,Oi as histogram_bin_edges,Ti as histogramdd,Lr as hsplit,Ur as hstack,Ke as hypot,Ae as i0,y as identity,wo as iindex,Yi as imag,di as in1d,io as indices,fr as inner,$o as insert,Ee as interp,ci as intersect1d,Gt as invert,Nt as isclose,it as iscomplex,ot as iscomplexobj,ct as isdtype,Za as isfinite,lt as isfortran,_i as isin,$a as isinf,et as isnan,rt as isnat,at as isneginf,tt as isposinf,st as isreal,nt as isrealobj,mt as isscalar,ss as issubdtype,Vo as item,dt as iterable,oo as ix_,cr as kron,ze as lcm,ke as ldexp,Jt as left_shift,ht as less,bt as less_equal,Dt as lexsort,qr as linalg,c as linspace,As as loadNpz,Cs as loadNpzSync,W as log,Z as log10,$ as log1p,X as log2,ee as logaddexp,re as logaddexp2,La as logical_and,Wa as logical_not,Qa as logical_or,Xa as logical_xor,_ as logspace,yo as mask_indices,lr as matmul,br as matrix_transpose,kr as matvec,ca as max,qe as maximum,bo as may_share_memory,ma as mean,ka as median,U as meshgrid,fa as min,os as min_scalar_type,Te as minimum,ps as mintypecode,se as mod,Ne as modf,Ar as moveaxis,B as multiply,Pe as nan_to_num,Ma as nanargmax,Ya as nanargmin,Ba as nancumprod,Va as nancumsum,Ua as nanmax,Ca as nanmean,Fa as nanmedian,Ia as nanmin,Ja as nanpercentile,Aa as nanprod,Ga as nanquantile,Ha as nanstd,ja as nansum,Ra as nanvar,Uo as ndim,ie as negative,yt as nextafter,jt as nonzero,kt as not_equal,n as ones,v as ones_like,ur as outer,Lt as packbits,es as pad,bs as parseNpy,ks as parseNpyData,zs as parseNpyHeader,Os as parseNpz,js as parseNpzSync,Ys as parseTxt,Pt as partition,Na as percentile,zr as permute_dims,ro as place,ms as poly,ds as polyadd,cs as polyder,_s as polydiv,fs as polyfit,us as polyint,ys as polymul,xs as polysub,gs as polyval,me as positive,J as pow,G as power,Po as printoptions,da as prod,_t as promote_types,ya as ptp,Qi as put,Xi as put_along_axis,ao as putmask,wa as quantile,Xe as rad2deg,Qe as radians,Pr as ravel,so as ravel_multi_index,Ui as real,pt as real_if_close,de as reciprocal,ye as remainder,Zr as repeat,T as require,Sr as reshape,ia as resize,is as result_type,Kt as right_shift,ni as rint,ta as roll,Cr as rollaxis,vs as roots,aa as rot90,ti as round,ai as round_,Br as row_stack,Ht as searchsorted,eo as select,ws as serializeNpy,Hs as serializeNpz,Is as serializeNpzSync,Fs as serializeTxt,So as set_printoptions,fi as setdiff1d,No as seterr,ui as setxor1d,Yo as shape,zo as shares_memory,oe as sign,ut as signbit,Re as sin,je as sinc,Ze as sinh,Mo as size,Tt as sort,Ot as sort_complex,xt as spacing,Gr as split,F as sqrt,ue as square,Er as squeeze,Hr as stack,za as std,V as subtract,la as sum,jr as swapaxes,Li as take,Wi as take_along_axis,Ie as tan,er as tanh,yr as tensordot,Xr as tile,Fo as tobytes,Ko as tofile,Bo as tolist,mr as trace,_r as transpose,ji as trapezoid,P as tri,E as tril,co as tril_indices,_o as tril_indices_from,xi as trim_zeros,O as triu,fo as triu_indices,uo as triu_indices_from,pe as true_divide,pi as trunc,ns as typename,yi as union1d,mi as unique,gi as unique_all,vi as unique_counts,hi as unique_inverse,bi as unique_values,Qt as unpackbits,no as unravel_index,Wr as unstack,Oe as unwrap,j as vander,ba as var,ha as var_,va as variance,vr as vdot,hr as vecdot,Nr as vecmat,Jo as view,Kr as vsplit,Ir as vstack,Rt as where,s as zeros,g as zeros_like};
+import{NDArrayCore as w,Complex as A,ArrayStorage as v}from"./types.js";import{zeros as k,ones as z,empty as N,full as q,array as I,arange as D,linspace as T,logspace as C,geomspace as O,eye as P,identity as E,asarray as L,zeros_like as j,ones_like as R,empty_like as G,full_like as H,copy as U,asanyarray as Y,ascontiguousarray as M,asfortranarray as V,asarray_chkfinite as B,require as F,diag as J,diagflat as K,tri as Q,tril as W,triu as X,vander as Z,frombuffer as $,fromfunction as rr,fromiter as er,fromstring as nr,fromfile as ar,meshgrid as tr}from"./creation.js";import{add as mr,subtract as ir,multiply as ur,sqrt as sr,power as br,pow as fr,exp as pr,exp2 as lr,expm1 as dr,log as cr,log2 as _r,log10 as gr,log1p as yr,logaddexp as xr,logaddexp2 as hr,absolute as wr,abs as Ar,negative as vr,sign as Sr,mod as kr,divide as zr,true_divide as Nr,floor_divide as qr,positive as Ir,reciprocal as Dr,cbrt as Tr,fabs as Cr,divmod as Or,square as Pr,remainder as Er,heaviside as Lr,float_power as jr,fmod as Rr,frexp as Gr,gcd as Hr,lcm as Ur,ldexp as Yr,modf as Mr,clip as Vr,maximum as Br,minimum as Fr,fmax as Jr,fmin as Kr,nan_to_num as Qr,interp as Wr,unwrap as Xr,sinc as Zr,i0 as $r}from"./arithmetic.js";import{sin as ee,cos as ne,tan as ae,arcsin as te,asin as oe,arccos as me,acos as ie,arctan as ue,atan as se,arctan2 as be,atan2 as fe,hypot as pe,degrees as le,radians as de,deg2rad as ce,rad2deg as _e,sinh as ge,cosh as ye,tanh as xe,arcsinh as he,asinh as we,arccosh as Ae,acosh as ve,arctanh as Se,atanh as ke}from"./trig.js";import{dot as Ne,matmul as qe,trace as Ie,diagonal as De,kron as Te,transpose as Ce,inner as Oe,outer as Pe,tensordot as Ee,einsum as Le,einsum_path as je,vdot as Re,vecdot as Ge,matrix_transpose as He,permute_dims as Ue,matvec as Ye,vecmat as Me,cross as Ve,linalg as Be}from"./linalg.js";import{reshape as Je,flatten as Ke,ravel as Qe,squeeze as We,expand_dims as Xe,swapaxes as Ze,moveaxis as $e,rollaxis as rn,concatenate as en,stack as nn,vstack as an,hstack as tn,dstack as on,concat as mn,column_stack as un,row_stack as sn,block as bn,split as fn,array_split as pn,vsplit as ln,hsplit as dn,dsplit as cn,unstack as _n,tile as gn,repeat as yn,flip as xn,fliplr as hn,flipud as wn,rot90 as An,roll as vn,resize as Sn,atleast_1d as kn,atleast_2d as zn,atleast_3d as Nn}from"./shape.js";import{sum as In,mean as Dn,prod as Tn,max as Cn,amax as On,min as Pn,amin as En,ptp as Ln,argmin as jn,argmax as Rn,variance as Gn,var_ as Hn,variance as Un,std as Yn,median as Mn,percentile as Vn,quantile as Bn,average as Fn,all as Jn,any as Kn,cumsum as Qn,cumulative_sum as Wn,cumprod as Xn,cumulative_prod as Zn,nansum as $n,nanprod as ra,nanmean as ea,nanvar as na,nanstd as aa,nanmin as ta,nanmax as oa,nanargmin as ma,nanargmax as ia,nancumsum as ua,nancumprod as sa,nanmedian as ba,nanquantile as fa,nanpercentile as pa}from"./reduction.js";import{logical_and as da,logical_or as ca,logical_not as _a,logical_xor as ga,isfinite as ya,isinf as xa,isnan as ha,isnat as wa,isneginf as Aa,isposinf as va,iscomplex as Sa,iscomplexobj as ka,isreal as za,isrealobj as Na,real_if_close as qa,isfortran as Ia,isscalar as Da,iterable as Ta,isdtype as Ca,promote_types as Oa,copysign as Pa,signbit as Ea,nextafter as La,spacing as ja,greater as Ra,greater_equal as Ga,less as Ha,less_equal as Ua,equal as Ya,not_equal as Ma,isclose as Va,allclose as Ba}from"./logic.js";import{sort as Ja,argsort as Ka,lexsort as Qa,partition as Wa,argpartition as Xa,sort_complex as Za,nonzero as $a,argwhere as rt,flatnonzero as et,where as nt,searchsorted as at,extract as tt,count_nonzero as ot}from"./sorting.js";import{bitwise_and as it,bitwise_or as ut,bitwise_xor as st,bitwise_not as bt,invert as ft,left_shift as pt,right_shift as lt,packbits as dt,unpackbits as ct,bitwise_count as _t,bitwise_invert as gt,bitwise_left_shift as yt,bitwise_right_shift as xt}from"./bitwise.js";import{around as wt,round as At,ceil as vt,fix as St,floor as kt,rint as zt,trunc as Nt}from"./rounding.js";import{unique as It,in1d as Dt,intersect1d as Tt,isin as Ct,setdiff1d as Ot,setxor1d as Pt,union1d as Et,trim_zeros as Lt,unique_all as jt,unique_counts as Rt,unique_inverse as Gt,unique_values as Ht}from"./sets.js";import{bincount as Yt,digitize as Mt,histogram as Vt,histogram2d as Bt,histogramdd as Ft,correlate as Jt,convolve as Kt,cov as Qt,corrcoef as Wt,histogram_bin_edges as Xt,trapezoid as Zt}from"./statistics.js";import{diff as ro,ediff1d as eo,gradient as no}from"./gradient.js";import{real as to,imag as oo,conj as mo,conjugate as io,angle as uo}from"./complex.js";import{broadcast_to as bo,broadcast_arrays as fo,broadcast_shapes as po,take as lo,put as co,take_along_axis as _o,put_along_axis as go,choose as yo,compress as xo,select as ho,place as wo,putmask as Ao,copyto as vo,indices as So,ix_ as ko,ravel_multi_index as zo,unravel_index as No,diag_indices as qo,diag_indices_from as Io,fill_diagonal as Do,tril_indices as To,tril_indices_from as Co,triu_indices as Oo,triu_indices_from as Po,mask_indices as Eo,array_equal as Lo,array_equiv as jo,apply_along_axis as Ro,apply_over_axes as Go,may_share_memory as Ho,shares_memory as Uo,geterr as Yo,seterr as Mo,iindex as Vo,bindex as Bo,vindex as Fo}from"./advanced.js";import{set_printoptions as Ko,get_printoptions as Qo,printoptions as Wo,format_float_positional as Xo,format_float_scientific as Zo,base_repr as $o,binary_repr as rm,array2string as em,array_repr as nm,array_str as am}from"./formatting.js";import{ndim as om,shape as mm,size as im,item as um,tolist as sm,tobytes as bm,byteswap as fm,view as pm,tofile as lm,fill as dm}from"./utility.js";import{append as _m,delete_ as gm,delete_ as ym,insert as xm,pad as hm}from"./shape-extra.js";import{can_cast as Am,common_type as vm,result_type as Sm,min_scalar_type as km,issubdtype as zm,typename as Nm,mintypecode as qm}from"./typechecking.js";import{poly as Dm,polyadd as Tm,polyder as Cm,polydiv as Om,polyfit as Pm,polyint as Em,polymul as Lm,polysub as jm,polyval as Rm,roots as Gm}from"./polynomial.js";import{parseNpy as Um,parseNpyHeader as Ym,parseNpyData as Mm}from"../io/npy/parser.js";import{serializeNpy as Bm}from"../io/npy/serializer.js";import{UnsupportedDTypeError as Jm,InvalidNpyError as Km,SUPPORTED_DTYPES as Qm,DTYPE_TO_DESCR as Wm}from"../io/npy/format.js";import{parseNpz as Zm,parseNpzSync as $m,loadNpz as ri,loadNpzSync as ei}from"../io/npz/parser.js";import{serializeNpz as ai,serializeNpzSync as ti}from"../io/npz/serializer.js";import{parseTxt as mi,genfromtxt as ii,fromregex as ui}from"../io/txt/parser.js";import{serializeTxt as bi}from"../io/txt/serializer.js";import*as o from"../common/ops/random.js";import{NDArrayCore as m}from"../common/ndarray-core.js";function t(r){return r&&typeof r=="object"&&"_data"in r&&"_shape"in r?m.fromStorage(r):r}const g={seed:o.seed,get_state:o.get_state,set_state:o.set_state,get_bit_generator:o.get_bit_generator,set_bit_generator:o.set_bit_generator,default_rng:r=>{const e=o.default_rng(r),n=a=>a instanceof m?a.storage:a;return{random:a=>t(e.random(a)),integers:(a,s,b)=>t(e.integers(a,s,b)),standard_normal:a=>t(e.standard_normal(a)),normal:(a,s,b)=>t(e.normal(a,s,b)),uniform:(a,s,b)=>t(e.uniform(a,s,b)),choice:(a,s,b,f)=>t(e.choice(n(a),s,b,f instanceof m?f.storage:f)),permutation:a=>t(e.permutation(n(a))),shuffle:a=>e.shuffle(a instanceof m?a.storage:a),exponential:(a,s)=>t(e.exponential(a,s)),poisson:(a,s)=>t(e.poisson(a,s)),binomial:(a,s,b)=>t(e.binomial(a,s,b))}},Generator:o.Generator,random:r=>t(o.random(r)),rand:(...r)=>t(o.rand(...r)),randn:(...r)=>t(o.randn(...r)),randint:(r,e,n,a)=>t(o.randint(r,e,n,a)),random_sample:r=>t(o.random_sample(r)),ranf:r=>t(o.ranf(r)),sample:r=>t(o.sample(r)),random_integers:(r,e,n)=>t(o.random_integers(r,e,n)),bytes:o.bytes,uniform:(r,e,n)=>t(o.uniform(r,e,n)),normal:(r,e,n)=>t(o.normal(r,e,n)),standard_normal:r=>t(o.standard_normal(r)),exponential:(r,e)=>t(o.exponential(r,e)),standard_exponential:r=>t(o.standard_exponential(r)),gamma:(r,e,n)=>t(o.gamma(r,e,n)),standard_gamma:(r,e)=>t(o.standard_gamma(r,e)),beta:(r,e,n)=>t(o.beta(r,e,n)),chisquare:(r,e)=>t(o.chisquare(r,e)),noncentral_chisquare:(r,e,n)=>t(o.noncentral_chisquare(r,e,n)),f:(r,e,n)=>t(o.f(r,e,n)),noncentral_f:(r,e,n,a)=>t(o.noncentral_f(r,e,n,a)),standard_cauchy:r=>t(o.standard_cauchy(r)),standard_t:(r,e)=>t(o.standard_t(r,e)),laplace:(r,e,n)=>t(o.laplace(r,e,n)),logistic:(r,e,n)=>t(o.logistic(r,e,n)),lognormal:(r,e,n)=>t(o.lognormal(r,e,n)),gumbel:(r,e,n)=>t(o.gumbel(r,e,n)),pareto:(r,e)=>t(o.pareto(r,e)),power:(r,e)=>t(o.power(r,e)),rayleigh:(r,e)=>t(o.rayleigh(r,e)),triangular:(r,e,n,a)=>t(o.triangular(r,e,n,a)),wald:(r,e,n)=>t(o.wald(r,e,n)),weibull:(r,e)=>t(o.weibull(r,e)),poisson:(r,e)=>t(o.poisson(r,e)),binomial:(r,e,n)=>t(o.binomial(r,e,n)),geometric:(r,e)=>t(o.geometric(r,e)),hypergeometric:(r,e,n,a)=>t(o.hypergeometric(r,e,n,a)),logseries:(r,e)=>t(o.logseries(r,e)),negative_binomial:(r,e,n)=>t(o.negative_binomial(r,e,n)),zipf:(r,e)=>t(o.zipf(r,e)),multinomial:(r,e,n)=>t(o.multinomial(r,e,n)),multivariate_normal:(r,e,n,a,s)=>t(o.multivariate_normal(r,e,n,a,s)),dirichlet:(r,e)=>t(o.dirichlet(r,e)),vonmises:(r,e,n)=>t(o.vonmises(r,e,n)),choice:(r,e,n,a)=>t(o.choice(r,e,n,a)),permutation:r=>t(o.permutation(r)),shuffle:r=>o.shuffle(r instanceof m?r.storage:r)};import*as i from"../common/ops/fft.js";function u(r){return r instanceof m?r.storage:r}const y={fft:(r,e,n,a)=>m.fromStorage(i.fft(u(r),e,n,a)),ifft:(r,e,n,a)=>m.fromStorage(i.ifft(u(r),e,n,a)),fft2:(r,e,n,a)=>m.fromStorage(i.fft2(u(r),e,n,a)),ifft2:(r,e,n,a)=>m.fromStorage(i.ifft2(u(r),e,n,a)),fftn:(r,e,n,a)=>m.fromStorage(i.fftn(u(r),e,n,a)),ifftn:(r,e,n,a)=>m.fromStorage(i.ifftn(u(r),e,n,a)),rfft:(r,e,n,a)=>m.fromStorage(i.rfft(u(r),e,n,a)),irfft:(r,e,n,a)=>m.fromStorage(i.irfft(u(r),e,n,a)),rfft2:(r,e,n,a)=>m.fromStorage(i.rfft2(u(r),e,n,a)),irfft2:(r,e,n,a)=>m.fromStorage(i.irfft2(u(r),e,n,a)),rfftn:(r,e,n,a)=>m.fromStorage(i.rfftn(u(r),e,n,a)),irfftn:(r,e,n,a)=>m.fromStorage(i.irfftn(u(r),e,n,a)),hfft:(r,e,n,a)=>m.fromStorage(i.hfft(u(r),e,n,a)),ihfft:(r,e,n,a)=>m.fromStorage(i.ihfft(u(r),e,n,a)),fftfreq:(r,e)=>m.fromStorage(i.fftfreq(r,e)),rfftfreq:(r,e)=>m.fromStorage(i.rfftfreq(r,e)),fftshift:(r,e)=>m.fromStorage(i.fftshift(u(r),e)),ifftshift:(r,e)=>m.fromStorage(i.ifftshift(u(r),e))};export{v as ArrayStorage,A as Complex,Wm as DTYPE_TO_DESCR,Km as InvalidNpyError,w as NDArrayCore,Qm as SUPPORTED_DTYPES,Jm as UnsupportedDTypeError,Ar as abs,wr as absolute,ie as acos,ve as acosh,mr as add,Jn as all,Ba as allclose,On as amax,En as amin,uo as angle,Kn as any,_m as append,Ro as apply_along_axis,Go as apply_over_axes,D as arange,me as arccos,Ae as arccosh,te as arcsin,he as arcsinh,ue as arctan,be as arctan2,Se as arctanh,Rn as argmax,jn as argmin,Xa as argpartition,Ka as argsort,rt as argwhere,wt as around,I as array,em as array2string,Lo as array_equal,jo as array_equiv,nm as array_repr,pn as array_split,am as array_str,Y as asanyarray,L as asarray,B as asarray_chkfinite,M as ascontiguousarray,V as asfortranarray,oe as asin,we as asinh,se as atan,fe as atan2,ke as atanh,kn as atleast_1d,zn as atleast_2d,Nn as atleast_3d,Fn as average,$o as base_repr,rm as binary_repr,Yt as bincount,Bo as bindex,it as bitwise_and,_t as bitwise_count,gt as bitwise_invert,yt as bitwise_left_shift,bt as bitwise_not,ut as bitwise_or,xt as bitwise_right_shift,st as bitwise_xor,bn as block,fo as broadcast_arrays,po as broadcast_shapes,bo as broadcast_to,fm as byteswap,Am as can_cast,Tr as cbrt,vt as ceil,yo as choose,Vr as clip,un as column_stack,vm as common_type,xo as compress,mn as concat,en as concatenate,mo as conj,io as conjugate,Kt as convolve,U as copy,Pa as copysign,vo as copyto,Wt as corrcoef,Jt as correlate,ne as cos,ye as cosh,ot as count_nonzero,Qt as cov,Ve as cross,Xn as cumprod,Qn as cumsum,Zn as cumulative_prod,Wn as cumulative_sum,ce as deg2rad,le as degrees,ym as delete,gm as delete_,J as diag,qo as diag_indices,Io as diag_indices_from,K as diagflat,De as diagonal,ro as diff,Mt as digitize,zr as divide,Or as divmod,Ne as dot,cn as dsplit,on as dstack,eo as ediff1d,Le as einsum,je as einsum_path,N as empty,G as empty_like,Ya as equal,pr as exp,lr as exp2,Xe as expand_dims,dr as expm1,tt as extract,P as eye,Cr as fabs,y as fft,dm as fill,Do as fill_diagonal,St as fix,et as flatnonzero,Ke as flatten,xn as flip,hn as fliplr,wn as flipud,jr as float_power,kt as floor,qr as floor_divide,Jr as fmax,Kr as fmin,Rr as fmod,Xo as format_float_positional,Zo as format_float_scientific,Gr as frexp,$ as frombuffer,ar as fromfile,rr as fromfunction,er as fromiter,ui as fromregex,nr as fromstring,q as full,H as full_like,Hr as gcd,ii as genfromtxt,O as geomspace,Qo as get_printoptions,Yo as geterr,no as gradient,Ra as greater,Ga as greater_equal,Lr as heaviside,Vt as histogram,Bt as histogram2d,Xt as histogram_bin_edges,Ft as histogramdd,dn as hsplit,tn as hstack,pe as hypot,$r as i0,E as identity,Vo as iindex,oo as imag,Dt as in1d,So as indices,Oe as inner,xm as insert,Wr as interp,Tt as intersect1d,ft as invert,Va as isclose,Sa as iscomplex,ka as iscomplexobj,Ca as isdtype,ya as isfinite,Ia as isfortran,Ct as isin,xa as isinf,ha as isnan,wa as isnat,Aa as isneginf,va as isposinf,za as isreal,Na as isrealobj,Da as isscalar,zm as issubdtype,um as item,Ta as iterable,ko as ix_,Te as kron,Ur as lcm,Yr as ldexp,pt as left_shift,Ha as less,Ua as less_equal,Qa as lexsort,Be as linalg,T as linspace,ri as loadNpz,ei as loadNpzSync,cr as log,gr as log10,yr as log1p,_r as log2,xr as logaddexp,hr as logaddexp2,da as logical_and,_a as logical_not,ca as logical_or,ga as logical_xor,C as logspace,Eo as mask_indices,qe as matmul,He as matrix_transpose,Ye as matvec,Cn as max,Br as maximum,Ho as may_share_memory,Dn as mean,Mn as median,tr as meshgrid,Pn as min,km as min_scalar_type,Fr as minimum,qm as mintypecode,kr as mod,Mr as modf,$e as moveaxis,ur as multiply,Qr as nan_to_num,ia as nanargmax,ma as nanargmin,sa as nancumprod,ua as nancumsum,oa as nanmax,ea as nanmean,ba as nanmedian,ta as nanmin,pa as nanpercentile,ra as nanprod,fa as nanquantile,aa as nanstd,$n as nansum,na as nanvar,om as ndim,vr as negative,La as nextafter,$a as nonzero,Ma as not_equal,z as ones,R as ones_like,Pe as outer,dt as packbits,hm as pad,Um as parseNpy,Mm as parseNpyData,Ym as parseNpyHeader,Zm as parseNpz,$m as parseNpzSync,mi as parseTxt,Wa as partition,Vn as percentile,Ue as permute_dims,wo as place,Dm as poly,Tm as polyadd,Cm as polyder,Om as polydiv,Pm as polyfit,Em as polyint,Lm as polymul,jm as polysub,Rm as polyval,Ir as positive,fr as pow,br as power,Wo as printoptions,Tn as prod,Oa as promote_types,Ln as ptp,co as put,go as put_along_axis,Ao as putmask,Bn as quantile,_e as rad2deg,de as radians,g as random,Qe as ravel,zo as ravel_multi_index,to as real,qa as real_if_close,Dr as reciprocal,Er as remainder,yn as repeat,F as require,Je as reshape,Sn as resize,Sm as result_type,lt as right_shift,zt as rint,vn as roll,rn as rollaxis,Gm as roots,An as rot90,At as round,sn as row_stack,at as searchsorted,ho as select,Bm as serializeNpy,ai as serializeNpz,ti as serializeNpzSync,bi as serializeTxt,Ko as set_printoptions,Ot as setdiff1d,Mo as seterr,Pt as setxor1d,mm as shape,Uo as shares_memory,Sr as sign,Ea as signbit,ee as sin,Zr as sinc,ge as sinh,im as size,Ja as sort,Za as sort_complex,ja as spacing,fn as split,sr as sqrt,Pr as square,We as squeeze,nn as stack,Yn as std,ir as subtract,In as sum,Ze as swapaxes,lo as take,_o as take_along_axis,ae as tan,xe as tanh,Ee as tensordot,gn as tile,bm as tobytes,lm as tofile,sm as tolist,Ie as trace,Ce as transpose,Zt as trapezoid,Q as tri,W as tril,To as tril_indices,Co as tril_indices_from,Lt as trim_zeros,X as triu,Oo as triu_indices,Po as triu_indices_from,Nr as true_divide,Nt as trunc,Nm as typename,Et as union1d,It as unique,jt as unique_all,Rt as unique_counts,Gt as unique_inverse,Ht as unique_values,ct as unpackbits,No as unravel_index,_n as unstack,Xr as unwrap,Z as vander,Un as var,Hn as var_,Gn as variance,Re as vdot,Ge as vecdot,Me as vecmat,pm as view,Fo as vindex,ln as vsplit,an as vstack,nt as where,k as zeros,j as zeros_like};

package/dist/esm/core/linalg.js

@@ -1 +1 @@
-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};
+import*as u from"../common/ops/linalg.js";import{toStorage as n,fromStorage as a,fromStorageView as A,Complex as C,ArrayStorage as s}from"./types.js";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 instanceof C?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/logic.js

@@ -1 +1 @@
-import*as n from"../common/ops/logic";import*as u from"../common/ops/comparison";import{toStorage as o,fromStorage as t}from"./types";function D(r,e){const a=typeof e=="number"?e:o(e);return t(n.logical_and(o(r),a))}function A(r,e){const a=typeof e=="number"?e:o(e);return t(n.logical_or(o(r),a))}function C(r){return t(n.logical_not(o(r)))}function N(r,e){const a=typeof e=="number"?e:o(e);return t(n.logical_xor(o(r),a))}function s(r){return t(n.isfinite(o(r)))}function c(r){return t(n.isinf(o(r)))}function f(r){return t(n.isnan(o(r)))}function l(r){return t(n.isnat(o(r)))}function m(r){return t(n.isneginf(o(r)))}function b(r){return t(n.isposinf(o(r)))}function g(r){return t(n.iscomplex(o(r)))}function x(r){return n.iscomplexobj(o(r))}function _(r){return t(n.isreal(o(r)))}function q(r){return n.isrealobj(o(r))}function S(r,e){return t(n.real_if_close(o(r),e))}function T(r){return n.isfortran(o(r))}function d(r){return n.isscalar(r)}function j(r){return n.iterable(r)}function E(r,e){return n.isdtype(r,e)}function k(r,e){return n.promote_types(r,e)}function w(r,e){const a=typeof e=="number"?e:o(e);return t(n.copysign(o(r),a))}function O(r){return t(n.signbit(o(r)))}function h(r,e){const a=typeof e=="number"?e:o(e);return t(n.nextafter(o(r),a))}function v(r){return t(n.spacing(o(r)))}function z(r,e){return t(u.greater(o(r),typeof e=="number"?e:o(e)))}function B(r,e){return t(u.greaterEqual(o(r),typeof e=="number"?e:o(e)))}function F(r,e){return t(u.less(o(r),typeof e=="number"?e:o(e)))}function G(r,e){return t(u.lessEqual(o(r),typeof e=="number"?e:o(e)))}function H(r,e){return t(u.equal(o(r),typeof e=="number"?e:o(e)))}function I(r,e){return t(u.notEqual(o(r),typeof e=="number"?e:o(e)))}function J(r,e,a=1e-5,i=1e-8){return t(u.isclose(o(r),typeof e=="number"?e:o(e),a,i))}function K(r,e,a=1e-5,i=1e-8){return u.allclose(o(r),typeof e=="number"?e:o(e),a,i)}export{K as allclose,w as copysign,H as equal,z as greater,B as greater_equal,J as isclose,g as iscomplex,x as iscomplexobj,E as isdtype,s as isfinite,T as isfortran,c as isinf,f as isnan,l as isnat,m as isneginf,b as isposinf,_ as isreal,q as isrealobj,d as isscalar,j as iterable,F as less,G as less_equal,D as logical_and,C as logical_not,A as logical_or,N as logical_xor,h as nextafter,I as not_equal,k as promote_types,S as real_if_close,O as signbit,v as spacing};
+import*as n from"../common/ops/logic.js";import*as u from"../common/ops/comparison.js";import{toStorage as o,fromStorage as t}from"./types.js";function D(r,e){const a=typeof e=="number"?e:o(e);return t(n.logical_and(o(r),a))}function A(r,e){const a=typeof e=="number"?e:o(e);return t(n.logical_or(o(r),a))}function C(r){return t(n.logical_not(o(r)))}function N(r,e){const a=typeof e=="number"?e:o(e);return t(n.logical_xor(o(r),a))}function s(r){return t(n.isfinite(o(r)))}function c(r){return t(n.isinf(o(r)))}function f(r){return t(n.isnan(o(r)))}function l(r){return t(n.isnat(o(r)))}function m(r){return t(n.isneginf(o(r)))}function b(r){return t(n.isposinf(o(r)))}function g(r){return t(n.iscomplex(o(r)))}function x(r){return n.iscomplexobj(o(r))}function _(r){return t(n.isreal(o(r)))}function q(r){return n.isrealobj(o(r))}function S(r,e){return t(n.real_if_close(o(r),e))}function T(r){return n.isfortran(o(r))}function d(r){return n.isscalar(r)}function j(r){return n.iterable(r)}function E(r,e){return n.isdtype(r,e)}function k(r,e){return n.promote_types(r,e)}function w(r,e){const a=typeof e=="number"?e:o(e);return t(n.copysign(o(r),a))}function O(r){return t(n.signbit(o(r)))}function h(r,e){const a=typeof e=="number"?e:o(e);return t(n.nextafter(o(r),a))}function v(r){return t(n.spacing(o(r)))}function z(r,e){return t(u.greater(o(r),typeof e=="number"?e:o(e)))}function B(r,e){return t(u.greaterEqual(o(r),typeof e=="number"?e:o(e)))}function F(r,e){return t(u.less(o(r),typeof e=="number"?e:o(e)))}function G(r,e){return t(u.lessEqual(o(r),typeof e=="number"?e:o(e)))}function H(r,e){return t(u.equal(o(r),typeof e=="number"?e:o(e)))}function I(r,e){return t(u.notEqual(o(r),typeof e=="number"?e:o(e)))}function J(r,e,a=1e-5,i=1e-8){return t(u.isclose(o(r),typeof e=="number"?e:o(e),a,i))}function K(r,e,a=1e-5,i=1e-8){return u.allclose(o(r),typeof e=="number"?e:o(e),a,i)}export{K as allclose,w as copysign,H as equal,z as greater,B as greater_equal,J as isclose,g as iscomplex,x as iscomplexobj,E as isdtype,s as isfinite,T as isfortran,c as isinf,f as isnan,l as isnat,m as isneginf,b as isposinf,_ as isreal,q as isrealobj,d as isscalar,j as iterable,F as less,G as less_equal,D as logical_and,C as logical_not,A as logical_or,N as logical_xor,h as nextafter,I as not_equal,k as promote_types,S as real_if_close,O as signbit,v as spacing};

package/dist/esm/core/polynomial.js

@@ -1 +1 @@
-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};
+import{NDArrayCore as T}from"../common/ndarray-core.js";import{ArrayStorage as q}from"../common/storage.js";import{isBigIntDType as F}from"../common/dtype.js";import{array as d}from"./creation.js";function C(e){return e instanceof T?e:d(e)}function R(e,i){const a=e[i];return typeof a=="bigint"?Number(a):a}function j(e){return e==="complex128"||e==="complex64"}function w(e){const i=e.data,a=e.size;if(j(e.dtype)){const t=[];for(let l=0;l<a;l++)t.push(Number(i[2*l]));return t}const s=[];for(let t=0;t<a;t++)s.push(R(i,t));return s}function _(e){if(e==="bool")throw new TypeError("numpy boolean subtract, the `-` operator, is not supported, use the bitwise_xor, the `^` operator, or the logical_xor function instead.")}function z(e){return e==="complex128"||e==="complex64"?"complex128":e==="float64"||e==="float32"||e==="float16"||e==="uint64"?"float64":"int64"}function E(e){return e==="complex128"||e==="complex64"?"complex128":"float64"}function B(e){return e==="complex128"||e==="complex64"?e:e==="float32"?"float32":e==="float16"?"float16":"float64"}function Z(e){const i=C(e),a=i.size,s=i.dtype==="float32"||i.dtype==="complex64"?"float32":"float64";if(a===0)return d([1],s);const t=w(i);let l=[1];for(let f=0;f<a;f++){const c=t[f],n=new Array(l.length+1).fill(0);for(let r=0;r<l.length;r++)n[r]+=l[r],n[r+1]-=l[r]*c;l=n}return d(l,s)}function G(e,i){const a=C(e),s=C(i),t=w(a),l=w(s),f=Math.max(t.length,l.length),c=new Array(f).fill(0);for(let r=0;r<t.length;r++)c[f-t.length+r]+=t[r];for(let r=0;r<l.length;r++)c[f-l.length+r]+=l[r];let n=0;for(;n<c.length-1&&c[n]===0;)n++;return d(c.slice(n),a.dtype)}function J(e,i=1){let a=C(e);const s=z(a.dtype);for(let t=0;t<i;t++){const l=w(a),f=l.length;if(f<=1)return d([0],s);const c=[];for(let n=0;n<f-1;n++){const r=f-1-n;c.push(l[n]*r)}a=d(c,s)}return a}function K(e,i){const a=C(e),s=C(i),t=w(a),l=w(s),f=B(a.dtype);if(l.length===0||l.length===1&&l[0]===0)throw new Error("Division by zero polynomial");for(;t.length>1&&t[0]===0;)t.shift();for(;l.length>1&&l[0]===0;)l.shift();if(t.length<l.length)return[d([0],f),d(t,f)];const c=[],n=[...t];for(;n.length>=l.length;){const r=n[0]/l[0];c.push(r);for(let g=0;g<l.length;g++)n[g]-=r*l[g];n.shift()}for(;n.length>1&&Math.abs(n[0])<1e-15;)n.shift();return[d(c.length>0?c:[0],f),d(n.length>0?n:[0],f)]}function U(e,i,a){if(e.dtype==="float16")throw new TypeError("array type float16 is unsupported in linalg");const s=e.size;if(a>=s)throw new Error("polyfit: degree must be less than number of points");const t=w(e),l=w(i),f=new Float64Array(s),c=new Float64Array(s);for(let o=0;o<s;o++)f[o]=t[o],c[o]=l[o];const n=[];for(let o=0;o<s;o++){const h=[];for(let m=a;m>=0;m--)h.push(Math.pow(f[o],m));n.push(h)}const r=[],g=[];for(let o=0;o<=a;o++){r.push([]);for(let m=0;m<=a;m++){let A=0;for(let D=0;D<s;D++)A+=n[D][o]*n[D][m];r[o].push(A)}let h=0;for(let m=0;m<s;m++)h+=n[m][o]*c[m];g.push(h)}const b=a+1,p=r.map((o,h)=>[...o,g[h]]);for(let o=0;o<b;o++){let h=o;for(let m=o+1;m<b;m++)Math.abs(p[m][o])>Math.abs(p[h][o])&&(h=m);[p[o],p[h]]=[p[h],p[o]];for(let m=o+1;m<b;m++){const A=p[m][o]/p[o][o];for(let D=o;D<=b;D++)p[m][D]-=A*p[o][D]}}const y=new Array(b).fill(0);for(let o=b-1;o>=0;o--){let h=p[o][b];for(let m=o+1;m<b;m++)h-=p[o][m]*y[m];y[o]=h/p[o][o]}const u=e.dtype==="complex128"||e.dtype==="complex64"?"complex128":"float64";return d(y,u)}function W(e,i=1,a=0){let s=C(e);const t=E(s.dtype),l=Array.isArray(a)?a:[a];for(let f=0;f<i;f++){const c=w(s),n=c.length,r=[];for(let b=0;b<n;b++){const p=n-b;r.push(c[b]/p)}const g=f<l.length?l[f]:0;r.push(g),s=d(r,t)}return s}function X(e,i){const a=C(e),s=C(i),t=w(a),l=w(s),f=t.length+l.length-1,c=new Array(f).fill(0);for(let n=0;n<t.length;n++)for(let r=0;r<l.length;r++)c[n+r]+=t[n]*l[r];return d(c,a.dtype)}function Y(e,i){const a=C(e),s=C(i);_(a.dtype);const t=w(a),l=w(s),f=Math.max(t.length,l.length),c=new Array(f).fill(0);for(let r=0;r<t.length;r++)c[f-t.length+r]+=t[r];for(let r=0;r<l.length;r++)c[f-l.length+r]-=l[r];let n=0;for(;n<c.length-1&&c[n]===0;)n++;return d(c.slice(n),a.dtype)}function $(e,i){const a=C(e),s=w(a);if(typeof i=="number"){let y=s[0];for(let u=1;u<s.length;u++)y=y*i+s[u];return y}const t=i instanceof T?i:d(i),l=w(t),f=t.size,c=s.length,n=a.dtype,r=q.empty(Array.from(t.shape),n),g=r.data,b=j(n),p=F(n);for(let y=0;y<f;y++){const u=l[y];let o=s[0];for(let h=1;h<c;h++)o=o*u+s[h];n==="bool"&&(o=o?1:0),b?(g[2*y]=o,g[2*y+1]=0):p?g[y]=BigInt(Math.round(o)):g[y]=o}return new T(r)}function H(e){const i=C(e);if(i.dtype==="float16")throw new TypeError("array type float16 is unsupported in linalg");const a=i.dtype==="float32"||i.dtype==="complex64"?"complex64":"complex128",s=w(i);for(;s.length>1&&s[0]===0;)s.shift();let t=0;for(;s.length>1&&s[s.length-1]===0;)s.pop(),t++;const l=s.length-1;if(l+t===0)return k([],[],a);let c=[],n=[];if(l===1)c.push(-s[1]/s[0]),n.push(0);else if(l===2){const p=s[0],y=s[1],u=s[2],o=y*y-4*p*u;if(o>=0){const h=Math.sqrt(o);c.push((-y+h)/(2*p),(-y-h)/(2*p)),n.push(0,0)}else{const h=Math.sqrt(-o);c.push(-y/(2*p),-y/(2*p)),n.push(h/(2*p),-h/(2*p))}}else if(l>=3){const p=S(s,l);for(const y of p)c.push(y.re),n.push(y.im)}for(let p=0;p<t;p++)c.push(0),n.push(0);const r=c.map((p,y)=>y);r.sort((p,y)=>{const u=Math.sqrt(c[p]**2+n[p]**2),o=Math.sqrt(c[y]**2+n[y]**2);return Math.abs(u-o)>1e-10?o-u:Math.abs(c[p]-c[y])>1e-10?c[y]-c[p]:n[y]-n[p]});const g=r.map(p=>c[p]),b=r.map(p=>n[p]);return k(g,b,a)}function k(e,i,a="complex128"){const s=e.length,t=q.empty([s],a),l=t.data;for(let f=0;f<s;f++)l[2*f]=e[f],l[2*f+1]=i[f];return new T(t)}function S(e,i){const a=Array.from({length:i},()=>new Array(i).fill(0)),s=e[0];for(let t=0;t<i;t++)a[0][t]=-e[t+1]/s;for(let t=1;t<i;t++)a[t][t-1]=1;return L(a,i)}function L(e,i){const a=[];let t=i,l=0;const f=100*i;let c=i,n=0;for(;t>0&&l<f;){if(l++,t===c?n++:(n=0,c=t),t===1){a.push({re:e[0][0],im:0}),t=0;break}if(t===2){a.push(...I(e[0][0],e[0][1],e[1][0],e[1][1])),t=0;break}let r=t-1;for(;r>0;){const u=Math.abs(e[r-1][r-1])+Math.abs(e[r][r]),o=222e-18*(u===0?1:u);if(Math.abs(e[r][r-1])<=o){e[r][r-1]=0;break}r--}const g=t-r;if(g===1){a.push({re:e[t-1][t-1],im:0}),t--;continue}if(g===2){a.push(...I(e[t-2][t-2],e[t-2][t-1],e[t-1][t-2],e[t-1][t-1])),t-=2;continue}let b;if(n>0&&n%10===0)b=Math.abs(e[t-1][t-2])+Math.abs(e[t-2][t-3]);else{const u=e[t-2][t-2],o=e[t-2][t-1],h=e[t-1][t-2],m=e[t-1][t-1],A=u+m,D=u*m-o*h,x=A*A-4*D;if(x>=0){const N=Math.sqrt(x),M=(A+N)/2,v=(A-N)/2;b=Math.abs(M-m)<Math.abs(v-m)?M:v}else b=m}for(let u=r;u<t;u++)e[u][u]=e[u][u]-b;const p=[],y=[];for(let u=r;u<t-1;u++){const o=Math.hypot(e[u][u],e[u+1][u]),h=o===0?1:e[u][u]/o,m=o===0?0:e[u+1][u]/o;p.push(h),y.push(m);for(let A=u;A<t;A++){const D=e[u][A],x=e[u+1][A];e[u][A]=h*D+m*x,e[u+1][A]=-m*D+h*x}}for(let u=0;u<p.length;u++){const o=r+u,h=p[u],m=y[u],A=Math.min(o+2,t-1);for(let D=r;D<=A;D++){const x=e[D][o],N=e[D][o+1];e[D][o]=h*x+m*N,e[D][o+1]=-m*x+h*N}}for(let u=r;u<t;u++)e[u][u]=e[u][u]+b}if(t>0)for(let r=0;r<t;r++)a.push({re:e[r][r],im:0});return a}function I(e,i,a,s){const t=e+s,l=e*s-i*a,f=t*t-4*l;if(f>=0){const c=Math.sqrt(f);return[{re:(t+c)/2,im:0},{re:(t-c)/2,im:0}]}else{const c=Math.sqrt(-f);return[{re:t/2,im:c/2},{re:t/2,im:-c/2}]}}export{Z as poly,G as polyadd,J as polyder,K as polydiv,U as polyfit,W as polyint,X as polymul,Y as polysub,$ as polyval,H as roots};

package/dist/esm/core/reduction.js

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

package/dist/esm/core/rounding.js

@@ -1 +1 @@
-import*as o from"../common/ops/rounding";import{toStorage as e,fromStorage as n}from"./types";function u(r,t=0){return n(o.around(e(r),t))}const y=u;function A(r,t=0){return n(o.round(e(r),t))}function C(r){return n(o.ceil(e(r)))}function D(r){return n(o.fix(e(r)))}function N(r){return n(o.floor(e(r)))}function f(r){return n(o.rint(e(r)))}function c(r){return n(o.trunc(e(r)))}export{u as around,C as ceil,D as fix,N as floor,f as rint,A as round,y as round_,c as trunc};
+import*as o from"../common/ops/rounding.js";import{toStorage as e,fromStorage as n}from"./types.js";function i(r,t=0){return n(o.around(e(r),t))}function y(r,t=0){return n(o.round(e(r),t))}function A(r){return n(o.ceil(e(r)))}function C(r){return n(o.fix(e(r)))}function D(r){return n(o.floor(e(r)))}function N(r){return n(o.rint(e(r)))}function f(r){return n(o.trunc(e(r)))}export{i as around,A as ceil,C as fix,D as floor,N as rint,y as round,f as trunc};

package/dist/esm/core/sets.js

@@ -1 +1 @@
-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};
+import*as u from"../common/ops/sets.js";import{toStorage as o,fromStorage as n}from"./types.js";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 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};
+import{NDArrayCore as z}from"../common/ndarray-core.js";import{ArrayStorage as k}from"../common/storage.js";import{Complex as I}from"../common/complex.js";import{array as N,zeros as F}from"./creation.js";import{concatenate as T,flatten as w}from"./shape.js";import{wasmPad2D as L}from"../common/wasm/pad.js";import{hasFloat16 as M,isComplexDType as P}from"../common/dtype.js";function G(u,a,c){const d=a instanceof z?a:N(Array.isArray(a)?a:[a]);if(c===void 0){const s=w(u),D=w(d);return T([s,D],0)}return T([u,d],c)}function v(u,a,c){const d=Array.isArray(a)?a:[a],s=[...u.shape],D=u.storage;if(c===void 0){const t=w(u),r=t.storage,n=[],p=new Set(d.map(l=>l<0?t.size+l:l));for(let l=0;l<t.size;l++)p.has(l)||n.push(r.iget(l));return N(n,u.dtype)}const m=c<0?s.length+c:c,g=s[m],x=new Set(d.map(t=>t<0?g+t:t)),S=g-x.size,b=[...s];b[m]=S;const h=F(b,u.dtype),C=h.storage,f=[];let o=1;for(let t=s.length-1;t>=0;t--)f.unshift(o),o*=s[t];const e=[];o=1;for(let t=b.length-1;t>=0;t--)e.unshift(o),o*=b[t];const y=s.reduce((t,r)=>t*r,1);for(let t=0;t<y;t++){const r=[];let n=t;for(let A=0;A<s.length;A++)r.push(Math.floor(n/f[A])),n=n%f[A];if(x.has(r[m]))continue;let p=0;for(const A of x)A<r[m]&&p++;const l=[...r];l[m]-=p;let i=0;for(let A=0;A<b.length;A++)i+=l[A]*e[A];C.iset(i,D.iget(t))}return h}function H(u,a,c,d){const s=Array.isArray(a)?a:[a],D=c instanceof z?c:N(Array.isArray(c)?c:[c]);if(d===void 0){const m=w(u),g=w(D),x=u.dtype;if(P(x)){const f=m.size,o=g.size,e=[];for(let r=0;r<f;r++)e.push(m.storage.iget(r));if(s.length===1){const r=s[0],n=r<0?f+r:r,p=[];for(let l=0;l<o;l++){const i=g.storage.iget(l);p.push(i instanceof I?i:new I(Number(i),0))}e.splice(n,0,...p)}else{const r=s.map((n,p)=>({idx:n<0?f+n:n,valIdx:p})).sort((n,p)=>n.idx-p.idx);for(let n=0;n<r.length;n++){const{idx:p,valIdx:l}=r[n],i=g.storage.iget(l%o),A=i instanceof I?i:new I(Number(i),0);e.splice(p+n,0,A)}}const y=k.empty([e.length],x),t=y.data;for(let r=0;r<e.length;r++)t[r*2]=e[r].re,t[r*2+1]=e[r].im;return new z(y)}const b=m.data,h=g.data,C=Array.from(b);if(s.length===1){const f=s[0],o=f<0?m.size+f:f,e=Array.from(h);C.splice(o,0,...e)}else{const f=s.map((o,e)=>({idx:o<0?m.size+o:o,valIdx:e})).sort((o,e)=>o.idx-e.idx);for(let o=0;o<f.length;o++){const{idx:e,valIdx:y}=f[o],t=h[y%h.length];C.splice(e+o,0,t)}}return N(C,x)}throw new Error("insert along axis not fully implemented in standalone")}function J(u,a,c="constant",d=0){const s=[...u.shape],D=s.length;let m;typeof a=="number"?m=Array(D).fill([a,a]):Array.isArray(a)&&typeof a[0]=="number"?m=Array(D).fill(a):m=a;const g=s.map((e,y)=>e+m[y][0]+m[y][1]);if(c==="constant"&&d===0&&D===2&&typeof a=="number"){const e=L(u.storage,a);if(e)return new z(e)}if(c!=="constant")throw new Error(`pad mode '${c}' not fully implemented in standalone`);const x=F(g,u.dtype),S=x.storage;if(d!==0)for(let e=0;e<x.size;e++)S.iset(e,d);const b=u.storage,h=[];let C=1;for(let e=s.length-1;e>=0;e--)h.unshift(C),C*=s[e];const f=[];C=1;for(let e=g.length-1;e>=0;e--)f.unshift(C),C*=g[e];const o=s.reduce((e,y)=>e*y,1);if(u.dtype==="float16"&&M&&b.isCContiguous&&S.isCContiguous){const e=new Float32Array(b.data.subarray(b.offset,b.offset+o)),y=S.size,t=new Float32Array(y);d!==0&&t.fill(d);for(let r=0;r<o;r++){const n=[];let p=r;for(let i=0;i<s.length;i++)n.push(Math.floor(p/h[i])),p=p%h[i];let l=0;for(let i=0;i<g.length;i++)l+=(n[i]+m[i][0])*f[i];t[l]=e[r]}return S.data.set(t),x}for(let e=0;e<o;e++){const y=[];let t=e;for(let n=0;n<s.length;n++)y.push(Math.floor(t/h[n])),t=t%h[n];let r=0;for(let n=0;n<g.length;n++)r+=(y[n]+m[n][0])*f[n];S.iset(r,b.iget(e))}return x}export{G as append,v as delete,v as delete_,H as insert,J as pad};

package/dist/esm/core/shape.js

@@ -1 +1 @@
-import*as o from"../common/ops/shape";import{toStorage as t,fromStorage as a,fromStorageView as u,fromStorageArray as A,fromStorageViewArray as y}from"./types";function m(r,e){return a(o.reshape(t(r),e))}function i(r){return a(o.flatten(t(r)))}function c(r){return a(o.ravel(t(r)))}function s(r,e){return u(o.squeeze(t(r),e),r)}function f(r,e){return u(o.expandDims(t(r),e),r)}function l(r,e,n){return u(o.swapaxes(t(r),e,n),r)}function b(r,e,n){return u(o.moveaxis(t(r),e,n),r)}function x(r,e,n=0){return u(o.rollaxis(t(r),e,n),r)}function k(r,e=0){return a(o.concatenate(r.map(n=>t(n)),e))}function d(r,e=0){return a(o.stack(r.map(n=>t(n)),e))}function D(r){return a(o.vstack(r.map(e=>t(e))))}function h(r){return a(o.hstack(r.map(e=>t(e))))}function g(r){return a(o.dstack(r.map(e=>t(e))))}function v(r,e=0){return a(o.concat(r.map(n=>t(n)),e))}function S(r){return a(o.columnStack(r.map(e=>t(e))))}const _=D;function w(r){return a(o.block(r.map(e=>t(e))))}function z(r,e,n=0){return y(o.split(t(r),e,n),r)}function q(r,e,n=0){return y(o.arraySplit(t(r),e,n),r)}function V(r,e){return y(o.vsplit(t(r),e),r)}function O(r,e){return y(o.hsplit(t(r),e),r)}function j(r,e){return y(o.dsplit(t(r),e),r)}function B(r,e=0){return A(o.unstack(t(r),e))}function E(r,e){return a(o.tile(t(r),e))}function F(r,e,n){return a(o.repeat(t(r),e,n))}function p(r,e){return a(o.flip(t(r),e))}function G(r){return p(r,1)}function H(r){return p(r,0)}function I(r,e=1,n=[0,1]){return a(o.rot90(t(r),e,n))}function J(r,e,n){return a(o.roll(t(r),e,n))}function K(r,e){return a(o.resize(t(r),e))}function L(...r){const e=o.atleast1d(r.map(n=>t(n)));return e.length===1?a(e[0]):A(e)}function M(...r){const e=o.atleast2d(r.map(n=>t(n)));return e.length===1?a(e[0]):A(e)}function P(...r){const e=o.atleast3d(r.map(n=>t(n)));return e.length===1?a(e[0]):A(e)}export{q as array_split,L as atleast_1d,M as atleast_2d,P as atleast_3d,w as block,S as column_stack,v as concat,k as concatenate,j as dsplit,g as dstack,f as expand_dims,i as flatten,p as flip,G as fliplr,H as flipud,O as hsplit,h as hstack,b as moveaxis,c as ravel,F as repeat,m as reshape,K as resize,J as roll,x as rollaxis,I as rot90,_ as row_stack,z as split,s as squeeze,d as stack,l as swapaxes,E as tile,B as unstack,V as vsplit,D as vstack};
+import*as o from"../common/ops/shape.js";import{toStorage as t,fromStorage as a,fromStorageView as u,fromStorageArray as A,fromStorageViewArray as y}from"./types.js";function m(r,e){return a(o.reshape(t(r),e))}function i(r){return a(o.flatten(t(r)))}function c(r){return a(o.ravel(t(r)))}function s(r,e){return u(o.squeeze(t(r),e),r)}function f(r,e){return u(o.expandDims(t(r),e),r)}function l(r,e,n){return u(o.swapaxes(t(r),e,n),r)}function b(r,e,n){return u(o.moveaxis(t(r),e,n),r)}function x(r,e,n=0){return u(o.rollaxis(t(r),e,n),r)}function k(r,e=0){return a(o.concatenate(r.map(n=>t(n)),e))}function d(r,e=0){return a(o.stack(r.map(n=>t(n)),e))}function D(r){return a(o.vstack(r.map(e=>t(e))))}function h(r){return a(o.hstack(r.map(e=>t(e))))}function g(r){return a(o.dstack(r.map(e=>t(e))))}function v(r,e=0){return a(o.concat(r.map(n=>t(n)),e))}function S(r){return a(o.columnStack(r.map(e=>t(e))))}const _=D;function w(r){return a(o.block(r.map(e=>t(e))))}function z(r,e,n=0){return y(o.split(t(r),e,n),r)}function q(r,e,n=0){return y(o.arraySplit(t(r),e,n),r)}function V(r,e){return y(o.vsplit(t(r),e),r)}function O(r,e){return y(o.hsplit(t(r),e),r)}function j(r,e){return y(o.dsplit(t(r),e),r)}function B(r,e=0){return A(o.unstack(t(r),e))}function E(r,e){return a(o.tile(t(r),e))}function F(r,e,n){return a(o.repeat(t(r),e,n))}function p(r,e){return a(o.flip(t(r),e))}function G(r){return p(r,1)}function H(r){return p(r,0)}function I(r,e=1,n=[0,1]){return a(o.rot90(t(r),e,n))}function J(r,e,n){return a(o.roll(t(r),e,n))}function K(r,e){return a(o.resize(t(r),e))}function L(...r){const e=o.atleast1d(r.map(n=>t(n)));return e.length===1?a(e[0]):A(e)}function M(...r){const e=o.atleast2d(r.map(n=>t(n)));return e.length===1?a(e[0]):A(e)}function P(...r){const e=o.atleast3d(r.map(n=>t(n)));return e.length===1?a(e[0]):A(e)}export{q as array_split,L as atleast_1d,M as atleast_2d,P as atleast_3d,w as block,S as column_stack,v as concat,k as concatenate,j as dsplit,g as dstack,f as expand_dims,i as flatten,p as flip,G as fliplr,H as flipud,O as hsplit,h as hstack,b as moveaxis,c as ravel,F as repeat,m as reshape,K as resize,J as roll,x as rollaxis,I as rot90,_ as row_stack,z as split,s as squeeze,d as stack,l as swapaxes,E as tile,B as unstack,V as vsplit,D as vstack};

package/dist/esm/core/sorting.js

@@ -1 +1 @@
-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};
+import*as t from"../common/ops/sorting.js";import{toStorage as o,fromStorage as n,fromStorageArray as D}from"./types.js";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/core/statistics.js

@@ -1 +1 @@
-import*as A from"../common/ops/statistics";import{NDArrayCore as d,toStorage as o,fromStorage as t}from"./types";function g(n,e,r){const a=e?o(e):void 0;return t(A.bincount(o(n),a,r))}function b(n,e,r){return t(A.digitize(o(n),o(e),r))}function f(n,e,r,a,u){const y=e instanceof d?o(e):e,i=u?o(u):void 0,m=A.histogram(o(n),y,r,a,i);return[t(m.hist),t(m.bin_edges)]}function l(n,e,r,a,u,y){let i;if(Array.isArray(r)&&r.length===2){const C=r[0]instanceof d?o(r[0]):r[0],D=r[1]instanceof d?o(r[1]):r[1];i=[C,D]}else i=r;const m=y?o(y):void 0,s=A.histogram2d(o(n),o(e),i,a,u,m);return[t(s.hist),t(s.x_edges),t(s.y_edges)]}function p(n,e,r,a,u){const y=u?o(u):void 0,i=A.histogramdd(o(n),e,r,a,y);return[t(i.hist),i.edges.map(m=>t(m))]}function S(n,e,r){return t(A.correlate(o(n),o(e),r))}function h(n,e,r){return t(A.convolve(o(n),o(e),r))}function x(n,e,r,a,u){return t(A.cov(o(n),e?o(e):void 0,r,a,u))}function v(n,e,r){return t(A.corrcoef(o(n),e?o(e):void 0,r))}function _(n,e,r,a){return t(A.histogram_bin_edges(o(n),e,r,a?o(a):void 0))}function z(n,e,r,a){const u=A.trapezoid(o(n),e?o(e):void 0,r,a);return typeof u=="number"?u:t(u)}export{g as bincount,h as convolve,v as corrcoef,S as correlate,x as cov,b as digitize,f as histogram,l as histogram2d,_ as histogram_bin_edges,p as histogramdd,z as trapezoid};
+import*as i from"../common/ops/statistics.js";import{Complex as D}from"../common/complex.js";import{NDArrayCore as C,toStorage as o,fromStorage as a}from"./types.js";function b(n,e,r){const u=e?o(e):void 0;return a(i.bincount(o(n),u,r))}function l(n,e,r){return a(i.digitize(o(n),o(e),r))}function p(n,e,r,u,t){const y=e instanceof C?o(e):e,A=t?o(t):void 0,m=i.histogram(o(n),y,r,u,A),s=m.hist;return[a(s),a(m.bin_edges)]}function S(n,e,r,u,t,y){let A;if(Array.isArray(r)&&r.length===2){const c=r[0]instanceof C?o(r[0]):r[0],d=r[1]instanceof C?o(r[1]):r[1];A=[c,d]}else A=r;const m=y?o(y):void 0,s=i.histogram2d(o(n),o(e),A,u,t,m);return[a(s.hist),a(s.x_edges),a(s.y_edges)]}function x(n,e,r,u,t){const y=t?o(t):void 0,A=i.histogramdd(o(n),e,r,u,y);return[a(A.hist),A.edges.map(m=>a(m))]}function h(n,e,r){return a(i.correlate(o(n),o(e),r))}function v(n,e,r){return a(i.convolve(o(n),o(e),r))}function _(n,e,r,u,t){return a(i.cov(o(n),e?o(e):void 0,r,u,t))}function z(n,e,r){return a(i.corrcoef(o(n),e?o(e):void 0,r))}function B(n,e,r,u){return a(i.histogram_bin_edges(o(n),e,r,u?o(u):void 0))}function q(n,e,r,u){const t=i.trapezoid(o(n),e?o(e):void 0,r,u);return typeof t=="number"||t instanceof D?t:a(t)}export{b as bincount,v as convolve,z as corrcoef,h as correlate,_ as cov,l as digitize,p as histogram,S as histogram2d,B as histogram_bin_edges,x as histogramdd,q as trapezoid};

package/dist/esm/core/trig.js

@@ -1 +1 @@
-import*as t from"../common/ops/trig";import*as a from"../common/ops/hyperbolic";import{toStorage as o,fromStorage as e}from"./types";function f(r){return e(t.sin(o(r)))}function x(r){return e(t.cos(o(r)))}function h(r){return e(t.tan(o(r)))}function y(r){return e(t.arcsin(o(r)))}const d=y;function u(r){return e(t.arccos(o(r)))}const g=u;function A(r){return e(t.arctan(o(r)))}const m=A;function C(r,n){const c=typeof n=="number"?n:o(n);return e(t.arctan2(o(r),c))}const b=C;function S(r,n){const c=typeof n=="number"?n:o(n);return e(t.hypot(o(r),c))}function O(r){return e(t.degrees(o(r)))}function l(r){return e(t.radians(o(r)))}function j(r){return e(t.deg2rad(o(r)))}function k(r){return e(t.rad2deg(o(r)))}function q(r){return e(a.sinh(o(r)))}function v(r){return e(a.cosh(o(r)))}function w(r){return e(a.tanh(o(r)))}function D(r){return e(a.arcsinh(o(r)))}const z=D;function N(r){return e(a.arccosh(o(r)))}const B=N;function s(r){return e(a.arctanh(o(r)))}const E=s;export{g as acos,B as acosh,u as arccos,N as arccosh,y as arcsin,D as arcsinh,A as arctan,C as arctan2,s as arctanh,d as asin,z as asinh,m as atan,b as atan2,E as atanh,x as cos,v as cosh,j as deg2rad,O as degrees,S as hypot,k as rad2deg,l as radians,f as sin,q as sinh,h as tan,w as tanh};
+import*as t from"../common/ops/trig.js";import*as a from"../common/ops/hyperbolic.js";import{toStorage as o,fromStorage as e}from"./types.js";function f(r){return e(t.sin(o(r)))}function x(r){return e(t.cos(o(r)))}function h(r){return e(t.tan(o(r)))}function y(r){return e(t.arcsin(o(r)))}const d=y;function u(r){return e(t.arccos(o(r)))}const g=u;function A(r){return e(t.arctan(o(r)))}const m=A;function C(r,n){const c=typeof n=="number"?n:o(n);return e(t.arctan2(o(r),c))}const b=C;function S(r,n){const c=typeof n=="number"?n:o(n);return e(t.hypot(o(r),c))}function O(r){return e(t.degrees(o(r)))}function l(r){return e(t.radians(o(r)))}function j(r){return e(t.deg2rad(o(r)))}function k(r){return e(t.rad2deg(o(r)))}function q(r){return e(a.sinh(o(r)))}function v(r){return e(a.cosh(o(r)))}function w(r){return e(a.tanh(o(r)))}function D(r){return e(a.arcsinh(o(r)))}const z=D;function N(r){return e(a.arccosh(o(r)))}const B=N;function s(r){return e(a.arctanh(o(r)))}const E=s;export{g as acos,B as acosh,u as arccos,N as arccosh,y as arcsin,D as arcsinh,A as arctan,C as arctan2,s as arctanh,d as asin,z as asinh,m as atan,b as atan2,E as atanh,x as cos,v as cosh,j as deg2rad,O as degrees,S as hypot,k as rad2deg,l as radians,f as sin,q as sinh,h as tan,w as tanh};

package/dist/esm/core/typechecking.js

@@ -1 +1 @@
-import{NDArrayCore as m}from"../common/ndarray-core";import{isComplexDType as p,isFloatDType as l,isIntegerDType as a}from"../common/dtype";const y=["int8","int16","int32","int64"],D=["uint8","uint16","uint32","uint64"],d=[...y,...D],T=["float32","float64"],x=["complex64","complex128"],R=[...d,...T,...x],h={bool:0,int8:1,uint8:2,int16:3,uint16:4,int32:5,uint32:6,int64:7,uint64:8,float32:9,float64:10,complex64:11,complex128:12};function C(t,n,e="safe"){const r=t instanceof m?t.dtype:t;if(e==="no"||e==="equiv")return r===n;if(e==="unsafe")return!0;const u=h[r]??-1,o=h[n]??-1;if(e==="safe"){if(a(r)&&l(n)){const c={int8:8,int16:16,int32:32,int64:64,uint8:8,uint16:16,uint32:32,uint64:64,bool:1,float32:0,float64:0,complex64:0,complex128:0},s={float32:24,float64:53,int8:0,int16:0,int32:0,int64:0,uint8:0,uint16:0,uint32:0,uint64:0,bool:0,complex64:24,complex128:53},i=c[r]??64,f=s[n]??0;return i<=f}return!!(a(r)&&p(n)||l(r)&&p(n)||u<=o)}return e==="same_kind"?!!(a(r)&&a(n)||l(r)&&l(n)||p(r)&&p(n)||r==="bool"&&a(n)||r==="bool"&&n==="bool"||l(r)&&p(n)):!1}function B(...t){if(t.length===0)return"float64";let n=!1,e=!1;for(const r of t){const u=r.dtype;p(u)?(n=!0,u==="complex128"&&(e=!0)):l(u)?u==="float64"&&(e=!0):a(u)&&(u==="int64"||u==="uint64"||u==="int32"||u==="uint32")&&(e=!0)}return n?e?"complex128":"complex64":e?"float64":"float32"}function _(...t){if(t.length===0)return"float64";const n=t.map(o=>o instanceof m?o.dtype:o),e=n.some(p),r=n.some(l),u=n.some(a);if(e){const o=n.some(i=>i==="float64"),c=n.some(i=>i==="complex128"),s=n.some(i=>i==="int64"||i==="uint64"||i==="int32"||i==="uint32");return c||o||s?"complex128":"complex64"}if(r){const o=n.some(c=>c==="float64");return u&&(n.some(s=>s==="int32"||s==="uint32"||s==="int64"||s==="uint64")||o)||o?"float64":"float32"}if(u){const o=n.some(f=>y.includes(f)),c=n.some(f=>D.includes(f)),s={bool:8,int8:8,uint8:8,int16:16,uint16:16,int32:32,uint32:32,int64:64,uint64:64,float32:0,float64:0,complex64:0,complex128:0},i=Math.max(...n.map(f=>s[f]??0));return o&&c?i>=32?"int64":i>=16?"int32":"int16":o?i>=64?"int64":i>=32?"int32":i>=16?"int16":"int8":i>=64?"uint64":i>=32?"uint32":i>=16?"uint16":"uint8"}return"bool"}function F(t){return typeof t=="boolean"?"bool":typeof t=="bigint"?t>=0n?t<=255n?"uint8":t<=65535n?"uint16":t<=4294967295n?"uint32":"uint64":t>=-128n&&t<=127n?"int8":t>=-32768n&&t<=32767n?"int16":t>=-2147483648n&&t<=2147483647n?"int32":"int64":Number.isInteger(t)?t>=0?t<=255?"uint8":t<=65535?"uint16":t<=4294967295?"uint32":"int64":t>=-128&&t<=127?"int8":t>=-32768&&t<=32767?"int16":t>=-2147483648&&t<=2147483647?"int32":"int64":"float64"}function S(t,n){const e=t instanceof m?t.dtype:t;if(typeof n=="string")switch(n){case"number":case"numeric":return R.includes(e);case"integer":case"int":return d.includes(e);case"signedinteger":return y.includes(e);case"unsignedinteger":return D.includes(e);case"floating":case"float":return T.includes(e);case"complexfloating":case"complex":return x.includes(e);case"inexact":return T.includes(e)||x.includes(e);default:return e===n}return e===n}function k(t){return t}function A(t,n="GDFgdf",e="d"){const r={b:"int8",B:"uint8",h:"int16",H:"uint16",i:"int32",I:"uint32",l:"int64",L:"uint64",f:"float32",d:"float64",F:"complex64",D:"complex128",g:"float64",G:"complex128"},u={b:0,B:1,h:2,H:3,i:4,I:5,l:6,L:7,f:8,g:9.1,d:9,F:10,G:11.1,D:11};if(t.length===0)return e;const o=[];for(const i of t){const f=r[i];f&&o.push(f)}if(o.length===0)return e;let c=e,s=1/0;for(const i of n){const f=r[i];if(!f)continue;const b=u[i]??1/0;let g=!0;for(const E of o)if(!C(E,f,"safe")){g=!1;break}g&&b<s&&(s=b,c=i)}return c}export{C as can_cast,B as common_type,S as issubdtype,F as min_scalar_type,A as mintypecode,_ as result_type,k as typename};
+import{NDArrayCore as m}from"../common/ndarray-core.js";import{isComplexDType as l,isFloatDType as a,isIntegerDType as p}from"../common/dtype.js";const y=["int8","int16","int32","int64"],D=["uint8","uint16","uint32","uint64"],d=[...y,...D],T=["float16","float32","float64"],x=["complex64","complex128"],I=[...d,...T,...x],h={bool:0,int8:1,uint8:2,int16:3,uint16:4,int32:5,uint32:6,int64:7,uint64:8,float16:9,float32:10,float64:11,complex64:12,complex128:13};function R(t,n,e="safe"){const o=t instanceof m?t.dtype:t;if(e==="no"||e==="equiv")return o===n;if(e==="unsafe")return!0;const c=h[o]??-1,u=h[n]??-1;if(e==="safe"){if(p(o)&&a(n)){const r={int8:8,int16:16,int32:32,int64:64,uint8:8,uint16:16,uint32:32,uint64:64,bool:1,float16:0,float32:0,float64:0,complex64:0,complex128:0},s={float16:11,float32:24,float64:53,int8:0,int16:0,int32:0,int64:0,uint8:0,uint16:0,uint32:0,uint64:0,bool:0,complex64:24,complex128:53},i=r[o]??64,f=s[n]??0;return i<=f}return!!(p(o)&&l(n)||a(o)&&l(n)||c<=u)}return e==="same_kind"?!!(p(o)&&p(n)||a(o)&&a(n)||l(o)&&l(n)||o==="bool"&&p(n)||o==="bool"&&n==="bool"||a(o)&&l(n)):!1}function N(...t){if(t.length===0)return"float64";let n=!1,e=!1,o=!1;for(const u of t){const r=u.dtype;l(r)?(n=!0,r==="complex128"&&(e=!0)):a(r)?r==="float64"?e=!0:r==="float32"&&(o=!0):p(r)&&(r==="int64"||r==="uint64"||r==="int32"||r==="uint32")&&(e=!0)}return n?e?"complex128":"complex64":e?"float64":o?"float32":t.every(u=>u.dtype==="float16")?"float16":"float32"}function B(...t){if(t.length===0)return"float64";const n=t.map(u=>u instanceof m?u.dtype:u),e=n.some(l),o=n.some(a),c=n.some(p);if(e){const u=n.some(i=>i==="float64"),r=n.some(i=>i==="complex128"),s=n.some(i=>i==="int64"||i==="uint64"||i==="int32"||i==="uint32");return r||u||s?"complex128":"complex64"}if(o){const u=n.some(s=>s==="float64"),r=n.some(s=>s==="float32");if(c){if(n.some(f=>f==="int32"||f==="uint32"||f==="int64"||f==="uint64")||u)return"float64";const i=n.some(f=>f==="int16"||f==="uint16");if(!r&&i)return"float32"}return u?"float64":r?"float32":"float16"}if(c){const u=n.some(f=>y.includes(f)),r=n.some(f=>D.includes(f)),s={bool:8,int8:8,uint8:8,int16:16,uint16:16,int32:32,uint32:32,int64:64,uint64:64,float16:0,float32:0,float64:0,complex64:0,complex128:0},i=Math.max(...n.map(f=>s[f]??0));return u&&r?i>=32?"int64":i>=16?"int32":"int16":u?i>=64?"int64":i>=32?"int32":i>=16?"int16":"int8":i>=64?"uint64":i>=32?"uint32":i>=16?"uint16":"uint8"}return"bool"}function _(t){return typeof t=="boolean"?"bool":typeof t=="bigint"?t>=0n?t<=255n?"uint8":t<=65535n?"uint16":t<=4294967295n?"uint32":"uint64":t>=-128n&&t<=127n?"int8":t>=-32768n&&t<=32767n?"int16":t>=-2147483648n&&t<=2147483647n?"int32":"int64":Number.isInteger(t)?t>=0?t<=255?"uint8":t<=65535?"uint16":t<=4294967295?"uint32":"int64":t>=-128&&t<=127?"int8":t>=-32768&&t<=32767?"int16":t>=-2147483648&&t<=2147483647?"int32":"int64":"float64"}function S(t,n){const e=t instanceof m?t.dtype:t;if(typeof n=="string")switch(n){case"number":case"numeric":return I.includes(e);case"integer":case"int":return d.includes(e);case"signedinteger":return y.includes(e);case"unsignedinteger":return D.includes(e);case"floating":case"float":return T.includes(e);case"complexfloating":case"complex":return x.includes(e);case"inexact":return T.includes(e)||x.includes(e);default:return e===n}return e===n}function k(t){return t}function A(t,n="GDFgdf",e="d"){const o={b:"int8",B:"uint8",h:"int16",H:"uint16",i:"int32",I:"uint32",l:"int64",L:"uint64",f:"float32",d:"float64",F:"complex64",D:"complex128",g:"float64",G:"complex128"},c={b:0,B:1,h:2,H:3,i:4,I:5,l:6,L:7,f:8,g:9.1,d:9,F:10,G:11.1,D:11};if(t.length===0)return e;const u=[];for(const i of t){const f=o[i];f&&u.push(f)}if(u.length===0)return e;let r=e,s=1/0;for(const i of n){const f=o[i];if(!f)continue;const b=c[i]??1/0;let g=!0;for(const E of u)if(!R(E,f,"safe")){g=!1;break}g&&b<s&&(s=b,r=i)}return r}export{R as can_cast,N as common_type,S as issubdtype,_ as min_scalar_type,A as mintypecode,B as result_type,k as typename};

package/dist/esm/core/types.js

@@ -1 +1 @@
-import{NDArrayCore as e}from"../common/ndarray-core";import{ArrayStorage as n}from"../common/storage";import{Complex as y}from"../common/complex";function p(r){return r instanceof e||r&&typeof r=="object"&&"storage"in r&&r.storage instanceof n?r.storage:r}function a(r,o){return e.fromStorage(r,o)}function A(r,o){const t=o.base??o;return e.fromStorage(r,t)}function u(r){return r.map(o=>a(o))}function i(r,o){return r.map(t=>A(t,o))}function S(r){return[a(r[0]),a(r[1])]}export{n as ArrayStorage,y as Complex,e as NDArrayCore,a as fromStorage,u as fromStorageArray,S as fromStorageTuple,A as fromStorageView,i as fromStorageViewArray,p as toStorage};
+import{NDArrayCore as e}from"../common/ndarray-core.js";import{ArrayStorage as n}from"../common/storage.js";import{Complex as y}from"../common/complex.js";function p(r){return r instanceof e||r&&typeof r=="object"&&"storage"in r&&r.storage instanceof n?r.storage:r}function a(r,o){return e.fromStorage(r,o)}function A(r,o){const t=o.base??o;return e.fromStorage(r,t)}function u(r){return r.map(o=>a(o))}function i(r,o){return r.map(t=>A(t,o))}function S(r){return[a(r[0]),a(r[1])]}export{n as ArrayStorage,y as Complex,e as NDArrayCore,a as fromStorage,u as fromStorageArray,S as fromStorageTuple,A as fromStorageView,i as fromStorageViewArray,p as toStorage};

package/dist/esm/core/utility.js

@@ -1 +1 @@
-import{NDArrayCore as g}from"../common/ndarray-core";import{Complex as C}from"../common/complex";function E(r){if(r instanceof g)return r.ndim;if(typeof r=="number"||typeof r=="bigint"||typeof r=="boolean")return 0;if(Array.isArray(r)){let t=0,e=r;for(;Array.isArray(e)&&e.length>0;)t++,e=e[0];return t}return 0}function d(r){if(r instanceof g)return r.shape;if(typeof r=="number"||typeof r=="bigint"||typeof r=="boolean")return[];if(Array.isArray(r)){const t=[];let e=r;for(;Array.isArray(e)&&e.length>0;)t.push(e.length),e=e[0];return t}return[]}function D(r){return r instanceof g?r.size:typeof r=="number"||typeof r=="bigint"||typeof r=="boolean"?1:Array.isArray(r)?d(r).reduce((e,o)=>e*o,1):1}function N(r,...t){const e=r.storage,o=r.shape;if(t.length===0){if(r.size!==1)throw new Error("can only convert an array of size 1 to a scalar");return e.isCContiguous?e.data[e.offset]:e.iget(0)}if(t.length===1)return e.isCContiguous?e.data[e.offset+t[0]]:e.iget(t[0]);if(t.length!==o.length)throw new Error("incorrect number of indices for array");return e.get(...t)}function k(r){const t=r.shape,e=r.storage,o=t.length;if(o===0)return e.iget(0);if(e.isCContiguous){let b=function(l,c){const y=t[c],w=s[c];if(c===o-1){const p=new Array(y);for(let m=0;m<y;m++)p[m]=n[i+l+m];return p}const A=new Array(y);for(let p=0;p<y;p++)A[p]=b(l+p*w,c+1);return A};var u=b;const n=e.data,i=e.offset;if(o===1){const l=t[0],c=new Array(l);for(let y=0;y<l;y++)c[y]=n[i+y];return c}const s=new Array(o);let a=1;for(let l=o-1;l>=0;l--)s[l]=a,a*=t[l];return b(0,0)}if(o===1){const n=[];for(let i=0;i<t[0];i++)n.push(e.iget(i));return n}function f(n,i){if(i===o)return e.get(...n);const s=[];for(let a=0;a<t[i];a++)n[i]=a,s.push(f(n,i+1));return s}return f(new Array(o),0)}function z(r,t="C"){const e=r.storage;t==="F"&&console.warn('tobytes with order="F" not fully implemented, returning C-order');const o=e.data,f=o.BYTES_PER_ELEMENT;if(e.isCContiguous){const i=o.byteOffset+e.offset*f,s=r.size*f;return new Uint8Array(o.buffer,i,s)}const n=r.copy().data;return new Uint8Array(n.buffer,n.byteOffset,r.size*f)}function T(r,t=!1){const o=r.data.BYTES_PER_ELEMENT;if(o===1)return t?r:r.copy();const f=t?r:r.copy(),u=f.data,n=new Uint8Array(u.buffer,u.byteOffset,u.byteLength);for(let i=0;i<u.length;i++){const s=i*o;for(let a=0;a<o/2;a++){const b=n[s+a];n[s+a]=n[s+o-1-a],n[s+o-1-a]=b}}return f}function _(r,t){if(!t||t===r.dtype)return r.copy();throw new Error("view with different dtype not fully implemented")}function U(r,t,e="",o=""){throw new Error("tofile requires Node.js file system access. Use serializeNpy for portable serialization.")}function B(r,t){const e=r.storage,o=r.dtype,f=r.size;if(t instanceof C){if(!(o==="complex64"||o==="complex128"))throw new Error("Cannot fill non-complex array with complex value");if(e.isCContiguous){const n=e.data,i=e.offset;for(let s=0;s<f;s++)n[(i+s)*2]=t.re,n[(i+s)*2+1]=t.im}else for(let n=0;n<f;n++)e.iset(n,t)}else if(typeof t=="bigint")if(e.isCContiguous)e.data.fill(t,e.offset,e.offset+f);else for(let u=0;u<f;u++)e.iset(u,t);else{const u=typeof t=="boolean"?t?1:0:t;if(e.isCContiguous)e.data.fill(u,e.offset,e.offset+f);else for(let n=0;n<f;n++)e.iset(n,u)}}export{T as byteswap,B as fill,N as item,E as ndim,d as shape,D as size,z as tobytes,U as tofile,k as tolist,_ as view};
+import{NDArrayCore as g}from"../common/ndarray-core.js";import{Complex as C}from"../common/complex.js";function E(r){if(r instanceof g)return r.ndim;if(typeof r=="number"||typeof r=="bigint"||typeof r=="boolean")return 0;if(Array.isArray(r)){let t=0,e=r;for(;Array.isArray(e)&&e.length>0;)t++,e=e[0];return t}return 0}function d(r){if(r instanceof g)return r.shape;if(typeof r=="number"||typeof r=="bigint"||typeof r=="boolean")return[];if(Array.isArray(r)){const t=[];let e=r;for(;Array.isArray(e)&&e.length>0;)t.push(e.length),e=e[0];return t}return[]}function D(r){return r instanceof g?r.size:typeof r=="number"||typeof r=="bigint"||typeof r=="boolean"?1:Array.isArray(r)?d(r).reduce((e,o)=>e*o,1):1}function N(r,...t){const e=r.storage,o=r.shape;if(t.length===0){if(r.size!==1)throw new Error("can only convert an array of size 1 to a scalar");return e.isCContiguous?e.data[e.offset]:e.iget(0)}if(t.length===1)return e.isCContiguous?e.data[e.offset+t[0]]:e.iget(t[0]);if(t.length!==o.length)throw new Error("incorrect number of indices for array");return e.get(...t)}function k(r){const t=r.shape,e=r.storage,o=t.length;if(o===0)return e.iget(0);if(e.isCContiguous){let b=function(l,c){const y=t[c],w=s[c];if(c===o-1){const p=new Array(y);for(let m=0;m<y;m++)p[m]=n[i+l+m];return p}const A=new Array(y);for(let p=0;p<y;p++)A[p]=b(l+p*w,c+1);return A};var u=b;const n=e.data,i=e.offset;if(o===1){const l=t[0],c=new Array(l);for(let y=0;y<l;y++)c[y]=n[i+y];return c}const s=new Array(o);let a=1;for(let l=o-1;l>=0;l--)s[l]=a,a*=t[l];return b(0,0)}if(o===1){const n=[];for(let i=0;i<t[0];i++)n.push(e.iget(i));return n}function f(n,i){if(i===o)return e.get(...n);const s=[];for(let a=0;a<t[i];a++)n[i]=a,s.push(f(n,i+1));return s}return f(new Array(o),0)}function z(r,t="C"){const e=r.storage;t==="F"&&console.warn('tobytes with order="F" not fully implemented, returning C-order');const o=e.data,f=o.BYTES_PER_ELEMENT;if(e.isCContiguous){const i=o.byteOffset+e.offset*f,s=r.size*f;return new Uint8Array(o.buffer,i,s)}const n=r.copy().data;return new Uint8Array(n.buffer,n.byteOffset,r.size*f)}function T(r,t=!1){const o=r.data.BYTES_PER_ELEMENT;if(o===1)return t?r:r.copy();const f=t?r:r.copy(),u=f.data,n=new Uint8Array(u.buffer,u.byteOffset,u.byteLength);for(let i=0;i<u.length;i++){const s=i*o;for(let a=0;a<o/2;a++){const b=n[s+a];n[s+a]=n[s+o-1-a],n[s+o-1-a]=b}}return f}function _(r,t){if(!t||t===r.dtype)return r.copy();throw new Error("view with different dtype not fully implemented")}function U(r,t,e="",o=""){throw new Error("tofile requires Node.js file system access. Use serializeNpy for portable serialization.")}function B(r,t){const e=r.storage,o=r.dtype,f=r.size;if(t instanceof C){if(!(o==="complex64"||o==="complex128"))throw new Error("Cannot fill non-complex array with complex value");if(e.isCContiguous){const n=e.data,i=e.offset;for(let s=0;s<f;s++)n[(i+s)*2]=t.re,n[(i+s)*2+1]=t.im}else for(let n=0;n<f;n++)e.iset(n,t)}else if(typeof t=="bigint")if(e.isCContiguous)e.data.fill(t,e.offset,e.offset+f);else for(let u=0;u<f;u++)e.iset(u,t);else{const u=typeof t=="boolean"?t?1:0:t;if(e.isCContiguous)e.data.fill(u,e.offset,e.offset+f);else for(let n=0;n<f;n++)e.iset(n,u)}}export{T as byteswap,B as fill,N as item,E as ndim,d as shape,D as size,z as tobytes,U as tofile,k as tolist,_ as view};

package/dist/esm/core.js

@@ -1 +1 @@
-export*from"./core/index";
+export*from"./core/index.js";

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

package/dist/esm/full/ndarray.js

@@ -1 +1 @@
-import{getTypedArrayConstructor as v,getDTypeSize as g,isBigIntDType as p,isComplexDType as z}from"../common/dtype";import{Complex as D}from"../common/complex";import{ArrayStorage as d}from"../common/storage";import{NDArrayCore as m}from"../common/ndarray-core";import*as t from"../core";const s=l=>{if(l instanceof b)return l;const r=l.base?s(l.base):void 0;return b.fromStorage(l.storage,r)};class b extends m{constructor(r,e){super(r,e),this._base=e}static fromStorage(r,e){return new b(r,e)}get base(){return this._base??null}get T(){return this.transpose()}get itemsize(){return g(this._storage.dtype)}get nbytes(){return this.size*this.itemsize}fill(r){const e=this._storage.dtype,n=this.size;if(p(e)){const a=typeof r=="bigint"?r:BigInt(Math.round(Number(r)));for(let o=0;o<n;o++)this._storage.iset(o,a)}else if(e==="bool"){const a=r?1:0;for(let o=0;o<n;o++)this._storage.iset(o,a)}else{const a=Number(r);for(let o=0;o<n;o++)this._storage.iset(o,a)}}*[Symbol.iterator](){if(this.ndim===0)yield this._storage.iget(0);else if(this.ndim===1)for(let r=0;r<this.shape[0];r++)yield this._storage.iget(r);else for(let r=0;r<this.shape[0];r++)yield this.slice(String(r))}get(r){if(r.length!==this.ndim)throw new Error(`Index has ${r.length} dimensions, but array has ${this.ndim} dimensions`);const e=r.map((n,a)=>{let o=n;if(o<0&&(o=this.shape[a]+o),o<0||o>=this.shape[a])throw new Error(`Index ${n} is out of bounds for axis ${a} with size ${this.shape[a]}`);return o});return this._storage.get(...e)}set(r,e){if(r.length!==this.ndim)throw new Error(`Index has ${r.length} dimensions, but array has ${this.ndim} dimensions`);const n=r.map((o,i)=>{let h=o;if(h<0&&(h=this.shape[i]+h),h<0||h>=this.shape[i])throw new Error(`Index ${o} is out of bounds for axis ${i} with size ${this.shape[i]}`);return h}),a=this.dtype;if(z(a))this._storage.set(n,e);else if(p(a)){const o=e instanceof D?e.re:Number(e),i=typeof e=="bigint"?e:BigInt(Math.round(o));this._storage.set(n,i)}else if(a==="bool"){const i=(e instanceof D?e.re:Number(e))?1:0;this._storage.set(n,i)}else{const o=e instanceof D?e.re:Number(e);this._storage.set(n,o)}}copy(){return new b(this._storage.copy())}astype(r,e=!0){const n=super.astype(r,e);return n instanceof b?n:new b(n._storage)}row(r){if(this.ndim<2)throw new Error("row() requires at least 2 dimensions");return this.slice(String(r),":")}col(r){if(this.ndim<2)throw new Error("col() requires at least 2 dimensions");return this.slice(":",String(r))}rows(r,e){if(this.ndim<2)throw new Error("rows() requires at least 2 dimensions");return this.slice(`${r}:${e}`,":")}cols(r,e){if(this.ndim<2)throw new Error("cols() requires at least 2 dimensions");return this.slice(":",`${r}:${e}`)}reshape(...r){const e=r.length===1&&Array.isArray(r[0])?r[0]:r,n=t.reshape(this,e).storage,o=n.data===this.data?this._base??this:void 0;return b.fromStorage(n,o)}ravel(){const r=t.ravel(this).storage,n=r.data===this.data?this._base??this:void 0;return b.fromStorage(r,n)}put(r,e){const n=e instanceof b?e._storage:e;t.put(this,r,n)}compress(r,e){const n=r instanceof b?r:b.fromStorage(d.fromData(new Uint8Array(r.map(a=>a?1:0)),[r.length],"bool"));return s(t.compress(n,this,e))}choose(r){return s(t.choose(this,r))}clip(r,e){return s(t.clip(this,r,e))}iindex(r,e=0){let n;if(r instanceof b){n=[];for(let a=0;a<r.size;a++){const o=r.storage.iget(a),i=typeof o=="bigint"?Number(o):o instanceof D?o.re:o;n.push(i)}}else Array.isArray(r)&&r.length>0&&Array.isArray(r[0])?n=r.flat():n=r;return this.take(n,e)}bindex(r,e){return s(t.compress(r,this,e))}toString(){return t.array_str(this)}toArray(){if(this.ndim===0)return this._storage.iget(0);const r=this.shape,e=r.length,n=(a,o)=>{if(o===e)return this._storage.get(...a);const i=[];for(let h=0;h<r[o];h++)a[o]=h,i.push(n(a,o+1));return i};return n(new Array(e),0)}tolist(){return this.toArray()}tobytes(){if(this._storage.isCContiguous){const n=this._storage.data,a=n.BYTES_PER_ELEMENT,o=this._storage.offset*a,i=this.size*a;return n.buffer.slice(o,o+i)}const e=this.copy()._storage.data;return e.buffer.slice(0,this.size*e.BYTES_PER_ELEMENT)}item(...r){if(r.length===0){if(this.size!==1)throw new Error("can only convert an array of size 1 to a Python scalar");return this._storage.iget(0)}if(r.length===1){const e=r[0];if(e<0||e>=this.size)throw new Error(`index ${e} is out of bounds for size ${this.size}`);return this._storage.iget(e)}return this.get(r)}byteswap(r=!1){const e=r?this:this.copy(),n=e._storage.data,a=n.BYTES_PER_ELEMENT;if(a===1)return e;const o=n.buffer,i=new DataView(o);for(let h=0;h<n.length;h++){const u=h*a;if(a===2){const c=i.getUint8(u),y=i.getUint8(u+1);i.setUint8(u,y),i.setUint8(u+1,c)}else if(a===4){const c=i.getUint8(u),y=i.getUint8(u+1),A=i.getUint8(u+2),f=i.getUint8(u+3);i.setUint8(u,f),i.setUint8(u+1,A),i.setUint8(u+2,y),i.setUint8(u+3,c)}else if(a===8){const c=i.getUint8(u),y=i.getUint8(u+1),A=i.getUint8(u+2),f=i.getUint8(u+3),N=i.getUint8(u+4),x=i.getUint8(u+5),w=i.getUint8(u+6),_=i.getUint8(u+7);i.setUint8(u,_),i.setUint8(u+1,w),i.setUint8(u+2,x),i.setUint8(u+3,N),i.setUint8(u+4,f),i.setUint8(u+5,A),i.setUint8(u+6,y),i.setUint8(u+7,c)}}return e}view(r){if(!r||r===this.dtype)return b.fromStorage(this._storage,this._base??this);const e=g(this.dtype),n=g(r);if(e!==n)throw new Error("When changing to a larger dtype, its size must be a divisor of the total size in bytes of the last axis of the array.");const a=v(r);if(!a)throw new Error(`Unsupported dtype: ${r}`);const o=this._storage.data,i=o.byteOffset+this._storage.offset*e,h=new a(o.buffer,i,this.size),u=d.fromData(h,[...this.shape],r,[...this._storage.strides],0);return b.fromStorage(u,this._base??this)}tofile(r,e="",n=""){throw new Error('tofile() requires file system access. Use the node module: import { save } from "numpy-ts/node"')}round(r=0){return this.around(r)}conjugate(){return this.conj()}around(r=0){return s(t.around(this,r))}allclose(r,e=1e-5,n=1e-8){return t.allclose(this,r,e,n)}isclose(r,e=1e-5,n=1e-8){return s(t.isclose(this,r,e,n))}average(r,e){const n=t.average(this,e,r);return n instanceof m?s(n):n}dot(r){const e=t.dot(this,r);return e instanceof m?s(e):e}trace(){const r=t.trace(this);return r instanceof m?s(r):r}inner(r){const e=t.inner(this,r);return e instanceof m?s(e):e}tensordot(r,e=2){const n=t.tensordot(this,r,e);return n instanceof m?s(n):n}divmod(r){const e=t.divmod(this,r);return[s(e[0]),s(e[1])]}searchsorted(r,e="left"){return s(t.searchsorted(this,r,e))}sqrt(){return s(t.sqrt(this))}exp(){return s(t.exp(this))}exp2(){return s(t.exp2(this))}expm1(){return s(t.expm1(this))}log(){return s(t.log(this))}log2(){return s(t.log2(this))}log10(){return s(t.log10(this))}log1p(){return s(t.log1p(this))}absolute(){return s(t.absolute(this))}negative(){return s(t.negative(this))}sign(){return s(t.sign(this))}positive(){return s(t.positive(this))}reciprocal(){return s(t.reciprocal(this))}ceil(){return s(t.ceil(this))}fix(){return s(t.fix(this))}floor(){return s(t.floor(this))}rint(){return s(t.rint(this))}trunc(){return s(t.trunc(this))}sin(){return s(t.sin(this))}cos(){return s(t.cos(this))}tan(){return s(t.tan(this))}arcsin(){return s(t.arcsin(this))}arccos(){return s(t.arccos(this))}arctan(){return s(t.arctan(this))}degrees(){return s(t.degrees(this))}radians(){return s(t.radians(this))}sinh(){return s(t.sinh(this))}cosh(){return s(t.cosh(this))}tanh(){return s(t.tanh(this))}arcsinh(){return s(t.arcsinh(this))}arccosh(){return s(t.arccosh(this))}arctanh(){return s(t.arctanh(this))}bitwise_not(){return s(t.bitwise_not(this))}invert(){return s(t.invert(this))}logical_not(){return s(t.logical_not(this))}isfinite(){return s(t.isfinite(this))}isinf(){return s(t.isinf(this))}isnan(){return s(t.isnan(this))}isnat(){return s(t.isnat(this))}signbit(){return s(t.signbit(this))}spacing(){return s(t.spacing(this))}cbrt(){return s(t.cbrt(this))}fabs(){return s(t.fabs(this))}square(){return s(t.square(this))}conj(){return s(t.conj(this))}flatten(){return s(t.flatten(this))}argwhere(){return s(t.argwhere(this))}add(r){return s(t.add(this,r))}subtract(r){return s(t.subtract(this,r))}multiply(r){return s(t.multiply(this,r))}divide(r){return s(t.divide(this,r))}mod(r){return s(t.mod(this,r))}floor_divide(r){return s(t.floor_divide(this,r))}power(r){return s(t.power(this,r))}logaddexp(r){return s(t.logaddexp(this,r))}logaddexp2(r){return s(t.logaddexp2(this,r))}arctan2(r){return s(t.arctan2(this,r))}hypot(r){return s(t.hypot(this,r))}greater(r){return s(t.greater(this,r))}greater_equal(r){return s(t.greater_equal(this,r))}less(r){return s(t.less(this,r))}less_equal(r){return s(t.less_equal(this,r))}equal(r){return s(t.equal(this,r))}not_equal(r){return s(t.not_equal(this,r))}bitwise_and(r){return s(t.bitwise_and(this,r))}bitwise_or(r){return s(t.bitwise_or(this,r))}bitwise_xor(r){return s(t.bitwise_xor(this,r))}left_shift(r){return s(t.left_shift(this,r))}right_shift(r){return s(t.right_shift(this,r))}logical_and(r){return s(t.logical_and(this,r))}logical_or(r){return s(t.logical_or(this,r))}logical_xor(r){return s(t.logical_xor(this,r))}copysign(r){return s(t.copysign(this,r))}nextafter(r){return s(t.nextafter(this,r))}remainder(r){return s(t.remainder(this,r))}heaviside(r){return s(t.heaviside(this,r))}matmul(r){return s(t.matmul(this,r))}outer(r){return s(t.outer(this,r))}sum(r,e=!1){const n=t.sum(this,r,e);return n instanceof m?s(n):n}mean(r,e=!1){const n=t.mean(this,r,e);return n instanceof m?s(n):n}prod(r,e=!1){const n=t.prod(this,r,e);return n instanceof m?s(n):n}max(r,e=!1){const n=t.max(this,r,e);return n instanceof m?s(n):n}min(r,e=!1){const n=t.min(this,r,e);return n instanceof m?s(n):n}ptp(r,e=!1){const n=t.ptp(this,r,e);return n instanceof m?s(n):n}nansum(r,e=!1){const n=t.nansum(this,r,e);return n instanceof m?s(n):n}nanprod(r,e=!1){const n=t.nanprod(this,r,e);return n instanceof m?s(n):n}nanmean(r,e=!1){const n=t.nanmean(this,r,e);return n instanceof m?s(n):n}nanmin(r,e=!1){const n=t.nanmin(this,r,e);return n instanceof m?s(n):n}nanmax(r,e=!1){const n=t.nanmax(this,r,e);return n instanceof m?s(n):n}argmin(r){const e=t.argmin(this,r);return e instanceof m?s(e):e}argmax(r){const e=t.argmax(this,r);return e instanceof m?s(e):e}nanargmin(r){const e=t.nanargmin(this,r);return e instanceof m?s(e):e}nanargmax(r){const e=t.nanargmax(this,r);return e instanceof m?s(e):e}var(r,e=0,n=!1){const a=t.variance(this,r,e,n);return a instanceof m?s(a):a}std(r,e=0,n=!1){const a=t.std(this,r,e,n);return a instanceof m?s(a):a}nanvar(r,e=0,n=!1){const a=t.nanvar(this,r,e,n);return a instanceof m?s(a):a}nanstd(r,e=0,n=!1){const a=t.nanstd(this,r,e,n);return a instanceof m?s(a):a}all(r,e=!1){const n=t.all(this,r,e);return n instanceof m?s(n):n}any(r,e=!1){const n=t.any(this,r,e);return n instanceof m?s(n):n}median(r,e=!1){const n=t.median(this,r,e);return n instanceof m?s(n):n}nanmedian(r,e=!1){const n=t.nanmedian(this,r,e);return n instanceof m?s(n):n}percentile(r,e,n=!1){const a=t.percentile(this,r,e,n);return a instanceof m?s(a):a}quantile(r,e,n=!1){const a=t.quantile(this,r,e,n);return a instanceof m?s(a):a}nanquantile(r,e,n=!1){const a=t.nanquantile(this,r,e,n);return a instanceof m?s(a):a}nanpercentile(r,e,n=!1){const a=t.nanpercentile(this,r,e,n);return a instanceof m?s(a):a}cumsum(r){return s(t.cumsum(this,r))}cumprod(r){return s(t.cumprod(this,r))}nancumsum(r){return s(t.nancumsum(this,r))}nancumprod(r){return s(t.nancumprod(this,r))}sort(r=-1){return s(t.sort(this,r))}argsort(r=-1){return s(t.argsort(this,r))}partition(r,e=-1){return s(t.partition(this,r,e))}argpartition(r,e=-1){return s(t.argpartition(this,r,e))}diagonal(r=0,e=0,n=1){return s(t.diagonal(this,r,e,n))}resize(r){return s(t.resize(this,r))}diff(r=1,e=-1){return s(t.diff(this,r,e))}take(r,e){return s(t.take(this,r,e))}repeat(r,e){return s(t.repeat(this,r,e))}transpose(r){return s(t.transpose(this,r))}squeeze(r){return s(t.squeeze(this,r))}expand_dims(r){return s(t.expand_dims(this,r))}swapaxes(r,e){return s(t.swapaxes(this,r,e))}moveaxis(r,e){return s(t.moveaxis(this,r,e))}nonzero(){return t.nonzero(this).map(s)}}function q(...l){let r=[],e="xy";for(const h of l)h instanceof b?r.push(h):typeof h=="object"&&"indexing"in h&&(e=h.indexing||"xy");if(r.length===0)return[];if(r.length===1)return[r[0].copy()];const n=r.map(h=>h.size);e==="xy"&&r.length>=2&&(r=[r[1],r[0],...r.slice(2)],[n[0],n[1]]=[n[1],n[0]]);const a=n,o=a.length,i=[];for(let h=0;h<r.length;h++){const u=r[h],c=u.size,y=new Array(o).fill(1);y[h]=c;const A=u.reshape(...y),f=t.broadcast_to(A,a),N=b.fromStorage(f.storage.copy());i.push(N)}return e==="xy"&&i.length>=2&&([i[0],i[1]]=[i[1],i[0]]),i}export{b as NDArray,q as meshgrid};
+import{getTypedArrayConstructor as z,getDTypeSize as g,isBigIntDType as p,isComplexDType as S}from"../common/dtype.js";import{Complex as N}from"../common/complex.js";import{ArrayStorage as d}from"../common/storage.js";import{NDArrayCore as u}from"../common/ndarray-core.js";import*as t from"../core/index.js";const s=c=>{if(c instanceof b)return c;const r=c.base?s(c.base):void 0;return b.fromStorage(c.storage,r)};class b extends u{constructor(r,e){super(r,e),this._base=e}static fromStorage(r,e){return new b(r,e)}get base(){return this._base??null}get T(){return this.transpose()}get itemsize(){return g(this._storage.dtype)}get nbytes(){return this.size*this.itemsize}fill(r){const e=this._storage.dtype,n=this.size;if(p(e)){const a=typeof r=="bigint"?r:BigInt(Math.round(Number(r)));for(let o=0;o<n;o++)this._storage.iset(o,a)}else if(e==="bool"){const a=r?1:0;for(let o=0;o<n;o++)this._storage.iset(o,a)}else{const a=Number(r);for(let o=0;o<n;o++)this._storage.iset(o,a)}}*[Symbol.iterator](){if(this.ndim===0)yield this._storage.iget(0);else if(this.ndim===1)for(let r=0;r<this.shape[0];r++)yield this._storage.iget(r);else for(let r=0;r<this.shape[0];r++)yield this.slice(String(r))}get(r){if(r.length!==this.ndim)throw new Error(`Index has ${r.length} dimensions, but array has ${this.ndim} dimensions`);const e=r.map((n,a)=>{let o=n;if(o<0&&(o=this.shape[a]+o),o<0||o>=this.shape[a])throw new Error(`Index ${n} is out of bounds for axis ${a} with size ${this.shape[a]}`);return o});return this._storage.get(...e)}set(r,e){if(r.length!==this.ndim)throw new Error(`Index has ${r.length} dimensions, but array has ${this.ndim} dimensions`);const n=r.map((o,h)=>{let i=o;if(i<0&&(i=this.shape[h]+i),i<0||i>=this.shape[h])throw new Error(`Index ${o} is out of bounds for axis ${h} with size ${this.shape[h]}`);return i}),a=this.dtype;if(S(a))this._storage.set(n,e);else if(p(a)){const o=e instanceof N?e.re:Number(e),h=typeof e=="bigint"?e:BigInt(Math.round(o));this._storage.set(n,h)}else if(a==="bool"){const h=(e instanceof N?e.re:Number(e))?1:0;this._storage.set(n,h)}else{const o=e instanceof N?e.re:Number(e);this._storage.set(n,o)}}copy(){return new b(this._storage.copy())}astype(r,e=!0){const n=super.astype(r,e);return n instanceof b?n:new b(n._storage)}row(r){if(this.ndim<2)throw new Error("row() requires at least 2 dimensions");return this.slice(String(r),":")}col(r){if(this.ndim<2)throw new Error("col() requires at least 2 dimensions");return this.slice(":",String(r))}rows(r,e){if(this.ndim<2)throw new Error("rows() requires at least 2 dimensions");return this.slice(`${r}:${e}`,":")}cols(r,e){if(this.ndim<2)throw new Error("cols() requires at least 2 dimensions");return this.slice(":",`${r}:${e}`)}reshape(...r){const e=r.length===1&&Array.isArray(r[0])?r[0]:r,n=t.reshape(this,e).storage,o=n.data===this.data?this._base??this:void 0;return b.fromStorage(n,o)}ravel(){const r=t.ravel(this).storage,n=r.data===this.data?this._base??this:void 0;return b.fromStorage(r,n)}put(r,e){const n=e instanceof b?e._storage:e;t.put(this,r,n)}compress(r,e){const n=r instanceof b?r:b.fromStorage(d.fromData(new Uint8Array(r.map(a=>a?1:0)),[r.length],"bool"));return s(t.compress(n,this,e))}choose(r){return s(t.choose(this,r))}clip(r,e){return s(t.clip(this,r,e))}iindex(r,e=0){let n;if(r instanceof b){n=[];for(let a=0;a<r.size;a++){const o=r.storage.iget(a),h=typeof o=="bigint"?Number(o):o instanceof N?o.re:o;n.push(h)}}else Array.isArray(r)&&r.length>0&&Array.isArray(r[0])?n=r.flat():n=r;return this.take(n,e)}bindex(r,e){return s(t.compress(r,this,e))}toString(){return t.array_str(this)}toArray(){if(this.ndim===0)return this._storage.iget(0);const r=this.shape,e=r.length,n=(a,o)=>{if(o===e)return this._storage.get(...a);const h=[];for(let i=0;i<r[o];i++)a[o]=i,h.push(n(a,o+1));return h};return n(new Array(e),0)}tolist(){return this.toArray()}tobytes(){if(this._storage.isCContiguous){const n=this._storage.data,a=n.BYTES_PER_ELEMENT,o=n.byteOffset+this._storage.offset*a,h=this.size*a;return n.buffer.slice(o,o+h)}const e=this.copy()._storage.data;return e.buffer.slice(e.byteOffset,e.byteOffset+this.size*e.BYTES_PER_ELEMENT)}item(...r){if(r.length===0){if(this.size!==1)throw new Error("can only convert an array of size 1 to a Python scalar");return this._storage.iget(0)}if(r.length===1){const e=r[0];if(e<0||e>=this.size)throw new Error(`index ${e} is out of bounds for size ${this.size}`);return this._storage.iget(e)}return this.get(r)}byteswap(r=!1){const e=r?this:this.copy(),n=e._storage.data,a=n.BYTES_PER_ELEMENT;if(a===1)return e;const o=n.buffer,h=n.byteOffset,i=new DataView(o);for(let y=0;y<n.length;y++){const m=h+y*a;if(a===2){const l=i.getUint8(m),f=i.getUint8(m+1);i.setUint8(m,f),i.setUint8(m+1,l)}else if(a===4){const l=i.getUint8(m),f=i.getUint8(m+1),A=i.getUint8(m+2),D=i.getUint8(m+3);i.setUint8(m,D),i.setUint8(m+1,A),i.setUint8(m+2,f),i.setUint8(m+3,l)}else if(a===8){const l=i.getUint8(m),f=i.getUint8(m+1),A=i.getUint8(m+2),D=i.getUint8(m+3),x=i.getUint8(m+4),w=i.getUint8(m+5),_=i.getUint8(m+6),v=i.getUint8(m+7);i.setUint8(m,v),i.setUint8(m+1,_),i.setUint8(m+2,w),i.setUint8(m+3,x),i.setUint8(m+4,D),i.setUint8(m+5,A),i.setUint8(m+6,f),i.setUint8(m+7,l)}}return e}view(r){if(!r||r===this.dtype)return b.fromStorage(this._storage,this._base??this);const e=g(this.dtype),n=g(r);if(e!==n)throw new Error("When changing to a larger dtype, its size must be a divisor of the total size in bytes of the last axis of the array.");const a=z(r);if(!a)throw new Error(`Unsupported dtype: ${r}`);const o=this._storage.data,h=o.byteOffset+this._storage.offset*e,i=new a(o.buffer,h,this.size),y=d.fromData(i,[...this.shape],r,[...this._storage.strides],0);return b.fromStorage(y,this._base??this)}tofile(r,e="",n=""){throw new Error('tofile() requires file system access. Use the node module: import { save } from "numpy-ts/node"')}round(r=0){return this.around(r)}conjugate(){return this.conj()}around(r=0){return s(t.around(this,r))}allclose(r,e=1e-5,n=1e-8){return t.allclose(this,r,e,n)}isclose(r,e=1e-5,n=1e-8){return s(t.isclose(this,r,e,n))}average(r,e){const n=t.average(this,e,r);return n instanceof u?s(n):n}dot(r){const e=t.dot(this,r);return e instanceof u?s(e):e}trace(){const r=t.trace(this);return r instanceof u?s(r):r}inner(r){const e=t.inner(this,r);return e instanceof u?s(e):e}tensordot(r,e=2){const n=t.tensordot(this,r,e);return n instanceof u?s(n):n}divmod(r){const e=t.divmod(this,r);return[s(e[0]),s(e[1])]}searchsorted(r,e="left"){return s(t.searchsorted(this,r,e))}sqrt(){return s(t.sqrt(this))}exp(){return s(t.exp(this))}exp2(){return s(t.exp2(this))}expm1(){return s(t.expm1(this))}log(){return s(t.log(this))}log2(){return s(t.log2(this))}log10(){return s(t.log10(this))}log1p(){return s(t.log1p(this))}absolute(){return s(t.absolute(this))}negative(){return s(t.negative(this))}sign(){return s(t.sign(this))}positive(){return s(t.positive(this))}reciprocal(){return s(t.reciprocal(this))}ceil(){return s(t.ceil(this))}fix(){return s(t.fix(this))}floor(){return s(t.floor(this))}rint(){return s(t.rint(this))}trunc(){return s(t.trunc(this))}sin(){return s(t.sin(this))}cos(){return s(t.cos(this))}tan(){return s(t.tan(this))}arcsin(){return s(t.arcsin(this))}arccos(){return s(t.arccos(this))}arctan(){return s(t.arctan(this))}degrees(){return s(t.degrees(this))}radians(){return s(t.radians(this))}sinh(){return s(t.sinh(this))}cosh(){return s(t.cosh(this))}tanh(){return s(t.tanh(this))}arcsinh(){return s(t.arcsinh(this))}arccosh(){return s(t.arccosh(this))}arctanh(){return s(t.arctanh(this))}bitwise_not(){return s(t.bitwise_not(this))}invert(){return s(t.invert(this))}logical_not(){return s(t.logical_not(this))}isfinite(){return s(t.isfinite(this))}isinf(){return s(t.isinf(this))}isnan(){return s(t.isnan(this))}isnat(){return s(t.isnat(this))}signbit(){return s(t.signbit(this))}spacing(){return s(t.spacing(this))}cbrt(){return s(t.cbrt(this))}fabs(){return s(t.fabs(this))}square(){return s(t.square(this))}conj(){return s(t.conj(this))}flatten(){return s(t.flatten(this))}argwhere(){return s(t.argwhere(this))}add(r){return s(t.add(this,r))}subtract(r){return s(t.subtract(this,r))}multiply(r){return s(t.multiply(this,r))}divide(r){return s(t.divide(this,r))}mod(r){return s(t.mod(this,r))}floor_divide(r){return s(t.floor_divide(this,r))}power(r){return s(t.power(this,r))}logaddexp(r){return s(t.logaddexp(this,r))}logaddexp2(r){return s(t.logaddexp2(this,r))}arctan2(r){return s(t.arctan2(this,r))}hypot(r){return s(t.hypot(this,r))}greater(r){return s(t.greater(this,r))}greater_equal(r){return s(t.greater_equal(this,r))}less(r){return s(t.less(this,r))}less_equal(r){return s(t.less_equal(this,r))}equal(r){return s(t.equal(this,r))}not_equal(r){return s(t.not_equal(this,r))}bitwise_and(r){return s(t.bitwise_and(this,r))}bitwise_or(r){return s(t.bitwise_or(this,r))}bitwise_xor(r){return s(t.bitwise_xor(this,r))}left_shift(r){return s(t.left_shift(this,r))}right_shift(r){return s(t.right_shift(this,r))}logical_and(r){return s(t.logical_and(this,r))}logical_or(r){return s(t.logical_or(this,r))}logical_xor(r){return s(t.logical_xor(this,r))}copysign(r){return s(t.copysign(this,r))}nextafter(r){return s(t.nextafter(this,r))}remainder(r){return s(t.remainder(this,r))}heaviside(r){return s(t.heaviside(this,r))}matmul(r){return s(t.matmul(this,r))}outer(r){return s(t.outer(this,r))}sum(r,e=!1){const n=t.sum(this,r,e);return n instanceof u?s(n):n}mean(r,e=!1){const n=t.mean(this,r,e);return n instanceof u?s(n):n}prod(r,e=!1){const n=t.prod(this,r,e);return n instanceof u?s(n):n}max(r,e=!1){const n=t.max(this,r,e);return n instanceof u?s(n):n}min(r,e=!1){const n=t.min(this,r,e);return n instanceof u?s(n):n}ptp(r,e=!1){const n=t.ptp(this,r,e);return n instanceof u?s(n):n}nansum(r,e=!1){const n=t.nansum(this,r,e);return n instanceof u?s(n):n}nanprod(r,e=!1){const n=t.nanprod(this,r,e);return n instanceof u?s(n):n}nanmean(r,e=!1){const n=t.nanmean(this,r,e);return n instanceof u?s(n):n}nanmin(r,e=!1){const n=t.nanmin(this,r,e);return n instanceof u?s(n):n}nanmax(r,e=!1){const n=t.nanmax(this,r,e);return n instanceof u?s(n):n}argmin(r){const e=t.argmin(this,r);return e instanceof u?s(e):e}argmax(r){const e=t.argmax(this,r);return e instanceof u?s(e):e}nanargmin(r){const e=t.nanargmin(this,r);return e instanceof u?s(e):e}nanargmax(r){const e=t.nanargmax(this,r);return e instanceof u?s(e):e}var(r,e=0,n=!1){const a=t.variance(this,r,e,n);return a instanceof u?s(a):a}std(r,e=0,n=!1){const a=t.std(this,r,e,n);return a instanceof u?s(a):a}nanvar(r,e=0,n=!1){const a=t.nanvar(this,r,e,n);return a instanceof u?s(a):a}nanstd(r,e=0,n=!1){const a=t.nanstd(this,r,e,n);return a instanceof u?s(a):a}all(r,e=!1){const n=t.all(this,r,e);return n instanceof u?s(n):n}any(r,e=!1){const n=t.any(this,r,e);return n instanceof u?s(n):n}median(r,e=!1){const n=t.median(this,r,e);return n instanceof u?s(n):n}nanmedian(r,e=!1){const n=t.nanmedian(this,r,e);return n instanceof u?s(n):n}percentile(r,e,n=!1){const a=t.percentile(this,r,e,n);return a instanceof u?s(a):a}quantile(r,e,n=!1){const a=t.quantile(this,r,e,n);return a instanceof u?s(a):a}nanquantile(r,e,n=!1){const a=t.nanquantile(this,r,e,n);return a instanceof u?s(a):a}nanpercentile(r,e,n=!1){const a=t.nanpercentile(this,r,e,n);return a instanceof u?s(a):a}cumsum(r){return s(t.cumsum(this,r))}cumprod(r){return s(t.cumprod(this,r))}nancumsum(r){return s(t.nancumsum(this,r))}nancumprod(r){return s(t.nancumprod(this,r))}sort(r=-1){return s(t.sort(this,r))}argsort(r=-1){return s(t.argsort(this,r))}partition(r,e=-1){return s(t.partition(this,r,e))}argpartition(r,e=-1){return s(t.argpartition(this,r,e))}diagonal(r=0,e=0,n=1){return s(t.diagonal(this,r,e,n))}resize(r){return s(t.resize(this,r))}diff(r=1,e=-1){return s(t.diff(this,r,e))}take(r,e){return s(t.take(this,r,e))}repeat(r,e){return s(t.repeat(this,r,e))}transpose(r){return s(t.transpose(this,r))}squeeze(r){return s(t.squeeze(this,r))}expand_dims(r){return s(t.expand_dims(this,r))}swapaxes(r,e){return s(t.swapaxes(this,r,e))}moveaxis(r,e){return s(t.moveaxis(this,r,e))}nonzero(){return t.nonzero(this).map(s)}}function T(...c){let r=[],e="xy";for(const i of c)i instanceof b?r.push(i):typeof i=="object"&&"indexing"in i&&(e=i.indexing||"xy");if(r.length===0)return[];if(r.length===1)return[r[0].copy()];const n=r.map(i=>i.size);e==="xy"&&r.length>=2&&(r=[r[1],r[0],...r.slice(2)],[n[0],n[1]]=[n[1],n[0]]);const a=n,o=a.length,h=[];for(let i=0;i<r.length;i++){const y=r[i],m=y.size,l=new Array(o).fill(1);l[i]=m;const f=y.reshape(...l),A=t.broadcast_to(f,a),D=b.fromStorage(A.storage.copy());h.push(D)}return e==="xy"&&h.length>=2&&([h[0],h[1]]=[h[1],h[0]]),h}export{b as NDArray,T as meshgrid};