common.h

Go to the documentation of this file.

 
#ifndef COMMON_H 
#define COMMON_H
 
#include <iomanip>
#include <iostream>
#include <string>
#include <cmath>
#include <cassert>
#include <vector>
#include <set>
#include <map>
#include <algorithm>
#include <sstream>
#include <chrono>
#include <functional> // has std::bind
 
#define BZ_HAVE_BOOST_SERIALIZATION
#include <blitz/array.h>
#include <boost/archive/text_oarchive.hpp>
#include <boost/archive/text_iarchive.hpp>
#include <boost/archive/binary_oarchive.hpp>
#include <boost/archive/binary_iarchive.hpp>
 
#include <boost/format.hpp>
#include <boost/ptr_container/ptr_vector.hpp>
#include <boost/ptr_container/ptr_map.hpp>
 
#include <boost/property_tree/ptree.hpp>
#include <boost/property_tree/xml_parser.hpp>
 
#include <boost/math/tools/minima.hpp> // find minima using Brent's method
 
/* Debugging librarys and definitions. There is only an effect if
 * the code is compiled with PIMC_DEBUG on.*/
#ifdef PIMC_DEBUG
 
#define BZ_DEBUG
 
#define PIMC_DEBUG_MESSAGE(X)\
{ std::cerr << "[pimc] " << __FILE__ << ": " << __LINE__ << " " << X << std::endl; }
 
#define PIMC_ASSERT(X) assert(X)
 
#else // !PIMC_DEBUG
 
#define PIMC_DEBUG_MESSAGE(X)
 
#define PIMC_ASSERT(X) 
 
#endif // ifdef PIMC_DEBUG 
 
/* We either specify the number of dimensions of the simulation at compile time,
 * or it defaults to 1D. */
#ifndef NDIM
#define NDIM 1 
#endif
 
/* Determine if we are performing a PIGS (T=0) or PIMC (T>0) simulation */
#ifndef PIGS
#define PIGS false
#endif
 
/* Used for getting the repo version number into the code */
#ifndef REPO_VERSION
#define REPO_VERSION "none"
#endif
 
/* We Default to turning on the NN lookup table. Comment this line out for testing
 * small systems where the lookup table might actually slow things down. */
#define NN_TABLE 
 
#include "MersenneTwister.h"
 
#define NPCFSEP 50      
#define NOBDMSEP 50     
#define NRADSEP 200     
#define NGRIDSEP 51     
#define EPS 1.0E-7      
#define DBL_EPS std::numeric_limits<double>::epsilon() //< Smallest double
#define BIG 1.0E30      
#define LBIG 69.07755279 
#define XXX -1          
 
using namespace std;
using namespace blitz;
using boost::format;
 
typedef unsigned long uint32;
 
typedef TinyMatrix<double,NDIM,NDIM> dMat;
 
typedef TinyVector<double,NDIM> dVec;
 
typedef TinyVector<int,NDIM> iVec;
 
typedef TinyVector<int,2> beadLocator;
 
typedef Array<int,1>::iterator intIter;
 
typedef Array<int,1>::const_iterator cintIter;
 
typedef Array<beadLocator,1>::iterator beadIter;
 
enum beadState {HEADTAIL,SPECIAL,NONE};
 
enum ensemble {DIAGONAL, OFFDIAGONAL, ANY};
 
inline int ipow (int base, int power) {
    return static_cast<int>(floor(pow(1.0*base,1.0*power) + EPS));
}
 
template<typename Ttype>
inline Ttype& min(Ttype& x, Ttype& y) { return (x < y ? x : y); } 
 
template<typename Ttype>
inline Ttype& max(Ttype& x, Ttype& y) { return (x > y ? x : y); } 
 
#endif