settings.h

Go to the documentation of this file.

/* ------------------------------------------------------------
 * This is the file "settings.h" of the H2Lib package.
 * All rights reserved, Steffen Boerm 2009
 * ------------------------------------------------------------ */

#ifndef SETTINGS_H
#define SETTINGS_H

#include <math.h>
#include <stddef.h>
#ifdef USE_COMPLEX
#include <complex.h>
#endif

#ifdef USE_SIMD
#include "simd.h"
#endif

/* ------------------------------------------------------------
 * Compilation settings
 * ------------------------------------------------------------ */

#ifdef __cplusplus
#define INLINE_PREFIX static extern "C"
#else
#define INLINE_PREFIX static
#endif

#ifdef __cplusplus
#define HEADER_PREFIX extern "C"
#else
#define HEADER_PREFIX
#endif

#ifdef __cplusplus
#define IMPORT_PREFIX extern "C"
#else
#define IMPORT_PREFIX
#endif

/* ------------------------------------------------------------
 * Types
 * ------------------------------------------------------------ */

typedef unsigned short bool;

extern const bool true;

extern const bool false;

typedef unsigned uint;

extern const int i_zero;

extern const int i_one;

extern const uint u_zero;

extern const uint u_one;

typedef size_t longindex;

#ifdef USE_FLOAT
typedef float real;
#else
typedef double real;
#endif

#ifndef VFLOAT
#define VFLOAT 1
#endif

#ifndef VDOUBLE
#define VDOUBLE 1
#endif

#ifndef VREAL
#define VREAL 1
#endif

#ifdef USE_FLOAT

#define H2_CHECK_TOLERANCE 1.0e-6

#define H2_ALMOST_ZERO 1e-30
#else

#define H2_CHECK_TOLERANCE 1.0e-12

#define H2_ALMOST_ZERO 1e-300
#endif

#ifdef USE_FLOAT
#define SCANF_PREFIX ""
#else
#define SCANF_PREFIX "l"
#endif

typedef real *preal;

typedef const real *pcreal;

extern const real r_zero;

extern const real r_one;

extern const real r_minusone;

#ifdef USE_FLOAT
#ifdef USE_COMPLEX
typedef float _Complex field;
#else
typedef float field;
#endif
#else
#ifdef USE_COMPLEX
typedef double _Complex field;
#else
typedef double field;
#endif
#endif

#ifndef VFIELD
#define VFIELD 1
#endif

typedef field *pfield;

typedef const field *pcfield;

extern const field f_zero;

extern const field f_one;

extern const field f_minusone;

#ifdef USE_COMPLEX

extern const field f_i;
#endif

#ifndef USE_COMPLEX
#ifdef I
#undef I
#endif
#define I 0.0
#endif

extern const char *_h2_ntrans;

extern const char *_h2_trans;

extern const char *_h2_adj;

extern const char *_h2_left;

extern const char *_h2_right;

extern const char *_h2_lower;

extern const char *_h2_upper;

extern const char *_h2_unit;

extern const char *_h2_nonunit;

extern const char *_h2_vectors;

extern const char *_h2_skinnyvectors;

extern const char *_h2_novectors;

/****************************************************
 * Define conversion specifier and argument macros
 * for different @ref field types.
 ****************************************************/

#ifdef USE_FLOAT
#ifdef USE_COMPLEX
#define FIELD_CS(format, cs) "%" #format #cs " + %" #format #cs "i"
#else
#define FIELD_CS(format, cs) "%" #format #cs
#endif
#else
#ifdef USE_COMPLEX
#define FIELD_CS(format, cs) "%" #format #cs " + %" #format #cs "i"
#else
#define FIELD_CS(format, cs) "%" #format #cs
#endif
#endif

#ifdef USE_FLOAT
#ifdef USE_COMPLEX
#define FIELD_SCANF_CS(format, cs) "%" #format #cs " + %" #format #cs "i"
#else
#define FIELD_SCANF_CS(format, cs) "%" #format #cs
#endif
#else
#ifdef USE_COMPLEX
#define FIELD_SCANF_CS(format, cs) "%l" #format #cs " + %l" #format #cs "i"
#else
#define FIELD_SCANF_CS(format, cs) "%l" #format #cs
#endif
#endif

#ifdef USE_FLOAT
#ifdef USE_COMPLEX
#define FIELD_ARG(z) crealf(z), cimagf(z)
#else
#define FIELD_ARG(z) z
#endif
#else
#ifdef USE_COMPLEX
#define FIELD_ARG(z) creal(z), cimag(z)
#else
#define FIELD_ARG(z) z
#endif
#endif

#ifdef USE_FLOAT
#ifdef USE_COMPLEX
#define FIELD_ADDR(z) (real*)(z), ((real*)(z)+1)
#else
#define FIELD_ADDR(z) (real*)(z)
#endif
#else
#ifdef USE_COMPLEX
#define FIELD_ADDR(z) (real*)(z), ((real*)(z)+1)
#else
#define FIELD_ADDR(z) (real*)(z)
#endif
#endif

typedef enum {
  AMATRIX = 0, 

  HMATRIX = 1, 

  H2MATRIX = 2, 

  SPARSEMATRIX = 3, 

  DH2MATRIX = 4 
} matrixtype;

#endif