1 /* cMarkovModelmodule.c
4 * Last modified: 11/26/02
6 * This optimizes some of the functions in MarkovModel.py.
13 /* Functions in this module. */
15 static char cMarkovModel__logadd__doc__[] =
16 "_logadd(logx, logy) -> log(x+y)\n";
18 static PyObject *cMarkovModel__logadd(PyObject *self, PyObject *args)
20 PyObject *py_logx, *py_logy;
21 double logx, logy, minxy;
24 if(!PyArg_ParseTuple(args, "OO", &py_logx, &py_logy))
26 logx = PyNumber_AsDouble(py_logx);
27 logy = PyNumber_AsDouble(py_logy);
31 if(logy-logx > 100.0) {
34 } else if (logx-logy > 100.0) {
38 minxy = (logx < logy) ? logx : logy;
39 sum = minxy + log(exp(logx-minxy) + exp(logy-minxy));
40 return PyFloat_FromDouble(sum);
44 /* Module definition stuff */
46 static PyMethodDef CMarkovModelMethods[] = {
47 {"_logadd", cMarkovModel__logadd, METH_VARARGS, cMarkovModel__logadd__doc__},
51 static char cMarkovModel__doc__[] =
52 "This module provides optimized replacement functions for MarkovModel.\n\
55 void initcMarkovModel(void)
57 Py_InitModule3("cMarkovModel", CMarkovModelMethods, cMarkovModel__doc__);