3
0
mirror of https://github.com/triqs/dft_tools synced 2024-12-25 05:43:40 +01:00

parameters : small details

- an include missing on OS X
- change serialization in the general case where hdf5 can be used,
  (version of hdf5 high enough) to use only one buffer,
  not a buffer per field.
- put back the mpi_param example.
This commit is contained in:
Olivier Parcollet 2014-05-21 17:47:54 +02:00
parent 27af80c04e
commit 4789dccda6
5 changed files with 43 additions and 25 deletions

View File

@ -26,12 +26,9 @@ int main(int argc, char* argv[]) {
P["d"] = 2.7; P["d"] = 2.7;
P["s"] = std::string("-14.3"); P["s"] = std::string("-14.3");
} }
/*
* WITH SERIALIZATION
std::cout << "Before bcast rank " << world.rank() << " : " << P << std::endl ; std::cout << "Before bcast rank " << world.rank() << " : " << P << std::endl ;
boost::mpi::broadcast(world,P,0); boost::mpi::broadcast(world,P,0);
std::cout << "After bcast rank " << world.rank() << " : " << P << std::endl ; std::cout << "After bcast rank " << world.rank() << " : " << P << std::endl ;
*/
return 0; return 0;
} }

View File

@ -48,7 +48,7 @@ namespace params {
#define TRIQS_UTIL_OPAQUE_OBJECT_PREDEFINED_CAST \ #define TRIQS_UTIL_OPAQUE_OBJECT_PREDEFINED_CAST \
(int)(long)(long long)(unsigned int)(unsigned long)(unsigned long long)(double)(bool)(std::string) (int)(long)(long long)(unsigned int)(unsigned long)(unsigned long long)(double)(bool)(std::string)
// the type actually stored // the type actually stored
template <typename T, typename Enable = void> struct storage_t_impl : std::conditional<std::is_integral<T>::value, long, T>{}; template <typename T, typename Enable = void> struct storage_t_impl : std::conditional<std::is_integral<T>::value, long, T>{};
template <typename T> using storage_t = typename storage_t_impl<std::c14::decay_t<T>>::type; template <typename T> using storage_t = typename storage_t_impl<std::c14::decay_t<T>>::type;
@ -109,14 +109,14 @@ namespace params {
std::string name_; // for_error_messages std::string name_; // for_error_messages
bool modified = false; bool modified = false;
bool no_default_value = false; bool no_default_value = false;
// only parameters will construct _field // only parameters will construct _field
template <typename T> template <typename T>
_field(T obj, std::string n, bool without_default_value) _field(T obj, std::string n, bool without_default_value)
: index(typeid(storage_t<T>)), p(new _data_impl<storage_t<T>>{std::move(obj)}), name_(n), no_default_value(without_default_value) { : index(typeid(storage_t<T>)), p(new _data_impl<storage_t<T>>{std::move(obj)}), name_(n), no_default_value(without_default_value) {
type_names.insert({index, get_triqs_hdf5_data_scheme(storage_t<T>{})}); type_names.insert({index, get_triqs_hdf5_data_scheme(storage_t<T>{})});
} }
std::string type_name() const { return type_name(index); } std::string type_name() const { return type_name(index); }
static std::string type_name(std::type_index i); static std::string type_name(std::type_index i);
@ -125,8 +125,8 @@ namespace params {
template <typename T> bool has_type() const { return index == typeid(T); } template <typename T> bool has_type() const { return index == typeid(T); }
public: public:
#ifdef TRIQS_SERIALIZATION_WITH_BOOST #ifdef TRIQS_SERIALIZATION_WITH_HDF5_IMPOSSIBLE
_field() : index(typeid(void)) {} // BREAKS invariant : only used for BOOST serialization. _field() : index(typeid(void)) {} // BREAKS invariant : only used for BOOST serialization failsafe mode.
#endif #endif
_field(_field &&c) = default; _field(_field &&c) = default;
_field(_field const &x) : index(x.index), p(x.p ? x.p->clone() : nullptr), name_(x.name_), modified(x.modified), no_default_value(x.no_default_value) {} // regular type _field(_field const &x) : index(x.index), p(x.p ? x.p->clone() : nullptr), name_(x.name_), modified(x.modified), no_default_value(x.no_default_value) {} // regular type
@ -137,7 +137,7 @@ namespace params {
template <typename RHS> _field &operator=(RHS rhs) { // pass by value (sink) template <typename RHS> _field &operator=(RHS rhs) { // pass by value (sink)
modified = true; modified = true;
using S_t = storage_t<RHS>; using S_t = storage_t<RHS>;
if (index != typeid(S_t)) { if (index != typeid(S_t)) {
TRIQS_RUNTIME_ERROR << "Field "<< name_<<" is of type "<< type_name(index)<<". I can not put a " << type_name(typeid(S_t)) << " into it."; TRIQS_RUNTIME_ERROR << "Field "<< name_<<" is of type "<< type_name(index)<<". I can not put a " << type_name(typeid(S_t)) << " into it.";
} }
p.reset(new _data_impl<S_t>{std::move(rhs)}); p.reset(new _data_impl<S_t>{std::move(rhs)});
@ -145,7 +145,7 @@ namespace params {
} }
// rewrite a few cases for practical convenience ... // rewrite a few cases for practical convenience ...
_field &operator=(int rhs) { // a special case where we can correct : int -> double _field &operator=(int rhs) { // a special case where we can correct : int -> double
if (index == typeid(double)) return operator=(double(rhs)); if (index == typeid(double)) return operator=(double(rhs));
return operator=(long(rhs));// beware infinite loop! Works because int != long ... Clean this ... return operator=(long(rhs));// beware infinite loop! Works because int != long ... Clean this ...
} }
@ -169,7 +169,7 @@ namespace params {
/// Is a modification required for this field /// Is a modification required for this field
bool is_modification_required() const { return no_default_value && (!modified); } bool is_modification_required() const { return no_default_value && (!modified); }
#ifdef TRIQS_WITH_PYTHON_SUPPORT #ifdef TRIQS_WITH_PYTHON_SUPPORT
/// Convertions python <-> C++ /// Convertions python <-> C++
bool from_python_convertible(PyObject *ob) const { return p->from_python_convertible(ob); } bool from_python_convertible(PyObject *ob) const { return p->from_python_convertible(ob); }
@ -182,6 +182,7 @@ namespace params {
BOOST_PP_SEQ_FOR_EACH(CAST_OPERATOR, nil, TRIQS_UTIL_OPAQUE_OBJECT_PREDEFINED_CAST); BOOST_PP_SEQ_FOR_EACH(CAST_OPERATOR, nil, TRIQS_UTIL_OPAQUE_OBJECT_PREDEFINED_CAST);
#undef CAST_OPERATOR #undef CAST_OPERATOR
#ifdef TRIQS_SERIALIZATION_WITH_HDF5_IMPOSSIBLE
// ----- Boost serialisation // ----- Boost serialisation
template <class Archive> void save(Archive &ar, const unsigned int version) const { template <class Archive> void save(Archive &ar, const unsigned int version) const {
@ -194,6 +195,7 @@ namespace params {
p->deserialize(s); p->deserialize(s);
} }
BOOST_SERIALIZATION_SPLIT_MEMBER(); BOOST_SERIALIZATION_SPLIT_MEMBER();
#endif
friend std::ostream &operator<<(std::ostream &out, _field const &ob) { return ob.p->print(out); } friend std::ostream &operator<<(std::ostream &out, _field const &ob) { return ob.p->print(out); }

View File

@ -4,6 +4,7 @@
#include <triqs/utility/serialization.hpp> #include <triqs/utility/serialization.hpp>
#include <algorithm> #include <algorithm>
#include <boost/algorithm/string.hpp> #include <boost/algorithm/string.hpp>
#include <iomanip>
namespace triqs { namespace triqs {
namespace params { namespace params {

View File

@ -31,14 +31,14 @@ namespace params {
* DOC to be written. * DOC to be written.
* *
* Provide a form for program parameters. * Provide a form for program parameters.
* *
* Fields of the form can be added with add_field method (only in C++), with a name, a type, a doc, and optionally a default value. * Fields of the form can be added with add_field method (only in C++), with a name, a type, a doc, and optionally a default value.
* *
* They can be retrieved from and put into the form, using the [] operator, in C++ and python. * They can be retrieved from and put into the form, using the [] operator, in C++ and python.
* *
* The data is stored with a strict type checking at runtime, with some type collapse (all integers into a long, * The data is stored with a strict type checking at runtime, with some type collapse (all integers into a long,
* all string, const char * into a std::string, etc). * all string, const char * into a std::string, etc).
* *
* Basic type cast to basic types are provided. * Basic type cast to basic types are provided.
* *
* The class is boost-serializable and implements hdf5 I/O operations. * The class is boost-serializable and implements hdf5 I/O operations.
@ -48,9 +48,11 @@ namespace params {
std::string key; std::string key;
_field f; _field f;
std::string doc; std::string doc;
#ifdef TRIQS_SERIALIZATION_WITH_HDF5_IMPOSSIBLE
template <class Archive> void serialize(Archive& ar, const unsigned int version) { template <class Archive> void serialize(Archive& ar, const unsigned int version) {
ar& TRIQS_MAKE_NVP("key", key) & TRIQS_MAKE_NVP("f", f) & TRIQS_MAKE_NVP("doc", doc); ar& TRIQS_MAKE_NVP("key", key) & TRIQS_MAKE_NVP("f", f) & TRIQS_MAKE_NVP("doc", doc);
} }
#endif
}; };
using _data_t = std::vector<_data_elem>; using _data_t = std::vector<_data_elem>;
_data_t _data; _data_t _data;
@ -59,7 +61,23 @@ namespace params {
_data_t::const_iterator find(std::string const& key) const; _data_t::const_iterator find(std::string const& key) const;
friend class boost::serialization::access; friend class boost::serialization::access;
template <class Archive> void serialize(Archive& ar, const unsigned int version) { ar& TRIQS_MAKE_NVP("_data", _data); } #ifdef TRIQS_SERIALIZATION_WITH_HDF5_IMPOSSIBLE
template <class Archive> void serialize(Archive& ar, const unsigned int version) {
ar& TRIQS_MAKE_NVP("_data", _data);
}
#else
// do it at once, with one buffer only.
template <class Archive> void save(Archive &ar, const unsigned int version) const {
std::string s = h5::serialize(*this);
ar << TRIQS_MAKE_NVP("seria_str", s);
}
template <class Archive> void load(Archive &ar, const unsigned int version) {
std::string s;
ar >> TRIQS_MAKE_NVP("seria_str", s);
*this = h5::deserialize<parameters>(s);
}
BOOST_SERIALIZATION_SPLIT_MEMBER();
#endif
public: public:
parameters(); parameters();
@ -88,13 +106,13 @@ namespace params {
/// Does the form have the key ? /// Does the form have the key ?
bool has_key(std::string const& key) const; bool has_key(std::string const& key) const;
/** /**
* Access the parameter key. * Access the parameter key.
* Key must be a valid key or a TRIQS_RUNTIME_ERROR is thrown. * Key must be a valid key or a TRIQS_RUNTIME_ERROR is thrown.
*/ */
_field& operator[](const char * key); _field& operator[](const char * key);
/** /**
* Access the parameter key. * Access the parameter key.
* Key must be a valid key or a TRIQS_RUNTIME_ERROR is thrown. * Key must be a valid key or a TRIQS_RUNTIME_ERROR is thrown.
*/ */
@ -103,14 +121,14 @@ namespace params {
/// generate help in form of a table of strings containing a list of required and optional parameters /// generate help in form of a table of strings containing a list of required and optional parameters
std::vector<std::vector<std::string>> generate_help(bool with_header = true) const; std::vector<std::vector<std::string>> generate_help(bool with_header = true) const;
/// help as a string /// help as a string
std::string help() const; std::string help() const;
/// hdf5 /// hdf5
friend std::string get_triqs_hdf5_data_scheme(parameters const &) { return ""; } friend std::string get_triqs_hdf5_data_scheme(parameters const &) { return ""; }
friend void h5_write(h5::group F, std::string const& subgroup_name, parameters const& p); friend void h5_write(h5::group F, std::string const& subgroup_name, parameters const& p);
friend void h5_read(h5::group F, std::string const& subgroup_name, parameters& p); friend void h5_read(h5::group F, std::string const& subgroup_name, parameters& p);
/// Ostream /// Ostream
friend std::ostream& operator<<(std::ostream& out, parameters const& p); friend std::ostream& operator<<(std::ostream& out, parameters const& p);

View File

@ -20,14 +20,14 @@
******************************************************************************/ ******************************************************************************/
#pragma once #pragma once
#if defined(TRIQS_SERIALIZATION_USE_BOOST) or not H5_VERSION_GE(1,8,9) #if defined(TRIQS_SERIALIZATION_DO_NOT_USE_HDF5) or not H5_VERSION_GE(1,8,9)
#define TRIQS_SERIALIZATION_WITH_HDF5_IMPOSSIBLE
#include "./boost_serialization.hpp" #include "./boost_serialization.hpp"
#define TRIQS_SERIALIZATION_WITH_BOOST
#else #else
#include "../h5/serialization.hpp"
#define TRIQS_SERIALIZATION_WITH_HDF5
#include "../h5/serialization.hpp"
namespace triqs { namespace triqs {
using h5::serialize; using h5::serialize;
using h5::deserialize; using h5::deserialize;