Substitution of pi for its numerical value moved from
useful-functions.hh to fix_c.r