How to use CoolProp from within MATLAB - nagimov/cyclotron-cryopanel GitHub Wiki

Intro

CoolProp is a thermophysical property database and wrappers for a selection of programming environments. It offers similar functionality to REFPROP, but CoolProp is open-source and free.

MATLAB Wrapper (easy syntax, slow)

Full instructions: http://www.coolprop.org/coolprop/wrappers/MATLAB/index.html#matlab

If in the end MATLAB gives Error using loadlibrary No supported compiler or SDK was found error, "MATLAB Support for the MinGW-w64 C/C++ Compiler from TDM-GCC" needs to be installed through Add-ons Explorer. Getting files

D:\CoolProp_wrapper
 |- SwigRef.m
 |- CoolPropsetup.m
 |- CoolPropMATLAB_wrap.mexw64
 |- +CoolProp
    |- AbstractState.m
    |- DmassHmass_INPUTS.m
    |- ...

Usage example

  • the library path must be included to MATLAB files using addpath command as follow: addpath('D:\CoolProp_wrapper\')

  • so the MWE of CoolProp usage would look as follows:

addpath('D:\CoolProp_wrapper\')
CoolProp.PropsSI('T','P',101325,'Q',0,'Water')

C-style Access (C-like syntax, fast)

Full instructions: http://www.coolprop.org/coolprop/wrappers/MATLAB/index.html#faster-properties

D:\CoolProp_wrapper_fast
 |- libCoolProp.dll
 |- libCoolProp.dll.a
 |- CoolProp.dll
 |- CoolProp.lib
 |- CoolPropLib.h
 |- exports.txt
  • some MWE is as follows:
%% Calling low-level interface from MATLAB through shared library (DLL)
%%
%% Originally developed by Edo Macchi, modified by Ian Bell
%% 
%% December 2015

path_to_lib = 'D:\CoolProp_wrapper_fast'; %specify path to coolprop shared library
path_to_include= 'D:\CoolProp_wrapper_fast'; %specify path to coolprop's include folder

% Loading shared library
if ~libisloaded('coolprop') %checking whether library is already loaded
    addpath(path_to_lib)
    addpath(path_to_include)
    libname = 'libCoolProp' % OSX and linux
    if ispc
        libname = 'CoolProp'
    end
    loadlibrary(libname,'CoolPropLib.h','includepath',path_to_include,'alias','coolprop'); % loading library with alias coolprop
    disp('loaded CoolProp shared library.')
    disp('loaded these functions: ')
    libfunctions coolprop
end

buffer_size = 1000;
ierr = 0;
b = (1:1:buffer_size);
herr= char(b);

%Selecting backend and fluid
backend = 'BICUBIC&HEOS';
fluid = 'Water';
[handle, sh] = calllib('coolprop','AbstractState_factory',backend,fluid,ierr,herr,buffer_size);

length=100000;
input1 = linspace(700000.0,1500000.0,length)';
input2 = linspace(2.8e6,3.0e6,length);
[input_pair,sip] = calllib('coolprop','get_input_pair_index','HmassP_INPUTS');

tic;

%Creating input and output pointers
input1Ptr = libpointer('doublePtr',input1);
input2Ptr = libpointer('doublePtr',input2);
TPtr = libpointer('doublePtr',zeros(length,1));
pPtr = libpointer('doublePtr',zeros(length,1));
rhomolarPtr = libpointer('doublePtr',zeros(length,1));
hmolarPtr = libpointer('doublePtr',zeros(length,1));
smolarPtr = libpointer('doublePtr',zeros(length,1));

calllib('coolprop','AbstractState_update_and_common_out',handle,input_pair,input1Ptr,input2Ptr,length,TPtr,pPtr,rhomolarPtr,hmolarPtr,smolarPtr,ierr,herr,buffer_size);

dt = toc;

T=get(TPtr,'Value');
p=get(pPtr,'Value');
rhomolar=get(rhomolarPtr,'Value');
hmolar=get(hmolarPtr,'Value');
smolar=get(smolarPtr,'Value');

fprintf('%d us/call\n', (dt*1e6/length));

outputs=zeros(5,1);
%Choosing parameters to compute
[outputs(1,1), so1] = calllib('coolprop','get_param_index','T');
[outputs(2,1), so2] = calllib('coolprop','get_param_index','P');
[outputs(3,1), so3] = calllib('coolprop','get_param_index','Dmolar');
[outputs(4,1), so4] = calllib('coolprop','get_param_index','Hmolar');
[outputs(5,1), so5] = calllib('coolprop','get_param_index','Smolar');

%Creating ouput pointers
out1Ptr = libpointer('doublePtr',zeros(length,1));
out2Ptr = libpointer('doublePtr',zeros(length,1));
out3Ptr = libpointer('doublePtr',zeros(length,1));
out4Ptr = libpointer('doublePtr',zeros(length,1));
out5Ptr = libpointer('doublePtr',zeros(length,1));

tic;

calllib('coolprop','AbstractState_update_and_5_out',handle,input_pair,input1Ptr,input2Ptr,length,outputs,out1Ptr,out2Ptr,out3Ptr,out4Ptr,out5Ptr,ierr,herr,buffer_size);

dt=toc;

%Saving computed values to array
out1=get(out1Ptr,'Value');
out2=get(out2Ptr,'Value');
out3=get(out3Ptr,'Value');
out4=get(out4Ptr,'Value');
out5=get(out5Ptr,'Value');

fprintf('%d us/call\n', (dt*1e6/length));