[ot][spam][random][crazy][random][crazy]
Undescribed Horrific Abuse, One Victim & Survivor of Many
gmkarl at gmail.com
Mon Nov 21 01:21:02 PST 2022
# this is after all summing (using multiplication by cos instead of addition):
(Pdb) p longvec[0]
8.797795339894686
# this is supposed to be the same, and isn't
(Pdb) p (inserting_spectrum @ extracting_ift)[0]
0.5488135039273254
0413
(Pdb) complex_shortspace_freqs = fftfreq(complex = True, dc_offset =
True, repetition_samples = len(randvec))
(Pdb) p complex_freq2time = create_freq2time(freqs=complex_shortspace_freqs)
(Pdb) p complex_freq2time[:,0]
array([0.0625+0.j, 0.0625+0.j, 0.0625+0.j, 0.0625+0.j, 0.0625+0.j,
0.0625+0.j, 0.0625+0.j, 0.0625+0.j, 0.0625+0.j, 0.0625+0.j,
0.0625+0.j, 0.0625+0.j, 0.0625+0.j, 0.0625+0.j, 0.0625+0.j,
0.0625+0.j])
Since the index is zero, the first sample is the average of all the magnitudes.
(Pdb) p extracting_ift[:,0]
array([ 0.0625, -0. , 0.125 , -0. , 0.125 , -0. , 0.125 ,
-0. , 0.125 , -0. , 0.125 , -0. , 0.125 , -0. ,
0.125 , -0. , 0.0625, -0. ])
extracting_ift is real-domain-only, so it is supposed to store some
values doubled. It looks like they cancel and amplify.
I made the wrong frequencies. These are the longspace ones. Oops. I
should probably check the shortspace ones before the longspace ones.
(Pdb) p randvec[0]
0.5488135039273248
Looks like the first value is right; so in theory something about my
cosine code is wrong.
Here's how it's generated properly:
(Pdb) p (inserting_spectrum * extracting_ift[:,0]).sum()
0.5488135039273254
I'm remembering that I'm ignoring the doubling performed by the
real-domain transform. That must be the next mistake.
More information about the cypherpunks
mailing list