from AS1_class import Asym_slab
import numpy as np
import matplotlib.pyplot as plt
import matplotlib
import pickle
def save():
with open(
'pickles/flow_c0={}_R1={}_R2={}_K={}_M_A={}.p'.format(
slab.c0, slab.R1, slab.R2, slab.K, slab.M_A), 'wb') as f:
pickle.dump(root_array, f)
slab = Asym_slab(c0=0.6, R1=1.4, R2=2.2, K=0.5, M_A=None)
x_range = np.linspace(0, 1.5, 201)
y_range = np.linspace(-1, 1, 201)
load = True
if load:
root_array = pickle.load(
open(
'pickles/flow_c0={}_R1={}_R2={}_K={}_M_A={}.p'.format(
slab.c0, slab.R1, slab.R2, slab.K, slab.M_A), 'rb'))
else:
root_array = sol.point_finder_sp(slab.disp_rel,
x_range,
y_range,
from AS1_class import Asym_slab
import numpy as np
import matplotlib.pyplot as plt
import matplotlib
import pickle
def save():
with open(
'pickles/wavenumber_c0={}_R1={}_R2={}_K={}_M_A={}.p'.format(
slab.c0, slab.R1, slab.R2, slab.K, slab.M_A), 'wb') as f:
pickle.dump(root_array, f)
slab = Asym_slab(c0=1.3, R1=5 / 9, R2=1, K=None, M_A=0)
x_range = np.linspace(0, 4, 101)
y_range = np.linspace(0, max(slab.c1, slab.c2, 1), 101)
load = False
if load:
root_array = pickle.load(
open(
'pickles/wavenumber_c0={}_R1={}_R2={}_K={}_M_A={}.p'.format(
slab.c0, slab.R1, slab.R2, slab.K, slab.M_A), 'rb'))
else:
root_array = sol.point_finder_sp(slab.disp_rel,
x_range,
y_range,
from AS1_class import Asym_slab
import numpy as np
import matplotlib.pyplot as plt
import matplotlib
import pickle
def save():
with open(
'pickles/flow_c0={}_R1={}_R2={}_K={}_M_A={}.p'.format(
slab.c0, slab.R1, slab.R2, slab.K, slab.M_A), 'wb') as f:
pickle.dump(root_array, f)
slab = Asym_slab(c0=0.6, R1=0.1, R2=0.1, K=10.0, M_A=None)
x_range = np.linspace(1, 3, 201)
y_range = np.linspace(0, slab.c1, 201)
load = False
if load:
root_array = pickle.load(
open(
'pickles/flow_c0={}_R1={}_R2={}_K={}_M_A={}.p'.format(
slab.c0, slab.R1, slab.R2, slab.K, slab.M_A), 'rb'))
else:
root_array = sol.point_finder_sp(slab.disp_rel,
x_range,
y_range,
from AS1_class import Asym_slab
import numpy as np
import matplotlib.pyplot as plt
import matplotlib
import pickle
def save():
with open(
'pickles/wavenumber_c0={}_R1={}_R2={}_K={}_M_A={}.p'.format(
slab.c0, slab.R1, slab.R2, slab.K, slab.M_A), 'wb') as f:
pickle.dump(root_array, f)
slab = Asym_slab(c0=1.3, R1=5 / 9, R2=5 / 9, K=None, M_A=1.0)
x_range = np.linspace(0, 4, 101)
y_range = np.linspace(-max(slab.c1, slab.c2, 1), max(slab.c1, slab.c2, 1), 101)
load = False
if load:
root_array = pickle.load(
open(
'pickles/wavenumber_c0={}_R1={}_R2={}_K={}_M_A={}.p'.format(
slab.c0, slab.R1, slab.R2, slab.K, slab.M_A), 'rb'))
else:
root_array = sol.point_finder_sp(slab.disp_rel,
x_range,
y_range,
from AS1_class import Asym_slab
import numpy as np
import matplotlib.pyplot as plt
import matplotlib
import pickle
def save():
with open(
'pickles/wavenumber_c0={}_R1={}_R2={}_K={}_M_A={}.p'.format(
slab.c0, slab.R1, slab.R2, slab.K, slab.M_A), 'wb') as f:
pickle.dump(root_array, f)
slab = Asym_slab(c0=0.6, R1=5 / 3, R2=5 / 3, K=None, M_A=0.9)
x_range = np.linspace(0, 4, 101)
y_range = np.linspace(-1, max(slab.c1, slab.c2, 1), 101)
load = False
if load:
root_array = pickle.load(
open(
'pickles/wavenumber_c0={}_R1={}_R2={}_K={}_M_A={}.p'.format(
slab.c0, slab.R1, slab.R2, slab.K, slab.M_A), 'rb'))
else:
root_array = sol.point_finder_sp(slab.disp_rel,
x_range,
y_range,
from AS1_class import Asym_slab
import numpy as np
import matplotlib.pyplot as plt
import matplotlib
import pickle
def save():
with open(
'pickles/wavenumber_c0={}_R1={}_R2={}_K={}_M_A={}.p'.format(
slab.c0, slab.R1, slab.R2, slab.K, slab.M_A), 'wb') as f:
pickle.dump(root_array, f)
slab = Asym_slab(c0=0.6, R1=1.25, R2=2.5, K=None, M_A=0)
x_range = np.linspace(0, 4, 101)
y_range = np.linspace(0, max(slab.c1, slab.c2, 1), 101)
load = False
if load:
root_array = pickle.load(
open(
'pickles/wavenumber_c0={}_R1={}_R2={}_K={}_M_A={}.p'.format(
slab.c0, slab.R1, slab.R2, slab.K, slab.M_A), 'rb'))
else:
root_array = sol.point_finder_sp(slab.disp_rel,
x_range,
y_range,
pickle.dump(threshold, f)
def load(K):
return pickle.load(
open(
'pickles/threshold_density_sym_c0={}_K={}_M_A={}.p'.format(
c0, K, slab.M_A), 'rb'))
#R_range = np.append([0.05], np.linspace(0.25, 10, 79))
R_range = np.linspace(7.125, 10, 24)
K = 10.0
c0 = 0.6
slab = Asym_slab(c0=c0, R1=R_range[0], R2=R_range[0], K=K, M_A=None)
"""First attempt at finding the first complex value for R_range[0].
Performed over a larger than usual linspace to make sure that it's found."""
x_range = np.linspace(0.7, 0.75, 200)
y_range = np.linspace(0, slab.c1, 200)
loc = sol.find_first_imag(slab.disp_rel,
x_range,
y_range,
args=(slab.K, slab.M_A))
loc = np.round(loc, 6)
"""Second attempt at finding the first complex value for R_range[0].
Performed over a smaller linspace in order to find an accurate value."""
x_range = np.linspace(loc[0] - 0.02, loc[0] + 0.02, 200)
from AS1_class import Asym_slab
import numpy as np
import matplotlib.pyplot as plt
import matplotlib
import pickle
def save():
with open(
'pickles/wavenumber_c0={}_R1={}_R2={}_K={}_M_A={}.p'.format(
slab.c0, slab.R1, slab.R2, slab.K, slab.M_A), 'wb') as f:
pickle.dump(root_array, f)
slab = Asym_slab(c0=0.6, R1=1.4, R2=1.6, K=None, M_A=1)
x_range = np.linspace(0, 2, 101)
y_range = np.linspace(0, 1, 101)
load = False
if load:
root_array = pickle.load(
open(
'pickles/wavenumber_c0={}_R1={}_R2={}_K={}_M_A={}.p'.format(
slab.c0, slab.R1, slab.R2, slab.K, slab.M_A), 'rb'))
else:
root_array = sol.point_finder_sp(slab.disp_rel,
x_range,
y_range,
from AS1_class import Asym_slab
import numpy as np
import matplotlib.pyplot as plt
import matplotlib
import pickle
def save():
with open(
'pickles/wavenumber_c0={}_R1={}_R2={}_K={}_M_A={}.p'.format(
slab.c0, slab.R1, slab.R2, slab.K, slab.M_A), 'wb') as f:
pickle.dump(root_array, f)
slab = Asym_slab(c0=1.3, R1=5 / 9, R2=5 / 7, K=None, M_A=0.4)
x_range = np.linspace(0, 4, 101)
y_range = np.linspace(-max(slab.c1, slab.c2, 1), max(slab.c1, slab.c2, 1), 101)
load = False
if load:
root_array = pickle.load(
open(
'pickles/wavenumber_c0={}_R1={}_R2={}_K={}_M_A={}.p'.format(
slab.c0, slab.R1, slab.R2, slab.K, slab.M_A), 'rb'))
else:
root_array = sol.point_finder_sp(slab.disp_rel,
x_range,
y_range,
from AS1_class import Asym_slab
import numpy as np
import matplotlib.pyplot as plt
import matplotlib
import pickle
def save():
with open(
'pickles/wavenumber_c0={}_R1={}_R2={}_K={}_M_A={}.p'.format(
slab.c0, slab.R1, slab.R2, slab.K, slab.M_A), 'wb') as f:
pickle.dump(root_array, f)
slab = Asym_slab(c0=1.3, R1=10 / 13, R2=10 / 13, K=None, M_A=0)
x_range = np.linspace(0, 4, 101)
y_range = np.linspace(0, max(slab.c1, slab.c2, 1), 101)
load = False
if load:
root_array = pickle.load(
open(
'pickles/wavenumber_c0={}_R1={}_R2={}_K={}_M_A={}.p'.format(
slab.c0, slab.R1, slab.R2, slab.K, slab.M_A), 'rb'))
else:
root_array = sol.point_finder_sp(slab.disp_rel,
x_range,
y_range,
import solvers as sol
from AS1_class import Asym_slab
import numpy as np
import matplotlib.pyplot as plt
import matplotlib
import pickle
def save():
with open('pickles/wavenumber_c0={}_R1={}_R2={}_K={}_M_A={}.p'.format(
slab.c0, slab.R1, slab.R2, slab.K, slab.M_A), 'wb') as f:
pickle.dump(root_array, f)
slab = Asym_slab(c0=0.6, R1=5/3, R2=2.5, K=None, M_A=0.4)
x_range = np.linspace(0, 4, 101)
y_range = np.linspace(-max(slab.c1, slab.c2, 1), max(slab.c1, slab.c2, 1), 101)
load = False
if load:
root_array = pickle.load(open('pickles/wavenumber_c0={}_R1={}_R2={}_K={}_M_A={}.p'.format(
slab.c0, slab.R1, slab.R2, slab.K, slab.M_A), 'rb'))
else:
root_array = sol.point_finder_sp(slab.disp_rel, x_range, y_range, args=(slab.K, slab.M_A))
font = {'size': 15}
matplotlib.rc('font', **font)
plt.figure(num=None, figsize=(8, 13), dpi=80, facecolor='w', edgecolor='k')
ax = plt.subplot()
