IsothermalSolver.cpp

Go to the documentation of this file.

 
#include "IsothermalSolver.hpp"
 
void IsothermalSolver::SetupMethod(Reservoir& rs, const OCPControl& ctrl)
{
    method = ctrl.GetMethod();
 
    switch (method) {
        case IMPEC:
            impec.Setup(rs, LSolver, ctrl);
            break;
        case AIMc:
            aimc.Setup(rs, LSolver, ctrl);
            break;
        case FIMn:
            fim_n.Setup(rs, LSolver, ctrl);
            break;
        case FIM:
        default:
            fim.Setup(rs, LSolver, ctrl);
            break;
    }
}
 
void IsothermalSolver::InitReservoir(Reservoir& rs) const
{
    switch (method) {
        case IMPEC:
            impec.InitReservoir(rs);
            break;
        case FIM:
            fim.InitReservoir(rs);
            break;
        case FIMn:
            fim_n.InitReservoir(rs);
            break;
        case AIMc:
            aimc.InitReservoir(rs);
            break;
        default:
            OCP_ABORT("Wrong method type!");
    }
}
 
void IsothermalSolver::Prepare(Reservoir& rs, OCPControl& ctrl)
{
    switch (method) {
        case IMPEC:
            impec.Prepare(rs, ctrl);
            break;
        case FIMn:
            fim_n.Prepare(rs, ctrl.GetCurDt());
            break;
        case FIM:
            fim.Prepare(rs, ctrl.GetCurDt());
            break;
        case AIMc:
            aimc.Prepare(rs, ctrl.GetCurDt());
            break;
        default:
            OCP_ABORT("Wrong method type!");
    }
}
 
void IsothermalSolver::AssembleMat(const Reservoir& rs, OCPControl& ctrl)
{
    const OCP_DBL dt = ctrl.GetCurDt();
 
    GetWallTime timer;
    timer.Start();
 
    switch (method) {
        case IMPEC:
            impec.AssembleMat(LSolver, rs, dt);
            break;
        case FIMn:
            fim_n.AssembleMat(LSolver, rs, dt);
            break;
        case FIM:
            fim.AssembleMat(LSolver, rs, dt);
            break;
        case AIMc:
            aimc.AssembleMat(LSolver, rs, dt);
            break;
        default:
            OCP_ABORT("Wrong method type!");
    }
 
    ctrl.RecordTimeAssembleMat(timer.Stop() / 1000);
}
 
void IsothermalSolver::SolveLinearSystem(Reservoir& rs, OCPControl& ctrl)
{
    switch (method) {
        case IMPEC:
            impec.SolveLinearSystem(LSolver, rs, ctrl);
            break;
        case FIMn:
            fim_n.SolveLinearSystem(LSolver, rs, ctrl);
            break;
        case FIM:
            fim.SolveLinearSystem(LSolver, rs, ctrl);
            break;
        case AIMc:
            aimc.SolveLinearSystem(LSolver, rs, ctrl);
            break;
        default:
            OCP_ABORT("Wrong method type!");
    }
}
 
OCP_BOOL IsothermalSolver::UpdateProperty(Reservoir& rs, OCPControl& ctrl)
{
    OCP_BOOL flag;
 
    GetWallTime timer;
    timer.Start();
 
    switch (method) {
        case IMPEC:
            flag = impec.UpdateProperty(rs, ctrl);
            break;
        case FIMn:
            flag = fim_n.UpdateProperty(rs, ctrl);
            break;
        case FIM:
            flag = fim.UpdateProperty(rs, ctrl);
            break;
        case AIMc:
            flag = aimc.UpdateProperty(rs, ctrl);
            break;
        default:
            OCP_ABORT("Wrong method type!");
    }
 
    ctrl.RecordTimeUpdateProperty(timer.Stop() / 1000);
 
    return flag;
}
 
OCP_BOOL IsothermalSolver::FinishNR(Reservoir& rs, OCPControl& ctrl)
{
    switch (method) {
        case IMPEC:
            return impec.FinishNR(rs);
        case FIMn:
            return fim_n.FinishNR(rs, ctrl);
        case FIM:
            return fim.FinishNR(rs, ctrl);
        case AIMc:
            return aimc.FinishNR(rs, ctrl);
        default:
            OCP_ABORT("Wrong method type!");
    }
}
 
void IsothermalSolver::FinishStep(Reservoir& rs, OCPControl& ctrl)
{
    switch (method) {
        case IMPEC:
            impec.FinishStep(rs, ctrl);
            break;
        case FIMn:
            fim_n.FinishStep(rs, ctrl);
            break;
        case FIM:
            fim.FinishStep(rs, ctrl);
            break;
        case AIMc:
            aimc.FinishStep(rs, ctrl);
            break;
        default:
            OCP_ABORT("Wrong method type!");
    }
    ctrl.UpdateIters();
}
 
/*----------------------------------------------------------------------------*/
/*  Brief Change History of This File                                         */
/*----------------------------------------------------------------------------*/
/*  Author              Date             Actions                              */
/*----------------------------------------------------------------------------*/
/*  Shizhe Li           Oct/21/2021      Create file                          */
/*  Chensong Zhang      Jan/08/2022      Update Doxygen                       */
/*----------------------------------------------------------------------------*/