TinyChan

var{abs,floor,cos,sin,cosh,sinh,tan,sqrt,atan2,exp,log,pow,PI}=Math;
TAU=PI*2;mod=(n,m)=>(n%m+m)%m;fract=x=>(x%1+1)%1;mix=(a,b,t)=>a+(b-a)*t;
cAdd=(a,b)=>[a[0]+b[0],a[1]+b[1]];cSub=(a,b)=>[a[0]-b[0],a[1]-b[1]];cMul=(a,b)=>[a[0]*b[0]-a[1]*b[1],a[0]*b[1]+a[1]*b[0]];cDiv=(a,b)=>[(a[0]*b[0]+a[1]*b[1])/(b[0]**2+b[1]**2),(a[1]*b[0]-a[0]*b[1])/(b[0]**2+b[1]**2)];
polar=(r,θ)=>[r*cos(θ),r*sin(θ)];
cAbs=c=>sqrt(c[0]**2+c[1]**2);cArg=c=>atan2(c[1],c[0]);
cExp=c=>[exp(c[0])*cos(c[1]),exp(c[0])*sin(c[1])];
cPow=(a,b)=>[pow(a[0]**2+a[1]**2,b[0]/2)*exp(-b[1]*cArg(a))*cos(.5*b[1]*log(a[0]**2+a[1]**2)+b[0]*cArg(a)),pow(a[0]**2+a[1]**2,b[0]/2)*exp(-b[1]*cArg(a))*sin(.5*b[1]*log(a[0]**2+a[1]**2)+b[0]*cArg(a))];
cLog=c=>[log(cAbs(c)),cArg(c)];
//(a+bi)/c = (a/c)+(b/c)i
cSgn=c=>{const r=cAbs(c);if(r===0)return[0,0];return[c[0]/r,c[1]/r];};
cEqual=(a,b,ε=1e-14)=>abs(a[0]-b[0])<ε&&abs(a[1]-b[1])<ε;
cSqrt=c=>polar(sqrt(cAbs(c)),cArg(c)/2);
cSin=c=>[sin(c[0])*cosh(c[1]),cos(c[0])*sinh(c[1])];cCos=c=>[cos(c[0])*cosh(c[1]),-sin(c[0])*sinh(c[1])];cTan=c=>cDiv(cSin(c),cCos(c));
cCot=c=>cDiv(cCos(c),cSin(c));cSec=c=>cDiv([1,0],cCos(c));cCsc=c=>cDiv([1,0],cSin(c));
cSinh=c=>cDiv(cSub(cExp(c),cExp([-c[0],-c[1]])),[2,0]);cCosh=c=>cDiv(cAdd(cExp(c),cExp([-c[0],-c[1]])),[2,0]);