zoulingqi

sorry,I confused the two functions. "bandE[i]=sqrt(sum[i])"is in the fuction compute_band_energy. Ex[i] and Ep[i] is already the ||x|| and ||p||,so Exp[i] should equal to Exp[i]/(1e-37+Ex[i]*Ep[i]) instead of Exp[i]/sqrt(1e-37+Ex[i]*Ep[i]). void compute_band_energy(float *bandE,...