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