Optimization changes to parRiemann
Minor changes to parRiemann.py
Added more GPU code to riemann theta, optimized existing code. Changes block structure in CUDA code so that more than 60,000 points can be computed at once. Limit is now 1,920,000 points, if enough memory is available.