lkpyramid.cpp:
sqrt(a)*sqrt(b) == sqrt(a*b),  fabs(c)<sqrt(d) == c*c<d