06dad4687f
* hb_slex.vc
* Minor corrections.
* include/hbapi.h
+ extern PHB_ITEM hb_paramError( int iParam ); /* Returns either the generic parameter or a NIL item if param not provided */
+ extern PHB_ITEM hb_arrayFromStack( USHORT uiLen ); /* Creates and returns an Array of n Elements from the Eval Stack - Does NOT pop the items. */
+ extern PHB_ITEM hb_arrayFromParams( void ); /* Creates and returns an Array of current Generic Parameters. */
* include/hbapierr.h
* extern USHORT hb_errRT_BASE ( ULONG ulGenCode, ULONG ulSubCode, char * szDescription, char * szOperation, USHORT uiArgCount, ... );
* extern USHORT hb_errRT_BASE_Ext1 ( ULONG ulGenCode, ULONG ulSubCode, char * szDescription, char * szOperation, USHORT uiOsCode, USHORT uiFlags, USHORT uiArgCount, ... );
* extern PHB_ITEM hb_errRT_BASE_Subst ( ULONG ulGenCode, ULONG ulSubCode, char * szDescription, char * szOperation, USHORT uiArgCount, ... );
* extern void hb_errRT_BASE_SubstR ( ULONG ulGenCode, ULONG ulSubCode, char * szDescription, char * szOperation, USHORT uiArgCount, ... );
+ Added support for variable number of Error Arguments.
* source/rtl/errorapi.c
* hb_errRT_BASE ( ULONG ulGenCode, ULONG ulSubCode, char * szDescription, char * szOperation, USHORT uiArgCount, ... );
* hb_errRT_BASE_Ext1 ( ULONG ulGenCode, ULONG ulSubCode, char * szDescription, char * szOperation, USHORT uiOsCode, USHORT uiFlags, USHORT uiArgCount, ... );
* hb_errRT_BASE_Subst ( ULONG ulGenCode, ULONG ulSubCode, char * szDescription, char * szOperation, USHORT uiArgCount, ... );
* hb_errRT_BASE_SubstR ( ULONG ulGenCode, ULONG ulSubCode, char * szDescription, char * szOperation, USHORT uiArgCount, ... );
+ Added support for variable number of Error Arguments.
* source/vm/arrays.c
+ PHB_ITEM hb_arrayFromStack( USHORT uiLen )
+ PHB_ITEM hb_arrayFromParams( void )
* source/vm/extend.c
+ PHB_ITEM hb_paramError( int iParam )
* source/rdd/dbcmd.c
* Added the Error Arguments as parameters to hb_errRT_BASE*().
* source/rtl/abs.c
* source/rtl/at.c
* source/rtl/chrasc.c
* source/rtl/copyfile.c
* source/rtl/datec.c
* source/rtl/dateshb.c
* source/rtl/diskspac.c
* source/rtl/disksphb.c
* source/rtl/do.c
* source/rtl/left.c
* source/rtl/len.c
* source/rtl/math.c
* source/rtl/minmax.c
* source/rtl/mod.c
* source/rtl/philes.c
* source/rtl/philesx.c
* source/rtl/replic.c
* source/rtl/round.c
* source/rtl/set.c
* source/rtl/space.c
* source/rtl/str.c
* source/rtl/strcase.c
* source/rtl/strtran.c
* source/rtl/strzero.c
* source/rtl/substr.c
* source/rtl/transfrm.c
* source/rtl/trim.c
* source/rtl/type.c
* source/rtl/val.c
* source/rtl/word.c
* Added the Error Arguments as parameters to hb_errRT_BASE*().
* source/vm/arrayshb.c
* source/vm/eval.c
* source/vm/hvm.c
* source/vm/macro.c
* source/vm/memvars.c
* source/vm/runner.c
* Added the Error Arguments as parameters to hb_errRT_BASE*().
* source/vm/classes.c
+ Added 0 arguments (J.F. can you please decide what should be the arguments! )
/* HBTEST now reports 173 failures ( down from 376 ) many of the failures are justified.
I have tested this as much as I could, but since it invloved lots of "Cust & Pasets" more testing should be done.
It would be great if someone could review each of HBTEST Reported Failures on a one-by-one basis.
BTW, HBTEST reports 1 UN-RELEASED Memory Block (from before my changes)!
*/
185 lines
5.1 KiB
C
185 lines
5.1 KiB
C
/*
|
|
* $Id$
|
|
*/
|
|
|
|
/*
|
|
* Harbour Project source code:
|
|
* Math functions
|
|
*
|
|
* Copyright 1999 Matthew Hamilton <mhamilton@bunge.com.au>
|
|
* www - http://www.harbour-project.org
|
|
*
|
|
* This program is free software; you can redistribute it and/or modify
|
|
* it under the terms of the GNU General Public License as published by
|
|
* the Free Software Foundation; either version 2 of the License, or
|
|
* (at your option) any later version, with one exception:
|
|
*
|
|
* The exception is that if you link the Harbour Runtime Library (HRL)
|
|
* and/or the Harbour Virtual Machine (HVM) with other files to produce
|
|
* an executable, this does not by itself cause the resulting executable
|
|
* to be covered by the GNU General Public License. Your use of that
|
|
* executable is in no way restricted on account of linking the HRL
|
|
* and/or HVM code into it.
|
|
*
|
|
* 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 (or visit
|
|
* their web site at http://www.gnu.org/).
|
|
*
|
|
*/
|
|
|
|
#include <math.h>
|
|
|
|
#include "hbapi.h"
|
|
#include "hbapiitm.h"
|
|
#include "hbapierr.h"
|
|
|
|
#if defined(__WATCOMC__)
|
|
#define HB_MATH_HANDLER
|
|
#define exception _exception
|
|
#elif defined(__BORLANDC__)
|
|
#if (__BORLANDC__ == 1328) && defined(__cplusplus)
|
|
/* NOTE: There seem to be a bug in Borland C++ 5.3 C++ mode which prevents
|
|
the redefinition of matherr, because nor "_exception" neither
|
|
"exception" will work. [vszakats] */
|
|
#else
|
|
#define HB_MATH_HANDLER
|
|
#define matherr _matherr
|
|
/* NOTE: This is needed for Borland C++ 5.5 in C++/STDC mode. [vszakats] */
|
|
#if (__BORLANDC__ >= 1360)
|
|
#define exception _exception
|
|
#endif
|
|
#endif
|
|
#elif defined(__MINGW32__)
|
|
#define HB_MATH_HANDLER
|
|
#define matherr _matherr
|
|
#define exception _exception
|
|
#endif
|
|
|
|
#if defined(HB_MATH_HANDLER)
|
|
|
|
static int s_internal_math_error = 0; /* TOFIX: This is not thread safe. */
|
|
|
|
/* define harbour specific error handler for math errors
|
|
*/
|
|
int matherr( struct exception * err )
|
|
{
|
|
HB_TRACE(HB_TR_DEBUG, ("matherr(%p)", err));
|
|
|
|
switch( err->type )
|
|
{
|
|
case DOMAIN:
|
|
/* a domain error has occured, such as sqrt( -1 ) */
|
|
s_internal_math_error = EG_ARG;
|
|
break;
|
|
case SING:
|
|
/* a singularity will result, such as pow( 0, -2 ) */
|
|
s_internal_math_error = EG_ARG;
|
|
break;
|
|
case OVERFLOW:
|
|
/* an overflow will result, such as pow( 10, 100 ) */
|
|
s_internal_math_error = EG_NUMOVERFLOW;
|
|
break;
|
|
case UNDERFLOW:
|
|
/* an underflow will result, such as pow( 10, -100 ) */
|
|
s_internal_math_error = EG_NUMOVERFLOW;
|
|
break;
|
|
case TLOSS:
|
|
/* total loss of significance will result, such as exp( 1000 ) */
|
|
s_internal_math_error = EG_NUMERR;
|
|
break;
|
|
case PLOSS:
|
|
/* partial loss of significance will result, such as sin( 10e70 ) */
|
|
s_internal_math_error = EG_NUMERR;
|
|
break;
|
|
default:
|
|
s_internal_math_error = EG_NUMERR;
|
|
break;
|
|
}
|
|
|
|
err->retval = 0.0;
|
|
|
|
return 1; /* don't print any message and don't set errno */
|
|
}
|
|
#endif
|
|
|
|
HB_FUNC( EXP )
|
|
{
|
|
if( ISNUM( 1 ) )
|
|
{
|
|
#if defined(HB_MATH_HANDLER)
|
|
double dResult = exp( hb_parnd( 1 ) );
|
|
|
|
if( s_internal_math_error )
|
|
{
|
|
hb_errRT_BASE_SubstR( s_internal_math_error, 1096, NULL, "EXP", 1, hb_paramError( 1 ) );
|
|
s_internal_math_error = 0;
|
|
}
|
|
else
|
|
hb_retnd( dResult );
|
|
#else
|
|
hb_retnd( exp( hb_parnd( 1 ) ) );
|
|
#endif
|
|
}
|
|
else
|
|
hb_errRT_BASE_SubstR( EG_ARG, 1096, NULL, "EXP", 1, hb_paramError( 1 ) );
|
|
}
|
|
|
|
HB_FUNC( LOG )
|
|
{
|
|
if( ISNUM( 1 ) )
|
|
{
|
|
#if defined(HB_MATH_HANDLER)
|
|
double dResult = log( hb_parnd( 1 ) );
|
|
|
|
if( s_internal_math_error )
|
|
{
|
|
hb_errRT_BASE_SubstR( s_internal_math_error, 1095, NULL, "LOG", 1, hb_paramError( 1 ) );
|
|
s_internal_math_error = 0;
|
|
}
|
|
else
|
|
hb_retnd( dResult );
|
|
#else
|
|
double dNumber = hb_parnd( 1 );
|
|
|
|
if( dNumber <= 0.0 )
|
|
/* Indicate overflow if called with an invalid argument */
|
|
hb_retndlen( log( dNumber ), 99, -1 );
|
|
else
|
|
hb_retnd( log( dNumber ) );
|
|
#endif
|
|
}
|
|
else
|
|
hb_errRT_BASE_SubstR( EG_ARG, 1095, NULL, "LOG", 1, hb_paramError( 1 ) );
|
|
}
|
|
|
|
HB_FUNC( SQRT )
|
|
{
|
|
if( ISNUM( 1 ) )
|
|
{
|
|
#if defined(HB_MATH_HANDLER)
|
|
double dResult = sqrt( hb_parnd( 1 ) );
|
|
|
|
if( s_internal_math_error )
|
|
{
|
|
hb_errRT_BASE_SubstR( s_internal_math_error, 1097, NULL, "SQRT", 1, hb_paramError( 1 ) );
|
|
s_internal_math_error = 0;
|
|
}
|
|
else
|
|
hb_retnd( dResult );
|
|
#else
|
|
double dNumber = hb_parnd( 1 );
|
|
|
|
hb_retnd( dNumber > 0 ? sqrt( dNumber ) : 0 ); /* Clipper doesn't error! */
|
|
#endif
|
|
}
|
|
else
|
|
hb_errRT_BASE_SubstR( EG_ARG, 1097, NULL, "SQRT", 1, hb_paramError( 1 ) );
|
|
}
|
|
|