github/CoolProp
1
1
Fork 0
mirror of https://github.com/CoolProp/CoolProp.git synced 2026-01-10 22:48:05 -05:00
Code Issues Packages Projects Releases Wiki Activity

Add test for #2582

Browse Source
This commit is contained in:
Ian Bell
2025-09-06 08:30:23 -04:00
parent 10ef4922d1
commit ae8886ac1e
1 changed files with 37 additions and 0 deletions
Download Patch File Download Diff File Expand all files Collapse all files

37
src/Tests/CoolProp-Tests.cpp
View File

@@ -2443,6 +2443,43 @@ TEST_CASE("Github issue #2491", "[2491]") {
CHECK(std::isfinite(AS->rhomolar()));
}
template <typename T>
std::vector<T> linspace(T start, T end, int num) {
std::vector<T> linspaced;
if (num <= 0) {
return linspaced; // Return empty vector for invalid num
}
if (num == 1) {
linspaced.push_back(start);
return linspaced;
}
T step = (end - start) / (num - 1);
for (int i = 0; i < num; ++i) {
linspaced.push_back(start + step * i);
}
return linspaced;
}
TEST_CASE("Github issue #2582", "[2582]") {
std::shared_ptr<CoolProp::AbstractState> AS(AbstractState::factory("HEOS", "CO2"));
double pc = AS->p_critical();
AS->update(PQ_INPUTS, 73.33e5, 0);
double hmass_liq = AS->saturated_liquid_keyed_output(iHmass);
double hmass_vap = AS->saturated_vapor_keyed_output(iHmass);
std::cout << pc << std::endl;
std::cout << hmass_liq << std::endl;
std::cout << hmass_vap << std::endl;
for (auto hmass: linspace(100e3, 700e3, 1000)){
CAPTURE(hmass);
CHECK_NOTHROW(AS->update(CoolProp::HmassP_INPUTS, hmass, 73.76e5));
}
for (auto hmass: linspace(100e3, 700e3, 1000)){
CAPTURE(hmass);
CHECK_NOTHROW(AS->update(CoolProp::HmassP_INPUTS, hmass, 73.33e5));
}
}
TEST_CASE("CoolProp.jl tests", "[2598]") {
// // Whoah, actually quite a few change meaningfully
// SECTION("Check pcrit doesn't change too much with SA on"){