/* SPDX-License-Identifier: BSD-3-Clause
 *
 * Copyright(c) 2021 Intel Corporation. All rights reserved.
 *
 * Author: Shriram Shastry <malladi.sastry@linux.intel.com>
 *
 */

#ifndef __SOF_MATH_CORDIC_H__
#define __SOF_MATH_CORDIC_H__

#include <stdint.h>

/**
 * cordic_atan2_lookup_table = atan(2.^-(0:N-1)) N = 31/16
 * CORDIC Gain is cordic_gain = prod(sqrt(1 + 2.^(-2*(0:31/16-1))))
 * Inverse CORDIC Gain,inverse_cordic_gain = 1 / cordic_gain
 */
static const int32_t cordic_lookup[CORDIC_31B_TABLE_SIZE] = { 843314857, 497837829,
	263043837, 133525159, 67021687, 33543516, 16775851, 8388437, 4194283, 2097149,
	1048576, 524288, 262144, 131072, 65536, 32768, 16384, 8192, 4096, 2048, 1024,
	512, 256, 128, 64, 32, 16, 8, 4, 2, 1 };
 /* calculate LUT = 2*atan(2.^(-1:-1:(nIters)+1)) */
static const int32_t cordic_ilookup[CORDIC_30B_ITABLE_SIZE] = {
	497837829, 263043836, 133525158, 67021686, 33543515, 16775850,
	8388437, 4194282, 2097149, 1048575, 524287, 262143,
	131071, 65535, 32767, 16383, 8191, 4095,
	2047, 1023, 511, 255, 127, 63,
	31, 15, 8, 4, 2, 1};

#endif /* __SOF_MATH_CORDIC_H__ */