jacob/FFT
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Compare commits

base: jacob:da394e3bda9570f7e3a4b56845f0a72cbe2fc76d
Branches Tags
jacob:master
...
compare: jacob:784aab8312b268e2b4709021a0d2923dcaeef796
Branches Tags
jacob:master

2 Commits
da394e3bda ... 784aab8312

Author SHA1 Message Date
Jacob
784aab8312 merge 2023-05-03 11:24:55 +02:00
Jacob
9152ad8e62 scc work 2023-05-03 09:51:45 +02:00
3 changed files with 1596 additions and 12 deletions
Show Stats Expand all files Collapse all files

1592
imgs/erklaerbaer_overlay.svg Normal file
View File

File diff suppressed because one or more lines are too long
Side by Side

After

Width:  |  Height:  |  Size: 508 KiB

12
software/ditact_pic.py
View File

@ -61,14 +61,6 @@ def norm(*intenss):
return max
def norm(*intenss):
max = 1e-10
for intens in intenss:
m = np.max(intens.intens)
max = np.maximum(max, m)
return max
def plot_all(intens_rutile, intens_mono, intens_mixed):
fig, axs = plt.subplots(3, 2)
fig.set_figheight(5.2)
@ -135,8 +127,8 @@ def load():
if __name__ == "__main__":
np.random.seed(1234)
# simulate()
# np.savez("intens.npz", r=r, mo=mo, mi=mi)
simulate()
np.savez("intens.npz", r=r, mo=mo, mi=mi)
r, mo, mi = load()
max = norm(r, mo, mi)
r.intens = r.intens/max

2
software/ising.py
View File

@ -60,7 +60,7 @@ def ising(file, num):
weighted_percentage = np.array(weighted_percentage)
percentage /= np.max(percentage)
np.savez(f"ising_rect_{num}.npz",
np.savez(f"ising_rect_{seed}.npz",
w_percentage=weighted_percentage, percentage=percentage, out_1=out_rect[0],
out_2=out_rect[1], out_3=out_rect[2], out_4=out_rect[3])