RGSpace.hermitian_decomposition() is incorrect for odd grid lengths

from nifty import *
import numpy as np

def test_hermitianity (length):
    arr=np.random.random(length)+ np.random.random(length)*1j
    s=RGSpace(length)
    x1,x2=s.hermitian_decomposition(arr)
    print x1

test_hermitianity(3)

The result should be three numbers; the first should only have a real part, and the second and third should be mutually complex conjugate.

This looks like the code implicitly (and incorrectly) assumes zerocenter=True.

Assignee Loading
Time tracking Loading