From: Ole Streicher Date: Sun, 13 Oct 2019 14:45:26 +0200 Subject: Check the values returned by the test --- ast_test.c | 26 ++++++++++++++++++++++++++ 1 file changed, 26 insertions(+) diff --git a/ast_test.c b/ast_test.c index 61e948e..5ee56d3 100644 --- a/ast_test.c +++ b/ast_test.c @@ -7,6 +7,8 @@ /* C header files. */ /* --------------- */ #include +#include +#include /* Main function. */ /* ============== */ @@ -100,8 +102,32 @@ int main( int argc, char *argv[] ) { astTran2( cvt, 10, xin, yin, 1, xout, yout ); printf( "\n"); + for ( i = 0; i < NCOORD; i++ ) { + // reference values taken from an example run on x86_64 + double x0[NCOORD] = { 0.019746, 0.191142, 0.363160, 0.536812, + 0.714043, 0.898732, 1.098905, 1.332477, + 1.644299, 2.152969 }; + double y0[NCOORD] = { 0.000018, 0.143869, 0.288499, 0.434597, + 0.582667, 0.732871, 0.884739, 1.036409, + 1.182151, 1.302710 }; + assert(fabs(xout[i] - x0[i]) < 1e-5); + assert(fabs(yout[i] - y0[i]) < 1e-5); + } + /* Perform the inverse transformation. */ astTran2( cvt, 10, xout, yout, 0, xin, yin ); + + for ( i = 0; i < NCOORD; i++ ) { + double x0 = sin(0.1*i) * cos(0.2*i); + double y0 = cos(0.1*i) * cos(0.2*i); + double z0 = sin(0.2*i); + double x1 = sin(xin[i]) * cos(yin[i]); + double y1 = cos(xin[i]) * cos(yin[i]); + double z1 = sin(yin[i]); + assert(fabs(x0-x1) < 1e-5); + assert(fabs(y0-y1) < 1e-5); + assert(fabs(z0-z1) < 1e-5); + } } /* End the AST context. */