mirror of
https://github.com/php/php-src.git
synced 2024-11-25 19:05:31 +08:00
1068 lines
24 KiB
C
1068 lines
24 KiB
C
/*
|
|
+----------------------------------------------------------------------+
|
|
| PHP version 4.0 |
|
|
+----------------------------------------------------------------------+
|
|
| Copyright (c) 1997-2001 The PHP Group |
|
|
+----------------------------------------------------------------------+
|
|
| This source file is subject to version 2.02 of the PHP license, |
|
|
| that is bundled with this package in the file LICENSE, and is |
|
|
| available at through the world-wide-web at |
|
|
| http://www.php.net/license/2_02.txt. |
|
|
| If you did not receive a copy of the PHP license and are unable to |
|
|
| obtain it through the world-wide-web, please send a note to |
|
|
| license@php.net so we can mail you a copy immediately. |
|
|
+----------------------------------------------------------------------+
|
|
| Authors: Jim Winstead (jimw@php.net) |
|
|
| Stig Sæther Bakken <ssb@guardian.no> |
|
|
| Zeev Suraski <zeev@zend.com> |
|
|
| PHP 4.0 patches by Thies C. Arntzen (thies@thieso.net) |
|
|
+----------------------------------------------------------------------+
|
|
*/
|
|
|
|
/* $Id$ */
|
|
|
|
#include "php.h"
|
|
#include "php_math.h"
|
|
|
|
#include <math.h>
|
|
#include <float.h>
|
|
|
|
#ifndef M_PI
|
|
#define M_PI 3.14159265358979323846
|
|
#endif
|
|
|
|
|
|
char *_php_math_number_format(double, int, char , char);
|
|
|
|
/* {{{ proto int abs(int number)
|
|
Return the absolute value of the number */
|
|
PHP_FUNCTION(abs)
|
|
{
|
|
zval **value;
|
|
|
|
if (ZEND_NUM_ARGS()!=1||zend_get_parameters_ex(1, &value)==FAILURE) {
|
|
WRONG_PARAM_COUNT;
|
|
}
|
|
|
|
convert_scalar_to_number_ex(value);
|
|
|
|
if (Z_TYPE_PP(value) == IS_DOUBLE) {
|
|
RETURN_DOUBLE(fabs(Z_DVAL_PP(value)));
|
|
} else if (Z_TYPE_PP(value) == IS_LONG) {
|
|
if (Z_LVAL_PP(value) == LONG_MIN) {
|
|
RETURN_DOUBLE(-(double)LONG_MIN);
|
|
} else {
|
|
RETURN_LONG(Z_LVAL_PP(value) < 0 ? -Z_LVAL_PP(value) : Z_LVAL_PP(value));
|
|
}
|
|
}
|
|
|
|
RETURN_FALSE;
|
|
}
|
|
|
|
/* }}} */
|
|
/* {{{ proto double ceil(double number)
|
|
Returns the next highest integer value of the number */
|
|
PHP_FUNCTION(ceil)
|
|
{
|
|
zval **value;
|
|
|
|
if (ZEND_NUM_ARGS()!=1||zend_get_parameters_ex(1, &value)==FAILURE) {
|
|
WRONG_PARAM_COUNT;
|
|
}
|
|
|
|
convert_scalar_to_number_ex(value);
|
|
|
|
if (Z_TYPE_PP(value) == IS_DOUBLE) {
|
|
RETURN_DOUBLE(ceil(Z_DVAL_PP(value)));
|
|
} else if (Z_TYPE_PP(value) == IS_LONG) {
|
|
RETURN_LONG(Z_LVAL_PP(value));
|
|
}
|
|
|
|
RETURN_FALSE;
|
|
}
|
|
|
|
/* }}} */
|
|
/* {{{ proto double floor(double number)
|
|
Returns the next lowest integer value from the number */
|
|
|
|
PHP_FUNCTION(floor)
|
|
{
|
|
zval **value;
|
|
|
|
if (ZEND_NUM_ARGS()!=1||zend_get_parameters_ex(1, &value)==FAILURE) {
|
|
WRONG_PARAM_COUNT;
|
|
}
|
|
|
|
convert_scalar_to_number_ex(value);
|
|
|
|
if (Z_TYPE_PP(value) == IS_DOUBLE) {
|
|
RETURN_DOUBLE(floor(Z_DVAL_PP(value)));
|
|
} else if (Z_TYPE_PP(value) == IS_LONG) {
|
|
RETURN_LONG(Z_LVAL_PP(value));
|
|
}
|
|
|
|
RETURN_FALSE;
|
|
}
|
|
|
|
/* }}} */
|
|
|
|
|
|
/* {{{ proto double round(double number [, int precision])
|
|
Returns the number rounded to specified precision. */
|
|
PHP_FUNCTION(round)
|
|
{
|
|
zval **value, **precision;
|
|
int places = 0;
|
|
double f, return_val;
|
|
|
|
if (ZEND_NUM_ARGS() < 1 || ZEND_NUM_ARGS() > 2 ||
|
|
zend_get_parameters_ex(ZEND_NUM_ARGS(), &value, &precision) == FAILURE) {
|
|
WRONG_PARAM_COUNT;
|
|
}
|
|
|
|
if (ZEND_NUM_ARGS() == 2) {
|
|
convert_to_long_ex(precision);
|
|
places = (int) Z_LVAL_PP(precision);
|
|
}
|
|
|
|
convert_scalar_to_number_ex(value);
|
|
|
|
switch (Z_TYPE_PP(value)) {
|
|
case IS_LONG:
|
|
/* Simple case - long that doesn't need to be rounded. */
|
|
if (places >= 0) {
|
|
RETURN_DOUBLE((double) Z_LVAL_PP(value));
|
|
}
|
|
/* break omitted intentionally */
|
|
|
|
case IS_DOUBLE:
|
|
return_val = (Z_TYPE_PP(value) == IS_LONG) ?
|
|
(double)Z_LVAL_PP(value) : Z_DVAL_PP(value);
|
|
|
|
f = pow(10.0, places);
|
|
|
|
return_val *= f;
|
|
if (return_val >= 0.0)
|
|
return_val = floor(return_val + 0.5);
|
|
else
|
|
return_val = ceil(return_val - 0.5);
|
|
return_val /= f;
|
|
|
|
RETURN_DOUBLE(return_val);
|
|
break;
|
|
|
|
default:
|
|
RETURN_FALSE;
|
|
break;
|
|
}
|
|
}
|
|
/* }}} */
|
|
/* {{{ proto double sin(double number)
|
|
Returns the sine of the number in radians */
|
|
|
|
PHP_FUNCTION(sin)
|
|
{
|
|
zval **num;
|
|
|
|
if (ZEND_NUM_ARGS() != 1 || zend_get_parameters_ex(1, &num) == FAILURE) {
|
|
WRONG_PARAM_COUNT;
|
|
}
|
|
convert_to_double_ex(num);
|
|
Z_DVAL_P(return_value) = sin(Z_DVAL_PP(num));
|
|
Z_TYPE_P(return_value) = IS_DOUBLE;
|
|
}
|
|
|
|
/* }}} */
|
|
/* {{{ proto double cos(double number)
|
|
Returns the cosine of the number in radians */
|
|
|
|
PHP_FUNCTION(cos)
|
|
{
|
|
zval **num;
|
|
|
|
if (ZEND_NUM_ARGS() != 1 || zend_get_parameters_ex(1, &num) == FAILURE) {
|
|
WRONG_PARAM_COUNT;
|
|
}
|
|
convert_to_double_ex(num);
|
|
Z_DVAL_P(return_value) = cos(Z_DVAL_PP(num));
|
|
Z_TYPE_P(return_value) = IS_DOUBLE;
|
|
}
|
|
/* }}} */
|
|
/* {{{ proto double tan(double number)
|
|
Returns the tangent of the number in radians */
|
|
PHP_FUNCTION(tan)
|
|
{
|
|
zval **num;
|
|
|
|
if (ZEND_NUM_ARGS() != 1 || zend_get_parameters_ex(1, &num) == FAILURE) {
|
|
WRONG_PARAM_COUNT;
|
|
}
|
|
convert_to_double_ex(num);
|
|
Z_DVAL_P(return_value) = tan(Z_DVAL_PP(num));
|
|
Z_TYPE_P(return_value) = IS_DOUBLE;
|
|
}
|
|
|
|
/* }}} */
|
|
/* {{{ proto double asin(double number)
|
|
Returns the arc sine of the number in radians */
|
|
|
|
PHP_FUNCTION(asin)
|
|
{
|
|
zval **num;
|
|
|
|
if (ZEND_NUM_ARGS() != 1 || zend_get_parameters_ex(1, &num) == FAILURE) {
|
|
WRONG_PARAM_COUNT;
|
|
}
|
|
convert_to_double_ex(num);
|
|
Z_DVAL_P(return_value) = asin(Z_DVAL_PP(num));
|
|
Z_TYPE_P(return_value) = IS_DOUBLE;
|
|
}
|
|
|
|
/* }}} */
|
|
/* {{{ proto double acos(double number)
|
|
Return the arc cosine of the number in radians */
|
|
|
|
PHP_FUNCTION(acos)
|
|
{
|
|
zval **num;
|
|
|
|
if (ZEND_NUM_ARGS() != 1 || zend_get_parameters_ex(1, &num) == FAILURE) {
|
|
WRONG_PARAM_COUNT;
|
|
}
|
|
convert_to_double_ex(num);
|
|
Z_DVAL_P(return_value) = acos(Z_DVAL_PP(num));
|
|
Z_TYPE_P(return_value) = IS_DOUBLE;
|
|
}
|
|
|
|
/* }}} */
|
|
/* {{{ proto double atan(double number)
|
|
Returns the arc tangent of the number in radians */
|
|
|
|
PHP_FUNCTION(atan)
|
|
{
|
|
zval **num;
|
|
|
|
if (ZEND_NUM_ARGS() != 1 || zend_get_parameters_ex(1, &num) == FAILURE) {
|
|
WRONG_PARAM_COUNT;
|
|
}
|
|
convert_to_double_ex(num);
|
|
Z_DVAL_P(return_value) = atan(Z_DVAL_PP(num));
|
|
Z_TYPE_P(return_value) = IS_DOUBLE;
|
|
}
|
|
|
|
/* }}} */
|
|
/* {{{ proto double atan2(double y, double x)
|
|
Returns the arc tangent of y/x, with the resulting quadrant determined by the signs of y and x */
|
|
|
|
PHP_FUNCTION(atan2)
|
|
{
|
|
zval **num1, **num2;
|
|
|
|
if (ZEND_NUM_ARGS() != 2 || zend_get_parameters_ex(2, &num1, &num2) == FAILURE) {
|
|
WRONG_PARAM_COUNT;
|
|
}
|
|
convert_to_double_ex(num1);
|
|
convert_to_double_ex(num2);
|
|
Z_DVAL_P(return_value) = atan2(Z_DVAL_PP(num1), Z_DVAL_PP(num2));
|
|
Z_TYPE_P(return_value) = IS_DOUBLE;
|
|
}
|
|
|
|
/* }}} */
|
|
/* {{{ proto double sinh(double number)
|
|
Returns the hyperbolic sine of the number,
|
|
defined as (exp(number) - exp(-number))/2 */
|
|
|
|
PHP_FUNCTION(sinh)
|
|
{
|
|
zval **num;
|
|
|
|
if (ZEND_NUM_ARGS() != 1 || zend_get_parameters_ex(1, &num) == FAILURE) {
|
|
WRONG_PARAM_COUNT;
|
|
}
|
|
convert_to_double_ex(num);
|
|
Z_DVAL_P(return_value) = sinh(Z_DVAL_PP(num));
|
|
Z_TYPE_P(return_value) = IS_DOUBLE;
|
|
}
|
|
|
|
/* }}} */
|
|
/* {{{ proto double cosh(double number)
|
|
Returns the hyperbolic cosine of the number,
|
|
defined as (exp(number) + exp(-number))/2 */
|
|
|
|
PHP_FUNCTION(cosh)
|
|
{
|
|
zval **num;
|
|
|
|
if (ZEND_NUM_ARGS() != 1 || zend_get_parameters_ex(1, &num) == FAILURE) {
|
|
WRONG_PARAM_COUNT;
|
|
}
|
|
convert_to_double_ex(num);
|
|
Z_DVAL_P(return_value) = cosh(Z_DVAL_PP(num));
|
|
Z_TYPE_P(return_value) = IS_DOUBLE;
|
|
}
|
|
/* }}} */
|
|
/* {{{ proto double tanh(double number)
|
|
Returns the hyperbolic tangent of the number,
|
|
defined as sinh(number)/cosh(number) */
|
|
PHP_FUNCTION(tanh)
|
|
{
|
|
zval **num;
|
|
|
|
if (ZEND_NUM_ARGS() != 1 || zend_get_parameters_ex(1, &num) == FAILURE) {
|
|
WRONG_PARAM_COUNT;
|
|
}
|
|
convert_to_double_ex(num);
|
|
Z_DVAL_P(return_value) = tanh(Z_DVAL_PP(num));
|
|
Z_TYPE_P(return_value) = IS_DOUBLE;
|
|
}
|
|
|
|
/* }}} */
|
|
|
|
#ifndef PHP_WIN32
|
|
/* {{{ proto double asinh(double number)
|
|
Returns the inverse hyperbolic sine of the number,
|
|
i.e. the value whose hyperbolic sine is number */
|
|
|
|
PHP_FUNCTION(asinh)
|
|
{
|
|
zval **num;
|
|
|
|
if (ZEND_NUM_ARGS() != 1 || zend_get_parameters_ex(1, &num) == FAILURE) {
|
|
WRONG_PARAM_COUNT;
|
|
}
|
|
convert_to_double_ex(num);
|
|
Z_DVAL_P(return_value) = asinh(Z_DVAL_PP(num));
|
|
Z_TYPE_P(return_value) = IS_DOUBLE;
|
|
}
|
|
|
|
/* }}} */
|
|
/* {{{ proto double acosh(double number)
|
|
Returns the inverse hyperbolic cosine of the number,
|
|
i.e. the value whose hyperbolic cosine is number */
|
|
|
|
PHP_FUNCTION(acosh)
|
|
{
|
|
zval **num;
|
|
|
|
if (ZEND_NUM_ARGS() != 1 || zend_get_parameters_ex(1, &num) == FAILURE) {
|
|
WRONG_PARAM_COUNT;
|
|
}
|
|
convert_to_double_ex(num);
|
|
Z_DVAL_P(return_value) = acosh(Z_DVAL_PP(num));
|
|
Z_TYPE_P(return_value) = IS_DOUBLE;
|
|
}
|
|
|
|
/* }}} */
|
|
/* {{{ proto double atanh(double number)
|
|
Returns the inverse hyperbolic tangent of the number,
|
|
i.e. the value whose hyperbolic tangent is number */
|
|
|
|
PHP_FUNCTION(atanh)
|
|
{
|
|
zval **num;
|
|
|
|
if (ZEND_NUM_ARGS() != 1 || zend_get_parameters_ex(1, &num) == FAILURE) {
|
|
WRONG_PARAM_COUNT;
|
|
}
|
|
convert_to_double_ex(num);
|
|
Z_DVAL_P(return_value) = atanh(Z_DVAL_PP(num));
|
|
Z_TYPE_P(return_value) = IS_DOUBLE;
|
|
}
|
|
|
|
/* }}} */
|
|
#endif
|
|
|
|
/* {{{ proto double pi(void)
|
|
Returns an approximation of pi */
|
|
|
|
PHP_FUNCTION(pi)
|
|
{
|
|
Z_DVAL_P(return_value) = M_PI;
|
|
Z_TYPE_P(return_value) = IS_DOUBLE;
|
|
}
|
|
|
|
/* }}} */
|
|
|
|
/* {{{ proto number pow(number base, number exponent)
|
|
Returns base raised to the power of exponent. Returns
|
|
integer result when possible. */
|
|
|
|
PHP_FUNCTION(pow)
|
|
{
|
|
/* FIXME: What is our policy on float-overflow? With pow, it's
|
|
* extremely easy to request results that won't fit in any double.
|
|
*/
|
|
|
|
zval **zbase, **zexp;
|
|
long lbase, lexp;
|
|
double dval;
|
|
|
|
if (ZEND_NUM_ARGS() != 2) {
|
|
WRONG_PARAM_COUNT;
|
|
}
|
|
zend_get_parameters_ex(ZEND_NUM_ARGS(), &zbase, &zexp);
|
|
convert_scalar_to_number_ex(zbase);
|
|
convert_scalar_to_number_ex(zexp);
|
|
if ((Z_TYPE_PP(zbase) != IS_LONG && Z_TYPE_PP(zbase) != IS_DOUBLE) ||
|
|
(Z_TYPE_PP(zexp ) != IS_LONG && Z_TYPE_PP(zexp ) != IS_DOUBLE)) {
|
|
php_error(E_WARNING, "Invalid argument(s) passed to pow()");
|
|
RETURN_FALSE;
|
|
}
|
|
|
|
if (Z_TYPE_PP(zexp) == IS_DOUBLE) {
|
|
/* pow(?, float), this is the ^^ case */
|
|
convert_to_double_ex(zbase);
|
|
|
|
if ( Z_DVAL_PP(zbase) <= 0.0 ) {
|
|
/* Note that with the old behaviour, php pow() returned bogus
|
|
results. Try pow(-1, 2.5) in PHP <= 4.0.6 ... */
|
|
php_error(E_WARNING, "Trying to raise a nonpositive value to a broken power");
|
|
RETURN_FALSE;
|
|
}
|
|
RETURN_DOUBLE(exp(log(Z_DVAL_PP(zbase)) * Z_DVAL_PP(zexp)));
|
|
}
|
|
|
|
/* pow(?, int), this is the ** case */
|
|
|
|
lexp = Z_LVAL_PP(zexp);
|
|
|
|
|
|
if (Z_TYPE_PP(zbase) == IS_DOUBLE) {
|
|
/* pow(float, int) */
|
|
if (lexp == 0) {
|
|
RETURN_DOUBLE(1.0);
|
|
}
|
|
if (Z_DVAL_PP(zbase) > 0.0) {
|
|
RETURN_DOUBLE(exp(log(Z_DVAL_PP(zbase)) * lexp));
|
|
} else if (Z_DVAL_PP(zbase) == 0.0) {
|
|
if (lexp < 0) {
|
|
php_error(E_WARNING,
|
|
"Division by zero: pow(0.0, [negative integer])");
|
|
RETURN_FALSE;
|
|
} else {
|
|
RETURN_DOUBLE(0.0);
|
|
}
|
|
} else { /* lbase < 0.0 */
|
|
dval = exp(log(-Z_DVAL_PP(zbase)) * (double)lexp);
|
|
RETURN_DOUBLE(lexp & 1 ? -dval : dval);
|
|
}
|
|
|
|
}
|
|
|
|
/* pow(int, int) */
|
|
if (lexp == 0) {
|
|
RETURN_LONG(1);
|
|
}
|
|
|
|
lbase = Z_LVAL_PP(zbase);
|
|
|
|
/* lexp != 0 */
|
|
switch (lbase) {
|
|
case -1:
|
|
RETURN_LONG( lexp & 1 ? -1 : 1 ); /* if lexp=odd ... */
|
|
case 0:
|
|
if (lexp < 0) {
|
|
php_error(E_WARNING,
|
|
"Division by zero: pow(0, [negative integer])");
|
|
RETURN_FALSE;
|
|
} else {
|
|
RETURN_LONG(0);
|
|
}
|
|
case 1:
|
|
RETURN_LONG(1);
|
|
default:
|
|
/* abs(lbase) > 1 */
|
|
dval = exp(log(lbase>0? (double)lbase : -(double)lbase ) *
|
|
(double) lexp);
|
|
if (lexp < 0 || dval > (double) LONG_MAX) {
|
|
/* 1/n ( abs(n) > 1 ) || overflow */
|
|
RETURN_DOUBLE(((lexp & 1) && lbase<0) ? -dval : dval);
|
|
}
|
|
|
|
Z_TYPE_P(return_value) = IS_LONG;
|
|
Z_LVAL_P(return_value) = 1;
|
|
|
|
/* loop runs at most log(log(LONG_MAX)) times, i.e. ~ 5 */
|
|
while (lexp > 0) {
|
|
if (lexp & 1) /* odd */
|
|
Z_LVAL_P(return_value) *= lbase;
|
|
lexp >>= 1;
|
|
lbase *= lbase;
|
|
}
|
|
/* return */
|
|
}
|
|
}
|
|
|
|
/* }}} */
|
|
/* {{{ proto double exp(double number)
|
|
Returns e raised to the power of the number */
|
|
|
|
PHP_FUNCTION(exp)
|
|
{
|
|
zval **num;
|
|
|
|
if (ZEND_NUM_ARGS() != 1 || zend_get_parameters_ex(1, &num) == FAILURE) {
|
|
WRONG_PARAM_COUNT;
|
|
}
|
|
convert_to_double_ex(num);
|
|
Z_DVAL_P(return_value) = exp(Z_DVAL_PP(num));
|
|
Z_TYPE_P(return_value) = IS_DOUBLE;
|
|
}
|
|
|
|
/* }}} */
|
|
/* {{{ proto double exp2(double number)
|
|
Returns 2 raised to the power of the number */
|
|
|
|
PHP_FUNCTION(exp2)
|
|
{
|
|
zval **num;
|
|
|
|
if (ZEND_NUM_ARGS() != 1 || zend_get_parameters_ex(1, &num) == FAILURE) {
|
|
WRONG_PARAM_COUNT;
|
|
}
|
|
convert_to_double_ex(num);
|
|
// libc function is broken in RH Linux 6.1, glibc 2.1.3
|
|
//Z_DVAL_P(return_value) = exp2(Z_DVAL_PP(num));
|
|
Z_DVAL_P(return_value) = pow(2.0, Z_DVAL_PP(num));
|
|
Z_TYPE_P(return_value) = IS_DOUBLE;
|
|
}
|
|
|
|
/* }}} */
|
|
/* {{{ proto double exp10(double number)
|
|
Returns 10 raised to the power of the number */
|
|
|
|
PHP_FUNCTION(exp10)
|
|
{
|
|
zval **num;
|
|
|
|
if (ZEND_NUM_ARGS() != 1 || zend_get_parameters_ex(1, &num) == FAILURE) {
|
|
WRONG_PARAM_COUNT;
|
|
}
|
|
convert_to_double_ex(num);
|
|
// libc function is broken in RH Linux 6.1, glibc 2.1.3
|
|
//Z_DVAL_P(return_value) = exp10(Z_DVAL_PP(num));
|
|
Z_DVAL_P(return_value) = pow(10.0, Z_DVAL_PP(num));
|
|
Z_TYPE_P(return_value) = IS_DOUBLE;
|
|
}
|
|
|
|
/* }}} */
|
|
|
|
#ifndef PHP_WIN32
|
|
/* {{{ proto double expm1(double number)
|
|
Returns exp(number) - 1, computed in a way that accurate even when
|
|
the value of number is close to zero */
|
|
|
|
PHP_FUNCTION(expm1)
|
|
{
|
|
zval **num;
|
|
|
|
if (ZEND_NUM_ARGS() != 1 || zend_get_parameters_ex(1, &num) == FAILURE) {
|
|
WRONG_PARAM_COUNT;
|
|
}
|
|
convert_to_double_ex(num);
|
|
Z_DVAL_P(return_value) = expm1(Z_DVAL_PP(num));
|
|
Z_TYPE_P(return_value) = IS_DOUBLE;
|
|
}
|
|
|
|
/* }}} */
|
|
/* {{{ proto double log1p(double number)
|
|
Returns log(1 + number), computed in a way that accurate even when
|
|
the value of number is close to zero */
|
|
|
|
PHP_FUNCTION(log1p)
|
|
{
|
|
zval **num;
|
|
|
|
if (ZEND_NUM_ARGS() != 1 || zend_get_parameters_ex(1, &num) == FAILURE) {
|
|
WRONG_PARAM_COUNT;
|
|
}
|
|
convert_to_double_ex(num);
|
|
Z_DVAL_P(return_value) = log1p(Z_DVAL_PP(num));
|
|
Z_TYPE_P(return_value) = IS_DOUBLE;
|
|
}
|
|
|
|
/* }}} */
|
|
|
|
#endif
|
|
/* {{{ proto double log(double number)
|
|
Returns the natural logarithm of the number */
|
|
|
|
PHP_FUNCTION(log)
|
|
{
|
|
zval **num;
|
|
|
|
if (ZEND_NUM_ARGS() != 1 || zend_get_parameters_ex(1, &num) == FAILURE) {
|
|
WRONG_PARAM_COUNT;
|
|
}
|
|
convert_to_double_ex(num);
|
|
Z_DVAL_P(return_value) = log(Z_DVAL_PP(num));
|
|
Z_TYPE_P(return_value) = IS_DOUBLE;
|
|
}
|
|
|
|
/* }}} */
|
|
/* {{{ proto double log2(double number)
|
|
Returns the base-2 logarithm of the number */
|
|
|
|
PHP_FUNCTION(log2)
|
|
{
|
|
zval **num;
|
|
|
|
if (ZEND_NUM_ARGS() != 1 || zend_get_parameters_ex(1, &num) == FAILURE) {
|
|
WRONG_PARAM_COUNT;
|
|
}
|
|
convert_to_double_ex(num);
|
|
// libc function is broken in RH Linux 6.1, glibc 2.1.3
|
|
//Z_DVAL_P(return_value) = log2(Z_DVAL_PP(num));
|
|
Z_DVAL_P(return_value) = log(Z_DVAL_PP(num))/log(2.0);
|
|
Z_TYPE_P(return_value) = IS_DOUBLE;
|
|
}
|
|
|
|
/* }}} */
|
|
/* {{{ proto double log10(double number)
|
|
Returns the base-10 logarithm of the number */
|
|
|
|
PHP_FUNCTION(log10)
|
|
{
|
|
zval **num;
|
|
|
|
if (ZEND_NUM_ARGS() != 1 || zend_get_parameters_ex(1, &num) == FAILURE) {
|
|
WRONG_PARAM_COUNT;
|
|
}
|
|
convert_to_double_ex(num);
|
|
Z_DVAL_P(return_value) = log10(Z_DVAL_PP(num));
|
|
Z_TYPE_P(return_value) = IS_DOUBLE;
|
|
}
|
|
|
|
/* }}} */
|
|
/* {{{ proto double sqrt(double number)
|
|
Returns the square root of the number */
|
|
|
|
PHP_FUNCTION(sqrt)
|
|
{
|
|
zval **num;
|
|
|
|
if (ZEND_NUM_ARGS() != 1 || zend_get_parameters_ex(1, &num) == FAILURE) {
|
|
WRONG_PARAM_COUNT;
|
|
}
|
|
convert_to_double_ex(num);
|
|
Z_DVAL_P(return_value) = sqrt(Z_DVAL_PP(num));
|
|
Z_TYPE_P(return_value) = IS_DOUBLE;
|
|
}
|
|
|
|
/* }}} */
|
|
|
|
#ifndef PHP_WIN32
|
|
/* {{{ proto double cbrt(double number)
|
|
Returns the cubic root of the number */
|
|
|
|
PHP_FUNCTION(cbrt)
|
|
{
|
|
zval **num;
|
|
|
|
if (ZEND_NUM_ARGS() != 1 || zend_get_parameters_ex(1, &num) == FAILURE) {
|
|
WRONG_PARAM_COUNT;
|
|
}
|
|
convert_to_double_ex(num);
|
|
Z_DVAL_P(return_value) = cbrt(Z_DVAL_PP(num));
|
|
Z_TYPE_P(return_value) = IS_DOUBLE;
|
|
}
|
|
|
|
/* }}} */
|
|
#endif
|
|
|
|
/* {{{ proto double hypot(double num1, double num2)
|
|
Returns sqrt( num1*num1 + num2*num2) */
|
|
|
|
PHP_FUNCTION(hypot)
|
|
{
|
|
zval **num1, **num2;
|
|
|
|
if (ZEND_NUM_ARGS() != 2 || zend_get_parameters_ex(2, &num1, &num2) == FAILURE) {
|
|
WRONG_PARAM_COUNT;
|
|
}
|
|
convert_to_double_ex(num1);
|
|
convert_to_double_ex(num2);
|
|
Z_DVAL_P(return_value) = hypot((*num1)->value.dval, (*num2)->value.dval);
|
|
Z_TYPE_P(return_value) = IS_DOUBLE;
|
|
}
|
|
|
|
/* }}} */
|
|
|
|
/* {{{ proto double deg2rad(double number)
|
|
Converts the number in degrees to the radian equivalent */
|
|
|
|
PHP_FUNCTION(deg2rad)
|
|
{
|
|
zval **deg;
|
|
|
|
if (ZEND_NUM_ARGS() != 1 || zend_get_parameters_ex(1, °) == FAILURE) {
|
|
WRONG_PARAM_COUNT;
|
|
}
|
|
convert_to_double_ex(deg);
|
|
RETVAL_DOUBLE((Z_DVAL_PP(deg) / 180.0) * M_PI);
|
|
}
|
|
|
|
/* }}} */
|
|
/* {{{ proto double rad2deg(double number)
|
|
Converts the radian number to the equivalent number in degrees */
|
|
|
|
PHP_FUNCTION(rad2deg)
|
|
{
|
|
zval **rad;
|
|
|
|
if (ZEND_NUM_ARGS() != 1 || zend_get_parameters_ex(1, &rad) == FAILURE) {
|
|
WRONG_PARAM_COUNT;
|
|
}
|
|
convert_to_double_ex(rad);
|
|
RETVAL_DOUBLE((Z_DVAL_PP(rad) / M_PI) * 180);
|
|
}
|
|
|
|
/* }}} */
|
|
/* {{{ _php_math_basetolong */
|
|
|
|
/*
|
|
* Convert a string representation of a base(2-36) number to a long.
|
|
*/
|
|
PHPAPI long
|
|
_php_math_basetolong(zval *arg, int base) {
|
|
long mult = 1, num = 0, digit;
|
|
int i;
|
|
char c, *s;
|
|
|
|
if (Z_TYPE_P(arg) != IS_STRING || base < 2 || base > 36) {
|
|
return 0;
|
|
}
|
|
|
|
s = Z_STRVAL_P(arg);
|
|
|
|
for (i = Z_STRLEN_P(arg) - 1; i >= 0; i--, mult *= base) {
|
|
c = toupper(s[i]);
|
|
if (c >= '0' && c <= '9') {
|
|
digit = (c - '0');
|
|
} else if (c >= 'A' && c <= 'Z') {
|
|
digit = (c - 'A' + 10);
|
|
} else {
|
|
continue;
|
|
}
|
|
if (digit >= base) {
|
|
continue;
|
|
}
|
|
if(!mult || digit > LONG_MAX/mult || num > LONG_MAX-mult*digit) {
|
|
php_error(E_WARNING, "base_to_long: number '%s' is too big to fit in long", s);
|
|
return LONG_MAX;
|
|
}
|
|
num += mult * digit;
|
|
}
|
|
|
|
return num;
|
|
}
|
|
|
|
/* }}} */
|
|
/* {{{ _php_math_longtobase */
|
|
|
|
/*
|
|
* Convert a long to a string containing a base(2-36) representation of
|
|
* the number.
|
|
*/
|
|
PHPAPI char *
|
|
_php_math_longtobase(zval *arg, int base)
|
|
{
|
|
static char digits[] = "0123456789abcdefghijklmnopqrstuvwxyz";
|
|
char *result, *ptr, *ret;
|
|
int len, digit;
|
|
unsigned long value;
|
|
|
|
if (Z_TYPE_P(arg) != IS_LONG || base < 2 || base > 36) {
|
|
return empty_string;
|
|
}
|
|
|
|
value = Z_LVAL_P(arg);
|
|
|
|
/* allocates space for the longest possible result with the lowest base */
|
|
len = (sizeof(Z_LVAL_P(arg)) * 8) + 1;
|
|
result = emalloc((sizeof(Z_LVAL_P(arg)) * 8) + 1);
|
|
|
|
ptr = result + len - 1;
|
|
*ptr-- = '\0';
|
|
|
|
do {
|
|
digit = value % base;
|
|
*ptr = digits[digit];
|
|
value /= base;
|
|
}
|
|
while (ptr-- > result && value);
|
|
ptr++;
|
|
ret = estrdup(ptr);
|
|
efree(result);
|
|
|
|
return ret;
|
|
}
|
|
|
|
/* }}} */
|
|
/* {{{ proto int bindec(string binary_number)
|
|
Returns the decimal equivalent of the binary number */
|
|
|
|
PHP_FUNCTION(bindec)
|
|
{
|
|
zval **arg;
|
|
long ret;
|
|
|
|
if (ZEND_NUM_ARGS() != 1 || zend_get_parameters_ex(1, &arg) == FAILURE) {
|
|
WRONG_PARAM_COUNT;
|
|
}
|
|
|
|
convert_to_string_ex(arg);
|
|
ret = _php_math_basetolong(*arg, 2);
|
|
|
|
RETVAL_LONG(ret);
|
|
}
|
|
|
|
/* }}} */
|
|
/* {{{ proto int hexdec(string hexadecimal_number)
|
|
Returns the decimal equivalent of the hexadecimal number */
|
|
|
|
PHP_FUNCTION(hexdec)
|
|
{
|
|
zval **arg;
|
|
long ret;
|
|
|
|
if (ZEND_NUM_ARGS() != 1 || zend_get_parameters_ex(1, &arg) == FAILURE) {
|
|
WRONG_PARAM_COUNT;
|
|
}
|
|
|
|
convert_to_string_ex(arg);
|
|
|
|
ret = _php_math_basetolong(*arg, 16);
|
|
RETVAL_LONG(ret);
|
|
}
|
|
|
|
/* }}} */
|
|
/* {{{ proto int octdec(string octal_number)
|
|
Returns the decimal equivalent of an octal string */
|
|
|
|
PHP_FUNCTION(octdec)
|
|
{
|
|
zval **arg;
|
|
long ret;
|
|
|
|
if (ZEND_NUM_ARGS() != 1 || zend_get_parameters_ex(1, &arg) == FAILURE) {
|
|
WRONG_PARAM_COUNT;
|
|
}
|
|
|
|
convert_to_string_ex(arg);
|
|
|
|
ret = _php_math_basetolong(*arg, 8);
|
|
RETVAL_LONG(ret);
|
|
}
|
|
|
|
/* }}} */
|
|
/* {{{ proto string decbin(int decimal_number)
|
|
Returns a string containing a binary representation of the number */
|
|
|
|
PHP_FUNCTION(decbin)
|
|
{
|
|
zval **arg;
|
|
char *result;
|
|
|
|
if (ZEND_NUM_ARGS() != 1 || zend_get_parameters_ex(1, &arg) == FAILURE) {
|
|
WRONG_PARAM_COUNT;
|
|
}
|
|
|
|
convert_to_long_ex(arg);
|
|
|
|
result = _php_math_longtobase(*arg, 2);
|
|
Z_TYPE_P(return_value) = IS_STRING;
|
|
Z_STRLEN_P(return_value) = strlen(result);
|
|
Z_STRVAL_P(return_value) = result;
|
|
}
|
|
|
|
/* }}} */
|
|
/* {{{ proto string decoct(int decimal_number)
|
|
Returns a string containing an octal representation of the given number */
|
|
|
|
PHP_FUNCTION(decoct)
|
|
{
|
|
zval **arg;
|
|
char *result;
|
|
|
|
if (ZEND_NUM_ARGS() != 1 || zend_get_parameters_ex(1, &arg) == FAILURE) {
|
|
WRONG_PARAM_COUNT;
|
|
}
|
|
|
|
convert_to_long_ex(arg);
|
|
|
|
result = _php_math_longtobase(*arg, 8);
|
|
Z_TYPE_P(return_value) = IS_STRING;
|
|
Z_STRLEN_P(return_value) = strlen(result);
|
|
Z_STRVAL_P(return_value) = result;
|
|
}
|
|
|
|
/* }}} */
|
|
/* {{{ proto string dechex(int decimal_number)
|
|
Returns a string containing a hexadecimal representation of the given number */
|
|
|
|
PHP_FUNCTION(dechex)
|
|
{
|
|
zval **arg;
|
|
char *result;
|
|
|
|
if (ZEND_NUM_ARGS() != 1 || zend_get_parameters_ex(1, &arg) == FAILURE) {
|
|
WRONG_PARAM_COUNT;
|
|
}
|
|
|
|
convert_to_long_ex(arg);
|
|
|
|
result = _php_math_longtobase(*arg, 16);
|
|
Z_TYPE_P(return_value) = IS_STRING;
|
|
Z_STRLEN_P(return_value) = strlen(result);
|
|
Z_STRVAL_P(return_value) = result;
|
|
}
|
|
|
|
/* }}} */
|
|
/* {{{ proto string base_convert(string number, int frombase, int tobase)
|
|
Converts a number in a string from any base <= 36 to any base <= 36.
|
|
*/
|
|
|
|
PHP_FUNCTION(base_convert)
|
|
{
|
|
zval **number, **frombase, **tobase, temp;
|
|
char *result;
|
|
|
|
if (ZEND_NUM_ARGS() != 3 || zend_get_parameters_ex(3, &number, &frombase, &tobase) == FAILURE) {
|
|
WRONG_PARAM_COUNT;
|
|
}
|
|
convert_to_string_ex(number);
|
|
convert_to_long_ex(frombase);
|
|
convert_to_long_ex(tobase);
|
|
if (Z_LVAL_PP(frombase) < 2 || Z_LVAL_PP(frombase) > 36) {
|
|
php_error(E_WARNING, "base_convert: invalid `from base' (%d)", Z_LVAL_PP(frombase));
|
|
RETURN_FALSE;
|
|
}
|
|
if (Z_LVAL_PP(tobase) < 2 || Z_LVAL_PP(tobase) > 36) {
|
|
php_error(E_WARNING, "base_convert: invalid `to base' (%d)", Z_LVAL_PP(tobase));
|
|
RETURN_FALSE;
|
|
}
|
|
Z_TYPE(temp) = IS_LONG;
|
|
Z_LVAL(temp) = _php_math_basetolong(*number, Z_LVAL_PP(frombase));
|
|
result = _php_math_longtobase(&temp, Z_LVAL_PP(tobase));
|
|
RETVAL_STRING(result, 0);
|
|
}
|
|
|
|
/* }}} */
|
|
/* {{{ _php_math_number_format */
|
|
|
|
char *_php_math_number_format(double d, int dec, char dec_point, char thousand_sep)
|
|
{
|
|
char *tmpbuf, *resbuf;
|
|
char *s, *t; /* source, target */
|
|
int tmplen, reslen=0;
|
|
int count=0;
|
|
int is_negative=0;
|
|
|
|
if (d<0) {
|
|
is_negative=1;
|
|
d = -d;
|
|
}
|
|
dec = MAX(0, dec);
|
|
tmpbuf = (char *) emalloc(1+DBL_MAX_10_EXP+1+dec+1);
|
|
|
|
tmplen=sprintf(tmpbuf, "%.*f", dec, d);
|
|
|
|
if (!isdigit((int)tmpbuf[0])) {
|
|
return tmpbuf;
|
|
}
|
|
|
|
if (dec) {
|
|
reslen = dec+1 + (tmplen-dec-1) + ((thousand_sep) ? (tmplen-1-dec-1)/3 : 0);
|
|
} else {
|
|
reslen = tmplen+((thousand_sep) ? (tmplen-1)/3 : 0);
|
|
}
|
|
if (is_negative) {
|
|
reslen++;
|
|
}
|
|
resbuf = (char *) emalloc(reslen+1);
|
|
|
|
s = tmpbuf+tmplen-1;
|
|
t = resbuf+reslen;
|
|
*t-- = 0;
|
|
|
|
if (dec) {
|
|
while (isdigit((int)*s)) {
|
|
*t-- = *s--;
|
|
}
|
|
*t-- = dec_point; /* copy that dot */
|
|
s--;
|
|
}
|
|
|
|
while(s>=tmpbuf) {
|
|
*t-- = *s--;
|
|
if (thousand_sep && (++count%3)==0 && s>=tmpbuf) {
|
|
*t-- = thousand_sep;
|
|
}
|
|
}
|
|
if (is_negative) {
|
|
*t-- = '-';
|
|
}
|
|
efree(tmpbuf);
|
|
return resbuf;
|
|
}
|
|
|
|
/* }}} */
|
|
/* {{{ proto string number_format(double number [, int num_decimal_places [, string dec_seperator, string thousands_seperator]])
|
|
Formats a number with grouped thousands */
|
|
|
|
PHP_FUNCTION(number_format)
|
|
{
|
|
zval **num, **dec, **t_s, **d_p;
|
|
char thousand_sep=',', dec_point='.';
|
|
|
|
switch(ZEND_NUM_ARGS()) {
|
|
case 1:
|
|
if (zend_get_parameters_ex(1, &num)==FAILURE) {
|
|
RETURN_FALSE;
|
|
}
|
|
convert_to_double_ex(num);
|
|
RETURN_STRING(_php_math_number_format(Z_DVAL_PP(num), 0, dec_point, thousand_sep), 0);
|
|
break;
|
|
case 2:
|
|
if (zend_get_parameters_ex(2, &num, &dec)==FAILURE) {
|
|
RETURN_FALSE;
|
|
}
|
|
convert_to_double_ex(num);
|
|
convert_to_long_ex(dec);
|
|
RETURN_STRING(_php_math_number_format(Z_DVAL_PP(num), Z_LVAL_PP(dec), dec_point, thousand_sep), 0);
|
|
break;
|
|
case 4:
|
|
if (zend_get_parameters_ex(4, &num, &dec, &d_p, &t_s)==FAILURE) {
|
|
RETURN_FALSE;
|
|
}
|
|
convert_to_double_ex(num);
|
|
convert_to_long_ex(dec);
|
|
convert_to_string_ex(d_p);
|
|
convert_to_string_ex(t_s);
|
|
if (Z_STRLEN_PP(d_p)==1) {
|
|
dec_point=Z_STRVAL_PP(d_p)[0];
|
|
}
|
|
if (Z_STRLEN_PP(t_s)==1) {
|
|
thousand_sep=Z_STRVAL_PP(t_s)[0];
|
|
} else if(Z_STRLEN_PP(t_s)==0) {
|
|
thousand_sep=0;
|
|
}
|
|
RETURN_STRING(_php_math_number_format(Z_DVAL_PP(num), Z_LVAL_PP(dec), dec_point, thousand_sep), 0);
|
|
break;
|
|
default:
|
|
WRONG_PARAM_COUNT;
|
|
break;
|
|
}
|
|
}
|
|
/* }}} */
|
|
|
|
/*
|
|
* Local variables:
|
|
* tab-width: 4
|
|
* c-basic-offset: 4
|
|
* End:
|
|
* vim600: sw=4 ts=4 tw=78 fdm=marker
|
|
* vim<600: sw=4 ts=4 tw=78
|
|
*/
|