33 |
|
register RAY *r; |
34 |
|
XF *bx; |
35 |
|
{ |
36 |
< |
extern double l_arg(); |
36 |
> |
extern double l_arg(), l_erf(), l_erfc(); |
37 |
|
extern long eclock; |
38 |
|
static char *initfile = "rayinit.cal"; |
39 |
|
static long lastrno = -1; |
52 |
|
scompile("Jx=$19;Jy=$20;Jz=$21;", NULL, 0); |
53 |
|
scompile("Kx=$22;Ky=$23;Kz=$24;", NULL, 0); |
54 |
|
funset("arg", 1, '=', l_arg); |
55 |
+ |
funset("erf", 1, ':', l_erf); |
56 |
+ |
funset("erfc", 1, ':', l_erfc); |
57 |
|
setnoisefuncs(); |
58 |
|
initfile = NULL; |
59 |
|
} |
135 |
|
objerror(fobj, USER, errmsg); |
136 |
|
} |
137 |
|
return(fobj->oargs.farg[n-1]); |
138 |
+ |
} |
139 |
+ |
|
140 |
+ |
|
141 |
+ |
double |
142 |
+ |
l_erf() /* error function */ |
143 |
+ |
{ |
144 |
+ |
extern double erf(); |
145 |
+ |
|
146 |
+ |
return(erf(argument(1))); |
147 |
+ |
} |
148 |
+ |
|
149 |
+ |
|
150 |
+ |
double |
151 |
+ |
l_erfc() /* cumulative error function */ |
152 |
+ |
{ |
153 |
+ |
extern double erfc(); |
154 |
+ |
|
155 |
+ |
return(erfc(argument(1))); |
156 |
|
} |
157 |
|
|
158 |
|
|