map01bf@620
|
1 // pyadbmodule.c
|
map01bf@620
|
2 //
|
map01bf@620
|
3 // the internal portion of the wrapper for audio
|
map01bf@620
|
4 // see pyadb.py for the public classes
|
map01bf@620
|
5 //
|
map01bf@620
|
6 // Created by Benjamin Fields on 2009-09-04.
|
map01bf@716
|
7 // Big update for direct data insertion 2010-June (Benjamin Fields)
|
map01bf@716
|
8 // Copyleft 2009, 2010 Goldsmith University of London.
|
map01bf@621
|
9 // Distributed and licensed under GPL2. See ../../license.txt for details.
|
map01bf@620
|
10 //
|
map01bf@620
|
11 #include <fcntl.h>
|
map01bf@621
|
12 #include <string.h>
|
map01bf@620
|
13 #include "Python.h"
|
map01bf@622
|
14 #include "structmember.h"
|
map01bf@620
|
15 #include "audioDB_API.h"
|
map01bf@620
|
16 #include "numpy/arrayobject.h"
|
map01bf@620
|
17
|
map01bf@620
|
18 static void _pyadb_close(void *ptr);
|
map01bf@620
|
19
|
map01bf@620
|
20
|
map01bf@620
|
21 /* create a new database */
|
map01bf@620
|
22 /* returns a struct or NULL on failure */
|
map01bf@620
|
23 /* api call: */
|
mas01cr@671
|
24 /* adb_t *audiodb_create(const char *path, unsigned datasize, unsigned ntracks, unsigned datadim);*/
|
map01bf@620
|
25 PyObject * _pyadb_create(PyObject *self, PyObject *args)
|
map01bf@620
|
26 {
|
map01bf@620
|
27 unsigned datasize, ntracks, datadim;
|
map01bf@620
|
28 const char *path;
|
map01bf@620
|
29 int ok;
|
mas01cr@671
|
30 adb_t *new_database;
|
map01bf@620
|
31 ok = PyArg_ParseTuple(args, "sIII", &path, &datasize, &ntracks, &datadim);
|
map01bf@620
|
32 if (!ok) return 0;
|
map01bf@620
|
33 new_database = audiodb_create(path, datasize, ntracks, datadim);
|
map01bf@620
|
34 if (!new_database) return 0;
|
map01bf@620
|
35
|
map01bf@620
|
36 return PyCObject_FromVoidPtr( new_database, _pyadb_close);
|
map01bf@620
|
37 }
|
map01bf@620
|
38
|
map01bf@620
|
39 /* open an existing database */
|
map01bf@620
|
40 /* returns a struct or NULL on failure */
|
map01bf@620
|
41 /* flags expects fcntl flags concerning the opening mode */
|
map01bf@620
|
42 /* api call: */
|
mas01cr@671
|
43 // adb_t *audiodb_open(const char *path, int flags);
|
map01bf@620
|
44 PyObject * _pyadb_open(PyObject *self, PyObject *args)
|
map01bf@620
|
45 {
|
map01bf@620
|
46 const char *path;
|
map01bf@620
|
47 char mode;
|
map01bf@620
|
48 int ok;//in python layer need to translate boolean flags to byte mask
|
mas01cr@671
|
49 adb_t *fresh_database;
|
map01bf@620
|
50 ok = PyArg_ParseTuple(args, "sc", &path, &mode);
|
map01bf@620
|
51 if (!ok) return 0;
|
map01bf@620
|
52 if (mode == 'r'){
|
map01bf@620
|
53 fresh_database = audiodb_open(path, O_RDONLY);
|
map01bf@620
|
54 }else if (mode == 'w'){
|
map01bf@620
|
55 fresh_database = audiodb_open(path, O_RDWR);
|
map01bf@620
|
56 }else{
|
map01bf@620
|
57 PyErr_SetString(PyExc_ValueError,
|
map01bf@620
|
58 "mode must be either \'r\' or \'w\'. It appears to be something else.");
|
map01bf@620
|
59 return 0;
|
map01bf@620
|
60 }
|
map01bf@620
|
61 if (!fresh_database) return 0;
|
map01bf@620
|
62
|
map01bf@620
|
63 return PyCObject_FromVoidPtr( fresh_database, _pyadb_close);
|
map01bf@620
|
64 }
|
map01bf@620
|
65
|
map01bf@620
|
66 /* database status */
|
map01bf@620
|
67 /* api call: */
|
mas01cr@671
|
68 // int audiodb_status(adb_t *mydb, adb_status_ptr status);
|
map01bf@620
|
69 PyObject * _pyadb_status(PyObject *self, PyObject *args)
|
map01bf@620
|
70 {
|
mas01cr@671
|
71 adb_t *check_db;
|
mas01cr@671
|
72 adb_status_t *status;
|
map01bf@620
|
73 int flags, ok;
|
map01bf@620
|
74 PyObject * incoming = 0;
|
mas01cr@671
|
75 status = (adb_status_t *)malloc(sizeof(adb_status_t));
|
map01bf@620
|
76
|
map01bf@620
|
77 ok = PyArg_ParseTuple(args, "O", &incoming);
|
map01bf@620
|
78 if (!ok) return 0;
|
mas01cr@671
|
79 check_db = (adb_t *)PyCObject_AsVoidPtr(incoming);
|
map01bf@620
|
80
|
map01bf@620
|
81
|
map01bf@620
|
82 flags = audiodb_status(check_db, status);
|
map01bf@620
|
83 return Py_BuildValue("IIIIILL", status->numFiles,
|
map01bf@620
|
84 status->dim,
|
map01bf@620
|
85 status->dudCount,
|
map01bf@620
|
86 status->nullCount,
|
map01bf@620
|
87 status->flags,
|
map01bf@620
|
88 status->length,
|
map01bf@620
|
89 status->data_region_size);
|
map01bf@620
|
90
|
map01bf@620
|
91 }
|
map01bf@620
|
92
|
map01bf@620
|
93 /*engage l2norm in the referenced db*/
|
map01bf@620
|
94 /*api call:*/
|
mas01cr@671
|
95 //int audiodb_l2norm(adb_t *mydb);
|
map01bf@620
|
96 PyObject * _pyadb_l2norm(PyObject *self, PyObject *args)
|
map01bf@620
|
97 {
|
mas01cr@671
|
98 adb_t *current_db;
|
map01bf@620
|
99 int ok;
|
map01bf@620
|
100 PyObject * incoming = 0;
|
map01bf@620
|
101
|
map01bf@620
|
102 ok = PyArg_ParseTuple(args, "O", &incoming);
|
map01bf@620
|
103 if (!ok) return 0;
|
mas01cr@671
|
104 current_db = (adb_t *)PyCObject_AsVoidPtr(incoming);
|
map01bf@620
|
105
|
map01bf@620
|
106 ok = audiodb_l2norm(current_db);
|
map01bf@622
|
107 return PyBool_FromLong(ok-1);
|
map01bf@620
|
108
|
map01bf@620
|
109 }
|
map01bf@620
|
110
|
map01bf@620
|
111 /*engage power thresholding in the referenced db*/
|
map01bf@620
|
112 /*api call:*/
|
mas01cr@671
|
113 // int audiodb_power(adb_t *mydb);
|
map01bf@620
|
114 PyObject * _pyadb_power(PyObject *self, PyObject *args)
|
map01bf@620
|
115 {
|
mas01cr@671
|
116 adb_t *current_db;
|
map01bf@620
|
117 int ok;
|
map01bf@620
|
118 PyObject * incoming = 0;
|
map01bf@620
|
119
|
map01bf@620
|
120 ok = PyArg_ParseTuple(args, "O", &incoming);
|
map01bf@620
|
121 if (!ok) return 0;
|
mas01cr@671
|
122 current_db = (adb_t *)PyCObject_AsVoidPtr(incoming);
|
map01bf@620
|
123
|
map01bf@620
|
124 ok = audiodb_power(current_db);
|
map01bf@622
|
125 return PyBool_FromLong(ok-1);
|
map01bf@620
|
126
|
map01bf@620
|
127 }
|
map01bf@716
|
128
|
map01bf@716
|
129 /* insert feature data from a numpy array */
|
map01bf@720
|
130 /* array given should have ndarray.shape = (numVectors, numDims)*/
|
map01bf@716
|
131 /* array datatype needs to be doubles (float may work...)*/
|
map01bf@716
|
132 /* if power is given, must be 1d array of length numVectors*/
|
map01bf@716
|
133 /* if times is given, must be 1d array of length 2*numVectors like this:*/
|
map01bf@716
|
134
|
map01bf@716
|
135 /* api call: */
|
map01bf@716
|
136 // typedef struct adb_datum {
|
map01bf@716
|
137 // uint32_t nvectors;
|
map01bf@716
|
138 // uint32_t dim;
|
map01bf@716
|
139 // const char *key;
|
map01bf@716
|
140 // double *data;
|
map01bf@716
|
141 // double *power;
|
map01bf@716
|
142 // double *times;
|
map01bf@716
|
143 // } adb_datum_t;
|
map01bf@716
|
144 // int audiodb_insert_datum(adb_t *, const adb_datum_t *);
|
map01bf@716
|
145 PyObject * _pyadb_insertFromArray(PyObject *self, PyObject *args, PyObject *keywds)
|
map01bf@716
|
146 {
|
map01bf@716
|
147 adb_t *current_db;
|
map01bf@716
|
148 adb_status_t *status;
|
map01bf@716
|
149 adb_datum_t *ins;
|
map01bf@716
|
150 int ok;
|
map01bf@716
|
151 npy_intp dims[1];
|
map01bf@716
|
152 unsigned int nDims = 0;
|
map01bf@716
|
153 unsigned int nVect = 0;
|
map01bf@716
|
154 PyObject *incoming = 0;
|
map01bf@718
|
155 PyArrayObject *features = 0;
|
map01bf@718
|
156 PyArrayObject *power = 0;
|
map01bf@718
|
157 PyArrayObject *times = 0;
|
map01bf@716
|
158 const char *key = NULL;
|
map01bf@716
|
159 PyArray_Descr *descr;
|
map01bf@716
|
160 static char *kwlist[] = { "db", "features", "nDim", "nVect", "power", "key", "times" , NULL};
|
map01bf@716
|
161
|
map01bf@718
|
162 ok = PyArg_ParseTupleAndKeywords(args, keywds, "OO!II|O!sO!", kwlist, &incoming, &PyArray_Type, &features, &nDims, &nVect, &PyArray_Type, &power, &key, &PyArray_Type, ×);
|
mas01mc@743
|
163 if (!ok){
|
mas01mc@743
|
164 PyErr_SetString(PyExc_TypeError, "Failed at PyArg_ParseTupleAndKeywords");
|
mas01mc@743
|
165 return NULL;
|
mas01mc@743
|
166 }
|
map01bf@716
|
167 //check our arrays
|
map01bf@718
|
168 // if (!PyArray_Check(features)){
|
map01bf@718
|
169 // PyErr_SetString(PyExc_TypeError, "features must be a numpy array (of floats or doubles)");
|
map01bf@718
|
170 // return NULL;
|
map01bf@718
|
171 // }
|
map01bf@716
|
172 if (!PyArray_ISFLOAT(features)){
|
map01bf@716
|
173 PyErr_SetString(PyExc_TypeError, "features numpy array must contain floats or doubles");
|
map01bf@716
|
174 return NULL;
|
map01bf@716
|
175 }
|
map01bf@716
|
176 if ((PyArray_NDIM(features) != 1) || (PyArray_DIMS(features)[0] != (nDims * nVect))){
|
map01bf@716
|
177 PyErr_SetString(PyExc_TypeError, "features numpy array must be flattened before call.");
|
map01bf@716
|
178 return NULL;
|
map01bf@716
|
179 }
|
map01bf@716
|
180 descr = PyArray_DescrFromType(NPY_DOUBLE);
|
map01bf@716
|
181
|
map01bf@716
|
182 if (power){
|
map01bf@716
|
183 if (!PyArray_Check(power)){
|
map01bf@716
|
184 PyErr_SetString(PyExc_TypeError, "power, if given, must be a numpy array (of floats or doubles)");
|
map01bf@716
|
185 return NULL;
|
map01bf@716
|
186 }
|
map01bf@716
|
187 if (!PyArray_ISFLOAT(power)){
|
map01bf@716
|
188 PyErr_SetString(PyExc_TypeError, "power numpy array, if given, must contain floats or doubles");
|
map01bf@716
|
189 return NULL;
|
map01bf@716
|
190 }
|
map01bf@716
|
191 // power = (PyArrayObject *)PyCObject_AsVoidPtr(incomingPow);
|
mas01mc@743
|
192 if (PyArray_NDIM(features) != 1 || PyArray_DIMS(power)[0] != nVect){
|
map01bf@716
|
193 PyErr_SetString(PyExc_ValueError, "power, if given must be a 1d numpy array with shape = (numVectors,)");
|
map01bf@716
|
194 return NULL;
|
map01bf@716
|
195 }
|
map01bf@716
|
196 }
|
map01bf@716
|
197 if (times){
|
map01bf@716
|
198 if (!PyArray_Check(times)){
|
map01bf@716
|
199 PyErr_SetString(PyExc_TypeError, "times, if given, must be a numpy array (of floats or doubles)");
|
map01bf@716
|
200 return NULL;
|
map01bf@716
|
201 }
|
map01bf@716
|
202 if (!PyArray_ISFLOAT(times)){
|
map01bf@716
|
203 PyErr_SetString(PyExc_TypeError, "times numpy array, if given, must contain floats or doubles");
|
map01bf@716
|
204 return NULL;
|
map01bf@716
|
205 }
|
map01bf@716
|
206 // times = (PyArrayObject *)PyCObject_AsVoidPtr(incomingTime);
|
map01bf@716
|
207 if (PyArray_NDIM(times) != 1 || PyArray_DIMS(times)[0] == (nVect*2)){
|
map01bf@716
|
208 PyErr_SetString(PyExc_ValueError, "times, if given must be a 1d numpy array with shape = (numVectors,)");
|
map01bf@716
|
209 return NULL;
|
map01bf@716
|
210 }
|
map01bf@716
|
211 }
|
map01bf@716
|
212 current_db = (adb_t *)PyCObject_AsVoidPtr(incoming);
|
map01bf@716
|
213 status = (adb_status_t *)malloc(sizeof(adb_status_t));
|
map01bf@716
|
214 //verify that the data to be inserted is the correct size for the database.
|
map01bf@716
|
215
|
map01bf@716
|
216 ins = (adb_datum_t *)malloc(sizeof(adb_datum_t));
|
map01bf@719
|
217 if (PyArray_AsCArray(&features, &(ins->data), dims, 1, descr)){
|
map01bf@716
|
218 PyErr_SetString(PyExc_RuntimeError, "Trouble expressing the feature np array as a C array.");
|
map01bf@716
|
219 return NULL;
|
map01bf@716
|
220 }
|
map01bf@716
|
221
|
map01bf@718
|
222 if (power){
|
map01bf@719
|
223 if (PyArray_AsCArray(&power, &(ins->power), dims, 1, descr)){
|
map01bf@716
|
224 PyErr_SetString(PyExc_RuntimeError, "Trouble expressing the power np array as a C array.");
|
map01bf@716
|
225 return NULL;
|
map01bf@716
|
226 }
|
map01bf@719
|
227 }else{
|
map01bf@719
|
228 ins->power=NULL;
|
map01bf@716
|
229 }
|
map01bf@716
|
230
|
map01bf@718
|
231 if (times){
|
map01bf@719
|
232 if (PyArray_AsCArray(×, &(ins->times), dims, 1, descr)){
|
map01bf@716
|
233 PyErr_SetString(PyExc_RuntimeError, "Trouble expressing the times np array as a C array.");
|
map01bf@716
|
234 return NULL;
|
map01bf@716
|
235 }
|
map01bf@719
|
236 }else{
|
map01bf@719
|
237 ins->times=NULL;
|
map01bf@716
|
238 }
|
map01bf@716
|
239 ins->key = key;
|
map01bf@716
|
240 ins->nvectors = (uint32_t)nVect;
|
map01bf@716
|
241 ins->dim = (uint32_t)nDims;
|
map01bf@716
|
242 //printf("features::%s\npower::%s\nkey::%s\ntimes::%s\n", ins->features, ins->power, ins->key, ins->times);
|
map01bf@716
|
243 ok = audiodb_insert_datum(current_db, ins);//(current_db, ins);
|
map01bf@719
|
244 return PyInt_FromLong(ok);
|
map01bf@716
|
245
|
map01bf@716
|
246 }
|
map01bf@716
|
247
|
map01bf@621
|
248 /* insert feature data stored in a file */
|
map01bf@621
|
249 /* this is a bit gross, */
|
map01bf@621
|
250 /* should be replaced eventually by a numpy based feature.*/
|
map01bf@621
|
251 /* api call: */
|
mas01cr@673
|
252 // struct adb_insert {
|
map01bf@621
|
253 // const char *features;
|
map01bf@621
|
254 // const char *power;
|
map01bf@621
|
255 // const char *key;
|
map01bf@621
|
256 // const char *times;
|
map01bf@621
|
257 // };
|
mas01cr@671
|
258 // int audiodb_insert(adb_t *mydb, adb_insert_t *ins);
|
map01bf@621
|
259 PyObject * _pyadb_insertFromFile(PyObject *self, PyObject *args, PyObject *keywds)
|
map01bf@621
|
260 {
|
mas01cr@671
|
261 adb_t *current_db;
|
mas01cr@671
|
262 adb_insert_t *ins;
|
map01bf@621
|
263 int ok;
|
map01bf@621
|
264 const char *features;
|
map01bf@621
|
265 const char *power = NULL;
|
map01bf@621
|
266 const char *key = NULL;
|
map01bf@621
|
267 const char *times = NULL;
|
map01bf@621
|
268 PyObject * incoming = 0;
|
map01bf@621
|
269 static char *kwlist[] = { "db", "features", "power", "key", "times" , NULL};
|
map01bf@621
|
270
|
map01bf@621
|
271 ok = PyArg_ParseTupleAndKeywords(args, keywds, "Os|sss", kwlist, &incoming, &features, &power, &key, ×);
|
map01bf@621
|
272 if (!ok){return NULL;}
|
map01bf@621
|
273
|
mas01cr@671
|
274 current_db = (adb_t *)PyCObject_AsVoidPtr(incoming);
|
mas01cr@671
|
275 ins = (adb_insert_t *)malloc(sizeof(adb_insert_t));
|
map01bf@621
|
276 ins->features = features;
|
map01bf@621
|
277 ins->power = power;
|
map01bf@621
|
278 ins->key = key;
|
map01bf@621
|
279 ins->times = times;
|
map01bf@622
|
280 //printf("features::%s\npower::%s\nkey::%s\ntimes::%s\n", ins->features, ins->power, ins->key, ins->times);
|
map01bf@621
|
281 ok = audiodb_insert(current_db, ins);
|
map01bf@622
|
282 return PyBool_FromLong(ok-1);
|
map01bf@621
|
283
|
map01bf@621
|
284 }
|
map01bf@621
|
285
|
map01bf@622
|
286 /* base query. The nomenclature here is about a far away as pythonic as is possible.
|
map01bf@622
|
287 * This should be taken care of via the higher level python structure
|
map01bf@622
|
288 * returns a dict that should be result ordered and key = result key
|
map01bf@622
|
289 * and value is a list of tuples one per result associated with that key, of the form:
|
map01bf@622
|
290 * (dist, qpos, ipos)
|
map01bf@622
|
291 * Note as well that this is by no means the most efficient way to cast from C, simply the most direct
|
map01bf@622
|
292 * and what it lacks in effeciency it gains in python side access. It remains to be seen if this is
|
map01bf@622
|
293 * a sensible trade.
|
map01bf@622
|
294 * api call:
|
map01bf@622
|
295 * adb_query_results_t *audiodb_query_spec(adb_t *, const adb_query_spec_t *);
|
map01bf@622
|
296 ***/
|
map01bf@622
|
297 PyObject * _pyadb_queryFromKey(PyObject *self, PyObject *args, PyObject *keywds)
|
map01bf@622
|
298 {
|
mas01cr@671
|
299 adb_t *current_db;
|
map01bf@622
|
300 adb_query_spec_t *spec;
|
map01bf@622
|
301 adb_query_results_t *result;
|
map01bf@622
|
302 int ok, exhaustive, falsePositives;
|
map01bf@622
|
303 uint32_t i;
|
map01bf@622
|
304 const char *key;
|
map01bf@622
|
305 const char *accuMode = "db";
|
map01bf@622
|
306 const char *distMode = "dot";
|
map01bf@624
|
307 const char *resFmt = "dict";
|
map01bf@622
|
308 uint32_t hop = 0;
|
map01bf@622
|
309 double radius = 0;
|
map01bf@622
|
310 double absThres = 0;
|
map01bf@622
|
311 double relThres = 0;
|
map01bf@622
|
312 double durRatio = 0;
|
map01bf@622
|
313 PyObject *includeKeys = NULL;
|
map01bf@622
|
314 PyObject *excludeKeys = NULL;
|
map01bf@622
|
315 PyObject *incoming = 0;
|
map01bf@622
|
316 PyObject *outgoing = NULL;
|
map01bf@622
|
317 PyObject *thisKey = NULL;
|
map01bf@622
|
318 PyObject *currentValue = 0;
|
map01bf@622
|
319 PyObject *newBits = 0;
|
map01bf@622
|
320 static char *kwlist[] = { "db", "key",
|
map01bf@622
|
321 "seqLength",
|
map01bf@622
|
322 "seqStart",
|
map01bf@622
|
323 "exhaustive",
|
map01bf@622
|
324 "falsePositives",
|
map01bf@622
|
325 "accumulation",
|
map01bf@622
|
326 "distance",
|
map01bf@622
|
327 "npoints",//nearest neighbor points per track
|
map01bf@624
|
328 "ntracks",
|
map01bf@622
|
329 "includeKeys",
|
map01bf@622
|
330 "excludeKeys",
|
map01bf@622
|
331 "radius",
|
map01bf@622
|
332 "absThres",
|
map01bf@622
|
333 "relThres",
|
map01bf@622
|
334 "durRatio",
|
map01bf@624
|
335 "hopSize",
|
map01bf@718
|
336 "resFmt",
|
map01bf@718
|
337 NULL
|
map01bf@622
|
338 };
|
map01bf@622
|
339 spec = (adb_query_spec_t *)malloc(sizeof(adb_query_spec_t));
|
map01bf@622
|
340 spec->qid.datum = (adb_datum_t *)malloc(sizeof(adb_datum_t));
|
map01bf@622
|
341 result = (adb_query_results_t *)malloc(sizeof(adb_query_results_t));
|
map01bf@622
|
342
|
map01bf@622
|
343 spec->qid.sequence_length = 16;
|
map01bf@622
|
344 spec->qid.sequence_start = 0;
|
map01bf@622
|
345 spec->qid.flags = 0;
|
map01bf@622
|
346 spec->params.npoints = 1;
|
map01bf@622
|
347 spec->params.ntracks = 100;//number of results returned in db mode
|
map01bf@622
|
348 spec->refine.flags = 0;
|
map01bf@622
|
349
|
map01bf@624
|
350 ok = PyArg_ParseTupleAndKeywords(args, keywds, "Os|iiiissIIOOddddIs", kwlist,
|
map01bf@622
|
351 &incoming, &key,
|
map01bf@622
|
352 &spec->qid.sequence_length,
|
map01bf@622
|
353 &spec->qid.sequence_start,
|
map01bf@622
|
354 &exhaustive, &falsePositives,
|
map01bf@622
|
355 &accuMode,&distMode,
|
map01bf@622
|
356 &spec->params.npoints,
|
map01bf@622
|
357 &spec->params.ntracks,
|
map01bf@622
|
358 &includeKeys, &excludeKeys,
|
map01bf@624
|
359 &radius, &absThres, &relThres, &durRatio, &hop,
|
map01bf@624
|
360 &resFmt
|
map01bf@622
|
361 );
|
map01bf@622
|
362
|
map01bf@622
|
363 if (!ok) {return NULL;}
|
mas01cr@671
|
364 current_db = (adb_t *)PyCObject_AsVoidPtr(incoming);
|
map01bf@622
|
365
|
map01bf@622
|
366 if (exhaustive){
|
map01bf@622
|
367 spec->qid.flags = spec->qid.flags | ADB_QID_FLAG_EXHAUSTIVE;
|
map01bf@622
|
368 }
|
map01bf@622
|
369 if (falsePositives){
|
map01bf@622
|
370 spec->qid.flags = spec->qid.flags | ADB_QID_FLAG_ALLOW_FALSE_POSITIVES;
|
map01bf@622
|
371 }
|
map01bf@622
|
372
|
map01bf@622
|
373 //set up spec->params
|
map01bf@622
|
374 if (strcmp(accuMode,"db")){
|
map01bf@622
|
375 spec->params.accumulation = ADB_ACCUMULATION_DB;
|
map01bf@622
|
376 } else if (strcmp(accuMode,"track")){
|
map01bf@622
|
377 spec->params.accumulation = ADB_ACCUMULATION_PER_TRACK;
|
map01bf@622
|
378 } else if (strcmp(accuMode,"one2one")){
|
map01bf@622
|
379 spec->params.accumulation = ADB_ACCUMULATION_ONE_TO_ONE;
|
map01bf@622
|
380 } else{
|
map01bf@624
|
381 PyErr_SetString(PyExc_ValueError,
|
map01bf@624
|
382 "Poorly specified distance mode. distance must either be \'db\', \'track\' or \'one2one\'.\n");
|
map01bf@622
|
383 return NULL;
|
map01bf@622
|
384 }
|
map01bf@622
|
385 if (strcmp(distMode, "dot")){
|
map01bf@622
|
386 spec->params.distance = ADB_DISTANCE_DOT_PRODUCT;
|
map01bf@622
|
387 }else if (strcmp(distMode, "eucNorm")){
|
map01bf@622
|
388 spec->params.distance = ADB_DISTANCE_EUCLIDEAN_NORMED;
|
map01bf@622
|
389 }else if (strcmp(distMode, "euclidean")){
|
map01bf@622
|
390 spec->params.distance = ADB_DISTANCE_EUCLIDEAN;
|
map01bf@622
|
391 }else{
|
map01bf@624
|
392 PyErr_SetString(PyExc_ValueError,
|
map01bf@624
|
393 "Poorly specified distance mode. distance must either be \'dot\', \'eucNorm\' or \'euclidean\'.\n");
|
map01bf@622
|
394 return NULL;
|
map01bf@622
|
395 }
|
map01bf@622
|
396
|
map01bf@622
|
397 //set up spec->refine
|
map01bf@622
|
398 //include/exclude keys
|
map01bf@622
|
399 if (includeKeys){
|
map01bf@622
|
400 if (!PyList_Check(includeKeys)){
|
map01bf@624
|
401 PyErr_SetString(PyExc_TypeError, "Include keys must be specified as a list of strings.\n");
|
map01bf@622
|
402 return NULL;
|
map01bf@622
|
403 }
|
map01bf@622
|
404 spec->refine.flags = spec->refine.flags | ADB_REFINE_INCLUDE_KEYLIST;
|
map01bf@622
|
405 spec->refine.include.nkeys = (uint32_t)PyList_Size(includeKeys);
|
map01bf@622
|
406 spec->refine.include.keys = (const char **)calloc(sizeof(const char *), spec->refine.include.nkeys);
|
map01bf@622
|
407 for (i=0;i<spec->refine.include.nkeys;i++){
|
map01bf@622
|
408 if (PyString_Check(PyList_GetItem(includeKeys, (Py_ssize_t)i))){
|
map01bf@622
|
409 spec->refine.include.keys[i] = PyString_AsString(PyList_GetItem(includeKeys, (Py_ssize_t)i));
|
map01bf@622
|
410 }else{
|
map01bf@624
|
411 PyErr_SetString(PyExc_TypeError, "Include keys must each be specified as a string.\nFound one that was not.\n");
|
map01bf@622
|
412 return NULL;
|
map01bf@622
|
413 }
|
map01bf@622
|
414 }
|
map01bf@622
|
415 }
|
map01bf@622
|
416 if (excludeKeys){
|
map01bf@622
|
417 if (!PyList_Check(excludeKeys)){
|
map01bf@624
|
418 PyErr_SetString(PyExc_TypeError, "Exclude keys must be specified as a list of strings.\n");
|
map01bf@622
|
419 return NULL;
|
map01bf@622
|
420 }
|
map01bf@622
|
421 spec->refine.flags = spec->refine.flags | ADB_REFINE_EXCLUDE_KEYLIST;
|
map01bf@622
|
422 spec->refine.exclude.nkeys = (uint32_t)PyList_Size(excludeKeys);
|
map01bf@622
|
423 spec->refine.exclude.keys = (const char **)calloc(sizeof(const char *), spec->refine.exclude.nkeys);
|
map01bf@622
|
424 for (i=0;i<spec->refine.exclude.nkeys;i++){
|
map01bf@622
|
425 if (PyString_Check(PyList_GetItem(excludeKeys, (Py_ssize_t)i))){
|
map01bf@622
|
426 spec->refine.exclude.keys[i] = PyString_AsString(PyList_GetItem(excludeKeys, (Py_ssize_t)i));
|
map01bf@622
|
427 }else{
|
map01bf@624
|
428 PyErr_SetString(PyExc_TypeError, "Exclude keys must each be specified as a string.\nFound one that was not.\n");
|
map01bf@622
|
429 return NULL;
|
map01bf@622
|
430 }
|
map01bf@622
|
431 }
|
map01bf@622
|
432 }
|
map01bf@622
|
433 //the rest of spec->refine
|
map01bf@622
|
434 if (radius){
|
map01bf@622
|
435 spec->refine.flags = spec->refine.flags | ADB_REFINE_RADIUS;
|
map01bf@622
|
436 spec->refine.radius = radius;
|
map01bf@622
|
437 }
|
map01bf@622
|
438 if (absThres){
|
map01bf@622
|
439 spec->refine.flags = spec->refine.flags | ADB_REFINE_ABSOLUTE_THRESHOLD;
|
map01bf@622
|
440 spec->refine.absolute_threshold = absThres;
|
map01bf@622
|
441 }
|
map01bf@622
|
442 if (relThres){
|
map01bf@622
|
443 spec->refine.flags = spec->refine.flags | ADB_REFINE_RELATIVE_THRESHOLD;
|
map01bf@622
|
444 spec->refine.relative_threshold = relThres;
|
map01bf@622
|
445 }
|
map01bf@622
|
446 if (durRatio){
|
map01bf@622
|
447 spec->refine.flags = spec->refine.flags | ADB_REFINE_DURATION_RATIO;
|
map01bf@622
|
448 spec->refine.duration_ratio = durRatio;
|
map01bf@622
|
449 }
|
map01bf@622
|
450 if (hop){
|
map01bf@622
|
451 spec->refine.flags = spec->refine.flags | ADB_REFINE_HOP_SIZE;
|
mas01cr@679
|
452 /* not ideal but a temporary bandage fix */
|
mas01cr@679
|
453 spec->refine.qhopsize = hop;
|
mas01cr@679
|
454 spec->refine.ihopsize = hop;
|
map01bf@622
|
455 }
|
map01bf@622
|
456 //setup the datum
|
map01bf@622
|
457 spec->qid.datum->data = NULL;
|
map01bf@622
|
458 spec->qid.datum->power = NULL;
|
map01bf@622
|
459 spec->qid.datum->times = NULL;
|
map01bf@622
|
460 //grab the datum from the key
|
map01bf@622
|
461 ok = audiodb_retrieve_datum(current_db, key, spec->qid.datum);
|
map01bf@622
|
462 if (ok != 0){
|
map01bf@622
|
463 PyErr_SetString(PyExc_RuntimeError, "Encountered an error while trying to retrieve the data associated with the passed key.\n");
|
map01bf@622
|
464 return NULL;
|
map01bf@622
|
465 }
|
map01bf@622
|
466 result = audiodb_query_spec(current_db, spec);
|
map01bf@622
|
467 if (result == NULL){
|
map01bf@718
|
468 PyErr_SetString(PyExc_RuntimeError, "Encountered an error while running the actual query, or there was nothing returned.\n");
|
map01bf@622
|
469 return NULL;
|
map01bf@622
|
470 }
|
map01bf@624
|
471 if(strcmp(resFmt, "dict")==0){
|
map01bf@624
|
472 outgoing = PyDict_New();
|
map01bf@624
|
473 for (i=0;i<result->nresults;i++){
|
mas01cr@672
|
474 thisKey = PyString_FromString(result->results[i].ikey);
|
map01bf@624
|
475 if (!PyDict_Contains(outgoing, thisKey)){
|
map01bf@624
|
476 newBits = Py_BuildValue("[(dII)]",
|
map01bf@624
|
477 result->results[i].dist,
|
map01bf@624
|
478 result->results[i].qpos,
|
map01bf@624
|
479 result->results[i].ipos);
|
map01bf@624
|
480 if (PyDict_SetItem(outgoing, thisKey,newBits)){
|
mas01cr@672
|
481 printf("key : %s\ndist : %f\nqpos : %i\nipos : %i\n", result->results[i].ikey, result->results[i].dist, result->results[i].qpos, result->results[i].ipos);
|
map01bf@624
|
482 PyErr_SetString(PyExc_AttributeError, "Error adding a tuple to the result dict\n");
|
map01bf@624
|
483 Py_XDECREF(newBits);
|
map01bf@624
|
484 return NULL;
|
map01bf@624
|
485 }
|
map01bf@624
|
486 Py_DECREF(newBits);
|
map01bf@624
|
487 }else {
|
map01bf@624
|
488 //the key already has a value, so we need to fetch the value, confirm it's a list and append another tuple to it.
|
map01bf@624
|
489 currentValue = PyDict_GetItem(outgoing, thisKey);
|
map01bf@624
|
490 if (!PyList_Check(currentValue)){
|
map01bf@624
|
491 PyErr_SetString(PyExc_TypeError, "The result dictionary appears to be malformed.\n");
|
map01bf@624
|
492 return NULL;
|
map01bf@624
|
493 }
|
map01bf@624
|
494 newBits = Py_BuildValue("dII",result->results[i].dist,
|
map01bf@624
|
495 result->results[i].qpos,
|
map01bf@624
|
496 result->results[i].ipos);
|
map01bf@624
|
497 if (PyList_Append(currentValue, newBits)){
|
map01bf@624
|
498 //error msg here
|
map01bf@624
|
499 Py_XDECREF(newBits);
|
map01bf@624
|
500 return NULL;
|
map01bf@624
|
501 }
|
map01bf@624
|
502 if (PyDict_SetItem(outgoing, thisKey, newBits)){
|
map01bf@624
|
503 PyErr_SetString(PyExc_AttributeError, "Error adding a tuple to the result dict\n");
|
map01bf@624
|
504 Py_XDECREF(newBits);
|
map01bf@624
|
505 return NULL;
|
map01bf@624
|
506 }
|
map01bf@624
|
507 Py_DECREF(newBits);
|
map01bf@624
|
508
|
map01bf@624
|
509 }
|
map01bf@624
|
510 }
|
map01bf@624
|
511 }else if(strcmp(resFmt, "list")==0){
|
map01bf@624
|
512 outgoing = PyList_New((Py_ssize_t)0);
|
map01bf@624
|
513 for (i=0;i<result->nresults;i++){
|
mas01cr@672
|
514 newBits = Py_BuildValue("sdII",result->results[i].ikey,
|
map01bf@622
|
515 result->results[i].dist,
|
map01bf@622
|
516 result->results[i].qpos,
|
map01bf@622
|
517 result->results[i].ipos);
|
map01bf@624
|
518 if (PyList_Append(outgoing, newBits)){
|
map01bf@624
|
519 //error msg here
|
map01bf@622
|
520 Py_XDECREF(newBits);
|
map01bf@622
|
521 return NULL;
|
map01bf@622
|
522 }
|
map01bf@622
|
523 Py_DECREF(newBits);
|
map01bf@624
|
524 }
|
map01bf@624
|
525 if(PyList_Reverse(outgoing)){//need to do this as things come off the accumulator backward.
|
map01bf@720
|
526 PyErr_SetString(PyExc_RuntimeError,
|
map01bf@720
|
527 "the reverse failed, hopefully a sensable error will follow.\nIf not, fix it.\n");
|
map01bf@624
|
528 return NULL;
|
map01bf@622
|
529 }
|
map01bf@624
|
530 }else{
|
map01bf@624
|
531 PyErr_SetString(PyExc_ValueError,
|
map01bf@624
|
532 "Poorly specified result mode. Result must be either \'dist\' or \'list\'.\n");
|
map01bf@624
|
533 return NULL;
|
map01bf@622
|
534 }
|
map01bf@628
|
535 if (audiodb_query_free_results(current_db, spec, result)){
|
map01bf@622
|
536 printf("bit of trouble freeing the result and spec...\ncheck for leaks.");
|
map01bf@622
|
537 }
|
map01bf@622
|
538
|
map01bf@622
|
539 return outgoing;
|
map01bf@622
|
540
|
map01bf@622
|
541
|
map01bf@622
|
542
|
map01bf@622
|
543 }
|
map01bf@622
|
544
|
map01bf@622
|
545
|
map01bf@622
|
546
|
map01bf@622
|
547
|
map01bf@620
|
548 /* close a database */
|
map01bf@620
|
549 /* api call: */
|
mas01cr@671
|
550 // void audiodb_close(adb_t *db);
|
map01bf@620
|
551 static void _pyadb_close(void *ptr)
|
map01bf@620
|
552 {
|
mas01cr@671
|
553 adb_t *stale_database;
|
mas01cr@671
|
554 stale_database = (adb_t *)ptr;
|
map01bf@620
|
555
|
map01bf@620
|
556 audiodb_close(stale_database);
|
map01bf@620
|
557 }
|
map01bf@620
|
558
|
map01bf@620
|
559 static PyMethodDef _pyadbMethods[] =
|
map01bf@620
|
560 {
|
map01bf@620
|
561 { "_pyadb_create", _pyadb_create, METH_VARARGS,
|
map01bf@620
|
562 "_pyadb_create(string path, unsigned datasize, unsigned ntracks, unsigned datadim)->adb object"},
|
map01bf@620
|
563 { "_pyadb_open", _pyadb_open, METH_VARARGS,
|
map01bf@620
|
564 "_pyadb_open(string path, [\'r\'|\'w\'])->adb object\nNote that specifing \'w\' opens the file in read and write mode. \
|
map01bf@620
|
565 There is currently no way to open in write only."},
|
map01bf@620
|
566 { "_pyadb_status", _pyadb_status, METH_VARARGS,
|
mas01cr@671
|
567 "_status(adb_t *)->(numFiles, dims, dudCount, nullCount, flags, length, data_region_size)"},
|
map01bf@620
|
568 { "_pyadb_l2norm", _pyadb_l2norm, METH_VARARGS,
|
mas01cr@671
|
569 "_pyadb_l2norm(adb_t *)->int return code (0 for sucess)"},
|
map01bf@620
|
570 { "_pyadb_power", _pyadb_power, METH_VARARGS,
|
mas01cr@671
|
571 "_pyadb_power(adb_t *)->int return code (0 for sucess)"},
|
map01bf@716
|
572 {"_pyadb_insertFromArray", (PyCFunction)_pyadb_insertFromArray, METH_VARARGS | METH_KEYWORDS,
|
map01bf@718
|
573 "_pyadb_insertFromArray(adb_t *, features=ndarray, [power=ndarray | key=keystring | times=ndarray])->\
|
map01bf@718
|
574 int return code (0 for sucess)\n\
|
map01bf@718
|
575 insert feature data from a numpy array\n\
|
map01bf@716
|
576 array given should have ndarray.shape = (numDims*numVectors,)\n\
|
map01bf@716
|
577 array datatype needs to be doubles (float may work...)\n\
|
map01bf@716
|
578 if power is given, must be 1d array of length numVectors\n\
|
map01bf@716
|
579 if times is given, must be 1d array of length 2*numVectors like this:\n\
|
map01bf@716
|
580 int audiodb_insert_datum(adb_t *, const adb_datum_t *);"},
|
map01bf@622
|
581 { "_pyadb_insertFromFile", (PyCFunction)_pyadb_insertFromFile, METH_VARARGS | METH_KEYWORDS,
|
mas01cr@671
|
582 "_pyadb_insertFromFile(adb_t *, features=featureFile, [power=powerfile | key=keystring | times=timingFile])->\
|
map01bf@621
|
583 int return code (0 for sucess)"},
|
map01bf@622
|
584 { "_pyadb_queryFromKey", (PyCFunction)_pyadb_queryFromKey, METH_VARARGS | METH_KEYWORDS,
|
map01bf@622
|
585 "base query. The nomenclature here is about a far away as pythonic as is possible.\n\
|
map01bf@622
|
586 This should be taken care of via the higher level python structure\n\
|
map01bf@622
|
587 returns a dict that should be result ordered and key = result key\n\
|
map01bf@622
|
588 and value is a list of tuples one per result associated with that key, of the form:\n\
|
map01bf@622
|
589 \t(dist, qpos, ipos)\n\
|
map01bf@622
|
590 Note as well that this is by no means the most efficient way to cast from C, simply the most direct\n\
|
map01bf@622
|
591 and what it lacks in effeciency it gains in python side access. It remains to be seen if this is\n\
|
map01bf@622
|
592 a sensible trade.\n\
|
mas01cr@671
|
593 _pyadb_queryFromKey(adb_t *, query key,\n\
|
map01bf@622
|
594 [seqLength = Int Sequence Length, \n\
|
map01bf@622
|
595 seqStart = Int offset from start for key, \n\
|
map01bf@622
|
596 exhaustive = boolean - True for exhaustive (false by default),\n\
|
map01bf@622
|
597 falsePositives= boolean - True to keep fps (false by defaults),\n\
|
map01bf@622
|
598 accumulation = [\"db\"|\"track\"|\"one2one\"] (\"db\" by default),\n\
|
map01bf@622
|
599 distance = [\"dot\"|\"eucNorm\"|\"euclidean\"] (\"dot\" by default),\n\
|
map01bf@622
|
600 npoints = int number of points per track,\n\
|
map01bf@622
|
601 ntracks = max number of results returned in db accu mode,\n\
|
map01bf@622
|
602 includeKeys = list of strings to include (use all by default),\n\
|
map01bf@622
|
603 excludeKeys = list of strings to exclude (none by default),\n\
|
map01bf@622
|
604 radius = double of nnRadius (1.0 default, overrides npoints if specified),\n\
|
map01bf@622
|
605 absThres = double absolute power threshold (db must have power),\n\
|
map01bf@622
|
606 relThres = double relative power threshold (db must have power),\n\
|
map01bf@622
|
607 durRatio = double time expansion/compresion ratio,\n\
|
map01bf@628
|
608 hopSize = int hopsize (1 by default)])->resultDict\n\
|
map01bf@628
|
609 resFmt = [\"list\"|\"dict\"](\"dict\" by default)"},
|
map01bf@620
|
610 {NULL,NULL, 0, NULL}
|
map01bf@620
|
611 };
|
map01bf@620
|
612
|
map01bf@620
|
613 void init_pyadb()
|
map01bf@620
|
614 {
|
map01bf@718
|
615 Py_InitModule3("_pyadb", _pyadbMethods, "internal c bindings for audioDB. Use pyadb for pythonic access to adb.");
|
map01bf@718
|
616 import_array();
|
map01bf@620
|
617 return;
|
map01bf@620
|
618 }
|
map01bf@620
|
619
|
map01bf@620
|
620
|