package.umd.incircle.min.js Maven / Gradle / Ivy

Go to download
Show more of this group Show more artifacts with this name
Show all versions of robust-predicates Show documentation
Fast robust predicates for computational geometry
The newest version!
!function(t,n){"object"==typeof exports&&"undefined"!=typeof module?n(exports):"function"==typeof define&&define.amd?define(["exports"],n):n((t="undefined"!=typeof globalThis?globalThis:t||self).predicates={})}(this,(function(t){"use strict";const n=134217729;function e(t,n,e,o,f){let i,r,c,s,u=n[0],a=o[0],l=0,d=0;a>u==a>-u?(i=u,u=n[++l]):(i=a,a=o[++d]);let b=0;if(lu==a>-u?(r=u+i,c=i-(r-u),u=n[++l]):(r=a+i,c=i-(r-a),a=o[++d]),i=r,0!==c&&(f[b++]=c);lu==a>-u?(r=i+u,s=r-i,c=i-(r-s)+(u-s),u=n[++l]):(r=i+a,s=r-i,c=i-(r-s)+(a-s),a=o[++d]),i=r,0!==c&&(f[b++]=c);for(;lit||-ot>it?ot:function(t,i,D,G,H,I,J,K,L){let N,Q,R,S,U,V,W,X,Y,Z,$,tt,nt,et,ot,ft,it,rt,ct,st,ut,at,lt,dt,bt,ht,pt,Mt,yt,xt,gt,jt,mt,Tt,_t;const vt=t-J,wt=D-J,At=H-J,Ft=i-K,Ot=G-K,Pt=I-K;gt=wt*Pt,lt=n*wt,dt=lt-(lt-wt),bt=wt-dt,lt=n*Pt,ht=lt-(lt-Pt),pt=Pt-ht,jt=bt*pt-(gt-dt*ht-bt*ht-dt*pt),mt=At*Ot,lt=n*At,dt=lt-(lt-At),bt=At-dt,lt=n*Ot,ht=lt-(lt-Ot),pt=Ot-ht,Tt=bt*pt-(mt-dt*ht-bt*ht-dt*pt),Mt=jt-Tt,at=jt-Mt,r[0]=jt-(Mt+at)+(at-Tt),yt=gt+Mt,at=yt-gt,xt=gt-(yt-at)+(Mt-at),Mt=xt-mt,at=xt-Mt,r[1]=xt-(Mt+at)+(at-mt),_t=yt+Mt,at=_t-yt,r[2]=yt-(_t-at)+(Mt-at),r[3]=_t,gt=At*Ft,lt=n*At,dt=lt-(lt-At),bt=At-dt,lt=n*Ft,ht=lt-(lt-Ft),pt=Ft-ht,jt=bt*pt-(gt-dt*ht-bt*ht-dt*pt),mt=vt*Pt,lt=n*vt,dt=lt-(lt-vt),bt=vt-dt,lt=n*Pt,ht=lt-(lt-Pt),pt=Pt-ht,Tt=bt*pt-(mt-dt*ht-bt*ht-dt*pt),Mt=jt-Tt,at=jt-Mt,c[0]=jt-(Mt+at)+(at-Tt),yt=gt+Mt,at=yt-gt,xt=gt-(yt-at)+(Mt-at),Mt=xt-mt,at=xt-Mt,c[1]=xt-(Mt+at)+(at-mt),_t=yt+Mt,at=_t-yt,c[2]=yt-(_t-at)+(Mt-at),c[3]=_t,gt=vt*Ot,lt=n*vt,dt=lt-(lt-vt),bt=vt-dt,lt=n*Ot,ht=lt-(lt-Ot),pt=Ot-ht,jt=bt*pt-(gt-dt*ht-bt*ht-dt*pt),mt=wt*Ft,lt=n*wt,dt=lt-(lt-wt),bt=wt-dt,lt=n*Ft,ht=lt-(lt-Ft),pt=Ft-ht,Tt=bt*pt-(mt-dt*ht-bt*ht-dt*pt),Mt=jt-Tt,at=jt-Mt,s[0]=jt-(Mt+at)+(at-Tt),yt=gt+Mt,at=yt-gt,xt=gt-(yt-at)+(Mt-at),Mt=xt-mt,at=xt-Mt,s[1]=xt-(Mt+at)+(at-mt),_t=yt+Mt,at=_t-yt,s[2]=yt-(_t-at)+(Mt-at),s[3]=_t,N=e(e(e(f(f(4,r,vt,A),A,vt,F),F,f(f(4,r,Ft,A),A,Ft,O),O,k),k,e(f(f(4,c,wt,A),A,wt,F),F,f(f(4,c,Ot,A),A,Ot,O),O,q),q,B),B,e(f(f(4,s,At,A),A,At,F),F,f(f(4,s,Pt,A),A,Pt,O),O,k),k,C);let kt=function(t,n){let e=n[0];for(let o=1;o=qt||-kt>=qt)return kt;if(at=t-vt,Q=t-(vt+at)+(at-J),at=i-Ft,U=i-(Ft+at)+(at-K),at=D-wt,R=D-(wt+at)+(at-J),at=G-Ot,V=G-(Ot+at)+(at-K),at=H-At,S=H-(At+at)+(at-J),at=I-Pt,W=I-(Pt+at)+(at-K),0===Q&&0===R&&0===S&&0===U&&0===V&&0===W)return kt;if(qt=5423418723394464e-46*L+33306690738754706e-32*Math.abs(kt),kt+=(vt*vt+Ft*Ft)*(wt*W+Pt*R-(Ot*S+At*V))+2*(vt*Q+Ft*U)*(wt*Pt-Ot*At)+((wt*wt+Ot*Ot)*(At*U+Ft*S-(Pt*Q+vt*W))+2*(wt*R+Ot*V)*(At*Ft-Pt*vt))+((At*At+Pt*Pt)*(vt*V+Ot*Q-(Ft*R+wt*U))+2*(At*S+Pt*W)*(vt*Ot-Ft*wt)),kt>=qt||-kt>=qt)return kt;if(0===R&&0===V&&0===S&&0===W||(gt=vt*vt,lt=n*vt,dt=lt-(lt-vt),bt=vt-dt,jt=bt*bt-(gt-dt*dt-(dt+dt)*bt),mt=Ft*Ft,lt=n*Ft,dt=lt-(lt-Ft),bt=Ft-dt,Tt=bt*bt-(mt-dt*dt-(dt+dt)*bt),Mt=jt+Tt,at=Mt-jt,u[0]=jt-(Mt-at)+(Tt-at),yt=gt+Mt,at=yt-gt,xt=gt-(yt-at)+(Mt-at),Mt=xt+mt,at=Mt-xt,u[1]=xt-(Mt-at)+(mt-at),_t=yt+Mt,at=_t-yt,u[2]=yt-(_t-at)+(Mt-at),u[3]=_t),0===S&&0===W&&0===Q&&0===U||(gt=wt*wt,lt=n*wt,dt=lt-(lt-wt),bt=wt-dt,jt=bt*bt-(gt-dt*dt-(dt+dt)*bt),mt=Ot*Ot,lt=n*Ot,dt=lt-(lt-Ot),bt=Ot-dt,Tt=bt*bt-(mt-dt*dt-(dt+dt)*bt),Mt=jt+Tt,at=Mt-jt,a[0]=jt-(Mt-at)+(Tt-at),yt=gt+Mt,at=yt-gt,xt=gt-(yt-at)+(Mt-at),Mt=xt+mt,at=Mt-xt,a[1]=xt-(Mt-at)+(mt-at),_t=yt+Mt,at=_t-yt,a[2]=yt-(_t-at)+(Mt-at),a[3]=_t),0===Q&&0===U&&0===R&&0===V||(gt=At*At,lt=n*At,dt=lt-(lt-At),bt=At-dt,jt=bt*bt-(gt-dt*dt-(dt+dt)*bt),mt=Pt*Pt,lt=n*Pt,dt=lt-(lt-Pt),bt=Pt-dt,Tt=bt*bt-(mt-dt*dt-(dt+dt)*bt),Mt=jt+Tt,at=Mt-jt,l[0]=jt-(Mt-at)+(Tt-at),yt=gt+Mt,at=yt-gt,xt=gt-(yt-at)+(Mt-at),Mt=xt+mt,at=Mt-xt,l[1]=xt-(Mt-at)+(mt-at),_t=yt+Mt,at=_t-yt,l[2]=yt-(_t-at)+(Mt-at),l[3]=_t),0!==Q&&(X=f(4,r,Q,h),N=E(N,o(f(X,h,2*vt,F),F,f(f(4,l,Q,A),A,Ot,O),O,f(f(4,a,Q,A),A,-Pt,P),P,k,z),z)),0!==U&&(Y=f(4,r,U,p),N=E(N,o(f(Y,p,2*Ft,F),F,f(f(4,a,U,A),A,At,O),O,f(f(4,l,U,A),A,-wt,P),P,k,z),z)),0!==R&&(Z=f(4,c,R,M),N=E(N,o(f(Z,M,2*wt,F),F,f(f(4,u,R,A),A,Pt,O),O,f(f(4,l,R,A),A,-Ft,P),P,k,z),z)),0!==V&&($=f(4,c,V,y),N=E(N,o(f($,y,2*Ot,F),F,f(f(4,l,V,A),A,vt,O),O,f(f(4,u,V,A),A,-At,P),P,k,z),z)),0!==S&&(tt=f(4,s,S,x),N=E(N,o(f(tt,x,2*At,F),F,f(f(4,a,S,A),A,Ft,O),O,f(f(4,u,S,A),A,-Ot,P),P,k,z),z)),0!==W&&(nt=f(4,s,W,g),N=E(N,o(f(nt,g,2*Pt,F),F,f(f(4,u,W,A),A,wt,O),O,f(f(4,a,W,A),A,-vt,P),P,k,z),z)),0!==Q||0!==U){if(0!==R||0!==V||0!==S||0!==W?(gt=R*Pt,lt=n*R,dt=lt-(lt-R),bt=R-dt,lt=n*Pt,ht=lt-(lt-Pt),pt=Pt-ht,jt=bt*pt-(gt-dt*ht-bt*ht-dt*pt),mt=wt*W,lt=n*wt,dt=lt-(lt-wt),bt=wt-dt,lt=n*W,ht=lt-(lt-W),pt=W-ht,Tt=bt*pt-(mt-dt*ht-bt*ht-dt*pt),Mt=jt+Tt,at=Mt-jt,d[0]=jt-(Mt-at)+(Tt-at),yt=gt+Mt,at=yt-gt,xt=gt-(yt-at)+(Mt-at),Mt=xt+mt,at=Mt-xt,d[1]=xt-(Mt-at)+(mt-at),_t=yt+Mt,at=_t-yt,d[2]=yt-(_t-at)+(Mt-at),d[3]=_t,gt=S*-Ot,lt=n*S,dt=lt-(lt-S),bt=S-dt,lt=n*-Ot,ht=lt-(lt- -Ot),pt=-Ot-ht,jt=bt*pt-(gt-dt*ht-bt*ht-dt*pt),mt=At*-V,lt=n*At,dt=lt-(lt-At),bt=At-dt,lt=n*-V,ht=lt-(lt- -V),pt=-V-ht,Tt=bt*pt-(mt-dt*ht-bt*ht-dt*pt),Mt=jt+Tt,at=Mt-jt,b[0]=jt-(Mt-at)+(Tt-at),yt=gt+Mt,at=yt-gt,xt=gt-(yt-at)+(Mt-at),Mt=xt+mt,at=Mt-xt,b[1]=xt-(Mt-at)+(mt-at),_t=yt+Mt,at=_t-yt,b[2]=yt-(_t-at)+(Mt-at),b[3]=_t,ot=e(4,d,4,b,m),gt=R*W,lt=n*R,dt=lt-(lt-R),bt=R-dt,lt=n*W,ht=lt-(lt-W),pt=W-ht,jt=bt*pt-(gt-dt*ht-bt*ht-dt*pt),mt=S*V,lt=n*S,dt=lt-(lt-S),bt=S-dt,lt=n*V,ht=lt-(lt-V),pt=V-ht,Tt=bt*pt-(mt-dt*ht-bt*ht-dt*pt),Mt=jt-Tt,at=jt-Mt,v[0]=jt-(Mt+at)+(at-Tt),yt=gt+Mt,at=yt-gt,xt=gt-(yt-at)+(Mt-at),Mt=xt-mt,at=xt-Mt,v[1]=xt-(Mt+at)+(at-mt),_t=yt+Mt,at=_t-yt,v[2]=yt-(_t-at)+(Mt-at),v[3]=_t,rt=4):(m[0]=0,ot=1,v[0]=0,rt=1),0!==Q){const t=f(ot,m,Q,P);N=E(N,e(f(X,h,Q,F),F,f(t,P,2*vt,k),k,z),z);const n=f(rt,v,Q,A);N=E(N,o(f(n,A,2*vt,F),F,f(n,A,Q,O),O,f(t,P,Q,k),k,q,B),B),0!==V&&(N=E(N,f(f(4,l,Q,A),A,V,F),F)),0!==W&&(N=E(N,f(f(4,a,-Q,A),A,W,F),F))}if(0!==U){const t=f(ot,m,U,P);N=E(N,e(f(Y,p,U,F),F,f(t,P,2*Ft,k),k,z),z);const n=f(rt,v,U,A);N=E(N,o(f(n,A,2*Ft,F),F,f(n,A,U,O),O,f(t,P,U,k),k,q,B),B)}}if(0!==R||0!==V){if(0!==S||0!==W||0!==Q||0!==U?(gt=S*Ft,lt=n*S,dt=lt-(lt-S),bt=S-dt,lt=n*Ft,ht=lt-(lt-Ft),pt=Ft-ht,jt=bt*pt-(gt-dt*ht-bt*ht-dt*pt),mt=At*U,lt=n*At,dt=lt-(lt-At),bt=At-dt,lt=n*U,ht=lt-(lt-U),pt=U-ht,Tt=bt*pt-(mt-dt*ht-bt*ht-dt*pt),Mt=jt+Tt,at=Mt-jt,d[0]=jt-(Mt-at)+(Tt-at),yt=gt+Mt,at=yt-gt,xt=gt-(yt-at)+(Mt-at),Mt=xt+mt,at=Mt-xt,d[1]=xt-(Mt-at)+(mt-at),_t=yt+Mt,at=_t-yt,d[2]=yt-(_t-at)+(Mt-at),d[3]=_t,st=-Pt,ut=-W,gt=Q*st,lt=n*Q,dt=lt-(lt-Q),bt=Q-dt,lt=n*st,ht=lt-(lt-st),pt=st-ht,jt=bt*pt-(gt-dt*ht-bt*ht-dt*pt),mt=vt*ut,lt=n*vt,dt=lt-(lt-vt),bt=vt-dt,lt=n*ut,ht=lt-(lt-ut),pt=ut-ht,Tt=bt*pt-(mt-dt*ht-bt*ht-dt*pt),Mt=jt+Tt,at=Mt-jt,b[0]=jt-(Mt-at)+(Tt-at),yt=gt+Mt,at=yt-gt,xt=gt-(yt-at)+(Mt-at),Mt=xt+mt,at=Mt-xt,b[1]=xt-(Mt-at)+(mt-at),_t=yt+Mt,at=_t-yt,b[2]=yt-(_t-at)+(Mt-at),b[3]=_t,ft=e(4,d,4,b,T),gt=S*U,lt=n*S,dt=lt-(lt-S),bt=S-dt,lt=n*U,ht=lt-(lt-U),pt=U-ht,jt=bt*pt-(gt-dt*ht-bt*ht-dt*pt),mt=Q*W,lt=n*Q,dt=lt-(lt-Q),bt=Q-dt,lt=n*W,ht=lt-(lt-W),pt=W-ht,Tt=bt*pt-(mt-dt*ht-bt*ht-dt*pt),Mt=jt-Tt,at=jt-Mt,w[0]=jt-(Mt+at)+(at-Tt),yt=gt+Mt,at=yt-gt,xt=gt-(yt-at)+(Mt-at),Mt=xt-mt,at=xt-Mt,w[1]=xt-(Mt+at)+(at-mt),_t=yt+Mt,at=_t-yt,w[2]=yt-(_t-at)+(Mt-at),w[3]=_t,ct=4):(T[0]=0,ft=1,w[0]=0,ct=1),0!==R){const t=f(ft,T,R,P);N=E(N,e(f(Z,M,R,F),F,f(t,P,2*wt,k),k,z),z);const n=f(ct,w,R,A);N=E(N,o(f(n,A,2*wt,F),F,f(n,A,R,O),O,f(t,P,R,k),k,q,B),B),0!==W&&(N=E(N,f(f(4,u,R,A),A,W,F),F)),0!==U&&(N=E(N,f(f(4,l,-R,A),A,U,F),F))}if(0!==V){const t=f(ft,T,V,P);N=E(N,e(f($,y,V,F),F,f(t,P,2*Ot,k),k,z),z);const n=f(ct,w,V,A);N=E(N,o(f(n,A,2*Ot,F),F,f(n,A,V,O),O,f(t,P,V,k),k,q,B),B)}}if(0!==S||0!==W){if(0!==Q||0!==U||0!==R||0!==V?(gt=Q*Ot,lt=n*Q,dt=lt-(lt-Q),bt=Q-dt,lt=n*Ot,ht=lt-(lt-Ot),pt=Ot-ht,jt=bt*pt-(gt-dt*ht-bt*ht-dt*pt),mt=vt*V,lt=n*vt,dt=lt-(lt-vt),bt=vt-dt,lt=n*V,ht=lt-(lt-V),pt=V-ht,Tt=bt*pt-(mt-dt*ht-bt*ht-dt*pt),Mt=jt+Tt,at=Mt-jt,d[0]=jt-(Mt-at)+(Tt-at),yt=gt+Mt,at=yt-gt,xt=gt-(yt-at)+(Mt-at),Mt=xt+mt,at=Mt-xt,d[1]=xt-(Mt-at)+(mt-at),_t=yt+Mt,at=_t-yt,d[2]=yt-(_t-at)+(Mt-at),d[3]=_t,st=-Ft,ut=-U,gt=R*st,lt=n*R,dt=lt-(lt-R),bt=R-dt,lt=n*st,ht=lt-(lt-st),pt=st-ht,jt=bt*pt-(gt-dt*ht-bt*ht-dt*pt),mt=wt*ut,lt=n*wt,dt=lt-(lt-wt),bt=wt-dt,lt=n*ut,ht=lt-(lt-ut),pt=ut-ht,Tt=bt*pt-(mt-dt*ht-bt*ht-dt*pt),Mt=jt+Tt,at=Mt-jt,b[0]=jt-(Mt-at)+(Tt-at),yt=gt+Mt,at=yt-gt,xt=gt-(yt-at)+(Mt-at),Mt=xt+mt,at=Mt-xt,b[1]=xt-(Mt-at)+(mt-at),_t=yt+Mt,at=_t-yt,b[2]=yt-(_t-at)+(Mt-at),b[3]=_t,et=e(4,d,4,b,j),gt=Q*V,lt=n*Q,dt=lt-(lt-Q),bt=Q-dt,lt=n*V,ht=lt-(lt-V),pt=V-ht,jt=bt*pt-(gt-dt*ht-bt*ht-dt*pt),mt=R*U,lt=n*R,dt=lt-(lt-R),bt=R-dt,lt=n*U,ht=lt-(lt-U),pt=U-ht,Tt=bt*pt-(mt-dt*ht-bt*ht-dt*pt),Mt=jt-Tt,at=jt-Mt,_[0]=jt-(Mt+at)+(at-Tt),yt=gt+Mt,at=yt-gt,xt=gt-(yt-at)+(Mt-at),Mt=xt-mt,at=xt-Mt,_[1]=xt-(Mt+at)+(at-mt),_t=yt+Mt,at=_t-yt,_[2]=yt-(_t-at)+(Mt-at),_[3]=_t,it=4):(j[0]=0,et=1,_[0]=0,it=1),0!==S){const t=f(et,j,S,P);N=E(N,e(f(tt,x,S,F),F,f(t,P,2*At,k),k,z),z);const n=f(it,_,S,A);N=E(N,o(f(n,A,2*At,F),F,f(n,A,S,O),O,f(t,P,S,k),k,q,B),B),0!==U&&(N=E(N,f(f(4,a,S,A),A,U,F),F)),0!==V&&(N=E(N,f(f(4,u,-S,A),A,V,F),F))}if(0!==W){const t=f(et,j,W,P);N=E(N,e(f(nt,g,W,F),F,f(t,P,2*Pt,k),k,z),z);const n=f(it,_,W,A);N=E(N,o(f(n,A,2*Pt,F),F,f(n,A,W,O),O,f(t,P,W,k),k,q,B),B)}}return C[N-1]}(t,i,D,G,H,I,J,K,ft)},t.incirclefast=function(t,n,e,o,f,i,r,c){const s=t-r,u=n-c,a=e-r,l=o-c,d=f-r,b=i-c;return(s*s+u*u)*(a*b-d*l)+(a*a+l*l)*(d*u-s*b)+(d*d+b*b)*(s*l-a*u)},Object.defineProperty(t,"__esModule",{value:!0})}));