inline Quad Sin( const Quad& alpha ) { return sinq(alpha); }
inline Complex<Quad> ComplexFromPolar( const Quad& r, const Quad& theta ) { const Quad realPart = r*cosq(theta); const Quad imagPart = r*sinq(theta); return Complex<Quad>(realPart,imagPart); }
__float128 h2(__float128 x, __float128 y) {return sinq(base(x,y));} //imaginary part
inline void eval_sin(float128_backend& result, const float128_backend& arg) { result.value() = sinq(arg.value()); }
//rapidly spinning harmonic! Should pick up on any and all distortions. __float128 base(__float128 x, __float128 y) {return 500.Q*(sinq(x)*sinq(x)+sinq(y)*sinq(y));} //ellipse parametrization