Function CP_PhaseSI

wrappers/MathCAD/CoolPropMathcad.cpp:540–583

This code executes the user function CP_PhaseSI, which is a wrapper for the CoolProp.PhaseSI() function, used to the fluid phase dependent on the state

// This code executes the user function CP_PhaseSI, which is a wrapper for
// the CoolProp.PhaseSI() function, used to the fluid phase dependent on the state
static LRESULT CP_PhaseSI(LPMCSTRING PhaseStr,          // output (CoolProp phase string)
                          LPCMCSTRING InputName1,       // CoolProp InputName1
                          LPCCOMPLEXSCALAR InputProp1,  // CoolProp InputProp1
                          LPCMCSTRING InputName2,       // CoolProp InputName2
                          LPCCOMPLEXSCALAR InputProp2,  // CoolProp InputProp2
                          LPCMCSTRING FluidName)        // CoolProp Fluid
{
    std::string ph;
    std::string Prop1Name(InputName1->str);
    std::string Prop2Name(InputName2->str);
    std::string FluidString = FluidName->str;

    // check that the first scalar argument is real
    LRESULT r = CheckRealOrError(InputProp1, 2);
    if (r) return r;

    // check that the second scalar argument is real
    r = CheckRealOrError(InputProp2, 4);
    if (r) return r;

    // pass the arguments to the CoolProp.phaseSI() function
    ph = CoolProp::PhaseSI(InputName1->str, InputProp1->real, InputName2->str, InputProp2->real, FluidName->str);

    // Note: PhaseSI does not throw exceptions, but instead appends the global parameter "errstring"
    // to the returned phase string if there is an error and returns "unknown" for the phase string.
    //
    // Use string search to see if PhaseSI failed with an error message (appended ": <error>")...
    if (ph.find("unknown:") != std::string::npos) {
        std::string emsg = ph.substr(ph.find(":") + 2, ph.length() - ph.find(":") + 2);
        CoolProp::set_error_string(emsg);  // reset error string so Mathcad can retrieve it
        // Use PropsSI error handling logic to determine the specific error message but with adjusted argument
        // positions for PhaseSI (which has no OutputName argument and thus shifts the positions of the InputName
        // and InputProp arguments by 1)
        return HandlePropsSIError(emsg, Prop1Name, 1u);
    }

    // Must use MathcadAllocate(size) so Mathcad can track and release, using Helper routine above
    char* c = AllocMathcadString(ph);
    // assign the string to the function's output parameter
    PhaseStr->str = c;

    // normal return
    return 0;
}

// Helper: centralize text-matching logic for HAPropsSI error handling
static LRESULT HandleHAPropsSIError(const std::string& emsg, const std::string& OutName, const std::string& Prop1Name, const std::string& Prop2Name,