Boost.Mathにある円周率を求める関数boost::math::constants::pi()ですが、単に定数を型Tにキャストして返すだけの実装だと思っていたら、実は違ったようです。定数で表現できる以上の桁数を要求された場合は、「ldexp(acos(T(0)), 1)」で計算するようになっ…