dv_pos.cc

Go to the documentation of this file.

 
#include "dv_pos.h"
#include "domain.h"
#include <cstdlib>
 
 
dv_pos::dv_pos(domain    *line,
               const      string s,
               const bool Lt,
               const bool Lo) {
 
    domn          = line;
    var_name      = s;
    L_transported = Lt;
    L_output      = Lo;
    d             = vector<double>(domn->ngrd, 0.0);
 
    double dx = domn->pram->domainLength / domn->ngrd;
    d.at(0) = domn->pram->xDomainCenter + 0.5*(-domn->pram->domainLength + dx);
    for(int i=1; i<domn->ngrd; i++)
        d.at(i) = d.at(i-1) + dx;
 
}
 
 
void dv_pos::splitCell(const int isplt,
                       const int nsplt,
                       const vector<double> &cellFaces) {
 
    d.insert( d.begin() +isplt+1, nsplt, 0.0);
    for(int i=isplt, j=0; i<=isplt+nsplt; i++,j++)
        d.at(i) = 0.5*(cellFaces.at(j)+cellFaces.at(j+1));
 
}
 
 
void dv_pos::merge2cells(const int    imrg,
                         const double m1,
                         const double m2,
                         const bool   LconstVolume) {
 
    if(LconstVolume || domn->pram->bcType=="WALL") {   
        d.erase(d.begin() + imrg+1);
        d.at(imrg) = 0.5*(domn->posf->d.at(imrg) + domn->posf->d.at(imrg+1));
    }
    else
        setVar();
}
 
 
void dv_pos::setVar(const int ipt){
 
    if(ipt != -1) {
        cout << endl << "ERROR in setVar: ipt must be = -1" << endl;
        exit(0);
    }
 
    d.resize(domn->posf->d.size()-1);
 
    for(int i=0; i<d.size(); i++)
        d.at(i) = 0.5*(domn->posf->d.at(i) + domn->posf->d.at(i+1));
 
}
 
 
void dv_pos::setDvFromRegion(const int i1, const int i2){
 
    // note, we are owned by the eddyline, so domn is eddl, so to get domn data, use domn->domn
    const vector<double> &domn_data = domn->domn->varMap.find(var_name)->second->d;
 
    if(i2 >= i1)
        d.assign(domn_data.begin()+i1, domn_data.begin()+i2+1  );
    else {           // wrap around (periodic assignment)
        d.assign(domn_data.begin()+i1, domn_data.end());
        d.insert(d.end(), domn_data.begin(), domn_data.begin()+i2+1 );
        int idmb = d.size();
        for(int i=idmb; i<d.size(); i++)
            d.at(i)+=domn->domn->Ldomain();
    }
 
}