## Define a mesh
Npts = 64
L = numpy.pi * 2.0
dim = 2
gamma = 1.4
problem = 'cylinder'
Lp = L * (Npts-1.0) / Npts
sys.path.append('../')
from meshTest import zoomMesh_solve
dxf = 4*Lp / float(Npts) * .3
xS = zoomMesh_solve(Npts,-2.*Lp,2.*Lp,-2.,2.,1.0,dxf)
def zoomMesh(i,j,k):
x = xS[i]
y = xS[j]
z = 0.0
return x,y,z
mesh_options = {}
mesh_options['coordsys'] = 3
mesh_options['function'] = zoomMesh
mesh_options['periodic'] = numpy.array([False, False, True])
mesh_options['gridPeriodic'] = numpy.array([False, False, False])
mesh_options['dim'] = 3
mesh_options['x1'] = [ -2*Lp , -2*Lp , 0.0 ]
mesh_options['xn'] = [ 2*Lp , 2*Lp , Lp ]
testName = None
## Define a mesh
#Npts = 32
L = numpy.pi * 2.0
dim = 2
gamma = 1.4
problem = 'scattering_test'
Lp = L * (Npts - 1.0) / Npts
from meshTest import zoomMesh_solve
dxf = 4 * Lp / float(Npts) * .5
xS = zoomMesh_solve(Npts, -2. * Lp, 2. * Lp, -5., 5., 3.0, dxf)
def zoomMesh(i, j, k):
x = xS[i]
y = xS[j]
z = 0.0
return x, y, z
def spec_1d(y, han=True):
if han:
ave = numpy.mean(y)
hw = numpy.hanning(y.shape[0])
y = (y - ave) * hw + ave
f1 = fftpack.fft(y)