 			cudaMemcpy(
 				damping_device,
 				damping,
-				hoppingAmplitudeSet->getBasisSize()*sizeof(complex<double>),
+				hoppingAmplitudeSet->getBasisSize()*sizeof(
+					complex<double>
+				),
 				cudaMemcpyHostToDevice
 			) == cudaSuccess,
 			"ChebyshevExpander::calculateCoefficientsGPU()",

	int energyResolution
){
	int e = blockIdx.x*blockDim.x + threadIdx.x;
	if(e < energyResolution)
		for(int n = 0; n < numCoefficients; n++)
			greensFunction[e] = cuCadd(greensFunction[e], cuCmul(lookupTable[n*energyResolution + e], coefficients[n]));
}

void ChebyshevExpander::loadLookupTableGPU(){


	//Calculate |j1>
	int block_size = 1024;
	int num_blocks = hoppingAmplitudeSet->getBasisSize()/block_size
		+ (hoppingAmplitudeSet->getBasisSize()%block_size == 0 ? 0:1);
	if(getGlobalVerbose() && getVerbose()){
		Streams::out << "\tCUDA Block size: " << block_size << "\n";
		Streams::out << "\tCUDA Num blocks: " << num_blocks << "\n";
