COSC-4P82-Final-Project/lib/lilgp/app/regression/function.c

115 lines
2.4 KiB
C
Raw Normal View History

2024-04-01 00:01:49 -04:00
/* lil-gp Genetic Programming System, version 1.0, 11 July 1995
* Copyright (C) 1995 Michigan State University
*
* This program is free software; you can redistribute it and/or modify
* it under the terms of version 2 of the GNU General Public License as
* published by the Free Software Foundation.
*
* This program is distributed in the hope that it will be useful,
* but WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
* GNU General Public License for more details.
*
* You should have received a copy of the GNU General Public License
* along with this program; if not, write to the Free Software
* Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
*
* Douglas Zongker (zongker@isl.cps.msu.edu)
* Dr. Bill Punch (punch@isl.cps.msu.edu)
*
* Computer Science Department
* A-714 Wells Hall
* Michigan State University
* East Lansing, Michigan 48824
* USA
*
*/
#include <math.h>
#include <stdio.h>
#include <lilgp.h>
DATATYPE f_multiply ( int tree, farg *args )
{
return args[0].d * args[1].d;
}
#ifdef TOLERANCE_ZERO
/* fixes the DEC Alpha "high performance arithmetic trap" error.
* provided by Glen Ropella.
*/
DATATYPE f_protdivide ( int tree, farg *args )
{
if (args[1].d >= -0.00000000001 || args[1].d <= 0.00000000001)
return 1.0;
else
return args[0].d / args[1].d;
}
#else
DATATYPE f_protdivide ( int tree, farg *args )
{
if ( args[1].d == 0.0 )
return 1.0;
else
return args[0].d / args[1].d;
}
#endif
DATATYPE f_add ( int tree, farg *args )
{
return args[0].d + args[1].d;
}
DATATYPE f_subtract ( int tree, farg *args )
{
return args[0].d - args[1].d;
}
DATATYPE f_sin ( int tree, farg *args )
{
return sin(args[0].d);
}
DATATYPE f_cos ( int tree, farg *args )
{
return cos(args[0].d);
}
DATATYPE f_exp ( int tree, farg *args )
{
return exp(args[0].d);
}
DATATYPE f_rlog ( int tree, farg *args )
{
if ( args[0].d == 0.0 )
return 0.0;
else
return log ( fabs ( args[0].d ) );
}
DATATYPE f_indepvar ( int tree, farg *args )
{
globaldata *g = get_globaldata();
return g->x;
}
void f_erc_gen ( DATATYPE *r )
{
*r = (random_double(&globrand)*2.0) - 1.0;
}
char *f_erc_print ( DATATYPE d )
{
static char buffer[20];
sprintf ( buffer, "%.5f", d );
return buffer;
}