Changeset 2bb0b26 in sasview for realSpaceModeling/pointsmodelpy/pointsmodelpymodule/misc.cc

Timestamp:

Sep 27, 2007 9:16:55 AM (17 years ago)

Author:

Mathieu Doucet <doucetm@…>

Branches:

master, ESS_GUI, ESS_GUI_Docs, ESS_GUI_batch_fitting, ESS_GUI_bumps_abstraction, ESS_GUI_iss1116, ESS_GUI_iss879, ESS_GUI_iss959, ESS_GUI_opencl, ESS_GUI_ordering, ESS_GUI_sync_sascalc, costrafo411, magnetic_scatt, release-4.1.1, release-4.1.2, release-4.2.2, release_4.0.1, ticket-1009, ticket-1094-headless, ticket-1242-2d-resolution, ticket-1243, ticket-1249, ticket885, unittest-saveload

Children:

02cb31a

Parents:

d1101bb

Message:

Added underlying 2D real-space simulation

File:

• realSpaceModeling/pointsmodelpy/pointsmodelpymodule/misc.cc (modified) (3 diffs)

realSpaceModeling/pointsmodelpy/pointsmodelpymodule/misc.cc

@@ -424 +424 @@
 }
 
+// method calculateI_Qxy(poitns, Qx,Qy)
+char pypointsmodelpy_calculateI_Qvxy__name__[] = "calculateI_Qvxy";
+char pypointsmodelpy_calculateI_Qvxy__doc__[] = "calculate scattering intensity on a 2D pixel";
+
+PyObject * pypointsmodelpy_calculateI_Qvxy(PyObject *, PyObject *args)
+{
+  PyObject *pylores = 0, *pypoint3dvec = 0;
+  double qx = 0, qy = 0;
+  double I = 0;
+
+  int ok = PyArg_ParseTuple(args, "OOdd", &pylores, &pypoint3dvec, &qx,&qy);
+  if(!ok) return NULL;
+
+  void *temp = PyCObject_AsVoidPtr(pylores);
+  LORESModel * thislores = static_cast<LORESModel *>(temp);
+  void *temp2 = PyCObject_AsVoidPtr(pypoint3dvec);
+  vector<Point3D> * thisvec = static_cast<vector<Point3D> *>(temp2);
+  
+  I = thislores->CalculateIQ_2D(*thisvec, qx,qy);
+
+  return Py_BuildValue("d",I);
+}
+
 // PDBModel method calculateIQ(iq)
 char pypointsmodelpy_get_pdb_iq__name__[] = "get_pdb_iq";
@@ -466 +489 @@
 }
 
+// PDBModel method calculateIQ_2Dv(points,qx,qy)
+char pypointsmodelpy_get_pdb_Iqvxy__name__[] = "get_pdb_Iqvxy";
+char pypointsmodelpy_get_pdb_Iqvxy__doc__[] = "calculate scattering intensity by a given (qx,qy) for PDB model";
+
+PyObject * pypointsmodelpy_get_pdb_Iqvxy(PyObject *, PyObject *args)
+{
+  PyObject *pypdb = 0, *pypoint3dvec = 0;
+  double qx = 0, qy = 0;
+  double I = 0;
+
+  int ok = PyArg_ParseTuple(args, "OOdd", &pypdb, &pypoint3dvec, &qx,&qy);
+  if(!ok) return NULL;
+
+  void *temp = PyCObject_AsVoidPtr(pypdb);
+  PDBModel * thispdb = static_cast<PDBModel *>(temp);
+  void *temp2 = PyCObject_AsVoidPtr(pypoint3dvec);
+  vector<Point3D> * thisvec = static_cast<vector<Point3D> *>(temp2);
+
+  I = thispdb->CalculateIQ_2D(*thisvec,qx,qy);
+
+  return Py_BuildValue("d",I);
+}
+
 // ComplexModel method calculateIQ(iq)
 char pypointsmodelpy_get_complex_iq__name__[] = "get_complex_iq";
@@ -485 +531 @@
 
   return Py_BuildValue("i",0);
+}
+
+//LORESModel method CalculateIQ_2D(points,qx,qy) 
+char pypointsmodelpy_get_complex_Iqxy__name__[] = "get_complex_iq_2D";
+char pypointsmodelpy_get_complex_Iqxy__doc__[] = "calculate averaged scattering intensity from a single q";
+
+PyObject * pypointsmodelpy_get_complex_Iqxy(PyObject *, PyObject *args)
+{
+  PyObject *pylores = 0, *pypoint3dvec = 0;
+  double qx = 0, qy = 0;
+  int ok = PyArg_ParseTuple(args, "OOdd", &pylores, &pypoint3dvec, &qx, &qy);
+  if(!ok) return NULL;
+
+  void *temp = PyCObject_AsVoidPtr(pylores);
+  ComplexModel * thiscomplex = static_cast<ComplexModel *>(temp);
+  void *temp2 = PyCObject_AsVoidPtr(pypoint3dvec);
+  vector<Point3D> * thisvec = static_cast<vector<Point3D> *>(temp2);
+
+  double I = thiscomplex->CalculateIQ_2D(*thisvec,qx,qy);
+
+  return Py_BuildValue("d",I);
 }