inline
std::ostream& operator<<(std::ostream& os,const measurement<Y>& val)
{
boost::io::ios_precision_saver precision_saver(os);
boost::io::ios_flags_saver flags_saver(os);
os << val.value() << "(+/-" << val.uncertainty() << ")";
return os;
}
inline
std::basic_ostream<Char,Traits>& operator<<(std::basic_ostream<Char,Traits>& os,const physical_constant<Y>& val)
{
boost::io::ios_precision_saver precision_saver(os);
//boost::io::ios_width_saver width_saver(os);
boost::io::ios_flags_saver flags_saver(os);
//os << std::setw(21);
typedef typename Y::value_type value_type;
if (val.uncertainty() > value_type())
{
const double relative_uncertainty = std::abs(val.uncertainty()/val.value());
const double exponent = std::log10(relative_uncertainty);
const long digits_of_precision = static_cast<long>(std::ceil(std::abs(exponent)))+3;
// should try to replicate NIST CODATA syntax
os << std::setprecision(digits_of_precision)
//<< std::setw(digits_of_precision+8)
//<< std::scientific
<< val.value();
// << long(10*(relative_uncertainty/std::pow(Y(10),Y(exponent))));
os << " (rel. unc. = "
<< std::setprecision(1)
//<< std::setw(7)
<< std::scientific
<< relative_uncertainty << ")";
}
else
{
os << val.value() << " (exact)";
}
return os;
}