link
T=t,t/=1.75,b=sqrt(t%4096)**.25*4e3,rb=((b^b*2)&255)*"1000010011000100"[t>>12&15],a=(t,p,T)=>t<0?0:((T*2**(([0,2,4,7,12,7,4,2][t*[1,1,2,[4,8][t>>13&1]][t>>14&3]>>12&7]-[0,2,4,5][t>>16&3])/12)*2&128)-64)*(1-t*[1,1,2,[4,8][t>>13&1]][t>>14&3]/4096%1),c=z=>(a(t,z,T)+a(t-(r=12288),z,T)/2+a(t-r*3,z,T)/3+a(t-r*4,z,T)/4+a(t-r*5,z,T)/5+a(t-r*6,z,T)/6)/1.5+128+(T*2**(-[0,2,4,5][t>>16&3]/12)&64)*(1*z(t*PI/32768)/4)-4,fd=(x,y)=>x%256*(y&2047)/2048,I=3,err=0/0,A=[err,err,1,1,I,I,1,I,1,I,err,err,I,I,1,I,err,err,err,err,I,I,1,I,1,I,err,err,I,I,1,I][(t*2>>12)%32],B=[err,err,err,err,I,I,1,I,1,I,err,err,I,I,1,I][(t*2>>12)%16],ABsnarehi=fd(t*tan(t>>[A,B][(t*2>>17)%2])&128*([[err,err,1,1,2,2,1,2,1,2,err,err,2,2,1,2,err,err,err,err,2,2,1,2,1,2,err,err,2,2,1,2][(t*2>>12)%32],[err,err,1,1,2,2,1,2,1,2,err,err,2,2,1,2,err,err,err,err,2,2,1,2,1,2,err,err,2,2,1,2][(t*2>>12)%32]][(t*2>>17)%2]>0),-t*2>>[2,2,2,2,2,2,1,2,1,2,2,2,2,2,1,2,2,2,2,2,2,2,1,2,1,2,2,2,2,2,1,2,2,2,2,2,2,2,1,2,1,2,2,2,2,2,1,2,2,2,2,2,2,2,1,2,1,2,2,2,2,2,1,2][(t*2>>12)%64])/256,[c(sin)+rb+ABsnarehi*128,c(cos)+rb+ABsnarehi*128].map(x=>x/128-1)