shells/applications/cylinderElastic/helpers/DetonationPressure1D.h

Go to the documentation of this file.

// -*- C++ -*-
//
// ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
//
//                                   Fehmi Cirak
//                        California Institute of Technology
//                           (C) 2004 All Rights Reserved
//
// <LicenseText>
//
// ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
//
#ifndef DETONATIONPRESSURE1D_H
#define DETONATIONPRESSURE1D_H
#include <vector>
#include <map>
#include <fstream>
#include <iostream>
#include <iterator>
#include <utility>
#include <cassert>


namespace cylExp {
    struct DetonationPressure1D;
}


struct DetonationPressure1D {
private:
    typedef std::map<double, double> DDMap;
    
    DDMap   _positionPressure;

public:

    DetonationPressure1D() {
        std::ifstream ifs("./Pressure_180.txt");
        assert(ifs.is_open());

        std::istream_iterator<double> begin(ifs);
        std::vector<double> data(begin, std::istream_iterator<double>());
        
        const int lines = static_cast<int>(data.size()/2.);
        
        for (int i=0; i<lines; ++i) {
            std::pair<double, double> posPres(data[2*i], data[2*i+1]);
            _positionPressure.insert(posPres);
        }
    }


    double operator()(double time, double x, double y, double z) const {

        const double vcj = 2365;
        double zshift = z-time*vcj;
        
        DDMap::const_iterator it = _positionPressure.upper_bound(zshift);

        if (it==_positionPressure.end()) {
          return 0.0;           
        } else {
          return (it->second-80000.0);
        }
    }   

};

#endif

---