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()