Added documentation
CUDA front-end print similar cpp program as CPU front-end
Linear Solver have choices solver = {CG, BiCGStab} and preconditioner = {none, diagonal}
Moved all threadIdx.x+... to a "__forceinline__ __device__ int myID()" function in the .cuh file
Moved reduction operations to be members of CollOfScalar