+++ /dev/null
-/* cMarkovModelmodule.c
- * jchang
- * Created: 1/13/01
- * Last modified: 11/26/02
- *
- * This optimizes some of the functions in MarkovModel.py.
- */
-
-#include "Python.h"
-#include "csupport.h"
-
-
-/* Functions in this module. */
-
-static char cMarkovModel__logadd__doc__[] =
-"_logadd(logx, logy) -> log(x+y)\n";
-
-static PyObject *cMarkovModel__logadd(PyObject *self, PyObject *args)
-{
- PyObject *py_logx, *py_logy;
- double logx, logy, minxy;
- double sum;
-
- if(!PyArg_ParseTuple(args, "OO", &py_logx, &py_logy))
- return NULL;
- logx = PyNumber_AsDouble(py_logx);
- logy = PyNumber_AsDouble(py_logy);
- if(PyErr_Occurred())
- return NULL;
-
- if(logy-logx > 100.0) {
- Py_INCREF(py_logy);
- return py_logy;
- } else if (logx-logy > 100.0) {
- Py_INCREF(py_logx);
- return py_logx;
- }
- minxy = (logx < logy) ? logx : logy;
- sum = minxy + log(exp(logx-minxy) + exp(logy-minxy));
- return PyFloat_FromDouble(sum);
-}
-
-
-/* Module definition stuff */
-
-static PyMethodDef CMarkovModelMethods[] = {
- {"_logadd", cMarkovModel__logadd, METH_VARARGS, cMarkovModel__logadd__doc__},
- {NULL, NULL}
-};
-
-static char cMarkovModel__doc__[] =
-"This module provides optimized replacement functions for MarkovModel.\n\
-";
-
-void initcMarkovModel(void)
-{
- Py_InitModule3("cMarkovModel", CMarkovModelMethods, cMarkovModel__doc__);
-}
-
-
-