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
|
/*
** -*- C -*-
**
** testSqrt.c
** Made by Bruno JOFRET <bruno.jofret@inria.fr>
**
** Started on Fri Dec 8 15:06:16 2006 jofret
** Last update Fri Apr 20 10:53:23 2007 jofret
**
** Copyright INRIA 2006
*/
#include "testSqrt.h"
void ssqrtsTest() {
printf(">> Float scalar\n");
printf("ssqrts(0) = %f\n", ssqrts((float) 0));
printf("ssqrts(PI) = %f\n", ssqrts(PI));
printf("ssqrts(PI/2) = %f\n", ssqrts(PI/2));
printf("ssqrts(PI/3) = %f\n", ssqrts(PI/3));
printf("ssqrts(PI/4) = %f\n", ssqrts(PI/4));
printf("ssqrts(PI/6) = %f\n", ssqrts(PI/6));
printf("ssqrts(-PI) = %f\n", ssqrts(-PI));
printf("ssqrts(-PI/2) = %f\n", ssqrts(-PI/2));
printf("ssqrts(-PI/3) = %f\n", ssqrts(-PI/3));
printf("ssqrts(-PI/4) = %f\n", ssqrts(-PI/4));
printf("ssqrts(-PI/6) = %f\n", ssqrts(-PI/6));
}
void dsqrtsTest() {
printf(">> Double scalar\n");
printf("dsqrts(0) = %e\n", dsqrts((double) 0));
printf("dsqrts(PI) = %e\n", dsqrts(PI));
printf("dsqrts(PI/2) = %e\n", dsqrts(PI/2));
printf("dsqrts(PI/3) = %e\n", dsqrts(PI/3));
printf("dsqrts(PI/4) = %e\n", dsqrts(PI/4));
printf("dsqrts(PI/6) = %e\n", dsqrts(PI/6));
printf("dsqrts(-PI) = %e\n", dsqrts(-PI));
printf("dsqrts(-PI/2) = %e\n", dsqrts(-PI/2));
printf("dsqrts(-PI/3) = %e\n", dsqrts(-PI/3));
printf("dsqrts(-PI/4) = %e\n", dsqrts(-PI/4));
printf("dsqrts(-PI/6) = %e\n", dsqrts(-PI/6));
}
int testSqrt() {
printf("\n>>>> Square Root Tests\n");
ssqrtsTest();
dsqrtsTest();
/* FIXME : Implement some test here ... */
/*
csqrtsTest();
zsqrtsTest();
ssqrtaTest();
dsqrtaTest();
csqrtaTest();
zsqrtaTest();
*/
return 0;
}
int main() {
assert(testSqrt() == 0);
return 0;
}
|
