/** This software was developed by the University of Tennessee as part of the Distributed Data Analysis of Neutron Scattering Experiments (DANSE) project funded by the US National Science Foundation. If you use DANSE applications to do scientific research that leads to publication, we ask that you acknowledge the use of the software with the following sentence: "This work benefited from DANSE software developed under NSF award DMR-0520547." copyright 2008, University of Tennessee */ /** CCore2ndMomentModel * * C extension * * WARNING: THIS FILE WAS GENERATED BY WRAPPERGENERATOR.PY * DO NOT MODIFY THIS FILE, MODIFY src\sans\models\include\coresecondmoment.h * AND RE-RUN THE GENERATOR SCRIPT * */ #define NO_IMPORT_ARRAY #define PY_ARRAY_UNIQUE_SYMBOL PyArray_API_sans extern "C" { #include #include #include "structmember.h" #include #include #include #include } #include "coresecondmoment.h" #include "dispersion_visitor.hh" /// Error object for raised exceptions static PyObject * CCore2ndMomentModelError = NULL; // Class definition typedef struct { PyObject_HEAD /// Parameters PyObject * params; /// Dispersion parameters PyObject * dispersion; /// Underlying model object Core2ndMomentModel * model; /// Log for unit testing PyObject * log; } CCore2ndMomentModel; static void CCore2ndMomentModel_dealloc(CCore2ndMomentModel* self) { Py_DECREF(self->params); Py_DECREF(self->dispersion); Py_DECREF(self->log); delete self->model; self->ob_type->tp_free((PyObject*)self); } static PyObject * CCore2ndMomentModel_new(PyTypeObject *type, PyObject *args, PyObject *kwds) { CCore2ndMomentModel *self; self = (CCore2ndMomentModel *)type->tp_alloc(type, 0); return (PyObject *)self; } static int CCore2ndMomentModel_init(CCore2ndMomentModel *self, PyObject *args, PyObject *kwds) { if (self != NULL) { // Create parameters self->params = PyDict_New(); self->dispersion = PyDict_New(); self->model = new Core2ndMomentModel(); // Initialize parameter dictionary PyDict_SetItemString(self->params,"scale",Py_BuildValue("d",1.000000000000)); PyDict_SetItemString(self->params,"second_moment",Py_BuildValue("d",23.000000000000)); PyDict_SetItemString(self->params,"density_poly",Py_BuildValue("d",0.700000000000)); PyDict_SetItemString(self->params,"ads_amount",Py_BuildValue("d",1.900000000000)); PyDict_SetItemString(self->params,"background",Py_BuildValue("d",0.000000000000)); PyDict_SetItemString(self->params,"volf_cores",Py_BuildValue("d",0.140000000000)); PyDict_SetItemString(self->params,"radius_core",Py_BuildValue("d",500.000000000000)); PyDict_SetItemString(self->params,"sld_poly",Py_BuildValue("d",0.000001500000)); PyDict_SetItemString(self->params,"sld_solv",Py_BuildValue("d",0.000006300000)); // Initialize dispersion / averaging parameter dict DispersionVisitor* visitor = new DispersionVisitor(); PyObject * disp_dict; disp_dict = PyDict_New(); self->model->radius_core.dispersion->accept_as_source(visitor, self->model->radius_core.dispersion, disp_dict); PyDict_SetItemString(self->dispersion, "radius_core", disp_dict); // Create empty log self->log = PyDict_New(); } return 0; } static char name_params[] = "params"; static char def_params[] = "Parameters"; static char name_dispersion[] = "dispersion"; static char def_dispersion[] = "Dispersion parameters"; static char name_log[] = "log"; static char def_log[] = "Log"; static PyMemberDef CCore2ndMomentModel_members[] = { {name_params, T_OBJECT, offsetof(CCore2ndMomentModel, params), 0, def_params}, {name_dispersion, T_OBJECT, offsetof(CCore2ndMomentModel, dispersion), 0, def_dispersion}, {name_log, T_OBJECT, offsetof(CCore2ndMomentModel, log), 0, def_log}, {NULL} /* Sentinel */ }; /** Read double from PyObject @param p PyObject @return double */ double CCore2ndMomentModel_readDouble(PyObject *p) { if (PyFloat_Check(p)==1) { return (double)(((PyFloatObject *)(p))->ob_fval); } else if (PyInt_Check(p)==1) { return (double)(((PyIntObject *)(p))->ob_ival); } else if (PyLong_Check(p)==1) { return (double)PyLong_AsLong(p); } else { return 0.0; } } /** * Function to call to evaluate model * @param args: input numpy array q[] * @return: numpy array object */ static PyObject *evaluateOneDim(Core2ndMomentModel* model, PyArrayObject *q){ PyArrayObject *result; // Check validity of array q , q must be of dimension 1, an array of double if (q->nd != 1 || q->descr->type_num != PyArray_DOUBLE) { //const char * message= "Invalid array: q->nd=%d,type_num=%d\n",q->nd,q->descr->type_num; //PyErr_SetString(PyExc_ValueError , message); return NULL; } result = (PyArrayObject *)PyArray_FromDims(q->nd, (int *)(q->dimensions), PyArray_DOUBLE); if (result == NULL) { const char * message= "Could not create result "; PyErr_SetString(PyExc_RuntimeError , message); return NULL; } #pragma omp parallel for for (int i = 0; i < q->dimensions[0]; i++){ double q_value = *(double *)(q->data + i*q->strides[0]); double *result_value = (double *)(result->data + i*result->strides[0]); *result_value =(*model)(q_value); } return PyArray_Return(result); } /** * Function to call to evaluate model * @param args: input numpy array [x[],y[]] * @return: numpy array object */ static PyObject * evaluateTwoDimXY( Core2ndMomentModel* model, PyArrayObject *x, PyArrayObject *y) { PyArrayObject *result; int x_len, y_len, dims[1]; //check validity of input vectors if (x->nd != 1 || x->descr->type_num != PyArray_DOUBLE || y->nd != 1 || y->descr->type_num != PyArray_DOUBLE || y->dimensions[0] != x->dimensions[0]){ const char * message= "evaluateTwoDimXY expect 2 numpy arrays"; PyErr_SetString(PyExc_ValueError , message); return NULL; } if (PyArray_Check(x) && PyArray_Check(y)) { x_len = dims[0]= x->dimensions[0]; y_len = dims[0]= y->dimensions[0]; // Make a new double matrix of same dims result=(PyArrayObject *) PyArray_FromDims(1,dims,NPY_DOUBLE); if (result == NULL){ const char * message= "Could not create result "; PyErr_SetString(PyExc_RuntimeError , message); return NULL; } /* Do the calculation. */ #pragma omp parallel for for (int i=0; i< x_len; i++) { double x_value = *(double *)(x->data + i*x->strides[0]); double y_value = *(double *)(y->data + i*y->strides[0]); double *result_value = (double *)(result->data + i*result->strides[0]); *result_value = (*model)(x_value, y_value); } return PyArray_Return(result); }else{ PyErr_SetString(CCore2ndMomentModelError, "CCore2ndMomentModel.evaluateTwoDimXY couldn't run."); return NULL; } } /** * evalDistribution function evaluate a model function with input vector * @param args: input q as vector or [qx, qy] where qx, qy are vectors * */ static PyObject * evalDistribution(CCore2ndMomentModel *self, PyObject *args){ PyObject *qx, *qy; PyArrayObject * pars; int npars ,mpars; // Get parameters // Reader parameter dictionary self->model->scale = PyFloat_AsDouble( PyDict_GetItemString(self->params, "scale") ); self->model->second_moment = PyFloat_AsDouble( PyDict_GetItemString(self->params, "second_moment") ); self->model->density_poly = PyFloat_AsDouble( PyDict_GetItemString(self->params, "density_poly") ); self->model->ads_amount = PyFloat_AsDouble( PyDict_GetItemString(self->params, "ads_amount") ); self->model->background = PyFloat_AsDouble( PyDict_GetItemString(self->params, "background") ); self->model->volf_cores = PyFloat_AsDouble( PyDict_GetItemString(self->params, "volf_cores") ); self->model->radius_core = PyFloat_AsDouble( PyDict_GetItemString(self->params, "radius_core") ); self->model->sld_poly = PyFloat_AsDouble( PyDict_GetItemString(self->params, "sld_poly") ); self->model->sld_solv = PyFloat_AsDouble( PyDict_GetItemString(self->params, "sld_solv") ); // Read in dispersion parameters PyObject* disp_dict; DispersionVisitor* visitor = new DispersionVisitor(); disp_dict = PyDict_GetItemString(self->dispersion, "radius_core"); self->model->radius_core.dispersion->accept_as_destination(visitor, self->model->radius_core.dispersion, disp_dict); // Get input and determine whether we have to supply a 1D or 2D return value. if ( !PyArg_ParseTuple(args,"O",&pars) ) { PyErr_SetString(CCore2ndMomentModelError, "CCore2ndMomentModel.evalDistribution expects a q value."); return NULL; } // Check params if(PyArray_Check(pars)==1) { // Length of list should 1 or 2 npars = pars->nd; if(npars==1) { // input is a numpy array if (PyArray_Check(pars)) { return evaluateOneDim(self->model, (PyArrayObject*)pars); } }else{ PyErr_SetString(CCore2ndMomentModelError, "CCore2ndMomentModel.evalDistribution expect numpy array of one dimension."); return NULL; } }else if( PyList_Check(pars)==1) { // Length of list should be 2 for I(qx,qy) mpars = PyList_GET_SIZE(pars); if(mpars!=2) { PyErr_SetString(CCore2ndMomentModelError, "CCore2ndMomentModel.evalDistribution expects a list of dimension 2."); return NULL; } qx = PyList_GET_ITEM(pars,0); qy = PyList_GET_ITEM(pars,1); if (PyArray_Check(qx) && PyArray_Check(qy)) { return evaluateTwoDimXY(self->model, (PyArrayObject*)qx, (PyArrayObject*)qy); }else{ PyErr_SetString(CCore2ndMomentModelError, "CCore2ndMomentModel.evalDistribution expect 2 numpy arrays in list."); return NULL; } } PyErr_SetString(CCore2ndMomentModelError, "CCore2ndMomentModel.evalDistribution couln't be run."); return NULL; } /** * Function to call to evaluate model * @param args: input q or [q,phi] * @return: function value */ static PyObject * run(CCore2ndMomentModel *self, PyObject *args) { double q_value, phi_value; PyObject* pars; int npars; // Get parameters // Reader parameter dictionary self->model->scale = PyFloat_AsDouble( PyDict_GetItemString(self->params, "scale") ); self->model->second_moment = PyFloat_AsDouble( PyDict_GetItemString(self->params, "second_moment") ); self->model->density_poly = PyFloat_AsDouble( PyDict_GetItemString(self->params, "density_poly") ); self->model->ads_amount = PyFloat_AsDouble( PyDict_GetItemString(self->params, "ads_amount") ); self->model->background = PyFloat_AsDouble( PyDict_GetItemString(self->params, "background") ); self->model->volf_cores = PyFloat_AsDouble( PyDict_GetItemString(self->params, "volf_cores") ); self->model->radius_core = PyFloat_AsDouble( PyDict_GetItemString(self->params, "radius_core") ); self->model->sld_poly = PyFloat_AsDouble( PyDict_GetItemString(self->params, "sld_poly") ); self->model->sld_solv = PyFloat_AsDouble( PyDict_GetItemString(self->params, "sld_solv") ); // Read in dispersion parameters PyObject* disp_dict; DispersionVisitor* visitor = new DispersionVisitor(); disp_dict = PyDict_GetItemString(self->dispersion, "radius_core"); self->model->radius_core.dispersion->accept_as_destination(visitor, self->model->radius_core.dispersion, disp_dict); // Get input and determine whether we have to supply a 1D or 2D return value. if ( !PyArg_ParseTuple(args,"O",&pars) ) { PyErr_SetString(CCore2ndMomentModelError, "CCore2ndMomentModel.run expects a q value."); return NULL; } // Check params if( PyList_Check(pars)==1) { // Length of list should be 2 for I(q,phi) npars = PyList_GET_SIZE(pars); if(npars!=2) { PyErr_SetString(CCore2ndMomentModelError, "CCore2ndMomentModel.run expects a double or a list of dimension 2."); return NULL; } // We have a vector q, get the q and phi values at which // to evaluate I(q,phi) q_value = CCore2ndMomentModel_readDouble(PyList_GET_ITEM(pars,0)); phi_value = CCore2ndMomentModel_readDouble(PyList_GET_ITEM(pars,1)); // Skip zero if (q_value==0) { return Py_BuildValue("d",0.0); } return Py_BuildValue("d",(*(self->model)).evaluate_rphi(q_value,phi_value)); } else { // We have a scalar q, we will evaluate I(q) q_value = CCore2ndMomentModel_readDouble(pars); return Py_BuildValue("d",(*(self->model))(q_value)); } } /** * Function to call to calculate_ER * @return: effective radius value */ static PyObject * calculate_ER(CCore2ndMomentModel *self) { // Get parameters // Reader parameter dictionary self->model->scale = PyFloat_AsDouble( PyDict_GetItemString(self->params, "scale") ); self->model->second_moment = PyFloat_AsDouble( PyDict_GetItemString(self->params, "second_moment") ); self->model->density_poly = PyFloat_AsDouble( PyDict_GetItemString(self->params, "density_poly") ); self->model->ads_amount = PyFloat_AsDouble( PyDict_GetItemString(self->params, "ads_amount") ); self->model->background = PyFloat_AsDouble( PyDict_GetItemString(self->params, "background") ); self->model->volf_cores = PyFloat_AsDouble( PyDict_GetItemString(self->params, "volf_cores") ); self->model->radius_core = PyFloat_AsDouble( PyDict_GetItemString(self->params, "radius_core") ); self->model->sld_poly = PyFloat_AsDouble( PyDict_GetItemString(self->params, "sld_poly") ); self->model->sld_solv = PyFloat_AsDouble( PyDict_GetItemString(self->params, "sld_solv") ); // Read in dispersion parameters PyObject* disp_dict; DispersionVisitor* visitor = new DispersionVisitor(); disp_dict = PyDict_GetItemString(self->dispersion, "radius_core"); self->model->radius_core.dispersion->accept_as_destination(visitor, self->model->radius_core.dispersion, disp_dict); return Py_BuildValue("d",(*(self->model)).calculate_ER()); } /** * Function to call to cal the ratio shell volume/ total volume * @return: the ratio shell volume/ total volume */ static PyObject * calculate_VR(CCore2ndMomentModel *self) { // Get parameters // Reader parameter dictionary self->model->scale = PyFloat_AsDouble( PyDict_GetItemString(self->params, "scale") ); self->model->second_moment = PyFloat_AsDouble( PyDict_GetItemString(self->params, "second_moment") ); self->model->density_poly = PyFloat_AsDouble( PyDict_GetItemString(self->params, "density_poly") ); self->model->ads_amount = PyFloat_AsDouble( PyDict_GetItemString(self->params, "ads_amount") ); self->model->background = PyFloat_AsDouble( PyDict_GetItemString(self->params, "background") ); self->model->volf_cores = PyFloat_AsDouble( PyDict_GetItemString(self->params, "volf_cores") ); self->model->radius_core = PyFloat_AsDouble( PyDict_GetItemString(self->params, "radius_core") ); self->model->sld_poly = PyFloat_AsDouble( PyDict_GetItemString(self->params, "sld_poly") ); self->model->sld_solv = PyFloat_AsDouble( PyDict_GetItemString(self->params, "sld_solv") ); // Read in dispersion parameters PyObject* disp_dict; DispersionVisitor* visitor = new DispersionVisitor(); disp_dict = PyDict_GetItemString(self->dispersion, "radius_core"); self->model->radius_core.dispersion->accept_as_destination(visitor, self->model->radius_core.dispersion, disp_dict); return Py_BuildValue("d",(*(self->model)).calculate_VR()); } /** * Function to call to evaluate model in cartesian coordinates * @param args: input q or [qx, qy]] * @return: function value */ static PyObject * runXY(CCore2ndMomentModel *self, PyObject *args) { double qx_value, qy_value; PyObject* pars; int npars; // Get parameters // Reader parameter dictionary self->model->scale = PyFloat_AsDouble( PyDict_GetItemString(self->params, "scale") ); self->model->second_moment = PyFloat_AsDouble( PyDict_GetItemString(self->params, "second_moment") ); self->model->density_poly = PyFloat_AsDouble( PyDict_GetItemString(self->params, "density_poly") ); self->model->ads_amount = PyFloat_AsDouble( PyDict_GetItemString(self->params, "ads_amount") ); self->model->background = PyFloat_AsDouble( PyDict_GetItemString(self->params, "background") ); self->model->volf_cores = PyFloat_AsDouble( PyDict_GetItemString(self->params, "volf_cores") ); self->model->radius_core = PyFloat_AsDouble( PyDict_GetItemString(self->params, "radius_core") ); self->model->sld_poly = PyFloat_AsDouble( PyDict_GetItemString(self->params, "sld_poly") ); self->model->sld_solv = PyFloat_AsDouble( PyDict_GetItemString(self->params, "sld_solv") ); // Read in dispersion parameters PyObject* disp_dict; DispersionVisitor* visitor = new DispersionVisitor(); disp_dict = PyDict_GetItemString(self->dispersion, "radius_core"); self->model->radius_core.dispersion->accept_as_destination(visitor, self->model->radius_core.dispersion, disp_dict); // Get input and determine whether we have to supply a 1D or 2D return value. if ( !PyArg_ParseTuple(args,"O",&pars) ) { PyErr_SetString(CCore2ndMomentModelError, "CCore2ndMomentModel.run expects a q value."); return NULL; } // Check params if( PyList_Check(pars)==1) { // Length of list should be 2 for I(qx, qy)) npars = PyList_GET_SIZE(pars); if(npars!=2) { PyErr_SetString(CCore2ndMomentModelError, "CCore2ndMomentModel.run expects a double or a list of dimension 2."); return NULL; } // We have a vector q, get the qx and qy values at which // to evaluate I(qx,qy) qx_value = CCore2ndMomentModel_readDouble(PyList_GET_ITEM(pars,0)); qy_value = CCore2ndMomentModel_readDouble(PyList_GET_ITEM(pars,1)); return Py_BuildValue("d",(*(self->model))(qx_value,qy_value)); } else { // We have a scalar q, we will evaluate I(q) qx_value = CCore2ndMomentModel_readDouble(pars); return Py_BuildValue("d",(*(self->model))(qx_value)); } } static PyObject * reset(CCore2ndMomentModel *self, PyObject *args) { return Py_BuildValue("d",0.0); } static PyObject * set_dispersion(CCore2ndMomentModel *self, PyObject *args) { PyObject * disp; const char * par_name; if ( !PyArg_ParseTuple(args,"sO", &par_name, &disp) ) { PyErr_SetString(CCore2ndMomentModelError, "CCore2ndMomentModel.set_dispersion expects a DispersionModel object."); return NULL; } void *temp = PyCObject_AsVoidPtr(disp); DispersionModel * dispersion = static_cast(temp); // Ugliness necessary to go from python to C // TODO: refactor this if (!strcmp(par_name, "radius_core")) { self->model->radius_core.dispersion = dispersion; } else { PyErr_SetString(CCore2ndMomentModelError, "CCore2ndMomentModel.set_dispersion expects a valid parameter name."); return NULL; } DispersionVisitor* visitor = new DispersionVisitor(); PyObject * disp_dict = PyDict_New(); dispersion->accept_as_source(visitor, dispersion, disp_dict); PyDict_SetItemString(self->dispersion, par_name, disp_dict); return Py_BuildValue("i",1); } static PyMethodDef CCore2ndMomentModel_methods[] = { {"run", (PyCFunction)run , METH_VARARGS, "Evaluate the model at a given Q or Q, phi"}, {"runXY", (PyCFunction)runXY , METH_VARARGS, "Evaluate the model at a given Q or Qx, Qy"}, {"calculate_ER", (PyCFunction)calculate_ER , METH_VARARGS, "Evaluate the model at a given Q or Q, phi"}, {"calculate_VR", (PyCFunction)calculate_VR , METH_VARARGS, "Evaluate VR"}, {"evalDistribution", (PyCFunction)evalDistribution , METH_VARARGS, "Evaluate the model at a given Q or Qx, Qy vector "}, {"reset", (PyCFunction)reset , METH_VARARGS, "Reset pair correlation"}, {"set_dispersion", (PyCFunction)set_dispersion , METH_VARARGS, "Set the dispersion model for a given parameter"}, {NULL} }; static PyTypeObject CCore2ndMomentModelType = { PyObject_HEAD_INIT(NULL) 0, /*ob_size*/ "CCore2ndMomentModel", /*tp_name*/ sizeof(CCore2ndMomentModel), /*tp_basicsize*/ 0, /*tp_itemsize*/ (destructor)CCore2ndMomentModel_dealloc, /*tp_dealloc*/ 0, /*tp_print*/ 0, /*tp_getattr*/ 0, /*tp_setattr*/ 0, /*tp_compare*/ 0, /*tp_repr*/ 0, /*tp_as_number*/ 0, /*tp_as_sequence*/ 0, /*tp_as_mapping*/ 0, /*tp_hash */ 0, /*tp_call*/ 0, /*tp_str*/ 0, /*tp_getattro*/ 0, /*tp_setattro*/ 0, /*tp_as_buffer*/ Py_TPFLAGS_DEFAULT | Py_TPFLAGS_BASETYPE, /*tp_flags*/ "CCore2ndMomentModel objects", /* tp_doc */ 0, /* tp_traverse */ 0, /* tp_clear */ 0, /* tp_richcompare */ 0, /* tp_weaklistoffset */ 0, /* tp_iter */ 0, /* tp_iternext */ CCore2ndMomentModel_methods, /* tp_methods */ CCore2ndMomentModel_members, /* tp_members */ 0, /* tp_getset */ 0, /* tp_base */ 0, /* tp_dict */ 0, /* tp_descr_get */ 0, /* tp_descr_set */ 0, /* tp_dictoffset */ (initproc)CCore2ndMomentModel_init, /* tp_init */ 0, /* tp_alloc */ CCore2ndMomentModel_new, /* tp_new */ }; //static PyMethodDef module_methods[] = { // {NULL} //}; /** * Function used to add the model class to a module * @param module: module to add the class to */ void addCCore2ndMomentModel(PyObject *module) { PyObject *d; if (PyType_Ready(&CCore2ndMomentModelType) < 0) return; Py_INCREF(&CCore2ndMomentModelType); PyModule_AddObject(module, "CCore2ndMomentModel", (PyObject *)&CCore2ndMomentModelType); d = PyModule_GetDict(module); static char error_name[] = "CCore2ndMomentModel.error"; CCore2ndMomentModelError = PyErr_NewException(error_name, NULL, NULL); PyDict_SetItemString(d, "CCore2ndMomentModelError", CCore2ndMomentModelError); }