When setting an odd number of neurons for the field, the dirac was not centered for x=0.
This was due to the default FE_TONEAREST rounding method used by nearbyint() function.
Setting FE_DOWNWARD method allows to get the floor int rather than ceiling for exact midway (.5) between 2 int
When setting an odd number of neurons for the field, the dirac was not centered for x=0. This was due to the default FE_TONEAREST rounding method used by nearbyint() function.
Setting FE_DOWNWARD method allows to get the floor int rather than ceiling for exact midway (.5) between 2 int