1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
|
#define _GNU_SOURCE 1
#ifdef HAVE_CONFIG_H
#include "config.h"
#endif
#include <stdio.h>
#include <stdlib.h>
#include <volume_io.h>
VIO_Real tolerance = 1e-8;
static int is_equal_real( VIO_Real e, VIO_Real a )
{
return fabs(e-a) < tolerance;
}
/* Args: expected, actual.
*/
static void assert_equal_point( VIO_Real ex, VIO_Real ey, VIO_Real ez,
VIO_Real ax, VIO_Real ay, VIO_Real az,
const char* msg )
{
if ( is_equal_real(ex,ax) &&
is_equal_real(ey,ay) &&
is_equal_real(ez,az) )
return;
printf( "%s failure.\n"
"Expected: %f %f %f\n"
" Actual: %f %f %f\n",
msg, ex,ey,ez, ax,ay,az );
exit(3);
}
int main( int ac, char* av[] )
{
int N;
VIO_General_transform xfm;
if ( ac != 3 && ac != 4 ) {
fprintf( stderr, "usage: %s N transform.xfm [tolerance]\n", av[0] );
return 1;
}
N = atoi( av[1] );
if ( input_transform_file( av[2], &xfm ) != VIO_OK ) {
fprintf( stderr, "Failed to load transform '%s'\n", av[2] );
return 2;
}
if ( ac == 4 ) {
tolerance = atof( av[3] );
printf( "Setting tolerance to %f.\n", tolerance );
}
while (N-- > 0) {
VIO_Real x = 500.0 * ( drand48() - 0.5 );
VIO_Real y = 500.0 * ( drand48() - 0.5 );
VIO_Real z = 500.0 * ( drand48() - 0.5 );
VIO_Real tx,ty,tz;
VIO_Real a,b,c;
general_transform_point( &xfm, x,y,z, &tx,&ty,&tz );
/* Check that general_inverse_transform_point() and
invert_general_transform() behave sensibly.
*/
general_inverse_transform_point( &xfm, tx,ty,tz, &a,&b,&c );
assert_equal_point( x,y,z, a,b,c,
"general_inverse_transform_point()" );
invert_general_transform( &xfm );
general_transform_point( &xfm, tx,ty,tz, &a,&b,&c );
assert_equal_point( x,y,z, a,b,c,
"general_transform_point() / inverted xfm" );
general_inverse_transform_point( &xfm, x,y,z, &a,&b,&c );
assert_equal_point( tx,ty,tz, a,b,c,
"general_inverse_transform_point() / inverted xfm" );
}
return 0;
}
|
