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

//		type == Property::GreensFunction::Type::Retarded,
		type == Type::Retarded,
		"ChebyshevExpander::generateGreensFunctionGPU()",
		"Only evaluation of retarded Green's function is implemented for GPU so far.",
		"Use CPU evaluation instead."
	);

/*	complex<double> *greensFunctionData = new complex<double>[lookupTableResolution];

	for(int e = 0; e < lookupTableResolution; e++)
		greensFunctionData[e] = 0.;*/

			coefficients[coefficientMap[n]][0] = jIn1[n];
//			coefficients[coefficientMap[n]*numCoefficients] = jIn1[n];

	const int numHoppingAmplitudes
		= hoppingAmplitudeSet->getNumMatrixElements();
	const int *cooHARowIndices_host
		= hoppingAmplitudeSet->getCOORowIndices();
	const int *cooHAColIndices_host
		= hoppingAmplitudeSet->getCOOColIndices();
	const complex<double> *cooHAValues_host
		= hoppingAmplitudeSet->getCOOValues();

	//Initialize GPU
	complex<double> *jIn1_device;
