Add bounding box calculation and other function in calculations.py

In some case _bsphere and _besphere give a different result.