/*
　荒い √(x*x+y*y) の計算
　　最も精度が悪い場所で1% 程度の誤差
　>>3 のような場所で荒い判定に使う

*/
double roughHypot1(double x, double y)
{
x=fabs(x);
y=fabs(y);
if (x<y){ double w=x;x=y;y=w;};
return x+ (y*y/x)* 0.429;
}