Python RiemannSolver.mx примеры использования

Пример #1

Файл: riemannSolverSerialTimer.py Проект: CodyCasteneda/pyclaw

rs.q = random.random((rs.meqn,rs.mx))
rs.aux = random.random((rs.maux,rs.mx))





 




j = 0
for i in range(max_mx/noOfTestSamples, max_mx,max_mx/noOfTestSamples):
	mxValues[j] = i
	rs.mx = i
	print
	print "Iteration",j+1,"of", noOfTestSamples
	print "mx is", i
	
	rs.q = random.random((rs.meqn,rs.mx))
	rs.aux = random.random((rs.maux,rs.mx))
	
	waves1, s1  = rs.solveVectorized(timer = True)
	timeResultsVectorized[j] = rs.elapsedTime
	
	
	
	waves2, s2  = rs.solvePointwize(timer = True)
	timeResultsPointwise[j] = rs.elapsedTime
	

Пример #2

timeResultsPointwise = empty((noOfTestSamples))
# error when max_mx is devisible by noOfTestSamples. size of arrays should be noOfTestSamples-1, use j to decide what to plot

rs = RiemannSolver(timeSteps=max_timeSteps,
                   mwaves=2,
                   mx=max_mx / noOfTestSamples,
                   meqn=2,
                   maux=2)

rs.q = random.random((rs.meqn, rs.mx))
rs.aux = random.random((rs.maux, rs.mx))

j = 0
for i in range(max_mx / noOfTestSamples, max_mx, max_mx / noOfTestSamples):
    mxValues[j] = i
    rs.mx = i
    print()
    print("Iteration", j + 1, "of", noOfTestSamples)
    print("mx is", i)

    rs.q = random.random((rs.meqn, rs.mx))
    rs.aux = random.random((rs.maux, rs.mx))

    waves1, s1 = rs.solveVectorized(timer=True)
    timeResultsVectorized[j] = rs.elapsedTime

    waves2, s2 = rs.solvePointwize(timer=True)
    timeResultsPointwise[j] = rs.elapsedTime

    j = j + 1