cstream.h

Go to the documentation of this file.

#ifndef _CSTREAM_H_

#include "tree.h"
#include "assert.h"
#include "stdlib.h"

#ifdef _OPENMP
# include "omp.h"
#endif

#ifndef TRUE
# define TRUE 1
#endif 

#ifndef FALSE
# define FALSE 0
#endif 


#ifdef DEBUG_LEVEL
# ifdef _OPENMP
#  define debug(level, ...) if (DEBUG_LEVEL >= level) {        \
    fprintf (stderr,"[%d] %s:%d: ", omp_get_thread_num(),    \
        __FILE__, __LINE__);              \
    fprintf (stderr,__VA_ARGS__);               \
  }
# else /*_OPENMP*/
#  define debug(level, ...) if (DEBUG_LEVEL >= level) {        \
    fprintf (stderr,"%s:%d: ", __FILE__, __LINE__);         \
    fprintf (stderr,__VA_ARGS__);               \
  }
# endif
#else
# define debug(level, ...)
# define NDEBUG  /* To supress the asserts */
#endif

#define warning(...) do{            \
    fprintf (stderr, "%s: Warning: ", invok_name); \
    fprintf (stderr, ## __VA_ARGS__);        \
  } while(0)

#define fatal(...) do{              \
    fprintf (stderr, "%s: Fatal error: ", invok_name);   \
    fprintf (stderr, ## __VA_ARGS__);        \
    exit(-1);                 \
  } while(0)
    
/* Useful to debug. */
#define show_double(VAR_) printf (#VAR_ " = %g\n", VAR_)
#define show_int(VAR_) printf (#VAR_ " = %d\n", VAR_)

/* Beware of side-effects! */
#define MYMAX(X_, Y_)  ((X_) > (Y_)? (X_): (Y_))
#define MYMIN(X_, Y_)  ((X_) < (Y_)? (X_): (Y_))
#define MAX_AT_LEVEL(X_, Y_, L_)  MYMAX(X_, (Y_) << (L_))
#define MIN_AT_LEVEL(X_, Y_, L_)  MYMIN(X_, (Y_) << (L_))
#define MAX_AT_LEVEL_WITH_SHIFT(X_, Y_, L_, S_)  \
  MYMAX(X_, ((Y_) << (L_)) + ((L_) > 0? ((S_) << (L_ - 1)): 0))
#define MIN_AT_LEVEL_WITH_SHIFT(X_, Y_, L_, S_)  \
  MYMIN(X_, ((Y_) << (L_)) + ((L_) > 0? ((S_) << (L_ - 1)): 0))

#define XCHG(X1_, X2_) do {         \
    typeof(X1_) TMP_;            \
    TMP_ = X1_;               \
    X1_ = X2_;             \
    X2_ = TMP_;               \
  } while(0)


#define SQ(X_)  ((X_) * (X_))

//typedef struct param_t param_t;
//struct param_t {
  //char *name;
  //char *desc;
  //char *type;
  //void *value;
//};

typedef char* string;
typedef double* doublep;

//#ifdef ALLOC_PARAMS
//# define decl_param(TYPE, NAME, DESC, DEFAULT)      \
       //TYPE NAME = DEFAULT;          \
       //param_t NAME ## _st = {          \
    //#NAME,                  \
    //DESC,                \
    //#TYPE,                  \
    //(void *) &NAME};
//# define decl_deprec_param(TYPE, NAME, DESC, DEFAULT)  \
       //TYPE NAME;              \
       //param_t NAME ## _st = {          \
    //#NAME,                  \
    //DESC,                \
    //"deprecated",              \
    //NULL};
//#else
//# define decl_param(TYPE, NAME, DESC, DEFAULT)   \
  //extern TYPE NAME;            \
  //extern param_t NAME ## _st;
//# define decl_deprec_param(TYPE, NAME, DESC, DEFAULT)  \
  //decl_param(TYPE, NAME, DESC, DEFAULT)
//#endif

typedef double REAL;

extern double *dr, *dz, dtheta;
extern double *w2k, *wk;
#define decl_field_comp(_DIR) \
  double (*ext_e_ ## _DIR) (double r, double z, double theta)

extern decl_field_comp(r);
extern decl_field_comp(z);
extern decl_field_comp(theta);

extern char *invok_name;

extern const double twopi;
extern const double invfourpi;
extern const double invpi32;

/**********
 * misc.c *
 **********/
void *xmalloc (size_t size);
void *xrealloc (void *ptr, size_t size);
void *xcalloc (size_t count, size_t size);

#define _CSTREAM_H_
#endif