Kernel Code (LLI) - 2022.2 English

#include <aie_api/aie_adf.hpp>

#include "kernel.hpp"

template<unsigned id>
void SecondOrderSection(
	input_window_float *idata,
	output_window_float *odata,iteration
	const float (&C_e)[48],			// run-time parameter: SIMD matrix of coefficients (even columns)
	const float (&C_o)[48]			// run-time parameter: SIMD matrix of coefficients (odd columns)
) {

	static v8float state_reg = null_v8float();

	for (auto i = 0; i < burst_cnt; i++) {

		v8float xreg_hi = window_readincr_v8(idata);
		v16float xreg = concat(state_reg, xreg_hi);

		v8float acc_e = null_v8float();
		v8float acc_o = null_v8float();

		v8float *ptr_coeff_e = (v8float *)(&C_e[0]);
		v8float *ptr_coeff_o = (v8float *)(&C_o[0]);

		for (auto j = 0; j < 6; j++)
		chess_flatten_loop
		{

			v8float coeff_e = *ptr_coeff_e++;
			acc_e = fpmac(acc_e, xreg, (2 * j + 4), 0, coeff_e, 0, 0x76543210);

			v8float coeff_o = *ptr_coeff_o++;
			acc_o = fpmac(acc_o, xreg, (2 * j + 5), 0, coeff_o, 0, 0x76543210);

		} // end for (auto j = 0; j < 6; j++)

		acc_o = fpadd(acc_o, acc_e);
		window_writeincr(odata, acc_o);

		// update states
		state_reg = xreg_hi;
		state_reg = upd_elem(state_reg, 4, ext_elem(acc_o, 6));
		state_reg = upd_elem(state_reg, 5, ext_elem(acc_o, 7));

	} // end for (auto i = 0; i < burst_cnt; i++)

} // end SecondOrderSection()