| 100 | R: `+this.R.toString()}}}function A(g){var w=v(g),b=w.R._data;if(g._data.length>0)for(var y=b[0][0].type==="Complex"?m(0):0,S=0;S<b.length;++S)for(var T=0;T<S&&T<(b[0]||[]).length;++T)b[S][T]=y;return w}function x(g){throw new Error("qr not implemented for sparse matrices yet")}});function Hbe(e,t,r,n){for(var i=e._values,a=e._index,o=e._ptr,s=e._size,l=e._datatype,u=s[0],c=s[1],f=n&&e._values?[]:null,d=[],p=[],h=0,m=0;m<c;m++){p[m]=h;for(var v=r?r[m]:m,A=o[v],x=o[v+1],g=A;g<x;g++){var w=t?t[a[g]]:a[g];d[h]=w,f&&(f[h]=i[g]),h++}}return p[c]=h,e.createSparseMatrix({values:f,index:d,ptr:p,size:[u,c],datatype:l})}function h7(e,t,r,n,i,a,o){var s=0;for(r[o]=e;s>=0;){var l=r[o+s],u=r[n+l];u===-1?(s--,a[t++]=l):(r[n+l]=r[i+u],++s,r[o+s]=u)}return t}function Vbe(e,t){if(!e)return null;var r=0,n,i=[],a=[],o=0,s=t,l=2*t;for(n=0;n<t;n++)a[o+n]=-1;for(n=t-1;n>=0;n--)e[n]!==-1&&(a[s+n]=a[o+e[n]],a[o+e[n]]=n);for(n=0;n<t;n++)e[n]===-1&&(r=h7(n,r,a,o,s,i,l));return i}function Gbe(e,t){if(!e)return null;var r=e._index,n=e._ptr,i=e._size,a=i[0],o=i[1],s=[],l=[],u=0,c=o,f,d;if(t)for(f=0;f<a;f++)l[c+f]=-1;for(var p=0;p<o;p++){s[p]=-1,l[u+p]=-1;for(var h=n[p],m=n[p+1],v=h;v<m;v++){var A=r[v];for(f=t?l[c+A]:A;f!==-1&&f<p;f=d)d=l[u+f],l[u+f]=p,d===-1&&(s[f]=p);t&&(l[c+A]=p)}}return s}function Wbe(e,t,r){for(var n=e._values,i=e._index,a=e._ptr,o=e._size,s=o[1],l=0,u=0;u<s;u++){var c=a[u];for(a[u]=l;c<a[u+1];c++)t(i[c],u,n?n[c]:1,r)&&(i[l]=i[c],n&&(n[l]=n[c]),l++)}return a[s]=l,i.splice(l,i.length-l),n&&n.splice(l,n.length-l),l}function ih(e){return-e-2}var Tut="csAmd",Nut=["add","multiply","transpose"],Qbe=X(Tut,Nut,e=>{var{add:t,multiply:r,transpose:n}=e;return function(c,f){if(!f||c<=0||c>3)return null;var d=f._size,p=d[0],h=d[1],m=0,v=Math.max(16,10*Math.sqrt(h));v=Math.min(h-2,v);var A=i(c,f,p,h,v);Wbe(A,l,null);for(var x=A._index,g=A._ptr,w=g[h],b=[],y=[],S=0,T=h+1,D=2*(h+1),k=3*(h+1),M=4*(h+1),E=5*(h+1),O=6*(h+1),P=7*(h+1),_=b,B=a(h,g,y,S,k,_,D,P,T,O,M,E),F=o(h,g,y,E,M,O,v,T,k,_,D),L=0,I,R,J,Q,q,ie,Ae,le,pe,ae,oe,ue,Se,H,ye,ne;F<h;){for(J=-1;L<h&&(J=y[k+L])===-1;L++);y[D+J]!==-1&&(_[y[D+J]]=-1),y[k+L]=y[D+J];var z=y[M+J],U=y[T+J];F+=U;var V=0;y[T+J]=-U;var Y=g[J],Z=z===0?Y:w,re=Z;for(Q=1;Q<=z+1;Q++){for(Q>z?(ie=J,Ae=Y,le=y[S+J]-z):(ie=x[Y++],Ae=g[ie],le=y[S+ie]),q=1;q<=le;q++)I=x[Ae++],!((pe=y[T+I])<=0)&&(V+=pe,y[T+I]=-pe,x[re++]=I,y[D+I]!==-1&&(_[y[D+I]]=_[I]),_[I]!==-1?y[D+_[I]]=y[D+I]:y[k+y[E+I]]=y[D+I]);ie!==J&&(g[ie]=ih(J),y[O+ie]=0)}for(z!==0&&(w=re),y[E+J]=V,g[J]=Z,y[S+J]=re-Z,y[M+J]=-2,B=s(B,m,y,O,h),ae=Z;ae<re;ae++)if(I=x[ae],!((oe=y[M+I])<=0)){pe=-y[T+I];var he=B-pe;for(Y=g[I],ue=g[I]+oe-1;Y<=ue;Y++)ie=x[Y],y[O+ie]>=B?y[O+ie]-=pe:y[O+ie]!==0&&(y[O+ie]=y[E+ie]+he)}for(ae=Z;ae<re;ae++){for(I=x[ae],ue=g[I],Se=ue+y[M+I]-1,H=ue,ye=0,ne=0,Y=ue;Y<=Se;Y++)if(ie=x[Y],y[O+ie]!==0){var fe=y[O+ie]-B;fe>0?(ne+=fe,x[H++]=ie,ye+=ie):(g[ie]=ih(J),y[O+ie]=0)}y[M+I]=H-ue+1;var ke=H,se=ue+y[S+I];for(Y=Se+1;Y<se;Y++){R=x[Y];var Le=y[T+R];Le<=0||(ne+=Le,x[H++]=R,ye+=R)}ne===0?(g[I]=ih(J),pe=-y[T+I],V-=pe,U+=pe,F+=pe,y[T+I]=0,y[M+I]=-1):(y[E+I]=Math.min(y[E+I],ne),x[H]=x[ke],x[ke]=x[ue],x[ue]=J,y[S+I]=H-ue+1,ye=(ye<0?-ye:ye)%h,y[D+I]=y[P+ye],y[P+ye]=I,_[I]=ye)}for(y[E+J]=V,m=Math.max(m,V),B=s(B+m,m,y,O,h),ae=Z;ae<re;ae++)if(I=x[ae],!(y[T+I]>=0))for(ye=_[I],I=y[P+ye],y[P+ye]=-1;I!==-1&&y[D+I]!==-1;I=y[D+I],B++){for(le=y[S+I],oe=y[M+I],Y=g[I]+1;Y<=g[I]+le-1;Y++)y[O+x[Y]]=B;var Ne=I;for(R=y[D+I];R!==-1;){var Ge=y[S+R]===le&&y[M+R]===oe;for(Y=g[R]+1;Ge&&Y<=g[R]+le-1;Y++)y[O+x[Y]]!==B&&(Ge=0);Ge?(g[R]=ih(I),y[T+I]+=y[T+R],y[T+R]=0,y[M+R]=-1,R=y[D+R],y[D+Ne]=R):(Ne=R,R=y[D+R])}}for(Y=Z,ae=Z;ae<re;ae++)I=x[ae],!((pe=-y[T+I])<=0)&&(y[T+I]=pe,ne=y[E+I]+V-pe,ne=Math.min(ne,h-F-pe),y[k+ne]!==-1&&(_[y[k+ne]]=I),y[D+I]=y[k+ne],_[I]=-1,y[k+ne]=I,L=Math.min(L,ne),y[E+I]=ne,x[Y++]=I);y[T+J]=U,(y[S+J]=Y-Z)===0&&(g[J]=-1,y[O+J]=0),z!==0&&(w=Y)}for(I=0;I<h;I++)g[I]=ih(g[I]);for(R=0;R<=h;R++)y[k+R]=-1;for(R=h;R>=0;R--)y[T+R]>0||(y[D+R]=y[k+g[R]],y[k+g[R]]=R);for(ie=h;ie>=0;ie--)y[T+ie]<=0||g[ie]!==-1&&(y[D+ie]=y[k+g[ie]],y[k+g[ie]]=ie);for(J=0,I=0;I<=h;I++)g[I]===-1&&(J=h7(I,J,y,k,D,b,O));return b.splice(b.length-1,1),b};function i(u,c,f,d,p){var h=n(c);if(u===1&&d===f)return t(c,h);if(u===2){for(var m=h._index,v=h._ptr,A=0,x=0;x<f;x++){var g=v[x];if(v[x]=A,!(v[x+1]-g>p))for(var w=v[x+1];g<w;g++)m[A++]=m[g]}return v[f]=A,c=n(h),r(h,c)}return r(h,c)}function a(u,c,f,d,p,h,m,v,A,x,g,w){for(var b=0;b<u;b++)f[d+b]=c[b+1]-c[b];f[d+u]=0;for(var y=0;y<=u;y++)f[p+y]=-1,h[y]=-1,f[m+y]=-1,f[v+y]=-1,f[A+y]=1,f[x+y]=1,f[g+y]=0,f[w+y]=f[d+y];var S=s(0,0,f,x,u);return f[g+u]=-2,c[u]=-1,f[x+u]=0,S}function o(u,c,f,d,p,h,m,v,A,x,g){for(var w=0,b=0;b<u;b++){var y=f[d+b];if(y===0)f[p+b]=-2,w++,c[b]=-1,f[h+b]=0;else if(y>m)f[v+b]=0,f[p+b]=-1,w++,c[b]=ih(u),f[v+u]++;else{var S=f[A+y];S!==-1&&(x[S]=b),f[g+b]=f[A+y],f[A+y]=b}}return w}function s(u,c,f,d,p){if(u<2||u+c<0){for(var h=0;h<p;h++)f[d+h]!==0&&(f[d+h]=1);u=2}return u}function l(u,c){return u!==c}});function Ybe(e,t,r,n,i,a,o){var s,l,u=0,c;if(e<=t||r[n+t]<=r[i+e])return-1;r[i+e]=r[n+t];var f=r[a+e];if(r[a+e]=t,f===-1)u=1,c=e;else{for(u=2,c=f;c!==r[o+c];c=r[o+c]);for(s=f;s!==c;s=l)l=r[o+s],r[o+s]=c}return{jleaf:u,q:c}}var Mut="csCounts",Out=["transpose"],Xbe=X(Mut,Out,e=>{var{transpose:t}=e;return function(r,n,i,a){if(!r||!n||!i)return null;var o=r._size,s=o[0],l=o[1],u,c,f,d,p,h,m,v=4*l+(a?l+s+1:0),A=[],x=0,g=l,w=2*l,b=3*l,y=4*l,S=5*l+1;for(f=0;f<v;f++)A[f]=-1;var T=[],D=t(r),k=D._index,M=D._ptr;for(f=0;f<l;f++)for(c=i[f],T[c]=A[b+c]===-1?1:0;c!==-1&&A[b+c]===-1;c=n[c])A[b+c]=f;if(a){for(f=0;f<l;f++)A[i[f]]=f;for(u=0;u<s;u++){for(f=l,h=M[u],m=M[u+1],p=h;p<m;p++)f=Math.min(f,A[k[p]]);A[S+u]=A[y+f],A[y+f]=u}}for(u=0;u<l;u++)A[x+u]=u;for(f=0;f<l;f++){for(c=i[f],n[c]!==-1&&T[n[c]]--,d=a?A[y+f]:c;d!==-1;d=a?A[S+d]:-1)for(p=M[d];p<M[d+1];p++){u=k[p];var E=Ybe(u,c,A,b,g,w,x);E.jleaf>=1&&T[c]++,E.jleaf===2&&T[E.q]--}n[c]!==-1&&(A[x+c]=n[c])}for(c=0;c<l;c++)n[c]!==-1&&(T[n[c]]+=T[c]);return T}});var Fut="csSqr",jut=["add","multiply","transpose"],Kbe=X(Fut,jut,e=>{var{add:t,multiply:r,transpose:n}=e,i=Qbe({add:t,multiply:r,transpose:n}),a=Xbe({transpose:n});return function(l,u,c){var f=u._ptr,d=u._size,p=d[1],h,m={};if(m.q=i(l,u),l&&!m.q)return null;if(c){var v=l?Hbe(u,null,m.q,0):u;m.parent=Gbe(v,1);var A=Vbe(m.parent,p);if(m.cp=a(v,m.parent,A,1),v&&m.parent&&m.cp&&o(v,m))for(m.unz=0,h=0;h<p;h++)m.unz+=m.cp[h]}else m.unz=4*f[p]+p,m.lnz=m.unz;return m};function o(s,l){var u=s._ptr,c=s._index,f=s._size,d=f[0],p=f[1];l.pinv=[],l.leftmost=[];var h=l.parent,m=l.pinv,v=l.leftmost,A=[],x=0,g=d,w=d+p,b=d+2*p,y,S,T,D,k;for(S=0;S<p;S++)A[g+S]=-1,A[w+S]=-1,A[b+S]=0;for(y=0;y<d;y++)v[y]=-1;for(S=p-1;S>=0;S--)for(D=u[S],k=u[S+1],T=D;T<k;T++)v[c[T]]=S;for(y=d-1;y>=0;y--)m[y]=-1,S=v[y],S!==-1&&(A[b+S]++===0&&(A[w+S]=y),A[x+y]=A[g+S],A[g+S]=y);for(l.lnz=0,l.m2=d,S=0;S<p;S++)if(y=A[g+S],l.lnz++,y<0&&(y=l.m2++),m[y]=S,!(--b[S]<=0)){l.lnz+=A[b+S];var M=h[S];M!==-1&&(A[b+M]===0&&(A[w+M]=A[w+S]),A[x+A[w+S]]=A[g+M],A[g+M]=A[x+y],A[b+M]+=A[b+S])}for(y=0;y<d;y++)m[y]<0&&(m[y]=S++);return!0}});function b3(e,t){return e[t]<0}function g7(e,t){e[t]=ih(e[t])}function _Z(e){return e<0?ih(e):e}function Jbe(e,t,r,n,i){var a=t._index,o=t._ptr,s=t._size,l=s[1],u,c,f,d=0;for(n[0]=e;d>=0;){e=n[d];var p=i?i[e]:e;b3(o,e)||(g7(o,e),n[l+d]=p<0?0:_Z(o[p]));var h=1;for(c=n[l+d],f=p<0?0:_Z(o[p+1]);c<f;c++)if(u=a[c],!b3(o,u)){n[l+d]=c,n[++d]=u,h=0;break}h&&(d--,n[--r]=e)}return r}function Zbe(e,t,r,n,i){var a=e._ptr,o=e._size,s=t._index,l=t._ptr,u=o[1],c,f,d,p=u;for(f=l[r],d=l[r+1],c=f;c<d;c++){var h=s[c];b3(a,h)||(p=Jbe(h,e,p,n,i))}for(c=p;c<u;c++)g7(a,n[c]);return p}var Put="csSpsolve",_ut=["divideScalar","multiply","subtract"],exe=X(Put,_ut,e=>{var{divideScalar:t,multiply:r,subtract:n}=e;return function(a,o,s,l,u,c,f){var d=a._values,p=a._index,h=a._ptr,m=a._size,v=m[1],A=o._values,x=o._index,g=o._ptr,w,b,y,S,T=Zbe(a,o,s,l,c);for(w=T;w<v;w++)u[l[w]]=0;for(b=g[s],y=g[s+1],w=b;w<y;w++)u[x[w]]=A[w];for(var D=T;D<v;D++){var k=l[D],M=c?c[k]:k;if(!(M<0))for(b=h[M],y=h[M+1],u[k]=t(u[k],d[f?b:y-1]),w=f?b+1:b,S=f?y:y-1;w<S;w++){var E=p[w];u[E]=n(u[E],r(d[w],u[k]))}}return T}});var Iut="csLu",But=["abs","divideScalar","multiply","subtract","larger","largerEq","SparseMatrix"],txe=X(Iut,But,e=>{var{abs:t,divideScalar:r,multiply:n,subtract:i,larger:a,largerEq:o,SparseMatrix:s}=e,l=exe({divideScalar:r,multiply:n,subtract:i});return function(c,f,d){if(!c)return null;var p=c._size,h=p[1],m,v=100,A=100;f&&(m=f.q,v=f.lnz||v,A=f.unz||A);var x=[],g=[],w=[],b=new s({values:x,index:g,ptr:w,size:[h,h]}),y=[],S=[],T=[],D=new s({values:y,index:S,ptr:T,size:[h,h]}),k=[],M,E,O=[],P=[];for(M=0;M<h;M++)O[M]=0,k[M]=-1,w[M+1]=0;v=0,A=0;for(var _=0;_<h;_++){w[_]=v,T[_]=A;var B=m?m[_]:_,F=l(b,c,B,P,O,k,1),L=-1,I=-1;for(E=F;E<h;E++)if(M=P[E],k[M]<0){var R=t(O[M]);a(R,I)&&(I=R,L=M)}else S[A]=k[M],y[A++]=O[M];if(L===-1||I<=0)return null;k[B]<0&&o(t(O[B]),n(I,d))&&(L=B);var J=O[L];for(S[A]=_,y[A++]=J,k[L]=_,g[v]=L,x[v++]=1,E=F;E<h;E++)M=P[E],k[M]<0&&(g[v]=M,x[v++]=r(O[M],J)),O[M]=0}for(w[h]=v,T[h]=A,E=0;E<v;E++)g[E]=k[g[E]];return x.splice(v,x.length-v),g.splice(v,g.length-v),y.splice(A,y.length-A),S.splice(A,S.length-A),{L:b,U:D,pinv:k}}});var rxe="slu",Lut=["typed","abs","add","multiply","transpose","divideScalar","subtract","larger","largerEq","SparseMatrix"],A7=X(rxe,Lut,e=>{var{typed:t,abs:r,add:n,multiply:i,transpose:a,divideScalar:o,subtract:s,larger:l,largerEq:u,SparseMatrix:c}=e,f=Kbe({add:n,multiply:i,transpose:a}),d=txe({abs:r,divideScalar:o,multiply:i,subtract:s,larger:l,largerEq:u,SparseMatrix:c});return t(rxe,{"SparseMatrix, number, number":function(h,m,v){if(!Ot(m)||m<0||m>3)throw new Error("Symbolic Ordering and Analysis order must be an integer number in the interval [0, 3]");if(v<0||v>1)throw new Error("Partial pivoting threshold must be a number from 0 to 1");var A=f(m,h,!1),x=d(h,A,v);return{L:x.L,U:x.U,p:x.pinv,q:A.q,toString:function(){return"L: "+this.L.toString()+` |