OCP.cpp

Go to the documentation of this file.

 
#include "OCP.hpp"
 
void OpenCAEPoro::ReadInputFile(const string& filename)
{
    ParamRead rp;
    rp.ReadInputFile(filename);
    // Read parameters from input file
    InputParam(rp);
}
 
void OpenCAEPoro::InputParam(ParamRead& param)
{
    reservoir.InputParam(param);
    control.InputParam(param.paramControl);
    output.InputParam(param.paramOutput);
}
 
void OpenCAEPoro::SetupSimulator(const USI& argc, const char* options[])
{
    GetWallTime timer;
    timer.Start();
 
    control.SetupFastControl(argc, options); // Read Fast control
 
    switch (control.GetModel()) {
        case ISOTHERMALMODEL:
            if (control.printLevel >= PRINT_MIN) {
                cout << endl << "Dynamic simulation for isothermal models" << endl;
            }
            break;
        case THERMALMODEL:
            if (control.printLevel >= PRINT_MIN) {
                cout << endl << "Dynamic simulation for thermal models" << endl;
            }
            break;
        default:
            OCP_ABORT("Wrong model type specified!");
    }
 
    solver.Setup(reservoir, control); // Setup static info for solver
 
    output.Setup(reservoir, control); // Setup output for dynamic simulation
 
    double finalTime = timer.Stop() / 1000;
    if (control.printLevel >= PRINT_MIN) {
        cout << endl
             << "Setup simulation done. Wall time : " << fixed << setprecision(3)
             << finalTime << " Sec" << endl;
    }
 
    control.RecordTotalTime(finalTime);
}
 
void OpenCAEPoro::InitReservoir()
{
    GetWallTime timer;
    timer.Start();
 
    solver.InitReservoir(reservoir);
 
    double finalTime = timer.Stop() / 1000;
    if (control.printLevel >= PRINT_MIN) {
        cout << endl
             << "Initialization done. Wall time : " << fixed << setprecision(3)
             << finalTime << " Sec" << endl;
    }
    control.RecordTotalTime(finalTime);
    control.initTime = finalTime;
}
 
void OpenCAEPoro::RunSimulation()
{
    switch (control.GetMethod()) {
        case IMPEC:
            if (control.printLevel >= PRINT_MIN) {
                cout << "\nDynamic simulation with IMPEC\n" << endl;
            }
            break;
        case FIM:
            if (control.printLevel >= PRINT_MIN) {
                cout << "\nDynamic simulation with FIM\n" << endl;
            }
            break;
        case FIMn:
            if (control.printLevel >= PRINT_MIN) {
                cout << "\nDynamic simulation with FIMn\n" << endl;
            }
            break;
        case AIMc:
            if (control.printLevel >= PRINT_MIN) {
                cout << "\nDynamic simulation with AIMc\n" << endl;
            }
            break;
        default:
            OCP_ABORT("Wrong method type is used!");
    }
 
    solver.RunSimulation(reservoir, control, output);
}
 
void OpenCAEPoro::OutputResults() const
{
    // find an appropriate size for printing times
    int fixWidth =
        MAX(log10(control.current_time), log10(MAX(control.totalSimTime, 1.0))) + 6;
 
    cout << "==================================================" << endl;
 
    // print numbers of steps
    cout << "Final time:                  " << right << fixed << setprecision(3)
         << setw(fixWidth) << control.current_time << " (Days)" << endl;
    cout << " - Avg time step size ......." << setw(fixWidth)
         << control.current_time / control.numTstep << " (" << control.numTstep
         << " steps)" << endl;
    cout << " - Avg Newton steps ........." << setw(fixWidth)
         << static_cast<double>(control.iterNR_total) / control.numTstep << " ("
         << control.iterNR_total << " succeeded + " << control.wastedIterNR
         << " wasted)" << endl;
    cout << " - Avg linear steps ........." << setw(fixWidth)
         << static_cast<double>(control.iterLS_total) / control.iterNR_total << " ("
         << control.iterLS_total << " succeeded + " << control.wastedIterLS
         << " wasted)" << endl;
 
    // print time usages
    cout << "Simulation time:             " << setw(fixWidth) << control.totalSimTime
         << " (Seconds)" << endl;
    cout << " - % Initialization ........." << setw(fixWidth)
         << 100.0 * control.initTime / control.totalSimTime << " (" << control.initTime
         << "s)" << endl;
    cout << " - % Assembling ............." << setw(fixWidth)
         << 100.0 * control.totalAssembleMatTime / control.totalSimTime << " ("
         << control.totalAssembleMatTime << "s)" << endl;
    cout << " - % Linear solver .........." << setw(fixWidth)
         << 100.0 * control.totalLStime / control.totalSimTime << " ("
         << control.totalLStime << "s)" << endl;
    cout << " - % Updating properties ...." << setw(fixWidth)
         << 100.0 * control.totalUpdatePropertyTime / control.totalSimTime << " ("
         << control.totalUpdatePropertyTime << "s)" << endl;
    cout << " - % Scheduled output ......." << setw(fixWidth)
         << 100.0 * output.outputTime / control.totalSimTime << " ("
         << output.outputTime << "s)" << endl;
 
    cout << "==================================================" << endl;
 
    output.PrintInfo();
}
 
/*----------------------------------------------------------------------------*/
/*  Brief Change History of This File                                         */
/*----------------------------------------------------------------------------*/
/*  Author              Date             Actions                              */
/*----------------------------------------------------------------------------*/
/*  Shizhe Li           Oct/01/2021      Create file                          */
/*  Chensong Zhang      Dec/05/2021      Format file                          */
/*----------------------------------------------------------------------------*/