feat: added a new function eexp_ for the exponential function.

- it can be numerically evaluated with the evaluate statement.
- Note: exp_ is already a FORM function, so the naming follows that of the symbol ee_.

Author: cbmarini

check/features.frm

@@ -1259,6 +1259,23 @@ assert result("LN") =~ expr("
        + 9.99999327347282003e-01*ln(2.71828e+00)
 ")
 *--#] evaluate_ln : 
+*--#[ evaluate_eexp :
+#Startfloat 84b
+CFunction exp;
+Local EXP = exp(0)*eexp_(0)+exp(1)*eexp_(1)+exp(-10)*eexp_(-10)+exp(2/11)*eexp_(2/11)+exp(36124.5)*eexp_(36124.5);
+Evaluate;
+Print +s;
+.end
+#pend_if wordsize == 2
+assert succeeded?
+assert result("EXP") =~ expr("
+         + 1.199396102035385909645662e+00*exp(2/11)
+       + 4.688258114961839863324201e+15688*exp(3.61245e+04)
+       + 4.539992976248485153559152e-05*exp(-10)
+       + 1.0e+00*exp(0)
+       + 2.718281828459045235360287e+00*exp(1)
+")
+*--#] evaluate_eexp : 
 *--#[ evaluate_li2 :
 #Startfloat 18d
 CFunction li2;

doc/manual/float.tex

@@ -124,9 +124,9 @@ \chapter{Floating point}
 the floating point numbers are printed in raw internal format. 
 \end{description}
 In addition to the above commands there are the following functions that 
-can be evaluated sqrt\_, ln\_, li2\_, gamma\_, agm\_, sin\_, cos\_, tan\_,
-asin\_, acos\_, atan\_, sinh\_, cosh\_, tanh\_, asinh\_, acosh\_, atanh\_.
-For the functions atan2\_ and lin\_ there is currently no code.
+can be evaluated sqrt\_, ln\_, eexp\_, li2\_, gamma\_, agm\_, sin\_, cos\_, tan\_,
+asin\_, acos\_, atan\_, atan2\_, sinh\_, cosh\_, tanh\_, asinh\_, acosh\_, atanh\_.
+For the function lin\_ there is currently no code.
 The agm\_ function is the arithmetic geometric mean of its two input 
 values.
 

doc/manual/functions.tex

@@ -1314,6 +1314,9 @@ \section{Extra reserved names}
 \leftvitem{3cm}{ln\_}\index{ln\_}\index{function!ln\_}
 \rightvitem{13cm}{The natural logarithm.}
 
+\leftvitem{3cm}{eexp\_}\index{eexp\_}\index{function!eexp\_}
+\rightvitem{13cm}{The exponential function.}
+
 \leftvitem{3cm}{sin\_}\index{sin\_}\index{function!sin\_}
 \rightvitem{13cm}{The sine function.}
 

sources/evaluate.c

@@ -531,6 +531,10 @@ int EvaluateFun(PHEAD WORD *term, WORD level, WORD *pars)
 					else mpfr_log(auxr3,auxr1,RND);
 					mpfr_mul(auxr2,auxr2,auxr3,RND);
 				break;
+				case EXPFUNCTION:
+					mpfr_exp(auxr3,auxr1,RND);
+					mpfr_mul(auxr2,auxr2,auxr3,RND);
+				break;
 				case LI2FUNCTION: /* should be between -1 and +1 */
 					if ( nsgn == 0 ) goto getout;
 					mpfr_abs(auxr3,auxr1,RND);
@@ -645,6 +649,7 @@ int EvaluateFun(PHEAD WORD *term, WORD level, WORD *pars)
 				case LNFUNCTION:
 				case LI2FUNCTION:
 				case LINFUNCTION:
+				case EXPFUNCTION:
 				case ASINFUNCTION:
 				case ACOSFUNCTION:
 				case ATANFUNCTION:

sources/float.c

@@ -2768,6 +2768,7 @@ int CoEvaluate(UBYTE *s)
 */
 			case SQRTFUNCTION:
 			case LNFUNCTION:
+			case EXPFUNCTION:
 			case SINFUNCTION:
 			case COSFUNCTION:
 			case TANFUNCTION:

sources/ftypes.h

@@ -477,7 +477,8 @@ typedef int (*TFUN1)();
 #define MZVHALF 121
 #define AGMFUNCTION 122
 #define GAMMAFUN 123
-#define MAXBUILTINFUNCTION 123
+#define EXPFUNCTION 124
+#define MAXBUILTINFUNCTION 124
 #else
 #define MAXBUILTINFUNCTION 115
 #endif

sources/inivar.h

@@ -252,6 +252,7 @@ static struct fixedfun {
 	,{"mzvhalf_"    ,0 ,0             ,0               ,0} /* MZVHALF */
 	,{"agm_"        ,0 ,0             ,0               ,0} /* AGMFUNCTION */
 	,{"gamma_"      ,0 ,0             ,0               ,0} /* GAMMAFUN */
+	,{"eexp_"       ,0 ,0             ,0               ,0} /* EXPFUNCTION */
 #endif
 };