coeffs = trexi.setup(
				function_name = function_name,
				M = N,
				h = h
			)

			# Convert to floating point
			coeffs = coeffs.toFloat()
			unique_id_string = trexi.getUniqueId()

			test_range = [-N*h*0.95, N*h*0.95]


		elif rexi_method == "cirexi":

			cirexi = CIREXI(efloat_mode=efloat_mode)

			coeffs = cirexi.setup(
				function_name = function_name,
				N = N,
				lambda_max_real = lambda_max_real,
				lambda_include_imag = lambda_include_imag
			)

			# Convert to floating point
			coeffs = coeffs.toFloat()

			unique_id_string = cirexi.getUniqueId()

			#test_range = [-lambda_include_imag*0.5, lambda_include_imag*0.5]
			test_range = [-1.0, 1.0]