benda/nonlinearbaseline2025
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

[lifsuscept] show negative frequencies

Browse Source
This commit is contained in:
Jan Benda 2025-07-04 11:14:17 +02:00
parent f9ace34d57
commit 7a6054d70e
2 changed files with 11 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

BIN
data/simulations/lif-mu011-D0010-chi2.npz
View File

Binary file not shown.

16
lifsuscept.py
View File

@ -84,7 +84,8 @@ def load_lifdata(mu, D, vT=1, vR=0, tau_ref=0):
file_path = sims_path / f'lif-mu{10*mu:03.0f}-D{10000*D:04.0f}-chi2.npz' file_path = sims_path / f'lif-mu{10*mu:03.0f}-D{10000*D:04.0f}-chi2.npz'
if not file_path.exists(): if not file_path.exists():
freqs1 = np.linspace(0.01, 1.0, 500) freqs1 = np.linspace(0.01, 1.0, 500)
freqs2 = np.linspace(0.01, 1.0, 200) #freqs2 = np.linspace(0.01, 1.0, 200)
freqs2 = np.linspace(-1.01, 1.0, 400)
r0, chi1, chi2 = susceptibilities(freqs1, freqs2, mu, D, r0, chi1, chi2 = susceptibilities(freqs1, freqs2, mu, D,
tau_ref, vR, vT) tau_ref, vR, vT)
np.savez(file_path, mu=mu, D=D, vT=vT, vR=vR, np.savez(file_path, mu=mu, D=D, vT=vT, vR=vR,
@ -120,13 +121,18 @@ def plot_chi2(ax, s, r0, freqs, chi2):
vmax = 300 vmax = 300
pc = ax.pcolormesh(freqs, freqs, chi2, vmin=0, vmax=vmax, pc = ax.pcolormesh(freqs, freqs, chi2, vmin=0, vmax=vmax,
rasterized=True) rasterized=True)
ax.plot([0, 0, 1], [1, 0, 0], **s.lsDiag)
ax.set_aspect('equal') ax.set_aspect('equal')
ax.set_xlabel('$f_1$') ax.set_xlabel('$f_1$')
ax.set_ylabel('$f_2$', labelpad=6) ax.set_ylabel('$f_2$', labelpad=6)
ax.set_xlim(0, 1) #ax.set_xlim(0, 1)
ax.set_ylim(0, 1) #ax.set_ylim(0, 1)
ax.set_xticks_delta(0.2) #ax.set_xticks_delta(0.2)
ax.set_yticks_delta(0.2) #ax.set_yticks_delta(0.2)
ax.set_xlim(-1, 1)
ax.set_ylim(-1, 1)
ax.set_xticks_delta(0.5)
ax.set_yticks_delta(0.5)
cax = ax.inset_axes([1.04, 0, 0.05, 1]) cax = ax.inset_axes([1.04, 0, 0.05, 1])
cax.set_spines_outward('lrbt', 0) cax.set_spines_outward('lrbt', 0)
cb = fig.colorbar(pc, cax=cax) cb = fig.colorbar(pc, cax=cax)