| 21 |
|
#include <stdio.h> |
| 22 |
|
#include <math.h> |
| 23 |
|
|
| 24 |
< |
#ifndef atof |
| 25 |
< |
extern double atof(); |
| 26 |
< |
#endif |
| 24 |
> |
#define PI 3.14159265358979323846 |
| 25 |
> |
#define DELTA 10. /* MINIMAL SUSTAINED ANGLE IN DEGREES */ |
| 26 |
|
|
| 28 |
– |
#define PI 3.141592653589793 |
| 29 |
– |
#define DELTA 5. /* MINIMAL SUSTAINED ANGLE IN DEGREES */ |
| 30 |
– |
|
| 27 |
|
double baseflat[4][3], baseblind[4][3][180]; |
| 28 |
|
double A[3],X[3]; |
| 29 |
|
char *material, *name; |
| 30 |
|
double height; |
| 31 |
|
int nslats, nsurf; |
| 32 |
|
|
| 33 |
+ |
#ifdef DCL_ATOF |
| 34 |
+ |
extern double atof(); |
| 35 |
+ |
#endif |
| 36 |
|
|
| 37 |
+ |
|
| 38 |
|
main(argc, argv) |
| 39 |
|
int argc; |
| 40 |
|
char *argv[]; |
| 41 |
|
{ |
| 42 |
– |
double fabs(); |
| 42 |
|
double width, delem, depth, rcurv = 0.0, angle; |
| 43 |
|
double beta, gamma, theta, chi; |
| 44 |
|
int i, j, k, l; |
| 55 |
|
angle = atof(argv[7]); |
| 56 |
|
if (argc == 10) |
| 57 |
|
if (!strcmp(argv[8], "-r")) |
| 58 |
< |
rcurv = atof(argv[8]); |
| 58 |
> |
rcurv = atof(argv[9]); |
| 59 |
|
else if (!strcmp(argv[8], "+r")) |
| 60 |
< |
rcurv = -atof(argv[8]); |
| 60 |
> |
rcurv = -atof(argv[9]); |
| 61 |
|
else |
| 62 |
|
goto userr; |
| 63 |
|
|
| 71 |
|
|
| 72 |
|
/* HOW MANY ELEMENTARY SURFACES SHOULD BE CALCULATED ? */ |
| 73 |
|
|
| 74 |
< |
nsurf = (theta / ((PI/180.)*DELTA)); |
| 74 |
> |
nsurf = (theta / ((PI/180.)*DELTA)) + 1; |
| 75 |
|
|
| 76 |
|
/* WHAT IS THE DEPTH OF THE ELEMENTARY SURFACES ? */ |
| 77 |
|
|
| 155 |
|
makeflat(w,d,a) |
| 156 |
|
double w, d, a; |
| 157 |
|
{ |
| 159 |
– |
double sin(), cos(); |
| 158 |
|
double h; |
| 159 |
|
|
| 160 |
|
h = d*sin(a); |
