johnpye/fprops/test.h

#ifndef HELMHOLTZ_TEST_H
#define HELMHOLTZ_TEST_H

#include "helmholtz.h"

/* Macros and type declarations for running simple test cases */

typedef struct{double T,p,rho,u,h,s,cv,cp,cp0,a;} TestData;

/**
    Run tests for p, u, h, s, a against values from a user-provided TestData array
    Tolerances are specified in the cdoe, in test.c.

    @return 1 if any failures occurred.
*/
int helm_run_test_cases(const HelmholtzData *d, unsigned ntd, const TestData *td);

/* a simple macro to actually do the testing */
#define ASSERT_TOL(FN,PARAM1,PARAM2,PARAM3,VAL,TOL) {\
        double cval; cval = FN(PARAM1,PARAM2,PARAM3);\
        double err; err = cval - (double)(VAL);\
        double relerrpc = (cval-(VAL))/(VAL)*100;\
        if(fabs(relerrpc)>maxerr)maxerr=fabs(relerrpc);\
        if(fabs(err)>fabs(TOL)){\
            fprintf(stderr,"ERROR in line %d: value of '%s(%f,%f,%s)' = %f,"\
                " should be %f, error is %f (%.2f%%)!\n"\
                , __LINE__, #FN,PARAM1,PARAM2,#PARAM3, cval, VAL,cval-(VAL)\
                ,relerrpc\
            );\
            exit(1);\
        }else{\
            fprintf(stderr,"    OK, %s(%f,%f,%s) = %8.2e with %.6f%% err.\n"\
                ,#FN,PARAM1,PARAM2,#PARAM3,VAL,relerrpc\
            );\
        }\
    }

#endif
1	#ifndef HELMHOLTZ_TEST_H
2	#define HELMHOLTZ_TEST_H
3
4	#include "helmholtz.h"
5
6	/* Macros and type declarations for running simple test cases */
7
8	typedef struct{double T,p,rho,u,h,s,cv,cp,cp0,a;} TestData;
9
10	/**
11	Run tests for p, u, h, s, a against values from a user-provided TestData array
12	Tolerances are specified in the cdoe, in test.c.
13
14	@return 1 if any failures occurred.
15	*/
16	int helm_run_test_cases(const HelmholtzData d, unsigned ntd, const TestData td);
17
18	/* a simple macro to actually do the testing */
19	#define ASSERT_TOL(FN,PARAM1,PARAM2,PARAM3,VAL,TOL) {\
20	double cval; cval = FN(PARAM1,PARAM2,PARAM3);\
21	double err; err = cval - (double)(VAL);\
22	double relerrpc = (cval-(VAL))/(VAL)*100;\
23	if(fabs(relerrpc)>maxerr)maxerr=fabs(relerrpc);\
24	if(fabs(err)>fabs(TOL)){\
25	fprintf(stderr,"ERROR in line %d: value of '%s(%f,%f,%s)' = %f,"\
26	" should be %f, error is %f (%.2f%%)!\n"\
27	, __LINE__, #FN,PARAM1,PARAM2,#PARAM3, cval, VAL,cval-(VAL)\
28	,relerrpc\
29	);\
30	exit(1);\
31	}else{\
32	fprintf(stderr," OK, %s(%f,%f,%s) = %8.2e with %.6f%% err.\n"\
33	,#FN,PARAM1,PARAM2,#PARAM3,VAL,relerrpc\
34	);\
35	}\
36	}
37
38	#endif