#include "mex.h"   /*--This one is required*/

/* Prototype for the Props function to be called.  Can't use the CoolProp.h header because there are a lot of
   c++ parts in the header that cannot be easily hidden when compiling */
double Props(char *Output, char Name1, double Prop1, char Name2, double Prop2, char * Ref);
double Props1(char *Output, char * Ref);
long get_global_param_string(char*, char*);
long get_fluid_param_string(char *fluid, char *param, char * Output);
long get_standard_unit_system(void);
void set_standard_unit_system(long);
#include "GlobalConstants.h"
#include "float.h"

bool ValidNumber(double x)
{
    // Idea from http://www.johndcook.com/IEEE_exceptions_in_cpp.html
    return (x <= DBL_MAX && x >= -DBL_MAX);
};

void mexFunction(int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[])
{
    int npol,Nel;
    int m,n,u,v,k,Output_len,Name1_len,Name2_len,Ref_len,Param_len;
    int *c,status;
    double *d,Prop1,Prop2,x,y;
	char *Output,*Name1,*Name2,*Ref,*Param,errstr[1000],errstr2[1000],fluidslist[10000];
    double val;
    mxArray *cMat[1];
    
    if (nrhs == 2 && mxIsChar (prhs[0]) && mxIsChar (prhs[1]))
    {
        /* Get the refrigerant (it is a string) (+1 for the NULL terminator)*/
        Ref_len=(mxGetM(prhs[0]) * mxGetN(prhs[0])) + 1;
        Ref = mxCalloc(Ref_len, sizeof(char));
        status = mxGetString(prhs[0], Ref, Ref_len);
        
        /* Get the output (it is a string) (+1 for the NULL terminator)*/
        Output_len=(mxGetM(prhs[1]) * mxGetN(prhs[1])) + 1;
        Output = mxCalloc(Output_len, sizeof(char));
        status = mxGetString(prhs[1], Output, Output_len);
        
        /* Try to shortcut to get the strings for the fluid */
        if (!strcmp(Output,"aliases") || !strcmp(Output,"CAS") || !strcmp(Output,"CAS_number") || !strcmp(Output,"ASHRAE34") || !strcmp(Output,"REFPROPName") || !strcmp(Output,"REFPROP_name") || !strcmp(Output,"TTSE_mode"))
        {
            get_fluid_param_string(Ref,Output,fluidslist);
            plhs[0] = mxCreateString(fluidslist);
            return;
        }
        else if (!strcmp(Output,"enable_TTSE"))
        {
            enable_TTSE_LUT(Ref);
            return;
        }
        else if (!strcmp(Output,"disable_TTSE"))
        {
            disable_TTSE_LUT(Ref);
            return;
        }
        
        /* Ok, now try to get the values for the fluid */
        
        /* Create matrix for the return argument. */
        plhs[0] = mxCreateDoubleMatrix(1,1, mxREAL);
        
        /* Get the value*/
        val = Props1(Ref,Output);
        /* If it is a good value, return it*/
        if (ValidNumber(val))
        {
            *mxGetPr(plhs[0]) = val; 
        }
        /* Otherwise there was an error, return the CoolProp error*/
        else
        {
            get_global_param_string("errstring",errstr);
            sprintf(errstr2,"CoolProp Error: %s",errstr);
            mexErrMsgTxt(errstr2);
        }
    }
    else if (nrhs == 6 && mxIsChar (prhs[0]) && mxIsChar (prhs[1])  && mxIsChar (prhs[3])  && mxIsChar (prhs[5])  )
    {
        /* Create matrix for the return argument. */
        plhs[0] = mxCreateDoubleMatrix(1,1, mxREAL);
        
        /* Get the refrigerant (it is a string) (+1 for the NULL terminator)*/
        Ref_len=(mxGetM(prhs[5]) * mxGetN(prhs[5])) + 1;
        Ref = mxCalloc(Ref_len, sizeof(char));
        
        /* Get the output (it is a string) (+1 for the NULL terminator)*/
        Output_len=(mxGetM(prhs[0]) * mxGetN(prhs[0])) + 1;
        Output = mxCalloc(Output_len, sizeof(char));
        
        status = mxGetString(prhs[5], Ref, Ref_len);
        status = mxGetString(prhs[0], Output, Output_len);
        
        Name1_len=(mxGetM(prhs[1]) * mxGetN(prhs[1])) + 1;
        Name1 = mxCalloc(Name1_len, sizeof(char));
        mxGetString(prhs[1], Name1, Name1_len);
        
        Name2_len=(mxGetM(prhs[3]) * mxGetN(prhs[3])) + 1;
        Name2 = mxCalloc(Name2_len, sizeof(char));
        mxGetString(prhs[3], Name2, Name2_len);
        
        Prop1 = mxGetScalar(prhs[2]);
        Prop2 = mxGetScalar(prhs[4]);

        val = Props(Output,Name1[0],Prop1,Name2[0],Prop2,Ref);
        
        /* If it is a good value, return it */
        if (ValidNumber(val))
        {
            *mxGetPr(plhs[0]) = val; 
        }
        /* Otherwise there was an error, return the CoolProp error*/
        else
        {
            get_global_param_string("errstring",errstr);
            sprintf(errstr2,"CoolProp Error: %s",errstr);
            mexErrMsgTxt(errstr2);
        }
    }
    else if (nrhs == 1 && mxIsChar (prhs[0]))
    {
        Ref_len=(mxGetM(prhs[0]) * mxGetN(prhs[0])) + 1;
        Ref = mxCalloc(Ref_len, sizeof(char));
        status = mxGetString(prhs[0], Ref, Ref_len);
        
        if (!strcmp(Ref,"FluidsList")){
            get_global_param_string("FluidsList",fluidslist);
        }
        else if (!strcmp(Ref,"version")){
            get_global_param_string("version",fluidslist);
        }
        else if (!strcmp(Ref,"gitrevision")){
            get_global_param_string("gitrevision", fluidslist);
        }
        else if (!strcmp(Ref,"set_UNIT_SYSTEM_SI"))
        {
            set_standard_unit_system(UNIT_SYSTEM_SI);
            mexPrintf("Unit system set to SI\n");
            return;
        }
        else if (!strcmp(Ref,"set_UNIT_SYSTEM_KSI"))
        {
            set_standard_unit_system(UNIT_SYSTEM_KSI);
            mexPrintf("Unit system set to KSI\n");
            return;
        }
        else if (!strcmp(Ref,"get_unit_system"))
        {
            switch (get_standard_unit_system())
            {
                case UNIT_SYSTEM_SI:
                    plhs[0] = mxCreateString("UNIT_SYSTEM_SI"); break;
                case UNIT_SYSTEM_KSI:
                    plhs[0] = mxCreateString("UNIT_SYSTEM_SI"); break;
            }
            return;
        }
        else
        {
            sprintf(errstr2,"single input is invalid: %s",Ref);
            mexErrMsgTxt(errstr2);
        }
        plhs[0] = mxCreateString(fluidslist);
        return;
    }
    else
    {
        mexErrMsgTxt("Props must either receive two strings or the signature Props(Output,Param1,Value1,Param2,Value2,FluidName)");
    }
}