In [1]:
import sys
sys.path.append('../')
In [2]:
import numpy as np
from zephyr.backend import MiniZephyr, Eurus, SparseKaiserSource
In [3]:
import matplotlib.pyplot as plt
import matplotlib.cm as cm
import matplotlib
%matplotlib inline
from IPython.display import set_matplotlib_formats
set_matplotlib_formats('png')
matplotlib.rcParams['savefig.dpi'] = 150 # Change this to adjust figure size
In [4]:
systemConfig = {
'dx': 1., # m
'dz': 1., # m
'c': 2500., # m/s
'rho': 1., # kg/m^3
'nx': 10, # count
'nz': 10, # count
'freq': 2e2, # Hz
'nPML': 2,
}
nx = systemConfig['nx']
nz = systemConfig['nz']
dx = systemConfig['dx']
dz = systemConfig['dz']
In [5]:
MZ = MiniZephyr(systemConfig)
EU = Eurus(systemConfig)
SKS = SparseKaiserSource(systemConfig)
In [6]:
import scipy
In [7]:
plt.figure()
plt.spy(MZ.A, markersize=0.1)
plt.title('MiniZephyr')
plt.figure()
plt.spy(EU.A.tolil()[:100,:100], markersize=0.1)
plt.title('Eurus')
Out[7]:
In [ ]: