From a98632df9d2d64b86968a9529b6c64eec5b52acf Mon Sep 17 00:00:00 2001 From: fedebenelli Date: Sun, 17 Nov 2024 21:56:05 -0300 Subject: [PATCH] test find hpll line --- test/test_hpll.f90 | 20 ++++++++++++++++++++ 1 file changed, 20 insertions(+) create mode 100644 test/test_hpll.f90 diff --git a/test/test_hpll.f90 b/test/test_hpll.f90 new file mode 100644 index 00000000..b88b9e9a --- /dev/null +++ b/test/test_hpll.f90 @@ -0,0 +1,20 @@ +program main + use yaeos + use fixtures_models, only: binary_PR76 + implicit none + + type(CubicEoS) :: model + type(PTEnvel2) :: env + real(pr) :: z(2), T, P + + model = binary_PR76() + + z = [0.5_pr, 0.5_pr] + + env = find_hpl(model, z, T0=500._pr, p0=150._pr) + + if(abs(env%points(1)%T - 152.55) > 0.01) error stop "HPLL Failed T1" + if(abs(env%points(1)%P - 150) > 0.1) error stop "HPLL Failed T1" + if (abs(env%points(size(env%points))%P) > 0.1) error stop "HPLL Failed Pend" + +end program \ No newline at end of file