LBflow > Reference
Home Research Publications LBflow LB theory Downloads Examples Reference LBview Collaborators Contact Links Site last updated 21st June 2007	lattice.h Go to the documentation of this file. /*************************************************************************** * * * LBflow - an extensible lattice Boltzmann simulation framework * * * * Copyright (C) 2007 by Ed Llewellin * * * * This file may not be copied or shared without permission * * Contact ed@lbflow.co.uk or visit www.lbflow.co.uk * * * ***************************************************************************/ #ifndef LATTICE_H #define LATTICE_H #include <vector> #include <set> #include <map> #include <iostream> #ifdef HAVE_CONFIG_H #include <config.h> #endif #include "simulation.h" //forward declarations class Parser; class TripleInt; class TripleDouble; class Node; using namespace std; //TODO: lattice should deal with all spatial information. Should have a where() function which returns a vector of indices // of node positions that conform to some positional descriptor. This would work with, e.g. "where(location)" in .sim files // (parsed by nodegroup), and any other place where position info is required, such as in outputs. class Lattice { protected: //data members Simulation *sim; map<int, Node*> m; //NOTE: starts from 0 //constructor Lattice(Simulation*); public: //destructor virtual ~Lattice(){}; virtual void read(Parser*)=0; virtual void writeToSim(ostream&)=0; //transform between coordinate and indexing space //for discrete space virtual TripleInt indexToLattCoords(int)=0; virtual TripleDouble indexToPhysCoords(int)=0; virtual TripleDouble lattCoordsToPhysCoords(TripleInt)=0; virtual TripleInt physCoordsToLattCoords(TripleDouble)=0; virtual int lattCoordsToIndex(TripleInt, bool)=0; virtual int physCoordsToIndex(TripleDouble, bool)=0; //for continuous space virtual TripleDouble physCoordsToLattCoordsCont(TripleDouble, bool)=0; virtual TripleDouble lattCoordsToPhysCoordsCont(TripleDouble)=0; //get bookkeeping information about lattice virtual int getDimensionality()=0; virtual string getLatticeType()=0; virtual TripleInt getLattDim()=0; virtual TripleDouble getPhysDim()=0; virtual TripleDouble getOrigin()=0; virtual TripleDouble getCentre()=0; virtual int getSize()=0; virtual string getDescription()=0; virtual string getDescriptorForDomainOutput()=0; //set nodes virtual void put(int, Node*)=0; //get nodes virtual Node* get(int)=0; virtual Node* get(int, TripleInt)=0; virtual vector<Node*>* getAdjacentNodes(int, vector<Node*>*)=0; virtual vector<Node*>* getAdjacentNodes(int, vector<Node*>*, bool)=0; #ifdef ENABLE_MPI protected: //data members map<int, int>* mpi_procmap; set<int>* mpi_mynodes; map<int, map<int, int> >* mpi_boundaryneighbours; public: void mpiPartitionDomain(int, Simulation::PartitionType, map<int, int>*); void mpiSetProcMap(map<int, int>*, int); map<int, int>* mpiGetProcMap(); set<int>* mpiGetMyNodes(); void mpiCreateShuttles(); map<int, map<int, int> >* mpiGetBoundaryNeighbours() {return mpi_boundaryneighbours;}; bool mpiHasIndex(int i) {return m.find(i) != m.end();}; virtual vector<int>* mpiGetBoundaryNeighbours(int, vector<int>*)=0; virtual void mpiInitBoundaryMap()=0; protected: virtual void mpiPartitionTrivial(int, map<int, int>*)=0; virtual void mpiCreateMetisAdjacencyStructure(vector<int>*, vector<int>*)=0; virtual vector<int>* mpiGetAdjacentIndices(int, vector<int>*)=0; private: void mpiPartitionMetis(int, map<int, int>*); #endif /*ENABLE_MPI*/ }; #endif

>