meshManager.h

Go to the documentation of this file.

 
#pragma once
 
#include "dv.h"
#include <vector>
 
class domain;
 
 
using namespace std;
 
class meshManager {
 
 
    public:
 
        domain                   *domn;     
 
        vector<dv*>              phi;       
 
        vector<double>           dx;        
        vector<int>              mark;      
 
        vector<double>           xf;        
        vector<vector<double> >  yf;        
        vector<double>           xnf;       
        vector<double>           X;         
 
        int                      ngrd;      
        int                      ngrdf;     
 
        int                      iLower;    
        int                      iUpper;    
        double                   posLower;  
        double                   posUpper;  
 
        vector<double>           lastDA;    
 
 
        void adaptGrid(int iLower, int iUpper);       // Public interface
 
        inline bool operator()(const int &a, const int &b) const {
            return dx[mark[a]] < dx[mark[b]];
        }                                                // Functor used for sorting
 
        void adaptAfterSufficientDiffTime(const double &time,
                                          double       &tLastDA,
                                          int          &cLastDA,
                                          double       &dtCUmax);
 
        void adaptEddyRegionOfMesh(const double &time,
                                   double &tLastDA,
                                   int &cLastDA);
 
        void splitCell(const int isplt,
                       const int nsplt,
                       const vector<double> &cellFaces);
 
        void merge2cells(const int imrg, const bool LconstVolume=false);
 
        void enforceDomainSize();
        void removeFaceNearZero();
        void makeCellWithZeroSymmetric();
        void splitCellWithZero();
 
        //--------------- domain operations
 
        void setGridDxc(const domain *line, vector<double> &dxc, double C);
        void setGridDx( const domain *line, vector<double>  &dx);
 
        void setGridFromDxc(const vector<double> &dxc2);
 
 
 
    private:
 
        void adaptGrid_details(const int iLower,
                               const int iUpper);        // called by adaptGrid
        void mergeSmallCells();
        void mergeSmallCellsMP();
        void impose2point5rule();
        void splitLargeCells();
 
        void fix2point5offender(const int mPos,
                                const int &iglobal);
        void setDxArray();
 
        int  findPos(const vector<double> &x,
                     const double         val,
                     const int            &istart);
 
        void interp1pt(const vector<double> &x,
                       const vector<double> &y,
                       const double         &xval,
                             double         &yval,
                             int            &istart);
 
        void interpVec(const vector<double> &x,
                       const vector<double> &y,
                       const vector<double> &xn,
                             vector<double> &yn);
 
        double calcDistance(const vector<double>          &x,
                            const vector<vector<double> > &y,
                                  vector<double>          &sDist);
        void set_iLower_iUpper();
 
        void updateDA(const double &time, double &tLastDA, int &cLastDA,
                      int iStart, int iEnd);
 
 
 
    public:
 
        void init(domain *p_domn, const vector<dv*> p_phi);
        meshManager(){};
        ~meshManager(){}
 
};