Tried to implement a 'continuous angular grid', but it either does not work or is too slow compared to the current method. Also added error messages whenever invalid option combinations are chosen for symmeterization, quadrature convergence study, and the type of Hv Calculator
Optimized some indices in Vv
Added reshuffling of indices in Vv for speed up purposes
Parallelized reshuffling of indices in Vv; but seems slower for a small basis, about to try a large one.
A larger basis seems slightly slower, but only about 1percent, which may be within the variance. There is a higher average CPU usage, however; so I will go with parallelized index shuffling in Vv
Optimized and parallelized Hv = Tv + Vv; now approx 10percent faster