Back to index

radiance  4R0+20100331
biggerlib.c
Go to the documentation of this file.
00001 #ifndef lint
00002 static const char    RCSid[] = "$Id: biggerlib.c,v 3.5 2003/08/04 22:37:53 greg Exp $";
00003 #endif
00004 /*
00005  *  biggerlib.c - functions for an even bigger library.
00006  *
00007  *   10/2/86
00008  */
00009 
00010 #include <stdio.h>
00011 #include <math.h>
00012 
00013 #include "calcomp.h"
00014 
00015 double  argument(int);
00016 static double  l_j0(char *), l_j1(char *), l_jn(char *);
00017 static double  l_y0(char *), l_y1(char *), l_yn(char *);
00018 static double  l_erf(char *), l_erfc(char *);
00019 
00020 
00021 void
00022 biggerlib()                 /* expand the library */
00023 {
00024                             /* the Bessel functions */
00025     funset("j0", 1, ':', l_j0);
00026     funset("j1", 1, ':', l_j1);
00027     funset("jn", 2, ':', l_jn);
00028     funset("y0", 1, ':', l_y0);
00029     funset("y1", 1, ':', l_y1);
00030     funset("yn", 2, ':', l_yn);
00031     funset("erf", 1, ':', l_erf);
00032     funset("erfc", 1, ':', l_erfc);
00033 }
00034 
00035 
00036 static double
00037 l_j0(char *nm)
00038 {
00039     return(j0(argument(1)));
00040 }
00041 
00042 
00043 static double
00044 l_j1(char *nm)
00045 {
00046     return(j1(argument(1)));
00047 }
00048 
00049 
00050 static double
00051 l_jn(char *nm)
00052 {
00053     return(jn((int)(argument(1)+.5), argument(2)));
00054 }
00055 
00056 
00057 static double
00058 l_y0(char *nm)
00059 {
00060     return(y0(argument(1)));
00061 }
00062 
00063 
00064 static double
00065 l_y1(char *nm)
00066 {
00067     return(y1(argument(1)));
00068 }
00069 
00070 
00071 static double
00072 l_yn(char *nm)
00073 {
00074     return(yn((int)(argument(1)+.5), argument(2)));
00075 }
00076 
00077 
00078 static double
00079 l_erf(char *nm)
00080 {
00081     extern double  erf();
00082 
00083     return(erf(argument(1)));
00084 }
00085 
00086 
00087 static double
00088 l_erfc(char *nm)
00089 {
00090     extern double  erfc();
00091 
00092     return(erfc(argument(1)));
00093 }