subroutine freeLMatFactors()
if (allocated(ycoulomb)) deallocate(ycoulomb)
if (allocated(qwprod)) deallocate(qwprod)
if (allocated(lMatCalcHKeys)) deallocate(lMatCalcHKeys)
if (allocated(lMatCalcHVals)) deallocate(lMatCalcHVals)
end subroutine freeLMatFactors