fixed idx mapping problem
fixed hermitian matrix multiplication
2D heatbath works, checked against fermiqcd
it works! it works! heathbath N-D SU(Nc) checked vs fermiqcd
cleanp