lb2d_prime.h

基于格子Boltzmann方法开源可视化软件的源代码 具有很高的实用价值。对学习LBM方法及其软件开发非常游泳

//##############################################################################
//
// Copyright (C), 2005, Michael Sukop and Danny Thorne
//
// lb2d_prime.h
//
//  - Header file for lb2d_prime.
//

#ifndef MCMP_PRIME_H
#define MCMP_PRIME_H

#include <stdio.h>
#include <stdlib.h>
//#include <malloc.h>
#include <memory.h>
#include <assert.h>
#include <math.h>
#include <string.h>
#include <time.h>

#define IJ2N(_i_,_j_) (_j_)*lattice->param.LX + (_i_)

#include "flags.h"
#include "bc_flags.h"
#include "ic_flags.h"
#include "lattice.h"
#ifdef PARALLEL
#include "lbmpi.h"
#endif /* PARALLEL */
#include "forward_declarations.h"
#ifdef PARALLEL
#include "lbmpi.c"
#endif /* PARALLEL */
#include "params.h"

int vx[9] = { 0, 1, 0, -1,  0, 1, -1, -1,  1};
int vy[9] = { 0, 0, 1,  0, -1, 1,  1, -1, -1};

#define EPS .0000000000001
#define PI  3.1415926535

#define HRULE0 "- - - - - - - - - - - - - - - - - - - - " \
               "- - - - - - - - - - - - - - - - - - - - "

#define HRULE1 "----------------------------------------" \
               "----------------------------------------"

#define HRULE2 "========================================" \
               "========================================"

#if DO_NOT_STORE_SOLIDS
// TODO: Incorporate version that omits storage of interior solids (which
// are not involved in flow).
//#include "min_nodes/compute.c"
//#include "min_nodes/stream.c"
//#include "min_nodes/bcs.c"
//#include "min_nodes/collide.c"
//#include "min_nodes/lbio.c"
//#include "min_nodes/latman.c"
#else /* !( DO_NOT_STORE_SOLIDS) */
#include "compute.c"
#include "stream.c"
#include "collide.c"
#include "bcs.c"
#include "lbio.h"
#include "lbio.c"
#include "latman.c"
#include "runman.c"
#endif /* DO_NOT_STORE_SOLIDS */

