python - Using scipy.weave and GSL ODE -


after previous question trying solve odes using gsl. complicating factor want written in python , trying scipy.weave help.

using standard gsl ode example have following based on this example

import scipy.weave weave  time =0.0 y1 = [0.0,0.0] # define function   support_code=r'''   int  func (double t, const double y[], double f[],        void *params)  {    double mu = *(double *)params;    f[0] = y[1];    f[1] = -y[0] - mu*y[1]*(y[0]*y[0] - 1);    return gsl_success;  }   int  jac (double t, const double y[], double *dfdy,        double dfdt[], void *params)  {    double mu = *(double *)params;    gsl_matrix_view dfdy_mat       = gsl_matrix_view_array (dfdy, 2, 2);    gsl_matrix * m = &dfdy_mat.matrix;     gsl_matrix_set (m, 0, 0, 0.0);    gsl_matrix_set (m, 0, 1, 1.0);    gsl_matrix_set (m, 1, 0, -2.0*mu*y[0]*y[1] - 1.0);    gsl_matrix_set (m, 1, 1, -mu*(y[0]*y[0] - 1.0));    dfdt[0] = 0.0;    dfdt[1] = 0.0;    return gsl_success;  } '''   # define main  c_main = r'''  int  main (void)  {    double mu = 10;    gsl_odeiv2_system sys = {func, jac, 2, &mu};     gsl_odeiv2_driver * d =       gsl_odeiv2_driver_alloc_y_new (&sys, gsl_odeiv2_step_rk8pd,                   1e-6, 1e-6, 0.0);    int i;    double t1 = 100.0;      (i = 1; <= 100; i++)      {        double ti = * t1 / 100.0;        int status = gsl_odeiv2_driver_apply (d, &t, ti, y);         if (status != gsl_success)     {       printf ("error, return value=%d\n", status);       break;     }         time = t;      }     gsl_odeiv2_driver_free (d);    return 0;  } '''   compiler = 'gcc' vars = ['time','y'] libs = ['gsl','gslcblas','m'] headers = ['<stdio.h>','<gsl/gsl_errno.h>','<gsl/gsl_odeiv2.h>','<gsl/gsl_matrix.h>']   res = weave.inline(c_main, vars, compiler=compiler,                    libraries=libs, headers=headers,                    support_code = support_code)  print time

which returns this horrible spew. don't have clue doing know should easy! trouble of examples of weave on web based on printf or short iterations.


Comments

Post a Comment

Popular posts from this blog

php - Why I am getting the Error "Commands out of sync; you can't run this command now" -

the documentation of error mentioned in title says if commands out of sync; can't run command now in client code, calling client functions in wrong order. this can happen, example, if using mysql_use_result() , try execute new query before have called mysql_free_result(). can happen if try execute 2 queries return data without calling mysql_use_result() or mysql_store_result() in between. from here: http://dev.mysql.com/doc/refman/5.0/en/commands-out-of-sync.html but in first query not fetching data mysql database, inserting. , in second query getting data database. here code $connection = mysqli_connect("localhost","username","password","tbl_msgs"); if(mysqli_connect_errno($connection)) { die("failed connect mysql: " . mysqli_connect_error()); } $query = "insert users (total_comments, total_views) values ({$total_comments}, {$total_views});"; $query .= "insert msgs (no...
Read more

linux - Does gcc have any options to add version info in ELF binary file? -

i mean whether gcc can insert source code version infor elf binary section or similar. not want change source file, add info gcc option in makefile. you can emit version info text file, turn text file object file statically link executable. the first step simple have write code: script or write version info in format plain text file. write makefile rule produce version.o version.txt, using objcopy . you'll have object file 2 useful symbols defined in it: beginning , end of textual version info. add generated object executable, , you'll able access version 2 ways: running strings on binary, or writing code in application print version string (you'll need declare start , end symbols variables in header file).
Read more

java - Are there any classes that implement javax.persistence.Parameter<T>? -

i'm using named queries fetch entities database using jpa , hibernate , came across puzzling me: how come there no default implementations of javax.persistence.parameter interface? i rather impressed how @staticmetamodel used criteriabuilder in jpa2 meaning type safety , don't have mess around setting string names. wanted apply similar of jpa1 ejb's i'm working on: dealergroup_.class: public dealergroup_ { public parameter<integer> id = new parameterimpl<integer>("id"); /// etc... } dealergroupfacade.class: // ... public dealergroup fetch(integer id) { typedquery<dealergroup> query = em.createnamedquery("dealergroup.fetchbyid", dealergroup.class); query.setparameter(dealergroup_.id, query); return query.getsingleresult(); } // ... i surprised find there no implementations of parameter interface - of course, implement myself absence makes me suspicious may mistake using in way? should implementing parameter...
Read more