Currently, the wave functions are stored in the minimal representation. However, when transformed to sample values with I(), as occurs in getE() and H(), the output lives in the full FFT box and consumes a huge amount of memory (sixteen times more than is needed). To avoid this, we should apply I() only to one column of the wave functions at a time. Make the following modifications to achieve this.