globals.h

Go to the documentation of this file.

#ifndef _GLOBALS__
#define _GLOBALS__



#include <fe/fe_values.h>
#include <lac/vector.h>
#include <assert.h>
#include <vector>
#include <string>
#include <iostream>
#include <base/point.h>
#include <grid/tria.h>
#include <lac/full_matrix.h>
#include <lac/full_matrix.h>
#include <lac/block_sparse_matrix.h>
#include <limits.h>
#define DIM 3
#define TOLERANCE 1E-08
//#define INVALID_INDEX -1

using namespace dealii;

using std::string;
using std::cout;
#define INDEX_MAX UINT_MAX
typedef unsigned Index;
typedef Point<2> Point2D;
typedef Point<1> Point1D;
typedef Point<3> Point3D;
typedef std::vector<bool> VecBool;
typedef std::map<unsigned,double> MapIntDouble;
typedef Triangulation<DIM>::active_cell_iterator Cell;
typedef Triangulation<2>::active_cell_iterator Cell2D;
typedef std::vector<std::string> ArrayString;
typedef Triangulation<1>::active_cell_iterator Cell1D;
typedef Triangulation<DIM>::active_face_iterator Face;
typedef DoFHandler<DIM>::active_cell_iterator DoFCell;
typedef DoFHandler<DIM>::active_face_iterator DoFFace;
typedef std::vector<int> VecInt;
typedef std::vector<Index> VecIndex;
typedef BlockVector<double> BlockVecDouble;
typedef Vector<double> DealVecDouble;
typedef BlockSparseMatrix<double> BlockMatrix;
typedef std::vector<char> VecTag; 
typedef  std::vector<std::string> StringList;
enum Axis {X,Y,Z};
enum CellDirection3D {LEFT_CELL,RIGHT_CELL,FRONT_CELL,BACK_CELL,BOTTOM_CELL,UP_CELL};
enum VertexDirection3D{VERTEX_000,VERTEX_100,VERTEX_010,VERTEX_110,VERTEX_001,VERTEX_101,VERTEX_011,VERTEX_111};


/*
 *       *-------*        *-------*
 *      /|       |       /       /|
 *     / |   3   |      /   5   / |
 *    /  |       |     /       /  |
 *   *   |       |    *-------*   |
 *   | 0 *-------*    |       | 1 *
 *   |  /       /     |       |  /
 *   | /   4   /      |   2   | /
 *   |/       /       |       |/
 *   *-------*        *-------*
 *
 */
enum FaceDirection3D{LEFT_FACE=0,RIGHT_FACE=1,FRONT_FACE=2,BACK_FACE=3,BOTTOM_FACE=4,UP_FACE=5};

enum VertexDirection2D {BL_VERTEX=0,BR_VERTEX=1,UL_VERTEX=2,UR_VERTEX=3,INVALID_VERTEX=4};


enum SolverType{CG_ITERATIVE,DIRECT};



#define EQUAL(X,Y) (fabs(X-Y) < 1e-08)


#endif

---