ohsl 0.12.0

A collection of numerical routines and mathematical types for use in scientific computing.
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72

#!/usr/bin/env python
import os
import subprocess
import numpy as np
import matplotlib.mlab as mlab
import matplotlib.pyplot as plt
from argparse import ArgumentParser

parser = ArgumentParser()
parser.add_argument('integers', metavar='INTS', type=int, nargs='+',
                    help='Values for J and N steps.')

args = parser.parse_args()

if len(args.integers)==2:
    J = args.integers[0]
    N = args.integers[1]
else:
    raise ValueError("J and N step values not specified.")

slowdown = 5

data = np.loadtxt( "./DATA/Heat_equation_J_"+ str(J) + "_N_" + str(N) + ".dat" )

x = data[:,0]
t = data[:,1]
u = data[:,2]

x_axis = x[0:J+1]
time = t[0:(J+1)*(N+1):J+1]
dt = time[1] - time[0]
Tmax = time[-1]

u.shape = (N+1,J+1)

files = []
print("Saving frames ...")

for i in range(N+1):
    plt.cla()
    axes = plt.gca()
    axes.set_xlim([0,1])
    axes.set_ylim([0,110])
    plt.xlabel('x')
    plt.ylabel('u(x)')
    title = 'Heat equation'
    plt.title(title)
    time_str = 't=%1.2f' %time[i]
    plt.text(0.07, 0.94, time_str,
         horizontalalignment='center',
         transform = axes.transAxes)
    plt.plot(x_axis, u[0,:], "k--", label="u(x,0)")
    plt.plot(x_axis, u[i,:], "k-", label="u(x,t)")
    plt.legend(loc='upper right')
    fname = '_tmp%03d.png' % i
    plt.savefig(fname)
    files.append(fname)

print("Saving complete.")
print("Making movie Heat_equation.avi ...")
fps = N / (Tmax * slowdown)
FNULL = open(os.devnull, 'w')  # For supressing subprocess output
subprocess.call("mencoder 'mf://_tmp*.png' -mf type=png:fps=" + str(fps) +
                " -ovc lavc -lavcopts vcodec=msmpeg4v2 -oac mp3lame -o " +
                "./DATA/Heat_equation.avi", shell=True, stdout=FNULL,
                stderr=subprocess.STDOUT)
#vcodec = wmv2 for .mpg and -oac copy -o
print("Movie complete.")

# cleanup
for fname in files:
    os.remove(fname)