16 |
|
* 4/3/87 |
17 |
|
*/ |
18 |
|
|
19 |
< |
#include <stdio.h> |
20 |
< |
#include "fvect.h" |
19 |
> |
#include "standard.h" |
20 |
|
|
21 |
|
#define XNAME "X_" /* x function name */ |
22 |
|
#define YNAME "Y_" /* y function name */ |
23 |
|
#define ZNAME "Z_" /* z function name */ |
24 |
|
|
26 |
– |
#define PI 3.14159265358979323846 |
27 |
– |
|
28 |
– |
#define FTINY 1e-7 |
29 |
– |
|
25 |
|
#define ABS(x) ((x)>=0 ? (x) : -(x)) |
26 |
|
|
27 |
|
#define pvect(p) printf(vformat, (p)[0], (p)[1], (p)[2]) |
329 |
|
register int i,j,k; |
330 |
|
register double temp; |
331 |
|
|
332 |
< |
bcopy(mat, m4tmp, sizeof(m4tmp)); |
332 |
> |
bcopy((char *)mat, (char *)m4tmp, sizeof(m4tmp)); |
333 |
|
/* set inverse to identity */ |
334 |
|
for (i = 0; i < 4; i++) |
335 |
|
for (j = 0; j < 4; j++) |