Changed three on the fly trig calculations from the Potential Pre-calculation to pre-computing them
Now also precomputing index in Potential Integration
More Comparison Calculations
Changed some 1/2 to 0.5 for speed purposes
Fixed grid spacing; had dx=x_max/n, now have dx=x_max/(n-1), where n is the number of points
Moved files around
Now precomputing the quadrature weights and spherical harmonics fully
Precomputing potential from point pieces
Parallelized potential calc with collapse 3
Parallelized spherical harmonics calc with collapse 2
Changed three parallel fors to a single parallel plus three fors
Made potential precalc a collapse 5 instead of 3 loop
Made spherical harmonics precalc a collapse 4 instead of 2 loop
Last for the day
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