lkpyramid.cpp,
const, a/b==a*(1/b), sqrt(a)*sqrt(b) == sqrt(a*b) // both faster