void report_flags()
{
  struct report_struct flags;
  report_open(          &flags, "./out/flags");
  report_integer_entry( &flags, "VERBOSITY_LEVEL", VERBOSITY_LEVEL, "");
  report_integer_entry( &flags, "SAY_HI", SAY_HI, "");
  report_integer_entry( &flags, "NUM_FLUID_COMPONENTS", NUM_FLUID_COMPONENTS, "");
  report_integer_entry( &flags, "INAMURO_SIGMA_COMPONENT", INAMURO_SIGMA_COMPONENT, "");
  report_integer_entry( &flags, "ZHANG_AND_CHEN_ENERGY_TRANSPORT", ZHANG_AND_CHEN_ENERGY_TRANSPORT, "");
  report_integer_entry( &flags, "POROUS_MEDIA", POROUS_MEDIA, "");
  report_integer_entry( &flags, "STORE_UEQ", STORE_UEQ, "");
  report_integer_entry( &flags, "DO_NOT_STORE_SOLIDS", DO_NOT_STORE_SOLIDS, "");
  report_integer_entry( &flags, "NON_LOCAL_FORCES", NON_LOCAL_FORCES, "");
  report_integer_entry( &flags, "MANAGE_BODY_FORCE", MANAGE_BODY_FORCE, "");
  report_integer_entry( &flags, "STORE_BTC", STORE_BTC, "");
  report_integer_entry( &flags, "DETERMINE_FLOW_DIRECTION", DETERMINE_FLOW_DIRECTION, "");
  report_integer_entry( &flags, "BC_SOLID_NODE", BC_SOLID_NODE, "");
  report_integer_entry( &flags, "BC_FLUID_NODE", BC_FLUID_NODE, "");
  report_integer_entry( &flags, "BC_SLIP_NODE", BC_SLIP_NODE, "");
  report_integer_entry( &flags, "BC_FILM_NODE", BC_FILM_NODE, "");
  report_integer_entry( &flags, "BC_PRESSURE_N_IN", BC_PRESSURE_N_IN, "");
  report_integer_entry( &flags, "BC_PRESSURE_S_IN", BC_PRESSURE_S_IN, "");
  report_integer_entry( &flags, "BC_PRESSURE_E_IN", BC_PRESSURE_E_IN, "");
  report_integer_entry( &flags, "BC_PRESSURE_W_IN", BC_PRESSURE_W_IN, "");
  report_integer_entry( &flags, "BC_PRESSURE_N_OUT", BC_PRESSURE_N_OUT, "");
  report_integer_entry( &flags, "BC_PRESSURE_S_OUT", BC_PRESSURE_S_OUT, "");
  report_integer_entry( &flags, "BC_PRESSURE_E_OUT", BC_PRESSURE_E_OUT, "");
  report_integer_entry( &flags, "BC_PRESSURE_W_OUT", BC_PRESSURE_W_OUT, "");
  report_integer_entry( &flags, "BC_VELOCITY_N_IN", BC_VELOCITY_N_IN, "");
  report_integer_entry( &flags, "BC_VELOCITY_S_IN", BC_VELOCITY_S_IN, "");
  report_integer_entry( &flags, "BC_VELOCITY_E_IN", BC_VELOCITY_E_IN, "");
  report_integer_entry( &flags, "BC_VELOCITY_W_IN", BC_VELOCITY_W_IN, "");
  report_integer_entry( &flags, "BC_VELOCITY_N_OUT", BC_VELOCITY_N_OUT, "");
  report_integer_entry( &flags, "BC_VELOCITY_S_OUT", BC_VELOCITY_S_OUT, "");
  report_integer_entry( &flags, "BC_VELOCITY_E_OUT", BC_VELOCITY_E_OUT, "");
  report_integer_entry( &flags, "BC_VELOCITY_W_OUT", BC_VELOCITY_W_OUT, "");
  report_integer_entry( &flags, "WRITE_MACRO_VAR_DAT_FILES", WRITE_MACRO_VAR_DAT_FILES, "");
  report_integer_entry( &flags, "WRITE_RHO_AND_U_TO_TXT", WRITE_RHO_AND_U_TO_TXT, "");
  report_integer_entry( &flags, "WRITE_PDF_DAT_FILES", WRITE_PDF_DAT_FILES, "");
  report_integer_entry( &flags, "WRITE_PDF_TO_TXT", WRITE_PDF_TO_TXT, "");
  report_integer_entry( &flags, "INACTIVE_NODE", INACTIVE_NODE, "");
  report_integer_entry( &flags, "PUKE_NEGATIVE_DENSITIES", PUKE_NEGATIVE_DENSITIES, "");
  report_integer_entry( &flags, "SOLID_COLOR_IS_CHECKERBOARD", SOLID_COLOR_IS_CHECKERBOARD, "");
  report_integer_entry( &flags, "SOLID_COLOR_IS_BLACK", SOLID_COLOR_IS_BLACK, "");
  report_integer_entry( &flags, "DELAY", DELAY, "");
  report_integer_entry( &flags, "END_GRAV", END_GRAV, "");
  report_integer_entry( &flags, "MARK_ORIGIN_FOR_REFERENCE", MARK_ORIGIN_FOR_REFERENCE, "");
  report_integer_entry( &flags, "PERTURBATIONS", PERTURBATIONS, "");
  report_integer_entry( &flags, "WRITE_CHEN_DAT_FILES", WRITE_CHEN_DAT_FILES, "");
  report_integer_entry( &flags, "IC_UNIFORM_RHO_A", IC_UNIFORM_RHO_A, "");
  report_integer_entry( &flags, "IC_UNIFORM_RHO_B", IC_UNIFORM_RHO_B, "");
  report_integer_entry( &flags, "IC_UNIFORM_RHO_IN", IC_UNIFORM_RHO_IN, "");
  report_integer_entry( &flags, "IC_BUBBLE", IC_BUBBLE, "");
  report_integer_entry( &flags, "IC_DIAGONAL", IC_DIAGONAL, "");
  report_integer_entry( &flags, "IC_2X2_CHECKERS", IC_2X2_CHECKERS, "");
  report_integer_entry( &flags, "IC_STATIC", IC_STATIC, "");
  report_integer_entry( &flags, "IC_RECTANGLE", IC_RECTANGLE, "");
  report_integer_entry( &flags, "IC_DOT", IC_DOT, "");
  report_integer_entry( &flags, "IC_WOLF_GLADROW_DIFFUSION", IC_WOLF_GLADROW_DIFFUSION, "");
  report_integer_entry( &flags, "IC_YIN_YANG", IC_YIN_YANG, "");
  report_integer_entry( &flags, "IC_HYDROSTATIC", IC_HYDROSTATIC, "");
  report_close(         &flags);
}

#endif /* MCMP_PRIME_H */