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