[EMAIL PROTECTED] wrote:
Hullo
I'm trying to use the permutation part of gsl through pygsl in a brute force
attack on my
Permutation (14) :
approx 87 billion cases, which will take about 10 days.
...
x = pygsl.permutation.Permutation (5)
f = open ("abc", "w")
pygsl.permutation.gslwrap.gsl_permutation_fwrite (f, x)
...
it bombs out on the last line:
p.gslwrap.gsl_permutation_fwrite (f, x)
TypeError: argument number 1: a 'FILE *' is expected, 'file(<open file 'abcd',
mode 'w' at 0xb7eba380>)' is received
I'm sure that I've missed something obvious, but would appreciate any pointers.
cheers
Tim
You have to use a the typemap.
%include typemaps/file_typemaps.i
then the py object should be automatically be converted to a FILE.
But one can also use swig to get methods for the permutation.
I attached the files
permutation.i and permutation.py
please copy them to
permutation.i -> src/gslwrap
permutation.py -> examples/
So be save to rebuild the whole module:
$touch src/gslwrap/gslwrap_wrap.i
then proceed with building and installing as usual.
The permutation.py illustrates how to use the read and write. Please
note that thes functions are not save. If you read less elements then
you have written it will not result in any error.
Sincerely yours
Pierre
--
+------------------------------------------------+
Pierre Schnizer <[EMAIL PROTECTED]>
Telephon : +49 6159 71 1557
Fax : +49 6159 71 2043
Adress :
Pierre Schnizer
BTE / MT
Gesellschaft fuer Schwerionenforschung
Planckstrasse 1
64291 Darmstadt
Germany
+-------------------------------------------------+
/* -*- C -*- */
/**
* Author : Fabian Jakobs
*
* 2. 2. 02007 Pierre Schnizer <[EMAIL PROTECTED]>
* Added methods for read and write.
* WARNING: These methods are not save and perhaps not portable between
* platforms
*
*/
/*
*
*/
%include file_typemaps.i
%{
#include <gsl/gsl_permutation.h>
/*
* Return a FILE descriptor for an approbriate object
*/
static FILE *
pygsl_object_to_file(PyObject * fob)
{
FILE * fp = NULL;
FUNC_MESS_BEGIN();
HANDLE_MINGW();
if (!PyFile_Check(fob)) {
PyErr_SetString(PyExc_TypeError, "Need a file object to read the
object!");
PyGSL_add_traceback(NULL, __FILE__, __FUNCTION__, __LINE__);
return NULL;
}
FUNC_MESS("Convert Python File to C File");
fp = PyFile_AsFile(fob);
DEBUG_MESS(2, "Using file at %p with filedes %d", (void *) fp, fileno(fp));
if(fp == NULL){
gsl_error("Object could not be converted to a file object!",
__FILE__, __LINE__, GSL_ESANITY);
return NULL;
}
FUNC_MESS_END();
return fp;
}
%}
%rename(Permutation) gsl_permutation_struct;
%include gsl/gsl_permutation.h
%include <typemaps/file_typemaps.i>
%extend gsl_permutation_struct{
%rename(__getitem__) get_item;
%rename(__len__) size;
%rename(__str__) printf;
gsl_permutation_struct(size_t n) {
return gsl_permutation_calloc(n);
}
~gls_permutation_struct() {
gsl_permutation_free(self);
}
gsl_error_flag_drop _linear_to_canonical(struct gsl_permutation_struct *q){
return gsl_permutation_linear_to_canonical(q, self);
}
gsl_error_flag_drop _canonical_to_linear(struct gsl_permutation_struct *q){
return gsl_permutation_canonical_to_linear(q, self);
}
gsl_error_flag_drop _mul(struct gsl_permutation_struct *res, struct
gsl_permutation_struct *m2){
return gsl_permutation_mul(res, self, m2);
}
size_t inversions(){
return gsl_permutation_inversions(self);
}
size_t linear_cycles(){
return gsl_permutation_linear_cycles(self);
}
size_t canonical_cycles(){
return gsl_permutation_canonical_cycles(self);
}
gsl_error_flag_drop _inverse(struct gsl_permutation_struct *inv){
return gsl_permutation_inverse(inv, self);
}
size_t get_item(const size_t i) {
return gsl_permutation_get(self, i);
}
gsl_error_flag_drop swap(const size_t i, const size_t j) {
return gsl_permutation_swap(self, i, j);
}
size_t size() {
return gsl_permutation_size(self);
}
gsl_error_flag_drop valid () {
return gsl_permutation_valid(self);
}
void reverse() {
gsl_permutation_reverse(self);
}
int next() {
return gsl_permutation_next(self);
}
int prev() {
return gsl_permutation_prev(self);
}
char *printf() {
/* FIXME this must return a string */
gsl_permutation_fprintf (stdout, self, " %u");
return NULL;
}
gsl_error_flag_drop read(PyObject *fob){
FILE *fp = pygsl_object_to_file(fob);
if(fp == NULL)
return GSL_EFAILED;
return gsl_permutation_fread(fp, self);
}
gsl_error_flag_drop write(PyObject *fob){
FILE *fp = pygsl_object_to_file(fob);
if(fp == NULL)
return GSL_EFAILED;
return gsl_permutation_fwrite(fp, self);
}
PyObject *tolist(){
PyObject *a_list = NULL, *a_int;
PyGSL_array_index_t size = 0, i;
size = (long) gsl_permutation_size(self);
a_list = PyList_New(size);
if (a_list == NULL)
return NULL;
for(i=0; i<size; i++){
a_int = PyInt_FromLong((long) gsl_permutation_get(self, i));
if (a_int == NULL){
Py_DECREF(a_list);
return NULL;
}
PyList_SET_ITEM(a_list, i, a_int);
}
return a_list;
}
PyObject *toarray(){
PyArrayObject * a_array = NULL;
long *data;
PyGSL_array_index_t size, i;
size = (int) gsl_permutation_size(self);
a_array = (PyArrayObject *) PyGSL_New_Array(1, &size, PyArray_LONG);
if(a_array == NULL)
return NULL;
data = (long *) a_array->data;
for(i=0; i<size; i++){
data[i] = (long) gsl_permutation_get(self, i);
}
return (PyObject *) a_array;
}
}
#!/usr/bin/env python
"""
Illustrate permutation use
"""
from pygsl import permutation
nelements = 5
testfile = "/tmp/permutation2"
p = permutation.Permutation(nelements)
print "Initial Permutation", p.toarray()
p.next()
p.next()
print "Iterated permutation", p.toarray()
p.write(open(testfile, "w"))
p2 = permutation.Permutation(nelements)
print "New permutation" , p2.toarray()
p2.read(open(testfile))
print "After read in", p2.toarray(), "Written was", p.toarray()
-------------------------------------------------------------------------
Using Tomcat but need to do more? Need to support web services, security?
Get stuff done quickly with pre-integrated technology to make your job easier.
Download IBM WebSphere Application Server v.1.0.1 based on Apache Geronimo
http://sel.as-us.falkag.net/sel?cmd=lnk&kid=120709&bid=263057&dat=121642
_______________________________________________
pygsl-discuss mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/pygsl-discuss
