import numpy as np

N = 20
vec = np.fft.fftfreq(N) * N
print(vec)
print("\n")

mat = np.zeros((N,N))
for i in range(0, N):
    for j in range(0, N):
        mat[i,j] = vec[i]**2 + vec[j]**2
print(mat)
print("\n")
np.savetxt("build/fftfreq_test_values_even.txt", mat)

N = 19
vec = np.fft.fftfreq(N) * N
print(vec)
print("\n")
mat = np.zeros((N,N))
for i in range(0, N):
    for j in range(0, N):
        mat[i,j] = vec[i]**2 + vec[j]**2
print(mat)
print("\n")
np.savetxt("build/fftfreq_test_values_odd.txt", mat)