Exemplo n.º 1
0
#!/usr/bin/env python
from Numeric import *
import os, sys

# plot NumPy arrays in Matlab, if pymat is present:
try:
    import pymat
except:
    print "pymat module is not available..."; sys.exit(1)

x = arrayrange(0,4*math.pi,0.1)
m = pymat.open()
pymat.put(m, 'x', x);
pymat.eval(m, 'y = sin(x)')
pymat.eval(m, 'plot(x,y)')
y = pymat.get(m, 'y')
import time; time.sleep(4)  # wait before killing the plot...
pymat.close(m)

# compare sin(x) in Matlab and Python:
print "Matlab: sin(%g)=%g. Python: sin(%g)=%g" % \
      (x[1], sin(x[1]), x[1], y[1])
Exemplo n.º 2
0
    def close(self):
	pymat.close(self.H)
Exemplo n.º 3
0
#!/usr/bin/python
# Simple test file for the PyMat interface

import pymat
from Numeric import *
x = array([1, 2, 3, 4, 5, 6])
H = pymat.open()
pymat.put(H, 'x', x)
pymat.eval(H, 'y = dct(x)')

y = pymat.get(H, 'y')
pymat.close(H)

expect = asarray([8.573214099741124e+000,
                  -4.162561795878958e+000,
                  -1.785561006679121e-015,
                  -4.082482904638622e-001,
                  -1.638955159695808e-015,
                  -8.007889124033019e-002])
err = sum((y-expect)*(y-expect)) 
print 'Norm of the error is:', err
if err < 1e-15:
  print "(That's pretty small -- the test passes)"
else:
  print "UMMM...THAT'S A BIT BIG! The test has failed."