| 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++) |
