OCPStructure.hpp

Go to the documentation of this file.

 
#ifndef __OCPSTRUCTURE_HEADER__
#define __OCPSTRUCTURE_HEADER__
 
// Standard header files
#include <vector>
 
// OpenCAEPoro header files
#include "OCPConst.hpp"
 
using namespace std;
 
class OCPRes
{
public:
    void Setup_IsoT(const OCP_USI& nb, const OCP_USI& nw, const OCP_USI& nc)
    {
        OCP_USI reslen = (nb + nw) * (nc + 1);
        resAbs.resize(reslen);
        resRelV.resize(nb);
        resRelN.resize(nb);
    }
    void SetupT(const OCP_USI& nb, const OCP_USI& nw, const OCP_USI& nc)
    {
        OCP_USI reslen = (nb + nw) * (nc + 2);
        resAbs.resize(reslen);
        resRelV.resize(nb);
        resRelN.resize(nb);
        resRelE.resize(nb);
    }
    void SetZero()
    {
        fill(resAbs.begin(), resAbs.end(), 0);
        fill(resRelV.begin(), resRelV.end(), 0);
        fill(resRelN.begin(), resRelN.end(), 0);
        fill(resRelE.begin(), resRelE.end(), 0);
        maxRelRes_V       = 0;
        maxRelRes_N       = 0;
        maxRelRes_E       = 0;
        maxWellRelRes_mol = 0;
        maxId_V           = 0;
        maxId_N           = 0;
        maxId_E           = 0;
    }
    void SetInitRes() { maxRelRes0_V = maxRelRes_V; }
 
    vector<OCP_DBL> resAbs;  
    vector<OCP_DBL> resRelV; 
    vector<OCP_DBL> resRelN; 
    vector<OCP_DBL> resRelE; 
    OCP_DBL maxRelRes0_V; 
    OCP_DBL maxRelRes_V;  
    OCP_DBL maxRelRes_N;  
    OCP_DBL maxRelRes_E;  
    OCP_DBL
        maxWellRelRes_mol; 
 
    // use negative number to represent well number (ToDo)
    OCP_INT maxId_V; 
    OCP_INT maxId_N; 
    OCP_INT maxId_E; 
};
 
#endif