removed hcross function, added ComputeHcross and AddCrossTerm function
redid changes from temporarily removed conditions Nc==1 and GetSize()==1, preparing for parallelization