src/expde/paramete.h

Go to the documentation of this file.

//    expde: expression templates for partial differential equations.
//    Copyright (C) 2001  Christoph Pflaum
//    This program is free software; you can redistribute it and/or modify
//    it under the terms of the GNU General Public License as published by
//    the Free Software Foundation; either version 2 of the License, or
//    (at your option) any later version.
//
//    This program is distributed in the hope that it will be useful,
//    but WITHOUT ANY WARRANTY; without even the implied warranty of
//    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
//    GNU General Public License for more details.
//
//                 SEE  Notice1.doc made by 
//                 LAWRENCE LIVERMORE NATIONAL LABORATORY
//

// ------------------------------------------------------------
//
//  First parameter file: parameter.h
//
// ------------------------------------------------------------

#ifndef PARAMETER_H_
#define PARAMETER_H_     

#ifdef COMP_GNUOLD

#else
 using std::cout;
 using std::endl;
#endif 
         

// For largest problems we need long long
// to index the grid points
#define long_int long long
// #define long_int int
              
// Some compilers need these definitions
//#define bool  char
//#define true  1
//#define false 0


// Choose one: MG restriction operator 
//             coefficients stored or not
// Similar implementation for prolongation operators is missing
#define mg_coefficients false
// #define mg_coefficients true


// maximale number of levels:
//----------------------------
#define Max_number_levels 20


// developer version with    error messages (slow) oder
// working   version without error messages (fast)
//----------------------------
//#define developer_version true
#define developer_version false


// multilevel variable? yes-no
//------------------------------
// no multilevel variable
// #define Parameter_MG_Variable 0
// with multilevel variable
#define Parameter_MG_Variable ebene


// size of hash tables
//--------------------------------
// start size of  hash table 0
#define hashtable0_lenght 300000

// factor for size of hash tables
#define Factor_hashtable0_lenght 2
#define Factor_hashtable1_lenght 2
#define Factor_hashtable2_lenght 2
#define Factor_hashtable3_lenght 2
#define Factor_hashtable4_lenght 1

#define Factor_hashtable_indices_lenght 2


// tolerance in domain.h
//-----------------------------
#define domain_eps 0.000000001

// parameters for boundary labels:
//------------------------------------
// a) shift in test of point is marked
#define shift_eps 0.001
// b) epsilon for v(p) > 0.0
#define eps_boundary_var 0.001

// factor for minimal and maximal distance to the boundary
//---------------------------------------------------------
#define factor_boundary_dis 0.3

// offset for surrounded square
//-----------------------------------
#define standard_offset_square 0.5132455
#define standard_stretch_square 1.00143
// #define standard_stretch_square 1.0

// standard parameter for parallelization
//----------------------------------------
#define standard_r_parallel 2

// #define standard_offset_square 0
// #define standard_stretch_square 1.0

// for boundary labels
//---------------------
#define TypLabel int           // type of storage
#define Max_label_number 30    // maximal number of boundary labels


// for finding bug in Tetraeder construction
//-------------------------------------------
#define Print_wrong_Tet false       
//#define Print_wrong_Tet true

// to make the bounding cube a bit larger in sca_dom.cc
//-----------------------------------------------------
#define Bounding_cube_enlargement 1.1


// Print status of grid generator
//-----------------------------------------------------
#define print_status_of_grid_generator false
//#define print_status_of_grid_generator true

// For grid generator
//-----------------------------------------------------
//#define pre_check_box_grid_generation false
#define pre_check_box_grid_generation true



// Tell about round off problems in subdivision of boundary cells
// such problems can only appear in the parallel version
//----------------------------------------------------------------
#define tell_round_off_problem false
//#define tell_round_off_problem true
#define eps_round_off_problem 0.000000000000000001

// Set periodic data structure for 
// periodic boundary conditions
// #define PERIODIC true
#endif
          

---