src/common/biggerlib.c

#ifndef lint
static const char       RCSid[] = "$Id: biggerlib.c,v 3.5 2003/08/04 22:37:53 greg Exp $";
#endif
/*
 *  biggerlib.c - functions for an even bigger library.
 *
 *   10/2/86
 */

#include <stdio.h>
#include <math.h>

#include "calcomp.h"

double  argument(int);
static double  l_j0(char *), l_j1(char *), l_jn(char *);
static double  l_y0(char *), l_y1(char *), l_yn(char *);
static double  l_erf(char *), l_erfc(char *);


void
biggerlib()                     /* expand the library */
{
                                /* the Bessel functions */
    funset("j0", 1, ':', l_j0);
    funset("j1", 1, ':', l_j1);
    funset("jn", 2, ':', l_jn);
    funset("y0", 1, ':', l_y0);
    funset("y1", 1, ':', l_y1);
    funset("yn", 2, ':', l_yn);
    funset("erf", 1, ':', l_erf);
    funset("erfc", 1, ':', l_erfc);
}


static double
l_j0(char *nm)
{
    return(j0(argument(1)));
}


static double
l_j1(char *nm)
{
    return(j1(argument(1)));
}


static double
l_jn(char *nm)
{
    return(jn((int)(argument(1)+.5), argument(2)));
}


static double
l_y0(char *nm)
{
    return(y0(argument(1)));
}


static double
l_y1(char *nm)
{
    return(y1(argument(1)));
}


static double
l_yn(char *nm)
{
    return(yn((int)(argument(1)+.5), argument(2)));
}


static double
l_erf(char *nm)
{
    extern double  erf(double x);

    return(erf(argument(1)));
}


static double
l_erfc(char *nm)
{
    extern double  erfc(double x);

    return(erfc(argument(1)));
}
Revision:	3.6
Committed:	Sat Apr 5 00:34:46 2025 UTC (6 weeks ago) by greg
Content type:	text/plain
Branch:	MAIN
CVS Tags:	HEAD
Changes since 3.5:	+3 -3 lines
Log Message:	fix: Correction to erf() and erfc() prototypes
#	User	Rev	Content
1	greg	3.1	#ifndef lint
2	greg	3.6	static const char RCSid[] = "$Id: biggerlib.c,v 3.5 2003/08/04 22:37:53 greg Exp $";
3	greg	3.1	#endif
4			/*
5			* biggerlib.c - functions for an even bigger library.
6			*
7			* 10/2/86
8			*/
9
10	greg	3.3	#include <stdio.h>
11	greg	3.1	#include <math.h>
12	schorsch	3.2
13			#include "calcomp.h"
14	greg	3.1
15	schorsch	3.4	double argument(int);
16	greg	3.5	static double l_j0(char ), l_j1(char ), l_jn(char *);
17			static double l_y0(char ), l_y1(char ), l_yn(char *);
18			static double l_erf(char ), l_erfc(char );
19	greg	3.1
20
21	schorsch	3.4	void
22	greg	3.1	biggerlib() /* expand the library */
23			{
24			/* the Bessel functions */
25			funset("j0", 1, ':', l_j0);
26			funset("j1", 1, ':', l_j1);
27			funset("jn", 2, ':', l_jn);
28			funset("y0", 1, ':', l_y0);
29			funset("y1", 1, ':', l_y1);
30			funset("yn", 2, ':', l_yn);
31			funset("erf", 1, ':', l_erf);
32			funset("erfc", 1, ':', l_erfc);
33			}
34
35
36			static double
37	greg	3.5	l_j0(char *nm)
38	greg	3.1	{
39			return(j0(argument(1)));
40			}
41
42
43			static double
44	greg	3.5	l_j1(char *nm)
45	greg	3.1	{
46			return(j1(argument(1)));
47			}
48
49
50			static double
51	greg	3.5	l_jn(char *nm)
52	greg	3.1	{
53			return(jn((int)(argument(1)+.5), argument(2)));
54			}
55
56
57			static double
58	greg	3.5	l_y0(char *nm)
59	greg	3.1	{
60			return(y0(argument(1)));
61			}
62
63
64			static double
65	greg	3.5	l_y1(char *nm)
66	greg	3.1	{
67			return(y1(argument(1)));
68			}
69
70
71			static double
72	greg	3.5	l_yn(char *nm)
73	greg	3.1	{
74			return(yn((int)(argument(1)+.5), argument(2)));
75			}
76
77
78			static double
79	greg	3.5	l_erf(char *nm)
80	greg	3.1	{
81	greg	3.6	extern double erf(double x);
82	greg	3.1
83			return(erf(argument(1)));
84			}
85
86
87			static double
88	greg	3.5	l_erfc(char *nm)
89	greg	3.1	{
90	greg	3.6	extern double erfc(double x);
91	greg	3.1
92			return(erfc(argument(1)));
93			}