/* Copyright 2000 by Jeffrey Chang.  All rights reserved.
 * This code is part of the Biopython distribution and governed by its
 * license.  Please see the LICENSE file that should have been included
 * as part of this package.
 *
 * cmathfnsmodule.c
 * Created 3 Jun 2000
 */

#include "Python.h"
#include <math.h>

#include "csupport.h"



/************************************** Exported Functions ***********/

static char cmathfns_intd__doc__[] = 
"intd(x[, digits_after_decimal]) -> int x, rounded\n\
\n\
Represent a floating point number with some digits after the\n\
decimal point as an integer.  This is useful when floating point\n\
comparisons are failing due to precision problems.  e.g.\n\
intd(5.35, 1) -> 54.\n\
\n\
";

static PyObject *cmathfns_intd(
    PyObject *self, PyObject *args, PyObject *keywds)
{
    PyObject *digits_after_decimal = Py_None;
    double x, digits;
    double precision;

    static char *kwlist[] = {"x", "digits_after_decimal", NULL};
    if(!PyArg_ParseTupleAndKeywords(args, keywds, "d|O", kwlist, 
				    &x, &digits_after_decimal))
	return NULL;

    if(digits_after_decimal == Py_None)
	digits = 0;
    else {
	digits = PyNumber_AsDouble(digits_after_decimal);
	if(PyErr_Occurred()) {
	    return NULL;
	}
    }
    precision = pow(10, digits);
    if(x >= 0)
	x = (int)(x * precision + 0.5);
    else
	x = (int)(x * precision - 0.5);
    return PyFloat_FromDouble(x);
}




static char cmathfns_fcmp__doc__[] = 
"fcmp(x, y, precision) -> -1, 0, or 1";

static PyObject *cmathfns_fcmp(
    PyObject *self, PyObject *args, PyObject *keywds)
{
    double x, y, precision;
    int result;

    static char *kwlist[] = {"x", "y", "precision", NULL};
    if(!PyArg_ParseTupleAndKeywords(args, keywds, "ddd", kwlist, 
				    &x, &y, &precision))
	return NULL;

    if(fabs(x-y) < precision)
	result = 0;
    else if(x < y)
	result = -1;
    else result = 1;
    return PyInt_FromLong(result);
}



static char cmathfns_safe_log__doc__[] = 
"safe_log(n, zero=None, neg=None) -> log(n)\n\
\n\
Calculate the log of n.  If n is 0, returns the value of zero.  If n is\n\
negative, returns the value of neg.\n\
\n\
";

static PyObject *cmathfns_safe_log(
    PyObject *self, PyObject *args, PyObject *keywds)
{
    PyObject *zero = Py_None,
	*neg = Py_None;
    double n;

    static char *kwlist[] = {"n", "zero", "neg", NULL};

    if(!PyArg_ParseTupleAndKeywords(args, keywds, "d|OO", kwlist, 
				    &n, &zero, &neg))
	return NULL;
    
    if(n < 0) {
	Py_INCREF(neg);
	return neg;
    } else if(n < 1E-100) {
	Py_INCREF(zero);
	return zero;
    }

    return PyFloat_FromDouble(log(n));
}




/************************************** Module definition stuff ******/

static PyMethodDef cmathfnsMethods[] = {
    {"fcmp", (PyCFunction)cmathfns_fcmp, METH_VARARGS|METH_KEYWORDS, 
     cmathfns_fcmp__doc__},
    {"intd", (PyCFunction)cmathfns_intd, METH_VARARGS|METH_KEYWORDS, 
     cmathfns_intd__doc__},
    {"safe_log", (PyCFunction)cmathfns_safe_log, METH_VARARGS|METH_KEYWORDS, 
     cmathfns_safe_log__doc__},
    {NULL, NULL}
};

static char cmathfns__doc__[] =
"This provides helper functions for the mathfns module.\n\
You should never import this module on its own.\n\
\n\
";

void initcmathfns(void)
{
    (void) Py_InitModule3("cmathfns", cmathfnsMethods, cmathfns__doc__);
}