Python Sundials - Revision history

Admin: /* NVECTOR.pxi: Taking Address of Non-Lvalue */

2014-03-26T16:49:05Z

NVECTOR.pxi: Taking Address of Non-Lvalue

← Older revision		Revision as of 16:49, 26 March 2014
Line 66:		Line 66:
	to this line:		to this line:

	<~~source lang="python"~~>		<pre>
	buffer.shape = <Py_ssize_t *>self.__shape		buffer.shape = <Py_ssize_t *>self.__shape
	</~~source~~>		</pre>

			removing the ampersand &.

	===Header Files Not Found===		===Header Files Not Found===

2014-01-30T20:17:35Z

2014-01-30T20:17:10Z

Usage

← Older revision		Revision as of 20:17, 30 January 2014
Line 156:		Line 156:

	=Usage=		=Usage=

			==Test Problem==

			The main python-sundials website gives a test problem:

			<source lang="python">
			import array
			import numpy as np
			from sundials import *

			class Ball:
			TINY = 1e-16
			def __init__(self, e = 0.7, g = 9.81):
			self.e = e
			self.g = g

			def f(self, t, y, sw):
			if sw[0]:
			return [0., 0.]
			else:
			return [-self.g, y[0]]

			def rootf(self, t, y, sw):
			return [y[1] + Ball.TINY]

			def solve(self, plot = True):
			t0 = 0
			y0 = array.array('d', [0.,1.])
			solver = CVodeSolver(RHS = self.f, ROOT = self.rootf, SW = [False],
			abstol = 1.0e-6, reltol = 1.0e-6)

			solver.init(t0,y0)

			dt = .01
			iter = solver.iter(t0, dt)
			tres = array.array('d', [t0])
			hres = array.array('d', [y0[1]])
			while True:
			try:
			t, y = next(iter)
			tres.append(t)
			hres.append(y[1])

			except CVodeRootException, info:
			if abs(info.y[0]) < 0.01:
			solver.SW[0] = True

			tres.append(info.t)
			hres.append(0.)

			solver.init(info.t, [-self.e*info.y[0], 0.])

			if t > 3.0:
			break

			if plot:
			import matplotlib
			import pylab

			pylab.plot(tres,hres, linewidth=1.0)
			pylab.xlabel('time (s)')
			pylab.ylabel('height (m)')
			pylab.title('Bouncing ball demo')
			pylab.grid(True)
			pylab.show()
			else:
			return tres, hres

			ball = Ball()
			ball.solve()
			</source>

2014-01-30T20:16:34Z

← Older revision		Revision as of 20:16, 30 January 2014
Line 1:		Line 1:
			This page details my installation procedure for the python-sundials library: https://code.google.com/p/python-sundials/

			=Installation=

			==Mac OS X==

			The usual setup.py installation procedure was followed, although I had to make several modifications to the source code and the compile line to get it working (see the [[#Errors]] section below). The final command that I used to build the library was:

	<pre>		<pre>
	sudo CC='~~clang -lsundials~~ -I/Users/charlesreid/pkg/sundials/std/include -L/Users/charlesreid/pkg/sundials/std/lib' LDFLAGS="~~-L/Users/charlesreid/pkg/sundials/std/lib~~" python setup.py build		sudo CFLAGS='-I/Users/charlesreid/pkg/sundials/std/include -L/Users/charlesreid/pkg/sundials/std/lib' LDFLAGS="$CFLAGS" python setup.py build
	</pre>		</pre>


	==Errors==		==Errors==
Line 149:		Line 154:
	sys.argv.append('bdist_wininst')		sys.argv.append('bdist_wininst')
	</pre>		</pre>

			=Usage=

2014-01-30T20:14:45Z

2014-01-30T20:03:40Z

New page

#REDIRECT [[PySundials]]