import RiemannSolver
from RiemannSolver import *
from pylab import plot, figure, suptitle
rs = RiemannSolver(timeSteps=1, mwaves=2, mx=800, meqn=2, maux=2)
rs.q = random.random((rs.meqn, rs.mx))
rs.aux = random.random((rs.maux, rs.mx))
waves1, s1 = rs.solveVectorized(timer=True)
waves2, s2 = rs.solvePointwize(timer=True)
figure(1)
waveIndex = 0
componentIndex = 0
suptitle(
"Figure 1: shows the curve of wave number {0}, for the component number {1}".format(
waveIndex + 1, componentIndex + 2
)
)
plot(waves1[componentIndex, waveIndex, :], "g")
from __future__ import absolute_import
import RiemannSolver
from RiemannSolver import *
from pylab import plot, figure, suptitle
rs = RiemannSolver(timeSteps=1, mwaves=2, mx=800, meqn=2, maux=2)
rs.q = random.random((rs.meqn, rs.mx))
rs.aux = random.random((rs.maux, rs.mx))
waves1, s1 = rs.solveVectorized(timer=True)
waves2, s2 = rs.solvePointwize(timer=True)
figure(1)
waveIndex = 0
componentIndex = 0
suptitle(
"Figure 1: shows the curve of wave number {0}, for the component number {1}"
.format(waveIndex + 1, componentIndex + 2))
plot(waves1[componentIndex, waveIndex, :], "g")
