123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323324325326327328329330331332333334335336337338339340341342343344345346347348349350351352353354355356357358359360361362363364365366367368369370371372373374375376377378379380381382383384385386387388389390391392393394395396397398399400401402403404405406407408409410411412413414415416417418419420421422423424425426427428429430431432433 |
- /** nonfinite_num_facet_serialization.cpp
- *
- * Copyright (c) 2011 Francois Mauger
- * Copyright (c) 2011 Paul A. Bristow
- *
- * Distributed under the Boost Software License, Version 1.0.
- * (See accompanying file LICENSE_1_0.txt
- * or copy at http://www.boost.org/LICENSE_1_0.txt)
- *
- * This sample program by Francois Mauger illustrates how to use the
- * `boost/math/nonfinite_num_facets.hpp' material from the original
- * Floating Point Utilities contribution by Johan Rade. Here it is
- * shown how non finite floating number can be serialized and
- * deserialized from I/O streams and/or Boost text/XML archives. It
- * produces two archives stored in `test.txt' and `test.xml' files.
- *
- * Tested with Boost 1.44, gcc 4.4.1, Linux/i686 (32bits).
- * Tested with Boost.1.46.1 MSVC 10.0 32 bit.
- */
- #ifdef _MSC_VER
- # pragma warning(push)
- //# pragma warning(disable : 4100) // unreferenced formal parameter.
- #endif
- #include <iostream>
- #include <sstream>
- #include <fstream>
- #include <limits>
- #include <boost/cstdint.hpp>
- #include <boost/serialization/nvp.hpp>
- #include <boost/archive/text_oarchive.hpp>
- #include <boost/archive/text_iarchive.hpp>
- #include <boost/archive/xml_oarchive.hpp>
- #include <boost/archive/xml_iarchive.hpp>
- #include <boost/archive/codecvt_null.hpp>
- // from the Floating Point Utilities :
- #include <boost/math/special_functions/nonfinite_num_facets.hpp>
- static const char sep = ','; // Separator of bracketed float and double values.
- // Use max_digits10 (or equivalent) to obtain
- // all potentially significant decimal digits for the floating-point types.
- #ifdef BOOST_NO_CXX11_NUMERIC_LIMITS
- std::streamsize max_digits10_float = 2 + std::numeric_limits<float>::digits * 30103UL / 100000UL;
- std::streamsize max_digits10_double = 2 + std::numeric_limits<double>::digits * 30103UL / 100000UL;
- #else
- // Can use new C++0X max_digits10 (the maximum potentially significant digits).
- std::streamsize max_digits10_float = std::numeric_limits<float>::max_digits10;
- std::streamsize max_digits10_double = std::numeric_limits<double>::max_digits10;
- #endif
- /* A class with a float and a double */
- struct foo
- {
- foo () : fvalue (3.1415927F), dvalue (3.1415926535897931)
- { // Construct using 32 and 64-bit max_digits10 decimal digits value of pi.
- }
- // Set the values at -infinity :
- void minus_infinity ()
- {
- fvalue = -std::numeric_limits<float>::infinity ();
- dvalue = -std::numeric_limits<double>::infinity ();
- return;
- }
- // Set the values at +infinity :
- void plus_infinity ()
- {
- fvalue = +std::numeric_limits<float>::infinity ();
- dvalue = +std::numeric_limits<double>::infinity ();
- return;
- }
- // Set the values at NaN :
- void nan ()
- {
- fvalue = +std::numeric_limits<float>::quiet_NaN ();
- dvalue = +std::numeric_limits<double>::quiet_NaN ();
- return;
- }
- // Print :
- void print (std::ostream & a_out, const std::string & a_title)
- {
- if (a_title.empty ()) a_out << "foo";
- else a_out << a_title;
- a_out << " : " << std::endl;
- a_out << "|-- " << "fvalue = ";
- a_out.precision (7);
- a_out << fvalue << std::endl;
- a_out << "`-- " << "dvalue = ";
- a_out.precision (15);
- a_out << dvalue << std::endl;
- return;
- }
- // I/O operators :
- friend std::ostream & operator<< (std::ostream & a_out, const foo & a_foo);
- friend std::istream & operator>> (std::istream & a_in, foo & a_foo);
- // Boost serialization :
- template <class Archive>
- void serialize (Archive & ar, int /*version*/)
- {
- ar & BOOST_SERIALIZATION_NVP (fvalue);
- ar & BOOST_SERIALIZATION_NVP (dvalue);
- return;
- }
- // Attributes :
- float fvalue; // Single precision floating-point number.
- double dvalue; // Double precision floating-point number.
- };
- std::ostream & operator<< (std::ostream & a_out, const foo & a_foo)
- { // Output bracketed FPs, for example "(3.1415927,3.1415926535897931)"
- a_out.precision (max_digits10_float);
- a_out << "(" << a_foo.fvalue << sep ;
- a_out.precision (max_digits10_double);
- a_out << a_foo.dvalue << ")";
- return a_out;
- }
- std::istream & operator>> (std::istream & a_in, foo & a_foo)
- { // Input bracketed floating-point values into a foo structure,
- // for example from "(3.1415927,3.1415926535897931)"
- char c = 0;
- a_in.get (c);
- if (c != '(')
- {
- std::cerr << "ERROR: operator>> No ( " << std::endl;
- a_in.setstate(std::ios::failbit);
- return a_in;
- }
- float f;
- a_in >> std::ws >> f;
- if (! a_in)
- {
- return a_in;
- }
- a_in >> std::ws;
- a_in.get (c);
- if (c != sep)
- {
- std::cerr << "ERROR: operator>> c='" << c << "'" << std::endl;
- std::cerr << "ERROR: operator>> No '" << sep << "'" << std::endl;
- a_in.setstate(std::ios::failbit);
- return a_in;
- }
- double d;
- a_in >> std::ws >> d;
- if (! a_in)
- {
- return a_in;
- }
- a_in >> std::ws;
- a_in.get (c);
- if (c != ')')
- {
- std::cerr << "ERROR: operator>> No ) " << std::endl;
- a_in.setstate(std::ios::failbit);
- return a_in;
- }
- a_foo.fvalue = f;
- a_foo.dvalue = d;
- return a_in;
- }
- int main (void)
- {
- std::clog << std::endl
- << "Nonfinite_serialization.cpp' example program." << std::endl;
- #ifdef BOOST_NO_CXX11_NUMERIC_LIMITS
- std::cout << "BOOST_NO_CXX11_NUMERIC_LIMITS is defined, so no max_digits10 available either,"
- "using our own version instead." << std::endl;
- #endif
- std::cout << "std::numeric_limits<float>::max_digits10 is " << max_digits10_float << std::endl;
- std::cout << "std::numeric_limits<double>::max_digits10 is " << max_digits10_double << std::endl;
- std::locale the_default_locale (std::locale::classic (),
- new boost::archive::codecvt_null<char>);
- // Demonstrate use of nonfinite facets with stringstreams.
- {
- std::clog << "Construct some foo structures with a finite and nonfinites." << std::endl;
- foo f0;
- foo f1; f1.minus_infinity ();
- foo f2; f2.plus_infinity ();
- foo f3; f3.nan ();
- // Display them.
- f0.print (std::clog, "f0");
- f1.print (std::clog, "f1");
- f2.print (std::clog, "f2");
- f3.print (std::clog, "f3");
- std::clog << " Write to a string buffer." << std::endl;
- std::ostringstream oss;
- std::locale the_out_locale (the_default_locale, new boost::math::nonfinite_num_put<char>);
- oss.imbue (the_out_locale);
- oss.precision (max_digits10_double);
- oss << f0 << f1 << f2 << f3;
- std::clog << "Output is: `" << oss.str () << "'" << std::endl;
- std::clog << "Done output to ostringstream." << std::endl;
- }
- {
- std::clog << "Read foo structures from a string buffer." << std::endl;
- std::string the_string = "(3.1415927,3.1415926535897931)(-inf,-inf)(inf,inf)(nan,nan)";
- std::clog << "Input is: `" << the_string << "'" << std::endl;
- std::locale the_in_locale (the_default_locale, new boost::math::nonfinite_num_get<char>);
- std::istringstream iss (the_string);
- iss.imbue (the_in_locale);
- foo f0, f1, f2, f3;
- iss >> f0 >> f1 >> f2 >> f3;
- if (! iss)
- {
- std::cerr << "Format error !" << std::endl;
- }
- else
- {
- std::cerr << "Read OK." << std::endl;
- f0.print (std::clog, "f0");
- f1.print (std::clog, "f1");
- f2.print (std::clog, "f2");
- f3.print (std::clog, "f3");
- }
- std::clog << "Done input from istringstream." << std::endl;
- }
- { // Demonstrate use of nonfinite facets for Serialization with Boost text archives.
- std::clog << "Serialize (using Boost text archive)." << std::endl;
- // Construct some foo structures with a finite and nonfinites.
- foo f0;
- foo f1; f1.minus_infinity ();
- foo f2; f2.plus_infinity ();
- foo f3; f3.nan ();
- // Display them.
- f0.print (std::clog, "f0");
- f1.print (std::clog, "f1");
- f2.print (std::clog, "f2");
- f3.print (std::clog, "f3");
- std::locale the_out_locale (the_default_locale, new boost::math::nonfinite_num_put<char>);
- // Use a temporary folder .temps (which contains "boost-no-inspect" so that it will not be inspected, and made 'hidden' too).
- std::ofstream fout ("./.temps/nonfinite_archive_test.txt");
- fout.imbue (the_out_locale);
- boost::archive::text_oarchive toar (fout, boost::archive::no_codecvt);
- // Write to archive.
- toar & f0;
- toar & f1;
- toar & f2;
- toar & f3;
- std::clog << "Done." << std::endl;
- }
- {
- std::clog << "Deserialize (Boost text archive)..." << std::endl;
- std::locale the_in_locale (the_default_locale, new boost::math::nonfinite_num_get<char>);
- // Use a temporary folder .temps (which contains "boost-no-inspect" so that it will not be inspected, and made 'hidden' too).
- std::ifstream fin ("./.temps/nonfinite_archive_test.txt");
- fin.imbue (the_in_locale);
- boost::archive::text_iarchive tiar (fin, boost::archive::no_codecvt);
- foo f0, f1, f2, f3;
- // Read from archive.
- tiar & f0;
- tiar & f1;
- tiar & f2;
- tiar & f3;
- // Display foos.
- f0.print (std::clog, "f0");
- f1.print (std::clog, "f1");
- f2.print (std::clog, "f2");
- f3.print (std::clog, "f3");
- std::clog << "Done." << std::endl;
- }
- { // Demonstrate use of nonfinite facets for Serialization with Boost XML Archive.
- std::clog << "Serialize (Boost XML archive)..." << std::endl;
- // Construct some foo structures with a finite and nonfinites.
- foo f0;
- foo f1; f1.minus_infinity ();
- foo f2; f2.plus_infinity ();
- foo f3; f3.nan ();
- // Display foos.
- f0.print (std::clog, "f0");
- f1.print (std::clog, "f1");
- f2.print (std::clog, "f2");
- f3.print (std::clog, "f3");
- std::locale the_out_locale (the_default_locale, new boost::math::nonfinite_num_put<char>);
- // Use a temporary folder /.temps (which contains "boost-no-inspect" so that it will not be inspected, and made 'hidden' too).
- std::ofstream fout ("./.temps/nonfinite_XML_archive_test.txt");
- fout.imbue (the_out_locale);
- boost::archive::xml_oarchive xoar (fout, boost::archive::no_codecvt);
- xoar & BOOST_SERIALIZATION_NVP (f0);
- xoar & BOOST_SERIALIZATION_NVP (f1);
- xoar & BOOST_SERIALIZATION_NVP (f2);
- xoar & BOOST_SERIALIZATION_NVP (f3);
- std::clog << "Done." << std::endl;
- }
- {
- std::clog << "Deserialize (Boost XML archive)..." << std::endl;
- std::locale the_in_locale (the_default_locale, new boost::math::nonfinite_num_get<char>);
- // Use a temporary folder /.temps (which contains "boost-no-inspect" so that it will not be inspected, and made 'hidden' too).
- std::ifstream fin ("./.temps/nonfinite_XML_archive_test.txt"); // Previously written above.
- fin.imbue (the_in_locale);
- boost::archive::xml_iarchive xiar (fin, boost::archive::no_codecvt);
- foo f0, f1, f2, f3;
- xiar & BOOST_SERIALIZATION_NVP (f0);
- xiar & BOOST_SERIALIZATION_NVP (f1);
- xiar & BOOST_SERIALIZATION_NVP (f2);
- xiar & BOOST_SERIALIZATION_NVP (f3);
- f0.print (std::clog, "f0");
- f1.print (std::clog, "f1");
- f2.print (std::clog, "f2");
- f3.print (std::clog, "f3");
- std::clog << "Done." << std::endl;
- }
- std::clog << "End nonfinite_serialization.cpp' example program." << std::endl;
- return 0;
- }
- /*
- Output:
- Nonfinite_serialization.cpp' example program.
- std::numeric_limits<float>::max_digits10 is 8
- std::numeric_limits<double>::max_digits10 is 17
- Construct some foo structures with a finite and nonfinites.
- f0 :
- |-- fvalue = 3.141593
- `-- dvalue = 3.14159265358979
- f1 :
- |-- fvalue = -1.#INF
- `-- dvalue = -1.#INF
- f2 :
- |-- fvalue = 1.#INF
- `-- dvalue = 1.#INF
- f3 :
- |-- fvalue = 1.#QNAN
- `-- dvalue = 1.#QNAN
- Write to a string buffer.
- Output is: `(3.1415927,3.1415926535897931)(-inf,-inf)(inf,inf)(nan,nan)'
- Done output to ostringstream.
- Read foo structures from a string buffer.
- Input is: `(3.1415927,3.1415926535897931)(-inf,-inf)(inf,inf)(nan,nan)'
- Read OK.
- f0 :
- |-- fvalue = 3.141593
- `-- dvalue = 3.14159265358979
- f1 :
- |-- fvalue = -1.#INF
- `-- dvalue = -1.#INF
- f2 :
- |-- fvalue = 1.#INF
- `-- dvalue = 1.#INF
- f3 :
- |-- fvalue = 1.#QNAN
- `-- dvalue = 1.#QNAN
- Done input from istringstream.
- Serialize (using Boost text archive).
- f0 :
- |-- fvalue = 3.141593
- `-- dvalue = 3.14159265358979
- f1 :
- |-- fvalue = -1.#INF
- `-- dvalue = -1.#INF
- f2 :
- |-- fvalue = 1.#INF
- `-- dvalue = 1.#INF
- f3 :
- |-- fvalue = 1.#QNAN
- `-- dvalue = 1.#QNAN
- Done.
- Deserialize (Boost text archive)...
- f0 :
- |-- fvalue = 3.141593
- `-- dvalue = 3.14159265358979
- f1 :
- |-- fvalue = -1.#INF
- `-- dvalue = -1.#INF
- f2 :
- |-- fvalue = 1.#INF
- `-- dvalue = 1.#INF
- f3 :
- |-- fvalue = 1.#QNAN
- `-- dvalue = 1.#QNAN
- Done.
- Serialize (Boost XML archive)...
- f0 :
- |-- fvalue = 3.141593
- `-- dvalue = 3.14159265358979
- f1 :
- |-- fvalue = -1.#INF
- `-- dvalue = -1.#INF
- f2 :
- |-- fvalue = 1.#INF
- `-- dvalue = 1.#INF
- f3 :
- |-- fvalue = 1.#QNAN
- `-- dvalue = 1.#QNAN
- Done.
- Deserialize (Boost XML archive)...
- f0 :
- |-- fvalue = 3.141593
- `-- dvalue = 3.14159265358979
- f1 :
- |-- fvalue = -1.#INF
- `-- dvalue = -1.#INF
- f2 :
- |-- fvalue = 1.#INF
- `-- dvalue = 1.#INF
- f3 :
- |-- fvalue = 1.#QNAN
- `-- dvalue = 1.#QNAN
- Done.
- End nonfinite_serialization.cpp' example program.
- */
|