 		c1 = -2 * exp(sigma * 1 / synth.T) * cos(omega * 1 / synth.T);
 		c0 = exp( 2 * sigma * 1 / synth.T);
 
-#if DEBUG == 1
-		printf("%d %d sigma %f\n", i, mu, sigma);
-		printf("      omega %f\n", omega);
-#endif
-
 		m_set(MatrixC, 0, 2*i  , 0);
 		m_set(MatrixC, 0, 2*i+1, a);
 


		m_swap(&pointer_MatrixAp_tmp, &pointer_MatrixAp);
	}

#if DEBUG == 3
	printf("MatrixA");
	m_print(MatrixA);
	printf("MatrixCA");
	m_print(MatrixCA);
	printf("state");
	m_print(state);
#endif
}



		cudaEventElapsedTime(&MatrixAp_time, MatrixAp_start, MatrixAp_stop);
		cudaEventElapsedTime(&Memcpy_time, Memcpy_start, Memcpy_stop);

		if(i >= 0) {
			m_swap(&pointer_device_state_read, &pointer_device_state_write);
			m_swap(&pointer_device_output_chunk_read, &pointer_device_output_chunk_write);
