0)if(++s>u)for(s-=u;s--;)l+="0";else s1&&((c[p+1]===null||c[p+1]===void 0)&&(c[p+1]=0),c[p+1]+=c[p]>>1,c[p]&=1)}return c.reverse()}function qB(t,e){if(t.isFinite()&&!t.isInteger()||e.isFinite()&&!e.isInteger())throw new Error("Integers expected in function bitXor");var r=t.constructor;if(t.isNaN()||e.isNaN())return new r(NaN);if(t.isZero())return e;if(e.isZero())return t;if(t.eq(e))return new r(0);var n=new r(-1);return t.eq(n)?Tf(e):e.eq(n)?Tf(t):!t.isFinite()||!e.isFinite()?!t.isFinite()&&!e.isFinite()?n:new r(t.isNegative()===e.isNegative()?1/0:-1/0):n1(t,e,function(i,a){return i^a})}function zB(t,e){if(t.isFinite()&&!t.isInteger()||e.isFinite()&&!e.isInteger())throw new Error("Integers expected in function leftShift");var r=t.constructor;return t.isNaN()||e.isNaN()||e.isNegative()&&!e.isZero()?new r(NaN):t.isZero()||e.isZero()?t:!t.isFinite()&&!e.isFinite()?new r(NaN):e.lt(55)?t.times(Math.pow(2,e.toNumber())+""):t.times(new r(2).pow(e))}function HB(t,e){if(t.isFinite()&&!t.isInteger()||e.isFinite()&&!e.isInteger())throw new Error("Integers expected in function rightArithShift");var r=t.constructor;return t.isNaN()||e.isNaN()||e.isNegative()&&!e.isZero()?new r(NaN):t.isZero()||e.isZero()?t:e.isFinite()?e.lt(55)?t.div(Math.pow(2,e.toNumber())+"").floor():t.div(new r(2).pow(e)).floor():t.isNegative()?new r(-1):t.isFinite()?new r(0):new r(NaN)}var GB="bitAnd",ZQ=["typed","matrix","equalScalar","concat"],sm=M(GB,ZQ,t=>{var{typed:e,matrix:r,equalScalar:n,concat:i}=t,a=In({typed:e,equalScalar:n}),o=Gu({typed:e,equalScalar:n}),s=Pr({typed:e,equalScalar:n}),l=Pt({typed:e,matrix:r,concat:i});return e(GB,{"number, number":Nw,"BigNumber, BigNumber":VB,"bigint, bigint":(u,c)=>u&c},l({SS:o,DS:a,Ss:s}))});var WB="bitNot",JQ=["typed"],i1=M(WB,JQ,t=>{var{typed:e}=t;return e(WB,{number:Fw,BigNumber:Tf,bigint:r=>~r,"Array | Matrix":e.referToSelf(r=>n=>nt(n,r))})});var YB="bitOr",XQ=["typed","matrix","equalScalar","DenseMatrix","concat"],um=M(YB,XQ,t=>{var{typed:e,matrix:r,equalScalar:n,DenseMatrix:i,concat:a}=t,o=Ri({typed:e}),s=Ef({typed:e,equalScalar:n}),l=Sa({typed:e,DenseMatrix:i}),u=Pt({typed:e,matrix:r,concat:a});return e(YB,{"number, number":Iw,"BigNumber, BigNumber":UB,"bigint, bigint":(c,f)=>c|f},u({SS:s,DS:o,Ss:l}))});var QQ="matAlgo07xSSf",KQ=["typed","SparseMatrix"],Gn=M(QQ,KQ,t=>{var{typed:e,SparseMatrix:r}=t;return function(a,o,s){var l=a._size,u=a._datatype||a._data===void 0?a._datatype:a.getDataType(),c=o._size,f=o._datatype||o._data===void 0?o._datatype:o.getDataType();if(l.length!==c.length)throw new Ct(l.length,c.length);if(l[0]!==c[0]||l[1]!==c[1])throw new RangeError("Dimension mismatch. Matrix A ("+l+") must match Matrix B ("+c+")");var m=l[0],p=l[1],d,h=0,x=s;typeof u=="string"&&u===f&&u!=="mixed"&&(d=u,h=e.convert(0,d),x=e.find(s,[d,d]));for(var b=[],T=[],v=new Array(p+1).fill(0),A=[],w=[],y=[],E=[],D=0;D{var{typed:e,matrix:r,DenseMatrix:n,concat:i,SparseMatrix:a}=t,o=jr({typed:e}),s=Gn({typed:e,SparseMatrix:a}),l=wr({typed:e,DenseMatrix:n}),u=Pt({typed:e,matrix:r,concat:i});return e(jB,{"number, number":Pw,"BigNumber, BigNumber":qB,"bigint, bigint":(c,f)=>c^f},u({SS:s,DS:o,Ss:l}))});var ZB="arg",tK=["typed"],o1=M(ZB,tK,t=>{var{typed:e}=t;return e(ZB,{number:function(n){return Math.atan2(0,n)},BigNumber:function(n){return n.constructor.atan2(0,n)},Complex:function(n){return n.arg()},"Array | Matrix":e.referToSelf(r=>n=>nt(n,r))})});var JB="conj",rK=["typed"],s1=M(JB,rK,t=>{var{typed:e}=t;return e(JB,{"number | BigNumber | Fraction":r=>r,Complex:r=>r.conjugate(),"Array | Matrix":e.referToSelf(r=>n=>nt(n,r))})});var XB="im",nK=["typed"],u1=M(XB,nK,t=>{var{typed:e}=t;return e(XB,{number:()=>0,"BigNumber | Fraction":r=>r.mul(0),Complex:r=>r.im,"Array | Matrix":e.referToSelf(r=>n=>nt(n,r))})});var QB="re",iK=["typed"],l1=M(QB,iK,t=>{var{typed:e}=t;return e(QB,{"number | BigNumber | Fraction":r=>r,Complex:r=>r.re,"Array | Matrix":e.referToSelf(r=>n=>nt(n,r))})});var KB="not",aK=["typed"],c1=M(KB,aK,t=>{var{typed:e}=t;return e(KB,{"null | undefined":()=>!0,number:Lw,Complex:function(n){return n.re===0&&n.im===0},BigNumber:function(n){return n.isZero()||n.isNaN()},bigint:r=>!r,Unit:e.referToSelf(r=>n=>e.find(r,n.valueType())(n.value)),"Array | Matrix":e.referToSelf(r=>n=>nt(n,r))})});var ek="or",oK=["typed","matrix","equalScalar","DenseMatrix","concat"],lm=M(ek,oK,t=>{var{typed:e,matrix:r,equalScalar:n,DenseMatrix:i,concat:a}=t,o=jr({typed:e}),s=Hu({typed:e,equalScalar:n}),l=wr({typed:e,DenseMatrix:i}),u=Pt({typed:e,matrix:r,concat:a});return e(ek,{"number, number":Jg,"Complex, Complex":function(f,m){return f.re!==0||f.im!==0||m.re!==0||m.im!==0},"BigNumber, BigNumber":function(f,m){return!f.isZero()&&!f.isNaN()||!m.isZero()&&!m.isNaN()},"bigint, bigint":Jg,"Unit, Unit":e.referToSelf(c=>(f,m)=>c(f.value||0,m.value||0))},u({SS:s,DS:o,Ss:l}))});var tk="xor",sK=["typed","matrix","DenseMatrix","concat","SparseMatrix"],f1=M(tk,sK,t=>{var{typed:e,matrix:r,DenseMatrix:n,concat:i,SparseMatrix:a}=t,o=jr({typed:e}),s=Gn({typed:e,SparseMatrix:a}),l=wr({typed:e,DenseMatrix:n}),u=Pt({typed:e,matrix:r,concat:i});return e(tk,{"number, number":Xg,"Complex, Complex":function(f,m){return(f.re!==0||f.im!==0)!=(m.re!==0||m.im!==0)},"bigint, bigint":Xg,"BigNumber, BigNumber":function(f,m){return(!f.isZero()&&!f.isNaN())!=(!m.isZero()&&!m.isNaN())},"Unit, Unit":e.referToSelf(c=>(f,m)=>c(f.value||0,m.value||0))},u({SS:s,DS:o,Ss:l}))});var rk="concat",uK=["typed","matrix","isInteger"],cm=M(rk,uK,t=>{var{typed:e,matrix:r,isInteger:n}=t;return e(rk,{"...Array | Matrix | number | BigNumber":function(a){var o,s=a.length,l=-1,u,c=!1,f=[];for(o=0;o0&&l>u)throw new Si(l,u+1)}else{var p=yt(m).valueOf(),d=pt(p);if(f[o]=p,u=l,l=d.length-1,o>0&&l!==u)throw new Ct(u+1,l+1)}}if(f.length===0)throw new SyntaxError("At least one matrix expected");for(var h=f.shift();f.length;)h=uw(h,f.shift(),l);return c?r(h):h},"...string":function(a){return a.join("")}})});var nk="column",lK=["typed","Index","matrix","range"],fm=M(nk,lK,t=>{var{typed:e,Index:r,matrix:n,range:i}=t;return e(nk,{"Matrix, number":a,"Array, number":function(s,l){return a(n(yt(s)),l).valueOf()}});function a(o,s){if(o.size().length!==2)throw new Error("Only two dimensional matrix is supported");Br(s,o.size()[1]);var l=i(0,o.size()[0]),u=new r(l,s),c=o.subset(u);return rt(c)?c:n([[c]])}});var ik="count",cK=["typed","size","prod"],p1=M(ik,cK,t=>{var{typed:e,size:r,prod:n}=t;return e(ik,{string:function(a){return a.length},"Matrix | Array":function(a){return n(r(a))}})});var ak="cross",fK=["typed","matrix","subtract","multiply"],m1=M(ak,fK,t=>{var{typed:e,matrix:r,subtract:n,multiply:i}=t;return e(ak,{"Matrix, Matrix":function(s,l){return r(a(s.toArray(),l.toArray()))},"Matrix, Array":function(s,l){return r(a(s.toArray(),l))},"Array, Matrix":function(s,l){return r(a(s,l.toArray()))},"Array, Array":a});function a(o,s){var l=Math.max(pt(o).length,pt(s).length);o=xf(o),s=xf(s);var u=pt(o),c=pt(s);if(u.length!==1||c.length!==1||u[0]!==3||c[0]!==3)throw new RangeError("Vectors with length 3 expected (Size A = ["+u.join(", ")+"], B = ["+c.join(", ")+"])");var f=[n(i(o[1],s[2]),i(o[2],s[1])),n(i(o[2],s[0]),i(o[0],s[2])),n(i(o[0],s[1]),i(o[1],s[0]))];return l>1?[f]:f}});var ok="diag",pK=["typed","matrix","DenseMatrix","SparseMatrix"],d1=M(ok,pK,t=>{var{typed:e,matrix:r,DenseMatrix:n,SparseMatrix:i}=t;return e(ok,{Array:function(u){return a(u,0,pt(u),null)},"Array, number":function(u,c){return a(u,c,pt(u),null)},"Array, BigNumber":function(u,c){return a(u,c.toNumber(),pt(u),null)},"Array, string":function(u,c){return a(u,0,pt(u),c)},"Array, number, string":function(u,c,f){return a(u,c,pt(u),f)},"Array, BigNumber, string":function(u,c,f){return a(u,c.toNumber(),pt(u),f)},Matrix:function(u){return a(u,0,u.size(),u.storage())},"Matrix, number":function(u,c){return a(u,c,u.size(),u.storage())},"Matrix, BigNumber":function(u,c){return a(u,c.toNumber(),u.size(),u.storage())},"Matrix, string":function(u,c){return a(u,0,u.size(),c)},"Matrix, number, string":function(u,c,f){return a(u,c,u.size(),f)},"Matrix, BigNumber, string":function(u,c,f){return a(u,c.toNumber(),u.size(),f)}});function a(l,u,c,f){if(!Qe(u))throw new TypeError("Second parameter in function diag must be an integer");var m=u>0?u:0,p=u<0?-u:0;switch(c.length){case 1:return o(l,u,f,c[0],p,m);case 2:return s(l,u,f,c,p,m)}throw new RangeError("Matrix for function diag must be 2 dimensional")}function o(l,u,c,f,m,p){var d=[f+m,f+p];if(c&&c!=="sparse"&&c!=="dense")throw new TypeError("Unknown matrix type ".concat(c,'"'));var h=c==="sparse"?i.diagonal(d,l,u):n.diagonal(d,l,u);return c!==null?h:h.valueOf()}function s(l,u,c,f,m,p){if(rt(l)){var d=l.diagonal(u);return c!==null?c!==d.storage()?r(d,c):d:d.valueOf()}for(var h=Math.min(f[0]-m,f[1]-p),x=[],b=0;b{var{typed:e}=t;return e("filter",{"Array, function":sk,"Matrix, function":function(n,i){return n.create(sk(n.valueOf(),i),n.datatype())},"Array, RegExp":ow,"Matrix, RegExp":function(n,i){return n.create(ow(n.valueOf(),i),n.datatype())}})});function sk(t,e){var r=bo(e,t,"filter");return qP(t,function(n,i,a){return r(n,[i],a)})}var uk="flatten",hK=["typed"],h1=M(uk,hK,t=>{var{typed:e}=t;return e(uk,{Array:function(n){return Rt(n)},Matrix:function(n){return n.create(Rt(n.toArray()),n.datatype())}})});var g1="forEach",gK=["typed"],mm=M(g1,gK,t=>{var{typed:e}=t;return e(g1,{"Array, function":vK,"Matrix, function":function(n,i){n.forEach(i)}})});function vK(t,e){tm(t,[],t,bo(e,t,g1))}var lk="getMatrixDataType",yK=["typed"],v1=M(lk,yK,t=>{var{typed:e}=t;return e(lk,{Array:function(n){return Yl(n,fr)},Matrix:function(n){return n.getDataType()}})});var ck="identity",xK=["typed","config","matrix","BigNumber","DenseMatrix","SparseMatrix"],y1=M(ck,xK,t=>{var{typed:e,config:r,matrix:n,BigNumber:i,DenseMatrix:a,SparseMatrix:o}=t;return e(ck,{"":function(){return r.matrix==="Matrix"?n([]):[]},string:function(c){return n(c)},"number | BigNumber":function(c){return l(c,c,r.matrix==="Matrix"?"dense":void 0)},"number | BigNumber, string":function(c,f){return l(c,c,f)},"number | BigNumber, number | BigNumber":function(c,f){return l(c,f,r.matrix==="Matrix"?"dense":void 0)},"number | BigNumber, number | BigNumber, string":function(c,f,m){return l(c,f,m)},Array:function(c){return s(c)},"Array, string":function(c,f){return s(c,f)},Matrix:function(c){return s(c.valueOf(),c.storage())},"Matrix, string":function(c,f){return s(c.valueOf(),f)}});function s(u,c){switch(u.length){case 0:return c?n(c):[];case 1:return l(u[0],u[0],c);case 2:return l(u[0],u[1],c);default:throw new Error("Vector containing two values expected")}}function l(u,c,f){var m=ft(u)||ft(c)?i:null;if(ft(u)&&(u=u.toNumber()),ft(c)&&(c=c.toNumber()),!Qe(u)||u<1)throw new Error("Parameters in function identity must be positive integers");if(!Qe(c)||c<1)throw new Error("Parameters in function identity must be positive integers");var p=m?new i(1):1,d=m?new m(0):0,h=[u,c];if(f){if(f==="sparse")return o.diagonal(h,p,0,d);if(f==="dense")return a.diagonal(h,p,0,d);throw new TypeError('Unknown matrix type "'.concat(f,'"'))}for(var x=Wo([],h,d),b=u{var{typed:e,matrix:r,multiplyScalar:n}=t;return e(fk,{"Matrix, Matrix":function(o,s){return r(i(o.toArray(),s.toArray()))},"Matrix, Array":function(o,s){return r(i(o.toArray(),s))},"Array, Matrix":function(o,s){return r(i(o,s.toArray()))},"Array, Array":i});function i(a,o){if(pt(a).length===1&&(a=[a]),pt(o).length===1&&(o=[o]),pt(a).length>2||pt(o).length>2)throw new RangeError("Vectors with dimensions greater then 2 are not supported expected (Size x = "+JSON.stringify(a.length)+", y = "+JSON.stringify(o.length)+")");var s=[],l=[];return a.map(function(u){return o.map(function(c){return l=[],s.push(l),u.map(function(f){return c.map(function(m){return l.push(n(f,m))})})})})&&s}});var b1="map",wK=["typed"],dm=M(b1,wK,t=>{var{typed:e}=t;return e(b1,{"Array, function":n,"Matrix, function":function(a,o){return a.map(o)},"Array|Matrix, Array|Matrix, ...Array|Matrix|function":(i,a,o)=>r([i,a,...o.slice(0,o.length-1)],o[o.length-1])});function r(i,a){if(typeof a!="function")throw new Error("Last argument must be a callback function");var o=i[0].isMatrix,s=em(...i.map(A=>A.isMatrix?A.size():pt(A))),l=o?(A,w)=>A.get(w):wf,u=o?i.map(A=>A.isMatrix?A.create(hs(A.toArray(),s),A.datatype()):i[0].create(hs(A.valueOf(),s))):i.map(A=>A.isMatrix?hs(A.toArray(),s):hs(A,s)),c;if(e.isTypedFunction(a)){var f=s.map(()=>0),m=u.map(A=>l(A,f)),p=v(a,m,f,u);c=b(p)}else{var d=i.length,h=T(a,d);c=b(h)}var x=(A,w)=>c([A,...u.slice(1).map(y=>l(y,w))],w);if(o)return u[0].map(x);return n(u[0],x);function b(A){switch(A){case 0:return w=>a(...w);case 1:return(w,y)=>a(...w,y);case 2:return(w,y)=>a(...w,y,...u)}}function T(A,w){return A.length>w+1?2:A.length===w+1?1:0}function v(A,w,y,E){return e.resolve(A,[...w,y,...E])!==null?2:e.resolve(A,[...w,y])!==null?1:(e.resolve(A,w)!==null,0)}}function n(i,a){return tm(i,[],i,bo(a,i,b1))}});var pk="diff",_K=["typed","matrix","subtract","number"],hm=M(pk,_K,t=>{var{typed:e,matrix:r,subtract:n,number:i}=t;return e(pk,{"Array | Matrix":function(c){return rt(c)?r(o(c.toArray())):o(c)},"Array | Matrix, number":function(c,f){if(!Qe(f))throw new RangeError("Dimension must be a whole number");return rt(c)?r(a(c.toArray(),f)):a(c,f)},"Array, BigNumber":e.referTo("Array,number",u=>(c,f)=>u(c,i(f))),"Matrix, BigNumber":e.referTo("Matrix,number",u=>(c,f)=>u(c,i(f)))});function a(u,c){if(rt(u)&&(u=u.toArray()),!Array.isArray(u))throw RangeError("Array/Matrix does not have that many dimensions");if(c>0){var f=[];return u.forEach(m=>{f.push(a(m,c-1))}),f}else{if(c===0)return o(u);throw RangeError("Cannot have negative dimension")}}function o(u){for(var c=[],f=u.length,m=1;m{var{typed:e,config:r,matrix:n,BigNumber:i}=t;return e("ones",{"":function(){return r.matrix==="Array"?a([]):a([],"default")},"...number | BigNumber | string":function(u){var c=u[u.length-1];if(typeof c=="string"){var f=u.pop();return a(u,f)}else return r.matrix==="Array"?a(u):a(u,"default")},Array:a,Matrix:function(u){var c=u.storage();return a(u.valueOf(),c)},"Array | Matrix, string":function(u,c){return a(u.valueOf(),c)}});function a(l,u){var c=o(l),f=c?new i(1):1;if(s(l),u){var m=n(u);return l.length>0?m.resize(l,f):m}else{var p=[];return l.length>0?Wo(p,l,f):p}}function o(l){var u=!1;return l.forEach(function(c,f,m){ft(c)&&(u=!0,m[f]=c.toNumber())}),u}function s(l){l.forEach(function(u){if(typeof u!="number"||!Qe(u)||u<0)throw new Error("Parameters in function ones must be positive integers")})}});function Df(){throw new Error('No "bignumber" implementation available')}function nv(){throw new Error('No "fraction" implementation available')}function iv(){throw new Error('No "matrix" implementation available')}var mk="range",EK=["typed","config","?matrix","?bignumber","smaller","smallerEq","larger","largerEq","add","isPositive"],gm=M(mk,EK,t=>{var{typed:e,config:r,matrix:n,bignumber:i,smaller:a,smallerEq:o,larger:s,largerEq:l,add:u,isPositive:c}=t;return e(mk,{string:m,"string, boolean":m,"number, number":function(x,b){return f(p(x,b,1,!1))},"number, number, number":function(x,b,T){return f(p(x,b,T,!1))},"number, number, boolean":function(x,b,T){return f(p(x,b,1,T))},"number, number, number, boolean":function(x,b,T,v){return f(p(x,b,T,v))},"BigNumber, BigNumber":function(x,b){var T=x.constructor;return f(p(x,b,new T(1),!1))},"BigNumber, BigNumber, BigNumber":function(x,b,T){return f(p(x,b,T,!1))},"BigNumber, BigNumber, boolean":function(x,b,T){var v=x.constructor;return f(p(x,b,new v(1),T))},"BigNumber, BigNumber, BigNumber, boolean":function(x,b,T,v){return f(p(x,b,T,v))},"Unit, Unit, Unit":function(x,b,T){return f(p(x,b,T,!1))},"Unit, Unit, Unit, boolean":function(x,b,T,v){return f(p(x,b,T,v))}});function f(h){return r.matrix==="Matrix"?n?n(h):iv():h}function m(h,x){var b=d(h);if(!b)throw new SyntaxError('String "'+h+'" is no valid range');return r.number==="BigNumber"?(i===void 0&&Df(),f(p(i(b.start),i(b.end),i(b.step)),x)):f(p(b.start,b.end,b.step,x))}function p(h,x,b,T){for(var v=[],A=c(b)?T?o:a:T?l:s,w=h;A(w,x);)v.push(w),w=u(w,b);return v}function d(h){var x=h.split(":"),b=x.map(function(v){return Number(v)}),T=b.some(function(v){return isNaN(v)});if(T)return null;switch(b.length){case 2:return{start:b[0],end:b[1],step:1};case 3:return{start:b[0],end:b[2],step:b[1]};default:return null}}});var dk="reshape",TK=["typed","isInteger","matrix"],_1=M(dk,TK,t=>{var{typed:e,isInteger:r}=t;return e(dk,{"Matrix, Array":function(i,a){return i.reshape(a,!0)},"Array, Array":function(i,a){return a.forEach(function(o){if(!r(o))throw new TypeError("Invalid size for dimension: "+o)}),Qp(i,a)}})});var DK="resize",MK=["config","matrix"],A1=M(DK,MK,t=>{var{config:e,matrix:r}=t;return function(a,o,s){if(arguments.length!==2&&arguments.length!==3)throw new tu("resize",arguments.length,2,3);if(rt(o)&&(o=o.valueOf()),ft(o[0])&&(o=o.map(function(c){return ft(c)?c.toNumber():c})),rt(a))return a.resize(o,s,!0);if(typeof a=="string")return n(a,o,s);var l=Array.isArray(a)?!1:e.matrix!=="Array";if(o.length===0){for(;Array.isArray(a);)a=a[0];return yt(a)}else{Array.isArray(a)||(a=[a]),a=yt(a);var u=Wo(a,o,s);return l?r(u):u}};function n(i,a,o){if(o!==void 0){if(typeof o!="string"||o.length!==1)throw new TypeError("Single character expected as defaultValue")}else o=" ";if(a.length!==1)throw new Ct(a.length,1);var s=a[0];if(typeof s!="number"||!Qe(s))throw new TypeError("Invalid size, must contain positive integers (size: "+Dt(a)+")");if(i.length>s)return i.substring(0,s);if(i.length{var{typed:e,multiply:r,rotationMatrix:n}=t;return e(hk,{"Array , number | BigNumber | Complex | Unit":function(o,s){i(o,2);var l=r(n(s),o);return l.toArray()},"Matrix , number | BigNumber | Complex | Unit":function(o,s){return i(o,2),r(n(s),o)},"Array, number | BigNumber | Complex | Unit, Array | Matrix":function(o,s,l){i(o,3);var u=r(n(s,l),o);return u},"Matrix, number | BigNumber | Complex | Unit, Array | Matrix":function(o,s,l){return i(o,3),r(n(s,l),o)}});function i(a,o){var s=Array.isArray(a)?pt(a):a.size();if(s.length>2)throw new RangeError("Vector must be of dimensions 1x".concat(o));if(s.length===2&&s[1]!==1)throw new RangeError("Vector must be of dimensions 1x".concat(o));if(s[0]!==o)throw new RangeError("Vector must be of dimensions 1x".concat(o))}});var gk="rotationMatrix",NK=["typed","config","multiplyScalar","addScalar","unaryMinus","norm","matrix","BigNumber","DenseMatrix","SparseMatrix","cos","sin"],E1=M(gk,NK,t=>{var{typed:e,config:r,multiplyScalar:n,addScalar:i,unaryMinus:a,norm:o,BigNumber:s,matrix:l,DenseMatrix:u,SparseMatrix:c,cos:f,sin:m}=t;return e(gk,{"":function(){return r.matrix==="Matrix"?l([]):[]},string:function(v){return l(v)},"number | BigNumber | Complex | Unit":function(v){return p(v,r.matrix==="Matrix"?"dense":void 0)},"number | BigNumber | Complex | Unit, string":function(v,A){return p(v,A)},"number | BigNumber | Complex | Unit, Array":function(v,A){var w=l(A);return d(w),b(v,w,void 0)},"number | BigNumber | Complex | Unit, Matrix":function(v,A){d(A);var w=A.storage()||(r.matrix==="Matrix"?"dense":void 0);return b(v,A,w)},"number | BigNumber | Complex | Unit, Array, string":function(v,A,w){var y=l(A);return d(y),b(v,y,w)},"number | BigNumber | Complex | Unit, Matrix, string":function(v,A,w){return d(A),b(v,A,w)}});function p(T,v){var A=ft(T),w=A?new s(-1):-1,y=f(T),E=m(T),D=[[y,n(w,E)],[E,y]];return x(D,v)}function d(T){var v=T.size();if(v.length<1||v[0]!==3)throw new RangeError("Vector must be of dimensions 1x3")}function h(T){return T.reduce((v,A)=>n(v,A))}function x(T,v){if(v){if(v==="sparse")return new c(T);if(v==="dense")return new u(T);throw new TypeError('Unknown matrix type "'.concat(v,'"'))}return T}function b(T,v,A){var w=o(v);if(w===0)throw new RangeError("Rotation around zero vector");var y=ft(T)?s:null,E=y?new y(1):1,D=y?new y(-1):-1,_=y?new y(v.get([0])/w):v.get([0])/w,N=y?new y(v.get([1])/w):v.get([1])/w,C=y?new y(v.get([2])/w):v.get([2])/w,B=f(T),R=i(E,a(B)),V=m(T),L=i(B,h([_,_,R])),P=i(h([_,N,R]),h([D,C,V])),k=i(h([_,C,R]),h([N,V])),W=i(h([_,N,R]),h([C,V])),O=i(B,h([N,N,R])),X=i(h([N,C,R]),h([D,_,V])),ae=i(h([_,C,R]),h([D,N,V])),J=i(h([N,C,R]),h([_,V])),U=i(B,h([C,C,R])),de=[[L,P,k],[W,O,X],[ae,J,U]];return x(de,A)}});var vk="row",FK=["typed","Index","matrix","range"],vm=M(vk,FK,t=>{var{typed:e,Index:r,matrix:n,range:i}=t;return e(vk,{"Matrix, number":a,"Array, number":function(s,l){return a(n(yt(s)),l).valueOf()}});function a(o,s){if(o.size().length!==2)throw new Error("Only two dimensional matrix is supported");Br(s,o.size()[0]);var l=i(0,o.size()[1]),u=new r(s,l),c=o.subset(u);return rt(c)?c:n([[c]])}});var yk="size",IK=["typed","config","?matrix"],T1=M(yk,IK,t=>{var{typed:e,config:r,matrix:n}=t;return e(yk,{Matrix:function(a){return a.create(a.size(),"number")},Array:pt,string:function(a){return r.matrix==="Array"?[a.length]:n([a.length],"dense","number")},"number | Complex | BigNumber | Unit | boolean | null":function(a){return r.matrix==="Array"?[]:n?n([],"dense","number"):iv()}})});var xk="squeeze",PK=["typed"],D1=M(xk,PK,t=>{var{typed:e}=t;return e(xk,{Array:function(n){return xf(yt(n))},Matrix:function(n){var i=xf(n.toArray());return Array.isArray(i)?n.create(i,n.datatype()):i},any:function(n){return yt(n)}})});var bk="subset",BK=["typed","matrix","zeros","add"],ym=M(bk,BK,t=>{var{typed:e,matrix:r,zeros:n,add:i}=t;return e(bk,{"Matrix, Index":function(s,l){return Wl(l)?r():(Xp(s,l),s.subset(l))},"Array, Index":e.referTo("Matrix, Index",function(o){return function(s,l){var u=o(r(s),l);return l.isScalar()?u:u.valueOf()}}),"Object, Index":$K,"string, Index":kK,"Matrix, Index, any, any":function(s,l,u,c){return Wl(l)?s:(Xp(s,l),s.clone().subset(l,a(u,l),c))},"Array, Index, any, any":e.referTo("Matrix, Index, any, any",function(o){return function(s,l,u,c){var f=o(r(s),l,u,c);return f.isMatrix?f.valueOf():f}}),"Array, Index, any":e.referTo("Matrix, Index, any, any",function(o){return function(s,l,u){return o(r(s),l,u,void 0).valueOf()}}),"Matrix, Index, any":e.referTo("Matrix, Index, any, any",function(o){return function(s,l,u){return o(s,l,u,void 0)}}),"string, Index, string":wk,"string, Index, string, string":wk,"Object, Index, any":OK});function a(o,s){if(typeof o=="string")throw new Error("can't boradcast a string");if(s._isScalar)return o;var l=s.size();if(l.every(u=>u>0))try{return i(o,n(l))}catch{return o}else return o}});function kK(t,e){if(!$u(e))throw new TypeError("Index expected");if(Wl(e))return"";if(Xp(Array.from(t),e),e.size().length!==1)throw new Ct(e.size().length,1);var r=t.length;Br(e.min()[0],r),Br(e.max()[0],r);var n=e.dimension(0),i="";return n.forEach(function(a){i+=t.charAt(a)}),i}function wk(t,e,r,n){if(!e||e.isIndex!==!0)throw new TypeError("Index expected");if(Wl(e))return t;if(Xp(Array.from(t),e),e.size().length!==1)throw new Ct(e.size().length,1);if(n!==void 0){if(typeof n!="string"||n.length!==1)throw new TypeError("Single character expected as defaultValue")}else n=" ";var i=e.dimension(0),a=i.size()[0];if(a!==r.length)throw new Ct(i.size()[0],r.length);var o=t.length;Br(e.min()[0]),Br(e.max()[0]);for(var s=[],l=0;lo)for(var u=o-1,c=s.length;u{var{typed:e,matrix:r}=t;return e(_k,{Array:o=>n(r(o)).valueOf(),Matrix:n,any:yt});function n(o){var s=o.size(),l;switch(s.length){case 1:l=o.clone();break;case 2:{var u=s[0],c=s[1];if(c===0)throw new RangeError("Cannot transpose a 2D matrix with no columns (size: "+Dt(s)+")");switch(o.storage()){case"dense":l=i(o,u,c);break;case"sparse":l=a(o,u,c);break}}break;default:throw new RangeError("Matrix must be a vector or two dimensional (size: "+Dt(s)+")")}return l}function i(o,s,l){for(var u=o._data,c=[],f,m=0;m{var{typed:e,transpose:r,conj:n}=t;return e(Ak,{any:function(a){return n(r(a))}})});var Sk="zeros",VK=["typed","config","matrix","BigNumber"],N1=M(Sk,VK,t=>{var{typed:e,config:r,matrix:n,BigNumber:i}=t;return e(Sk,{"":function(){return r.matrix==="Array"?a([]):a([],"default")},"...number | BigNumber | string":function(u){var c=u[u.length-1];if(typeof c=="string"){var f=u.pop();return a(u,f)}else return r.matrix==="Array"?a(u):a(u,"default")},Array:a,Matrix:function(u){var c=u.storage();return a(u.valueOf(),c)},"Array | Matrix, string":function(u,c){return a(u.valueOf(),c)}});function a(l,u){var c=o(l),f=c?new i(0):0;if(s(l),u){var m=n(u);return l.length>0?m.resize(l,f):m}else{var p=[];return l.length>0?Wo(p,l,f):p}}function o(l){var u=!1;return l.forEach(function(c,f,m){ft(c)&&(u=!0,m[f]=c.toNumber())}),u}function s(l){l.forEach(function(u){if(typeof u!="number"||!Qe(u)||u<0)throw new Error("Parameters in function zeros must be positive integers")})}});var Ek="fft",UK=["typed","matrix","addScalar","multiplyScalar","divideScalar","exp","tau","i","dotDivide","conj","pow","ceil","log2"],F1=M(Ek,UK,t=>{var{typed:e,matrix:r,addScalar:n,multiplyScalar:i,divideScalar:a,exp:o,tau:s,i:l,dotDivide:u,conj:c,pow:f,ceil:m,log2:p}=t;return e(Ek,{Array:d,Matrix:function(v){return v.create(d(v.valueOf()),v.datatype())}});function d(T){var v=pt(T);return v.length===1?b(T,v[0]):h(T.map(A=>d(A,v.slice(1))),0)}function h(T,v){var A=pt(T);if(v!==0)return new Array(A[0]).fill(0).map((y,E)=>h(T[E],v-1));if(A.length===1)return b(T);function w(y){var E=pt(y);return new Array(E[1]).fill(0).map((D,_)=>new Array(E[0]).fill(0).map((N,C)=>y[C][_]))}return w(h(w(T),1))}function x(T){for(var v=T.length,A=o(a(i(-1,i(l,s)),v)),w=[],y=1-v;yi(T[k],w[v-1+k])),...new Array(E-v).fill(0)],_=[...new Array(v+v-1).fill(0).map((P,k)=>a(1,w[k])),...new Array(E-(v+v-1)).fill(0)],N=b(D),C=b(_),B=new Array(E).fill(0).map((P,k)=>i(N[k],C[k])),R=u(c(d(c(B))),E),V=[],L=v-1;L_%2===0),v/2),...b(T.filter((D,_)=>_%2===1),v/2)],w=0;w{var{typed:e,fft:r,dotDivide:n,conj:i}=t;return e(Tk,{"Array | Matrix":function(o){var s=rt(o)?o.size():pt(o);return n(i(r(i(o))),s.reduce((l,u)=>l*u,1))}})});function Jl(t){"@babel/helpers - typeof";return Jl=typeof Symbol=="function"&&typeof Symbol.iterator=="symbol"?function(e){return typeof e}:function(e){return e&&typeof Symbol=="function"&&e.constructor===Symbol&&e!==Symbol.prototype?"symbol":typeof e},Jl(t)}function Dk(t,e){if(Jl(t)!="object"||!t)return t;var r=t[Symbol.toPrimitive];if(r!==void 0){var n=r.call(t,e||"default");if(Jl(n)!="object")return n;throw new TypeError("@@toPrimitive must return a primitive value.")}return(e==="string"?String:Number)(t)}function Mk(t){var e=Dk(t,"string");return Jl(e)=="symbol"?e:e+""}function Er(t,e,r){return(e=Mk(e))in t?Object.defineProperty(t,e,{value:r,enumerable:!0,configurable:!0,writable:!0}):t[e]=r,t}function Ck(t,e){var r=Object.keys(t);if(Object.getOwnPropertySymbols){var n=Object.getOwnPropertySymbols(t);e&&(n=n.filter(function(i){return Object.getOwnPropertyDescriptor(t,i).enumerable})),r.push.apply(r,n)}return r}function zK(t){for(var e=1;e{var{typed:e,add:r,subtract:n,multiply:i,divide:a,max:o,map:s,abs:l,isPositive:u,isNegative:c,larger:f,smaller:m,matrix:p,bignumber:d,unaryMinus:h}=t;function x(D){return function(_,N,C,B){var R=!(N.length===2&&(N.every(y)||N.every(bi)));if(R)throw new Error('"tspan" must be an Array of two numeric values or two units [tStart, tEnd]');var V=N[0],L=N[1],P=f(L,V),k=B.firstStep;if(k!==void 0&&!u(k))throw new Error('"firstStep" must be positive');var W=B.maxStep;if(W!==void 0&&!u(W))throw new Error('"maxStep" must be positive');var O=B.minStep;if(O&&c(O))throw new Error('"minStep" must be positive or zero');var X=[V,L,k,O,W].filter(St=>St!==void 0);if(!(X.every(y)||X.every(bi)))throw new Error('Inconsistent type of "t" dependant variables');for(var ae=1,J=B.tol?B.tol:1e-4,U=B.minDelta?B.minDelta:.2,de=B.maxDelta?B.maxDelta:5,we=B.maxIter?B.maxIter:1e4,se=[V,L,...C,W,O].some(ft),[be,K,j,ne]=se?[d(D.a),d(D.c),d(D.b),d(D.bp)]:[D.a,D.c,D.b,D.bp],pe=k?P?k:h(k):a(n(L,V),ae),fe=[V],Ae=[C],_e=n(j,ne),G=0,Z=0,ge=A(P),xe=w(P);ge(fe[G],L);){var Ee=[];pe=xe(fe[G],L,pe),Ee.push(_(fe[G],Ae[G]));for(var Be=1;Bebi(St)?St.value:St)));Ve1/4&&(fe.push(r(fe[G],pe)),Ae.push(r(Ae[G],i(pe,j,Ee))),G++);var et=.84*(J/Ve)**(1/5);if(m(et,U)?et=U:f(et,de)&&(et=de),et=se?d(et):et,pe=i(pe,et),W&&f(l(pe),W)?pe=P?W:h(W):O&&m(l(pe),O)&&(pe=P?O:h(O)),Z++,Z>we)throw new Error("Maximum number of iterations reached, try changing options")}return{t:fe,y:Ae}}}function b(D,_,N,C){var B=[[],[.5],[0,.75],[.2222222222222222,.3333333333333333,.4444444444444444]],R=[null,1/2,3/4,1],V=[2/9,1/3,4/9,0],L=[7/24,1/4,1/3,1/8],P={a:B,c:R,b:V,bp:L};return x(P)(D,_,N,C)}function T(D,_,N,C){var B=[[],[.2],[.075,.225],[.9777777777777777,-3.7333333333333334,3.5555555555555554],[2.9525986892242035,-11.595793324188385,9.822892851699436,-.2908093278463649],[2.8462752525252526,-10.757575757575758,8.906422717743473,.2784090909090909,-.2735313036020583],[.09114583333333333,0,.44923629829290207,.6510416666666666,-.322376179245283,.13095238095238096]],R=[null,1/5,3/10,4/5,8/9,1,1],V=[35/384,0,500/1113,125/192,-2187/6784,11/84,0],L=[5179/57600,0,7571/16695,393/640,-92097/339200,187/2100,1/40],P={a:B,c:R,b:V,bp:L};return x(P)(D,_,N,C)}function v(D,_,N,C){var B=C.method?C.method:"RK45",R={RK23:b,RK45:T};if(B.toUpperCase()in R){var V=zK({},C);return delete V.method,R[B.toUpperCase()](D,_,N,V)}else{var L=Object.keys(R).map(k=>'"'.concat(k,'"')),P="".concat(L.slice(0,-1).join(", ")," and ").concat(L.slice(-1));throw new Error('Unavailable method "'.concat(B,'". Available methods are ').concat(P))}}function A(D){return D?m:f}function w(D){var _=D?f:m;return function(N,C,B){var R=r(N,B);return _(R,C)?n(C,N):B}}function y(D){return ft(D)||At(D)}function E(D,_,N,C){var B=v(D,_.toArray(),N.toArray(),C);return{t:p(B.t),y:p(B.y)}}return e("solveODE",{"function, Array, Array, Object":v,"function, Matrix, Matrix, Object":E,"function, Array, Array":(D,_,N)=>v(D,_,N,{}),"function, Matrix, Matrix":(D,_,N)=>E(D,_,N,{}),"function, Array, number | BigNumber | Unit":(D,_,N)=>{var C=v(D,_,[N],{});return{t:C.t,y:C.y.map(B=>B[0])}},"function, Matrix, number | BigNumber | Unit":(D,_,N)=>{var C=v(D,_.toArray(),[N],{});return{t:p(C.t),y:p(C.y.map(B=>B[0]))}},"function, Array, number | BigNumber | Unit, Object":(D,_,N,C)=>{var B=v(D,_,[N],C);return{t:B.t,y:B.y.map(R=>R[0])}},"function, Matrix, number | BigNumber | Unit, Object":(D,_,N,C)=>{var B=v(D,_.toArray(),[N],C);return{t:p(B.t),y:p(B.y.map(R=>R[0]))}}})});var WK="erf",YK=["typed"],B1=M(WK,YK,t=>{var{typed:e}=t;return e("name",{number:function(o){var s=Math.abs(o);return s>=JK?vo(o):s<=jK?vo(o)*r(s):s<=4?vo(o)*(1-n(s)):vo(o)*(1-i(s))},"Array | Matrix":e.referToSelf(a=>o=>nt(o,a))});function r(a){var o=a*a,s=ru[0][4]*o,l=o,u;for(u=0;u<3;u+=1)s=(s+ru[0][u])*o,l=(l+Mf[0][u])*o;return a*(s+ru[0][3])/(l+Mf[0][3])}function n(a){var o=ru[1][8]*a,s=a,l;for(l=0;l<7;l+=1)o=(o+ru[1][l])*a,s=(s+Mf[1][l])*a;var u=(o+ru[1][7])/(s+Mf[1][7]),c=parseInt(a*16)/16,f=(a-c)*(a+c);return Math.exp(-c*c)*Math.exp(-f)*u}function i(a){var o=1/(a*a),s=ru[2][5]*o,l=o,u;for(u=0;u<4;u+=1)s=(s+ru[2][u])*o,l=(l+Mf[2][u])*o;var c=o*(s+ru[2][4])/(l+Mf[2][4]);c=(ZK-c)/a,o=parseInt(a*16)/16;var f=(a-o)*(a+o);return Math.exp(-o*o)*Math.exp(-f)*c}}),jK=.46875,ZK=.5641895835477563,ru=[[3.1611237438705655,113.86415415105016,377.485237685302,3209.3775891384694,.18577770618460315],[.5641884969886701,8.883149794388377,66.11919063714163,298.6351381974001,881.952221241769,1712.0476126340707,2051.0783778260716,1230.3393547979972,21531153547440383e-24],[.30532663496123236,.36034489994980445,.12578172611122926,.016083785148742275,.0006587491615298378,.016315387137302097]],Mf=[[23.601290952344122,244.02463793444417,1282.6165260773723,2844.236833439171],[15.744926110709835,117.6939508913125,537.1811018620099,1621.3895745666903,3290.7992357334597,4362.619090143247,3439.3676741437216,1230.3393548037495],[2.568520192289822,1.8729528499234604,.5279051029514285,.06051834131244132,.0023352049762686918]],JK=Math.pow(2,53);var Nk="zeta",XK=["typed","config","multiply","pow","divide","factorial","equal","smallerEq","isNegative","gamma","sin","subtract","add","?Complex","?BigNumber","pi"],k1=M(Nk,XK,t=>{var{typed:e,config:r,multiply:n,pow:i,divide:a,factorial:o,equal:s,smallerEq:l,isNegative:u,gamma:c,sin:f,subtract:m,add:p,Complex:d,BigNumber:h,pi:x}=t;return e(Nk,{number:y=>b(y,E=>E,()=>20),BigNumber:y=>b(y,E=>new h(E),()=>Math.abs(Math.log10(r.relTol))),Complex:T});function b(y,E,D){return s(y,0)?E(-.5):s(y,1)?E(NaN):isFinite(y)?v(y,E,D,_=>_):u(y)?E(NaN):E(1)}function T(y){return y.re===0&&y.im===0?new d(-.5):y.re===1?new d(NaN,NaN):y.re===1/0&&y.im===0?new d(1):y.im===1/0||y.re===-1/0?new d(NaN,NaN):v(y,E=>E,E=>Math.round(1.3*15+.9*Math.abs(E.im)),E=>E.re)}function v(y,E,D,_){var N=D(y);if(_(y)>-(N-1)/2)return w(y,E(N),E);var C=n(i(2,y),i(E(x),m(y,1)));return C=n(C,f(n(a(E(x),2),y))),C=n(C,c(m(1,y))),n(C,v(m(1,y),E,D,_))}function A(y,E){for(var D=y,_=y;l(_,E);_=p(_,1)){var N=a(n(o(p(E,m(_,1))),i(4,_)),n(o(m(E,_)),o(n(2,_))));D=p(D,N)}return n(E,D)}function w(y,E,D){for(var _=a(1,n(A(D(0),E),m(1,i(2,m(1,y))))),N=D(0),C=D(1);l(C,E);C=p(C,1))N=p(N,a(n((-1)**(C-1),A(C,E)),i(C,y)));return n(_,N)}});var Fk="mode",QK=["typed","isNaN","isNumeric"],$1=M(Fk,QK,t=>{var{typed:e,isNaN:r,isNumeric:n}=t;return e(Fk,{"Array | Matrix":i,"...":function(o){return i(o)}});function i(a){a=Rt(a.valueOf());var o=a.length;if(o===0)throw new Error("Cannot calculate mode of an empty array");for(var s={},l=[],u=0,c=0;cu&&(u=s[f],l=[f])}return l}});function An(t,e,r){var n;return String(t).includes("Unexpected type")?(n=arguments.length>2?" (type: "+fr(r)+", value: "+JSON.stringify(r)+")":" (type: "+t.data.actual+")",new TypeError("Cannot calculate "+e+", unexpected type of argument"+n)):String(t).includes("complex numbers")?(n=arguments.length>2?" (type: "+fr(r)+", value: "+JSON.stringify(r)+")":"",new TypeError("Cannot calculate "+e+", no ordering relation is defined for complex numbers"+n)):t}var Ik="prod",KK=["typed","config","multiplyScalar","numeric"],O1=M(Ik,KK,t=>{var{typed:e,config:r,multiplyScalar:n,numeric:i}=t;return e(Ik,{"Array | Matrix":a,"Array | Matrix, number | BigNumber":function(s,l){throw new Error("prod(A, dim) is not yet supported")},"...":function(s){return a(s)}});function a(o){var s;if(Oi(o,function(l){try{s=s===void 0?l:n(s,l)}catch(u){throw An(u,"prod",l)}}),typeof s=="string"&&(s=i(s,Bi(s,r))),s===void 0)throw new Error("Cannot calculate prod of an empty array");return s}});var Pk="format",eee=["typed"],R1=M(Pk,eee,t=>{var{typed:e}=t;return e(Pk,{any:Dt,"any, Object | function | number | BigNumber":Dt})});var Bk="bin",tee=["typed","format"],L1=M(Bk,tee,t=>{var{typed:e,format:r}=t;return e(Bk,{"number | BigNumber":function(i){return r(i,{notation:"bin"})},"number | BigNumber, number | BigNumber":function(i,a){return r(i,{notation:"bin",wordSize:a})}})});var kk="oct",ree=["typed","format"],V1=M(kk,ree,t=>{var{typed:e,format:r}=t;return e(kk,{"number | BigNumber":function(i){return r(i,{notation:"oct"})},"number | BigNumber, number | BigNumber":function(i,a){return r(i,{notation:"oct",wordSize:a})}})});var $k="hex",nee=["typed","format"],U1=M($k,nee,t=>{var{typed:e,format:r}=t;return e($k,{"number | BigNumber":function(i){return r(i,{notation:"hex"})},"number | BigNumber, number | BigNumber":function(i,a){return r(i,{notation:"hex",wordSize:a})}})});var av=/\$([\w.]+)/g;var Ok="print",iee=["typed"],xm=M(Ok,iee,t=>{var{typed:e}=t;return e(Ok,{"string, Object | Array":Rk,"string, Object | Array, number | Object":Rk})});function Rk(t,e,r){return t.replace(av,function(n,i){var a=i.split("."),o=e[a.shift()];for(o!==void 0&&o.isMatrix&&(o=o.toArray());a.length&&o!==void 0;){var s=a.shift();o=s?o[s]:o+"."}return o!==void 0?bn(o)?o:Dt(o,r):n})}var Lk="to",aee=["typed","matrix","concat"],q1=M(Lk,aee,t=>{var{typed:e,matrix:r,concat:n}=t,i=Pt({typed:e,matrix:r,concat:n});return e(Lk,{"Unit, Unit | string":(a,o)=>a.to(o)},i({Ds:!0}))});var Vk="isPrime",oee=["typed"],z1=M(Vk,oee,t=>{var{typed:e}=t;return e(Vk,{number:function(n){if(n<=3)return n>1;if(n%2===0||n%3===0)return!1;for(var i=5;i*i<=n;i+=6)if(n%i===0||n%(i+2)===0)return!1;return!0},bigint:function(n){if(n<=3n)return n>1n;if(n%2n===0n||n%3n===0n)return!1;for(var i=5n;i*i<=n;i+=6n)if(n%i===0n||n%(i+2n)===0n)return!1;return!0},BigNumber:function(n){if(n.lte(3))return n.gt(1);if(n.mod(2).eq(0)||n.mod(3).eq(0))return!1;if(n.lt(Math.pow(2,32))){for(var i=n.toNumber(),a=5;a*a<=i;a+=6)if(i%a===0||i%(a+2)===0)return!1;return!0}function o(T,v,A){for(var w=1;!v.eq(0);)v.mod(2).eq(0)?(v=v.div(2),T=T.mul(T).mod(A)):(v=v.sub(1),w=T.mul(w).mod(A));return w}var s=n.constructor.clone({precision:n.toFixed(0).length*2});n=new s(n);for(var l=0,u=n.sub(1);u.mod(2).eq(0);)u=u.div(2),l+=1;var c=null;if(n.lt("3317044064679887385961981"))c=[2,3,5,7,11,13,17,19,23,29,31,37,41].filter(T=>Tn=>nt(n,r))})});var see="numeric",uee=["number","?bignumber","?fraction"],H1=M(see,uee,t=>{var{number:e,bignumber:r,fraction:n}=t,i={string:!0,number:!0,BigNumber:!0,Fraction:!0},a={number:o=>e(o),BigNumber:r?o=>r(o):Df,bigint:o=>BigInt(o),Fraction:n?o=>n(o):nv};return function(s){var l=arguments.length>1&&arguments[1]!==void 0?arguments[1]:"number",u=arguments.length>2?arguments[2]:void 0;if(u!==void 0)throw new SyntaxError("numeric() takes one or two arguments");var c=fr(s);if(!(c in i))throw new TypeError("Cannot convert "+s+' of type "'+c+'"; valid input types are '+Object.keys(i).join(", "));if(!(l in a))throw new TypeError("Cannot convert "+s+' to type "'+l+'"; valid output types are '+Object.keys(a).join(", "));return l===c?s:a[l](s)}});var Uk="divideScalar",lee=["typed","numeric"],G1=M(Uk,lee,t=>{var{typed:e,numeric:r}=t;return e(Uk,{"number, number":function(i,a){return i/a},"Complex, Complex":function(i,a){return i.div(a)},"BigNumber, BigNumber":function(i,a){return i.div(a)},"bigint, bigint":function(i,a){return i/a},"Fraction, Fraction":function(i,a){return i.div(a)},"Unit, number | Complex | Fraction | BigNumber | Unit":(n,i)=>n.divide(i),"number | Fraction | Complex | BigNumber, Unit":(n,i)=>i.divideInto(n)})});var qk="pow",cee=["typed","config","identity","multiply","matrix","inv","fraction","number","Complex"],W1=M(qk,cee,t=>{var{typed:e,config:r,identity:n,multiply:i,matrix:a,inv:o,number:s,fraction:l,Complex:u}=t;return e(qk,{"number, number":c,"Complex, Complex":function(d,h){return d.pow(h)},"BigNumber, BigNumber":function(d,h){return h.isInteger()||d>=0||r.predictable?d.pow(h):new u(d.toNumber(),0).pow(h.toNumber(),0)},"bigint, bigint":(p,d)=>p**d,"Fraction, Fraction":function(d,h){var x=d.pow(h);if(x!=null)return x;if(r.predictable)throw new Error("Result of pow is non-rational and cannot be expressed as a fraction");return c(d.valueOf(),h.valueOf())},"Array, number":f,"Array, BigNumber":function(d,h){return f(d,h.toNumber())},"Matrix, number":m,"Matrix, BigNumber":function(d,h){return m(d,h.toNumber())},"Unit, number | BigNumber":function(d,h){return d.pow(h)}});function c(p,d){if(r.predictable&&!Qe(d)&&p<0)try{var h=l(d),x=s(h);if((d===x||Math.abs((d-x)/d)<1e-14)&&h.d%2n===1n)return(h.n%2n===0n?1:-1)*Math.pow(-p,d)}catch{}return r.predictable&&(p<-1&&d===1/0||p>-1&&p<0&&d===-1/0)?NaN:Qe(d)||p>=0||r.predictable?Cw(p,d):p*p<1&&d===1/0||p*p>1&&d===-1/0?0:new u(p,0).pow(d,0)}function f(p,d){if(!Qe(d))throw new TypeError("For A^b, b must be an integer (value is "+d+")");var h=pt(p);if(h.length!==2)throw new Error("For A^b, A must be 2 dimensional (A has "+h.length+" dimensions)");if(h[0]!==h[1])throw new Error("For A^b, A must be square (size is "+h[0]+"x"+h[1]+")");if(d<0)try{return f(o(p),-d)}catch(T){throw T.message==="Cannot calculate inverse, determinant is zero"?new TypeError("For A^b, when A is not invertible, b must be a positive integer (value is "+d+")"):T}for(var x=n(h[0]).valueOf(),b=p;d>=1;)(d&1)===1&&(x=i(b,x)),d>>=1,b=i(b,b);return x}function m(p,d){return a(f(p.valueOf(),d))}});var Cf="Number of decimals in function round must be an integer",zk="round",fee=["typed","config","matrix","equalScalar","zeros","BigNumber","DenseMatrix"],Y1=M(zk,fee,t=>{var{typed:e,config:r,matrix:n,equalScalar:i,zeros:a,BigNumber:o,DenseMatrix:s}=t,l=Pr({typed:e,equalScalar:i}),u=wr({typed:e,DenseMatrix:s}),c=Fn({typed:e});function f(m){return Math.abs(df(m).exponent)}return e(zk,{number:function(p){var d=Af(p,f(r.relTol)),h=Kr(p,d,r.relTol,r.absTol)?d:p;return Af(h)},"number, number":function(p,d){var h=f(r.relTol);if(d>=h)return Af(p,d);var x=Af(p,h),b=Kr(p,x,r.relTol,r.absTol)?x:p;return Af(b,d)},"number, BigNumber":function(p,d){if(!d.isInteger())throw new TypeError(Cf);return new o(p).toDecimalPlaces(d.toNumber())},Complex:function(p){return p.round()},"Complex, number":function(p,d){if(d%1)throw new TypeError(Cf);return p.round(d)},"Complex, BigNumber":function(p,d){if(!d.isInteger())throw new TypeError(Cf);var h=d.toNumber();return p.round(h)},BigNumber:function(p){var d=new o(p).toDecimalPlaces(f(r.relTol)),h=mn(p,d,r.relTol,r.absTol)?d:p;return h.toDecimalPlaces(0)},"BigNumber, BigNumber":function(p,d){if(!d.isInteger())throw new TypeError(Cf);var h=f(r.relTol);if(d>=h)return p.toDecimalPlaces(d.toNumber());var x=p.toDecimalPlaces(h),b=mn(p,x,r.relTol,r.absTol)?x:p;return b.toDecimalPlaces(d.toNumber())},Fraction:function(p){return p.round()},"Fraction, number":function(p,d){if(d%1)throw new TypeError(Cf);return p.round(d)},"Fraction, BigNumber":function(p,d){if(!d.isInteger())throw new TypeError(Cf);return p.round(d.toNumber())},"Unit, number, Unit":e.referToSelf(m=>function(p,d,h){var x=p.toNumeric(h);return h.multiply(m(x,d))}),"Unit, BigNumber, Unit":e.referToSelf(m=>(p,d,h)=>m(p,d.toNumber(),h)),"Array | Matrix, number | BigNumber, Unit":e.referToSelf(m=>(p,d,h)=>nt(p,x=>m(x,d,h),!0)),"Array | Matrix | Unit, Unit":e.referToSelf(m=>(p,d)=>m(p,0,d)),"Array | Matrix":e.referToSelf(m=>p=>nt(p,m,!0)),"SparseMatrix, number | BigNumber":e.referToSelf(m=>(p,d)=>l(p,d,m,!1)),"DenseMatrix, number | BigNumber":e.referToSelf(m=>(p,d)=>c(p,d,m,!1)),"Array, number | BigNumber":e.referToSelf(m=>(p,d)=>c(n(p),d,m,!1).valueOf()),"number | Complex | BigNumber | Fraction, SparseMatrix":e.referToSelf(m=>(p,d)=>i(p,0)?a(d.size(),d.storage()):u(d,p,m,!0)),"number | Complex | BigNumber | Fraction, DenseMatrix":e.referToSelf(m=>(p,d)=>i(p,0)?a(d.size(),d.storage()):c(d,p,m,!0)),"number | Complex | BigNumber | Fraction, Array":e.referToSelf(m=>(p,d)=>c(n(d),p,m,!0).valueOf())})});var Hk="log",pee=["config","typed","typeOf","divideScalar","Complex"],j1=M(Hk,pee,t=>{var{typed:e,typeOf:r,config:n,divideScalar:i,Complex:a}=t;return e(Hk,{number:function(s){return s>=0||n.predictable?jP(s):new a(s,0).log()},Complex:function(s){return s.log()},BigNumber:function(s){return!s.isNegative()||n.predictable?s.ln():new a(s.toNumber(),0).log()},"any, any":e.referToSelf(o=>(s,l)=>{if(r(s)==="Fraction"&&r(l)==="Fraction"){var u=s.log(l);if(u!==null)return u}return i(o(s),o(l))})})});var Gk="log1p",mee=["typed","config","divideScalar","log","Complex"],Z1=M(Gk,mee,t=>{var{typed:e,config:r,divideScalar:n,log:i,Complex:a}=t;return e(Gk,{number:function(l){return l>=-1||r.predictable?Ng(l):o(new a(l,0))},Complex:o,BigNumber:function(l){var u=l.plus(1);return!u.isNegative()||r.predictable?u.ln():o(new a(l.toNumber(),0))},"Array | Matrix":e.referToSelf(s=>l=>nt(l,s)),"any, any":e.referToSelf(s=>(l,u)=>n(s(l),i(u)))});function o(s){var l=s.re+1;return new a(Math.log(Math.sqrt(l*l+s.im*s.im)),Math.atan2(s.im,l))}});var Wk="nthRoots",dee=["config","typed","divideScalar","Complex"],J1=M(Wk,dee,t=>{var{typed:e,config:r,divideScalar:n,Complex:i}=t,a=[function(l){return new i(l,0)},function(l){return new i(0,l)},function(l){return new i(-l,0)},function(l){return new i(0,-l)}];function o(s,l){if(l<0)throw new Error("Root must be greater than zero");if(l===0)throw new Error("Root must be non-zero");if(l%1!==0)throw new Error("Root must be an integer");if(s===0||s.abs()===0)return[new i(0,0)];var u=typeof s=="number",c;(u||s.re===0||s.im===0)&&(u?c=2*+(s<0):s.im===0?c=2*+(s.re<0):c=2*+(s.im<0)+1);for(var f=s.arg(),m=s.abs(),p=[],d=Math.pow(m,1/l),h=0;h{var{typed:e,equalScalar:r,matrix:n,pow:i,DenseMatrix:a,concat:o,SparseMatrix:s}=t,l=jr({typed:e}),u=Gn({typed:e,SparseMatrix:s}),c=Pr({typed:e,equalScalar:r}),f=wr({typed:e,DenseMatrix:a}),m=Pt({typed:e,matrix:n,concat:o}),p={};for(var d in i.signatures)Object.prototype.hasOwnProperty.call(i.signatures,d)&&!d.includes("Matrix")&&!d.includes("Array")&&(p[d]=i.signatures[d]);var h=e(p);return e(Yk,m({elop:h,SS:u,DS:l,Ss:c,sS:f}))});var jk="dotDivide",gee=["typed","matrix","equalScalar","divideScalar","DenseMatrix","concat","SparseMatrix"],Q1=M(jk,gee,t=>{var{typed:e,matrix:r,equalScalar:n,divideScalar:i,DenseMatrix:a,concat:o,SparseMatrix:s}=t,l=In({typed:e,equalScalar:n}),u=jr({typed:e}),c=Gn({typed:e,SparseMatrix:s}),f=Pr({typed:e,equalScalar:n}),m=wr({typed:e,DenseMatrix:a}),p=Pt({typed:e,matrix:r,concat:o});return e(jk,p({elop:i,SS:c,DS:u,SD:l,Ss:f,sS:m}))});function gs(t){var{DenseMatrix:e}=t;return function(n,i,a){var o=n.size();if(o.length!==2)throw new RangeError("Matrix must be two dimensional (size: "+Dt(o)+")");var s=o[0],l=o[1];if(s!==l)throw new RangeError("Matrix must be square (size: "+Dt(o)+")");var u=[];if(rt(i)){var c=i.size(),f=i._data;if(c.length===1){if(c[0]!==s)throw new RangeError("Dimension mismatch. Matrix columns must match vector length.");for(var m=0;m{var{typed:e,matrix:r,divideScalar:n,multiplyScalar:i,subtractScalar:a,equalScalar:o,DenseMatrix:s}=t,l=gs({DenseMatrix:s});return e(Zk,{"SparseMatrix, Array | Matrix":function(m,p){return c(m,p)},"DenseMatrix, Array | Matrix":function(m,p){return u(m,p)},"Array, Array | Matrix":function(m,p){var d=r(m),h=u(d,p);return h.valueOf()}});function u(f,m){m=l(f,m,!0);for(var p=m._data,d=f._size[0],h=f._size[1],x=[],b=f._data,T=0;TA&&(E.push(x[C]),D.push(B))}if(o(y,0))throw new Error("Linear system cannot be solved since matrix is singular");for(var R=n(w,y),V=0,L=D.length;V{var{typed:e,matrix:r,divideScalar:n,multiplyScalar:i,subtractScalar:a,equalScalar:o,DenseMatrix:s}=t,l=gs({DenseMatrix:s});return e(Jk,{"SparseMatrix, Array | Matrix":function(m,p){return c(m,p)},"DenseMatrix, Array | Matrix":function(m,p){return u(m,p)},"Array, Array | Matrix":function(m,p){var d=r(m),h=u(d,p);return h.valueOf()}});function u(f,m){m=l(f,m,!0);for(var p=m._data,d=f._size[0],h=f._size[1],x=[],b=f._data,T=h-1;T>=0;T--){var v=p[T][0]||0,A=void 0;if(o(v,0))A=0;else{var w=b[T][T];if(o(w,0))throw new Error("Linear system cannot be solved since matrix is singular");A=n(v,w);for(var y=T-1;y>=0;y--)p[y]=[a(p[y][0]||0,i(A,b[y][T]))]}x[T]=[A]}return new s({data:x,size:[d,1]})}function c(f,m){m=l(f,m,!0);for(var p=m._data,d=f._size[0],h=f._size[1],x=f._values,b=f._index,T=f._ptr,v=[],A=h-1;A>=0;A--){var w=p[A][0]||0;if(o(w,0))v[A]=[0];else{for(var y=0,E=[],D=[],_=T[A],N=T[A+1],C=N-1;C>=_;C--){var B=b[C];B===A?y=x[C]:B{var{typed:e,matrix:r,divideScalar:n,multiplyScalar:i,subtractScalar:a,equalScalar:o,DenseMatrix:s}=t,l=gs({DenseMatrix:s});return e(Xk,{"SparseMatrix, Array | Matrix":function(m,p){return c(m,p)},"DenseMatrix, Array | Matrix":function(m,p){return u(m,p)},"Array, Array | Matrix":function(m,p){var d=r(m),h=u(d,p);return h.map(x=>x.valueOf())}});function u(f,m){for(var p=[l(f,m,!0)._data.map(D=>D[0])],d=f._data,h=f._size[0],x=f._size[1],b=0;bnew s({data:D.map(_=>[_]),size:[h,1]}))}function c(f,m){for(var p=[l(f,m,!0)._data.map(ae=>ae[0])],d=f._size[0],h=f._size[1],x=f._values,b=f._index,T=f._ptr,v=0;vv&&(E.push(x[B]),D.push(R))}if(o(C,0))if(o(y[v],0)){if(w===0){var k=[...y];k[v]=1;for(var W=0,O=D.length;Wnew s({data:ae.map(J=>[J]),size:[d,1]}))}});var Qk="usolveAll",bee=["typed","matrix","divideScalar","multiplyScalar","subtractScalar","equalScalar","DenseMatrix"],rA=M(Qk,bee,t=>{var{typed:e,matrix:r,divideScalar:n,multiplyScalar:i,subtractScalar:a,equalScalar:o,DenseMatrix:s}=t,l=gs({DenseMatrix:s});return e(Qk,{"SparseMatrix, Array | Matrix":function(m,p){return c(m,p)},"DenseMatrix, Array | Matrix":function(m,p){return u(m,p)},"Array, Array | Matrix":function(m,p){var d=r(m),h=u(d,p);return h.map(x=>x.valueOf())}});function u(f,m){for(var p=[l(f,m,!0)._data.map(D=>D[0])],d=f._data,h=f._size[0],x=f._size[1],b=x-1;b>=0;b--)for(var T=p.length,v=0;v=0;E--)y[E]=a(y[E],d[E][b]);p.push(y)}}else{if(v===0)return[];p.splice(v,1),v-=1,T-=1}else{A[b]=n(A[b],d[b][b]);for(var w=b-1;w>=0;w--)A[w]=a(A[w],i(A[b],d[w][b]))}}return p.map(D=>new s({data:D.map(_=>[_]),size:[h,1]}))}function c(f,m){for(var p=[l(f,m,!0)._data.map(ae=>ae[0])],d=f._size[0],h=f._size[1],x=f._values,b=f._index,T=f._ptr,v=h-1;v>=0;v--)for(var A=p.length,w=0;w=_;B--){var R=b[B];R===v?C=x[B]:Rnew s({data:ae.map(J=>[J]),size:[d,1]}))}});var wee="matAlgo08xS0Sid",_ee=["typed","equalScalar"],Nf=M(wee,_ee,t=>{var{typed:e,equalScalar:r}=t;return function(i,a,o){var s=i._values,l=i._index,u=i._ptr,c=i._size,f=i._datatype||i._data===void 0?i._datatype:i.getDataType(),m=a._values,p=a._index,d=a._ptr,h=a._size,x=a._datatype||a._data===void 0?a._datatype:a.getDataType();if(c.length!==h.length)throw new Ct(c.length,h.length);if(c[0]!==h[0]||c[1]!==h[1])throw new RangeError("Dimension mismatch. Matrix A ("+c+") must match Matrix B ("+h+")");if(!s||!m)throw new Error("Cannot perform operation on Pattern Sparse Matrices");var b=c[0],T=c[1],v,A=r,w=0,y=o;typeof f=="string"&&f===x&&f!=="mixed"&&(v=f,A=e.find(r,[v,v]),w=e.convert(0,v),y=e.find(o,[v,v]));for(var E=[],D=[],_=[],N=[],C=[],B,R,V,L,P=0;P{var{typed:e,matrix:r}=t;return{"Array, number":e.referTo("DenseMatrix, number",n=>(i,a)=>n(r(i),a).valueOf()),"Array, BigNumber":e.referTo("DenseMatrix, BigNumber",n=>(i,a)=>n(r(i),a).valueOf()),"number, Array":e.referTo("number, DenseMatrix",n=>(i,a)=>n(i,r(a)).valueOf()),"BigNumber, Array":e.referTo("BigNumber, DenseMatrix",n=>(i,a)=>n(i,r(a)).valueOf())}});var Kk="leftShift",Aee=["typed","matrix","equalScalar","zeros","DenseMatrix","concat"],nA=M(Kk,Aee,t=>{var{typed:e,matrix:r,equalScalar:n,zeros:i,DenseMatrix:a,concat:o}=t,s=Ri({typed:e}),l=In({typed:e,equalScalar:n}),u=Nf({typed:e,equalScalar:n}),c=Sa({typed:e,DenseMatrix:a}),f=Pr({typed:e,equalScalar:n}),m=Fn({typed:e}),p=Pt({typed:e,matrix:r,concat:o}),d=Ff({typed:e,matrix:r});return e(Kk,{"number, number":Bw,"BigNumber, BigNumber":zB,"bigint, bigint":(h,x)=>h<(x,b)=>n(b,0)?x.clone():f(x,b,h,!1)),"DenseMatrix, number | BigNumber":e.referToSelf(h=>(x,b)=>n(b,0)?x.clone():m(x,b,h,!1)),"number | BigNumber, SparseMatrix":e.referToSelf(h=>(x,b)=>n(x,0)?i(b.size(),b.storage()):c(b,x,h,!0)),"number | BigNumber, DenseMatrix":e.referToSelf(h=>(x,b)=>n(x,0)?i(b.size(),b.storage()):m(b,x,h,!0))},d,p({SS:u,DS:s,SD:l}))});var e$="rightArithShift",See=["typed","matrix","equalScalar","zeros","DenseMatrix","concat"],iA=M(e$,See,t=>{var{typed:e,matrix:r,equalScalar:n,zeros:i,DenseMatrix:a,concat:o}=t,s=Ri({typed:e}),l=In({typed:e,equalScalar:n}),u=Nf({typed:e,equalScalar:n}),c=Sa({typed:e,DenseMatrix:a}),f=Pr({typed:e,equalScalar:n}),m=Fn({typed:e}),p=Pt({typed:e,matrix:r,concat:o}),d=Ff({typed:e,matrix:r});return e(e$,{"number, number":kw,"BigNumber, BigNumber":HB,"bigint, bigint":(h,x)=>h>>x,"SparseMatrix, number | BigNumber":e.referToSelf(h=>(x,b)=>n(b,0)?x.clone():f(x,b,h,!1)),"DenseMatrix, number | BigNumber":e.referToSelf(h=>(x,b)=>n(b,0)?x.clone():m(x,b,h,!1)),"number | BigNumber, SparseMatrix":e.referToSelf(h=>(x,b)=>n(x,0)?i(b.size(),b.storage()):c(b,x,h,!0)),"number | BigNumber, DenseMatrix":e.referToSelf(h=>(x,b)=>n(x,0)?i(b.size(),b.storage()):m(b,x,h,!0))},d,p({SS:u,DS:s,SD:l}))});var t$="rightLogShift",Eee=["typed","matrix","equalScalar","zeros","DenseMatrix","concat"],aA=M(t$,Eee,t=>{var{typed:e,matrix:r,equalScalar:n,zeros:i,DenseMatrix:a,concat:o}=t,s=Ri({typed:e}),l=In({typed:e,equalScalar:n}),u=Nf({typed:e,equalScalar:n}),c=Sa({typed:e,DenseMatrix:a}),f=Pr({typed:e,equalScalar:n}),m=Fn({typed:e}),p=Pt({typed:e,matrix:r,concat:o}),d=Ff({typed:e,matrix:r});return e(t$,{"number, number":$w,"SparseMatrix, number | BigNumber":e.referToSelf(h=>(x,b)=>n(b,0)?x.clone():f(x,b,h,!1)),"DenseMatrix, number | BigNumber":e.referToSelf(h=>(x,b)=>n(b,0)?x.clone():m(x,b,h,!1)),"number | BigNumber, SparseMatrix":e.referToSelf(h=>(x,b)=>n(x,0)?i(b.size(),b.storage()):c(b,x,h,!0)),"number | BigNumber, DenseMatrix":e.referToSelf(h=>(x,b)=>n(x,0)?i(b.size(),b.storage()):m(b,x,h,!0))},d,p({SS:u,DS:s,SD:l}))});var r$="and",Tee=["typed","matrix","equalScalar","zeros","not","concat"],bm=M(r$,Tee,t=>{var{typed:e,matrix:r,equalScalar:n,zeros:i,not:a,concat:o}=t,s=In({typed:e,equalScalar:n}),l=Gu({typed:e,equalScalar:n}),u=Pr({typed:e,equalScalar:n}),c=Fn({typed:e}),f=Pt({typed:e,matrix:r,concat:o});return e(r$,{"number, number":Qg,"Complex, Complex":function(p,d){return(p.re!==0||p.im!==0)&&(d.re!==0||d.im!==0)},"BigNumber, BigNumber":function(p,d){return!p.isZero()&&!d.isZero()&&!p.isNaN()&&!d.isNaN()},"bigint, bigint":Qg,"Unit, Unit":e.referToSelf(m=>(p,d)=>m(p.value||0,d.value||0)),"SparseMatrix, any":e.referToSelf(m=>(p,d)=>a(d)?i(p.size(),p.storage()):u(p,d,m,!1)),"DenseMatrix, any":e.referToSelf(m=>(p,d)=>a(d)?i(p.size(),p.storage()):c(p,d,m,!1)),"any, SparseMatrix":e.referToSelf(m=>(p,d)=>a(p)?i(p.size(),p.storage()):u(d,p,m,!0)),"any, DenseMatrix":e.referToSelf(m=>(p,d)=>a(p)?i(p.size(),p.storage()):c(d,p,m,!0)),"Array, any":e.referToSelf(m=>(p,d)=>m(r(p),d).valueOf()),"any, Array":e.referToSelf(m=>(p,d)=>m(p,r(d)).valueOf())},f({SS:l,DS:s}))});var ov="compare",Dee=["typed","config","matrix","equalScalar","BigNumber","Fraction","DenseMatrix","concat"],oA=M(ov,Dee,t=>{var{typed:e,config:r,equalScalar:n,matrix:i,BigNumber:a,Fraction:o,DenseMatrix:s,concat:l}=t,u=jr({typed:e}),c=Hu({typed:e,equalScalar:n}),f=wr({typed:e,DenseMatrix:s}),m=Pt({typed:e,matrix:i,concat:l}),p=_o({typed:e});return e(ov,Mee({typed:e,config:r}),{"boolean, boolean":function(h,x){return h===x?0:h>x?1:-1},"BigNumber, BigNumber":function(h,x){return mn(h,x,r.relTol,r.absTol)?new a(0):new a(h.cmp(x))},"bigint, bigint":function(h,x){return h===x?0n:h>x?1n:-1n},"Fraction, Fraction":function(h,x){return new o(h.compare(x))},"Complex, Complex":function(){throw new TypeError("No ordering relation is defined for complex numbers")}},p,m({SS:c,DS:u,Ss:f}))}),Mee=M(ov,["typed","config"],t=>{var{typed:e,config:r}=t;return e(ov,{"number, number":function(i,a){return Kr(i,a,r.relTol,r.absTol)?0:i>a?1:-1}})});var Xl=zr(i$(),1);var a$="compareNatural",Cee=["typed","compare"],sA=M(a$,Cee,t=>{var{typed:e,compare:r}=t,n=r.signatures["boolean,boolean"];return e(a$,{"any, any":i});function i(l,u){var c=fr(l),f=fr(u),m;if((c==="number"||c==="BigNumber"||c==="Fraction")&&(f==="number"||f==="BigNumber"||f==="Fraction"))return m=r(l,u),m.toString()!=="0"?m>0?1:-1:(0,Xl.default)(c,f);var p=["Array","DenseMatrix","SparseMatrix"];if(p.includes(c)||p.includes(f))return m=a(i,l,u),m!==0?m:(0,Xl.default)(c,f);if(c!==f)return(0,Xl.default)(c,f);if(c==="Complex")return Nee(l,u);if(c==="Unit")return l.equalBase(u)?i(l.value,u.value):o(i,l.formatUnits(),u.formatUnits());if(c==="boolean")return n(l,u);if(c==="string")return(0,Xl.default)(l,u);if(c==="Object")return s(i,l,u);if(c==="null"||c==="undefined")return 0;throw new TypeError('Unsupported type of value "'+c+'"')}function a(l,u,c){return Uo(u)&&Uo(c)?o(l,u.toJSON().values,c.toJSON().values):Uo(u)?a(l,u.toArray(),c):Uo(c)?a(l,u,c.toArray()):Ul(u)?a(l,u.toJSON().data,c):Ul(c)?a(l,u,c.toJSON().data):Array.isArray(u)?Array.isArray(c)?o(l,u,c):a(l,u,[c]):a(l,[u],c)}function o(l,u,c){for(var f=0,m=Math.min(u.length,c.length);fc.length?1:u.lengthe.re?1:t.ree.im?1:t.im{var{typed:e,matrix:r,concat:n}=t,i=Pt({typed:e,matrix:r,concat:n});return e(o$,Hg,i({elop:Hg,Ds:!0}))});var sv="equal",Iee=["typed","matrix","equalScalar","DenseMatrix","concat","SparseMatrix"],lA=M(sv,Iee,t=>{var{typed:e,matrix:r,equalScalar:n,DenseMatrix:i,concat:a,SparseMatrix:o}=t,s=jr({typed:e}),l=Gn({typed:e,SparseMatrix:o}),u=wr({typed:e,DenseMatrix:i}),c=Pt({typed:e,matrix:r,concat:a});return e(sv,Pee({typed:e,equalScalar:n}),c({elop:n,SS:l,DS:s,Ss:u}))}),Pee=M(sv,["typed","equalScalar"],t=>{var{typed:e,equalScalar:r}=t;return e(sv,{"any, any":function(i,a){return i===null?a===null:a===null?i===null:i===void 0?a===void 0:a===void 0?i===void 0:r(i,a)}})});var s$="equalText",Bee=["typed","compareText","isZero"],cA=M(s$,Bee,t=>{var{typed:e,compareText:r,isZero:n}=t;return e(s$,{"any, any":function(a,o){return n(r(a,o))}})});var uv="smaller",kee=["typed","config","matrix","DenseMatrix","concat","SparseMatrix"],fA=M(uv,kee,t=>{var{typed:e,config:r,matrix:n,DenseMatrix:i,concat:a,SparseMatrix:o}=t,s=jr({typed:e}),l=Gn({typed:e,SparseMatrix:o}),u=wr({typed:e,DenseMatrix:i}),c=Pt({typed:e,matrix:n,concat:a}),f=_o({typed:e});return e(uv,$ee({typed:e,config:r}),{"boolean, boolean":(m,p)=>mm
m.compare(p)===-1,"Complex, Complex":function(p,d){throw new TypeError("No ordering relation is defined for complex numbers")}},f,c({SS:l,DS:s,Ss:u}))}),$ee=M(uv,["typed","config"],t=>{var{typed:e,config:r}=t;return e(uv,{"number, number":function(i,a){return i{var{typed:e,config:r,matrix:n,DenseMatrix:i,concat:a,SparseMatrix:o}=t,s=jr({typed:e}),l=Gn({typed:e,SparseMatrix:o}),u=wr({typed:e,DenseMatrix:i}),c=Pt({typed:e,matrix:n,concat:a}),f=_o({typed:e});return e(lv,Ree({typed:e,config:r}),{"boolean, boolean":(m,p)=>m<=p,"BigNumber, BigNumber":function(p,d){return p.lte(d)||mn(p,d,r.relTol,r.absTol)},"bigint, bigint":(m,p)=>m<=p,"Fraction, Fraction":(m,p)=>m.compare(p)!==1,"Complex, Complex":function(){throw new TypeError("No ordering relation is defined for complex numbers")}},f,c({SS:l,DS:s,Ss:u}))}),Ree=M(lv,["typed","config"],t=>{var{typed:e,config:r}=t;return e(lv,{"number, number":function(i,a){return i<=a||Kr(i,a,r.relTol,r.absTol)}})});var cv="larger",Lee=["typed","config","matrix","DenseMatrix","concat","SparseMatrix"],mA=M(cv,Lee,t=>{var{typed:e,config:r,matrix:n,DenseMatrix:i,concat:a,SparseMatrix:o}=t,s=jr({typed:e}),l=Gn({typed:e,SparseMatrix:o}),u=wr({typed:e,DenseMatrix:i}),c=Pt({typed:e,matrix:n,concat:a}),f=_o({typed:e});return e(cv,Vee({typed:e,config:r}),{"boolean, boolean":(m,p)=>m>p,"BigNumber, BigNumber":function(p,d){return p.gt(d)&&!mn(p,d,r.relTol,r.absTol)},"bigint, bigint":(m,p)=>m>p,"Fraction, Fraction":(m,p)=>m.compare(p)===1,"Complex, Complex":function(){throw new TypeError("No ordering relation is defined for complex numbers")}},f,c({SS:l,DS:s,Ss:u}))}),Vee=M(cv,["typed","config"],t=>{var{typed:e,config:r}=t;return e(cv,{"number, number":function(i,a){return i>a&&!Kr(i,a,r.relTol,r.absTol)}})});var fv="largerEq",Uee=["typed","config","matrix","DenseMatrix","concat","SparseMatrix"],dA=M(fv,Uee,t=>{var{typed:e,config:r,matrix:n,DenseMatrix:i,concat:a,SparseMatrix:o}=t,s=jr({typed:e}),l=Gn({typed:e,SparseMatrix:o}),u=wr({typed:e,DenseMatrix:i}),c=Pt({typed:e,matrix:n,concat:a}),f=_o({typed:e});return e(fv,qee({typed:e,config:r}),{"boolean, boolean":(m,p)=>m>=p,"BigNumber, BigNumber":function(p,d){return p.gte(d)||mn(p,d,r.relTol,r.absTol)},"bigint, bigint":function(p,d){return p>=d},"Fraction, Fraction":(m,p)=>m.compare(p)!==-1,"Complex, Complex":function(){throw new TypeError("No ordering relation is defined for complex numbers")}},f,c({SS:l,DS:s,Ss:u}))}),qee=M(fv,["typed","config"],t=>{var{typed:e,config:r}=t;return e(fv,{"number, number":function(i,a){return i>=a||Kr(i,a,r.relTol,r.absTol)}})});var u$="deepEqual",zee=["typed","equal"],hA=M(u$,zee,t=>{var{typed:e,equal:r}=t;return e(u$,{"any, any":function(a,o){return n(a.valueOf(),o.valueOf())}});function n(i,a){if(Array.isArray(i))if(Array.isArray(a)){var o=i.length;if(o!==a.length)return!1;for(var s=0;s{var{typed:e,config:r,equalScalar:n,matrix:i,DenseMatrix:a,concat:o,SparseMatrix:s}=t,l=jr({typed:e}),u=Gn({typed:e,SparseMatrix:s}),c=wr({typed:e,DenseMatrix:a}),f=Pt({typed:e,matrix:i,concat:o});return e(pv,Gee({typed:e,equalScalar:n}),f({elop:m,SS:u,DS:l,Ss:c}));function m(p,d){return!n(p,d)}}),Gee=M(pv,["typed","equalScalar"],t=>{var{typed:e,equalScalar:r}=t;return e(pv,{"any, any":function(i,a){return i===null?a!==null:a===null?i!==null:i===void 0?a!==void 0:a===void 0?i!==void 0:!r(i,a)}})});var l$="partitionSelect",Wee=["typed","isNumeric","isNaN","compare"],vA=M(l$,Wee,t=>{var{typed:e,isNumeric:r,isNaN:n,compare:i}=t,a=i,o=(u,c)=>-i(u,c);return e(l$,{"Array | Matrix, number":function(c,f){return s(c,f,a)},"Array | Matrix, number, string":function(c,f,m){if(m==="asc")return s(c,f,a);if(m==="desc")return s(c,f,o);throw new Error('Compare string must be "asc" or "desc"')},"Array | Matrix, number, function":s});function s(u,c,f){if(!Qe(c)||c<0)throw new Error("k must be a non-negative integer");if(rt(u)){var m=u.size();if(m.length>1)throw new Error("Only one dimensional matrices supported");return l(u.valueOf(),c,f)}if(Array.isArray(u))return l(u,c,f)}function l(u,c,f){if(c>=u.length)throw new Error("k out of bounds");for(var m=0;m