/* USER CODE BEGIN Header */
/**
  ******************************************************************************
  * @file           : main.c
  * @brief          : Main program body
  ******************************************************************************
  * @attention
  *
  * Copyright (c) 2024 STMicroelectronics.
  * All rights reserved.
  *
  * This software is licensed under terms that can be found in the LICENSE file
  * in the root directory of this software component.
  * If no LICENSE file comes with this software, it is provided AS-IS.
  *
  ******************************************************************************
  */
/* USER CODE END Header */
/* Includes ------------------------------------------------------------------*/
#include "main.h"

/* Private includes ----------------------------------------------------------*/
/* USER CODE BEGIN Includes */
#include "stdio.h"
#include "string.h"
#include "ai_platform.h"
#include "network.h"
#include "network_data.h"
#include "arm_math.h"
/* USER CODE END Includes */

/* Private typedef -----------------------------------------------------------*/
/* USER CODE BEGIN PTD */

/* USER CODE END PTD */

/* Private define ------------------------------------------------------------*/
/* USER CODE BEGIN PD */

/* USER CODE END PD */

/* Private macro -------------------------------------------------------------*/
/* USER CODE BEGIN PM */
#define DATA_USER 512

/* Input signal sampling rate */
#define SAMPLE_RATE  16000U

#define FRAME_SIZE 512//in typical mfcc 1024,512 and 256

#define FRAME_SHIFT 256 //Always 50% of the frame_size

#define NUM_MFCC_COEFFS 13

#define NUM_FILTER_BANKS 26//always double of the num_filter_bank
/* USER CODE END PM */

/* Private variables ---------------------------------------------------------*/
ADC_HandleTypeDef hadc1;

CRC_HandleTypeDef hcrc;

UART_HandleTypeDef huart2;

/* USER CODE BEGIN PV */
int n = 0;
uint16_t ecg_x = 0.0;
float ecg_buffer[DATA_USER] = {0};
arm_rfft_fast_instance_f32 fft_instance;
arm_dct4_instance_f32 dct_instance;
arm_cfft_radix4_instance_f32 cfft_instance;
arm_rfft_instance_f32 rfft_instance;
float mfcc_features[NUM_MFCC_COEFFS * (DATA_USER - FRAME_SIZE) / FRAME_SHIFT + 1];
uint32_t num_frames = (DATA_USER - FRAME_SIZE) / FRAME_SHIFT + 1;
float32_t frame[FRAME_SIZE];
float32_t fft_output[FRAME_SIZE];
float32_t power_spectrum[FRAME_SIZE / 2 + 1];
float32_t mel_energies[NUM_FILTER_BANKS];
float32_t log_mel_energies[NUM_FILTER_BANKS];
float32_t dct_output[NUM_FILTER_BANKS];
float32_t pre_emphasized[DATA_USER];
//for cube ai
ai_handle network;
float aiInData[AI_NETWORK_IN_1_SIZE];
float aiOutData[AI_NETWORK_OUT_1_SIZE];
ai_u8 activations[AI_NETWORK_DATA_ACTIVATIONS_SIZE];
const char* activities[AI_NETWORK_OUT_1_SIZE] = {
  "1"
};
ai_buffer * ai_input;
ai_buffer * ai_output;
int write_index = 0;
/* USER CODE END PV */

/* Private function prototypes -----------------------------------------------*/
void SystemClock_Config(void);
static void MX_GPIO_Init(void);
static void MX_USART2_UART_Init(void);
static void MX_ADC1_Init(void);
static void MX_CRC_Init(void);
/* USER CODE BEGIN PFP */
void fill_ecg_buffer(void);
int __io_putchar(int ch);
void pre_emphasis(float32_t *input, float32_t *output, uint32_t length, float32_t alpha);
void apply_hamming_window(float32_t *frame, uint32_t length);
void compute_mel_filter_banks(float32_t *power_spectrum, float32_t *mel_energies);
void compute_dct(float32_t *input, float32_t *output, uint32_t length);
void compute_mfcc(float32_t *input, uint32_t length, float32_t *mfcc);
static void AI_Init(void);
static void AI_Run(float *pIn, float *pOut);
static uint32_t argmax(const float * values, uint32_t len);
/* USER CODE END PFP */

/* Private user code ---------------------------------------------------------*/
/* USER CODE BEGIN 0 */

/* USER CODE END 0 */

/**
  * @brief  The application entry point.
  * @retval int
  */
int main(void)
{

  /* USER CODE BEGIN 1 */
  uint32_t write_index = 0;
  /* USER CODE END 1 */

  /* MCU Configuration--------------------------------------------------------*/

  /* Reset of all peripherals, Initializes the Flash interface and the Systick. */
  HAL_Init();

  /* USER CODE BEGIN Init */

  /* USER CODE END Init */

  /* Configure the system clock */
  SystemClock_Config();

  /* USER CODE BEGIN SysInit */

  /* USER CODE END SysInit */

  /* Initialize all configured peripherals */
  MX_GPIO_Init();
  MX_USART2_UART_Init();
  MX_ADC1_Init();
  MX_CRC_Init();
  /* USER CODE BEGIN 2 */
  printf("Initializing the system");
  AI_Init();
  HAL_ADC_Start(&hadc1);
  arm_rfft_fast_init_f32(&fft_instance, FRAME_SIZE);
  arm_dct4_init_f32(&dct_instance, &rfft_instance, &cfft_instance, NUM_FILTER_BANKS, NUM_FILTER_BANKS / 2, 0.5);
  /* USER CODE END 2 */

  /* Infinite loop */
  /* USER CODE BEGIN WHILE */
  while (1)
  {
    /* USER CODE END WHILE */

    /* USER CODE BEGIN 3 */
	ecg_x = HAL_ADC_GetValue(&hadc1);
	ecg_buffer[n] = ecg_x;
	printf("%d ",ecg_x);
	HAL_ADC_Start(&hadc1);
	n++;
	if(n == DATA_USER) {
		n = 0;
	    compute_mfcc(ecg_buffer, DATA_USER, mfcc_features);
	    for(int i =0;i<AI_NETWORK_IN_1_SIZE;i++){
	    	aiInData[i] = mfcc_features[i];
	     }
	     if (write_index == AI_NETWORK_IN_1_SIZE) {
	    	 write_index = 0;
	         printf("Running inference\r\n");
	         AI_Run(aiInData, aiOutData);

	         /* Output results */
	         for (uint32_t i = 0; i < AI_NETWORK_OUT_1_SIZE; i++) {
	        	 printf("%8.6f ", aiOutData[i]);
	         }
	         uint32_t class = argmax(aiOutData, AI_NETWORK_OUT_1_SIZE);
	         printf(": %d - %s\r\n", (int) class, activities[class]);
	          	 }
	         }
	        // finish of the function
  }
  /* USER CODE END 3 */
}

/**
  * @brief System Clock Configuration
  * @retval None
  */
void SystemClock_Config(void)
{
  RCC_OscInitTypeDef RCC_OscInitStruct = {0};
  RCC_ClkInitTypeDef RCC_ClkInitStruct = {0};

  /** Configure the main internal regulator output voltage
  */
  __HAL_RCC_PWR_CLK_ENABLE();
  __HAL_PWR_VOLTAGESCALING_CONFIG(PWR_REGULATOR_VOLTAGE_SCALE1);

  /** Initializes the RCC Oscillators according to the specified parameters
  * in the RCC_OscInitTypeDef structure.
  */
  RCC_OscInitStruct.OscillatorType = RCC_OSCILLATORTYPE_HSI;
  RCC_OscInitStruct.HSIState = RCC_HSI_ON;
  RCC_OscInitStruct.HSICalibrationValue = RCC_HSICALIBRATION_DEFAULT;
  RCC_OscInitStruct.PLL.PLLState = RCC_PLL_ON;
  RCC_OscInitStruct.PLL.PLLSource = RCC_PLLSOURCE_HSI;
  RCC_OscInitStruct.PLL.PLLM = 8;
  RCC_OscInitStruct.PLL.PLLN = 100;
  RCC_OscInitStruct.PLL.PLLP = RCC_PLLP_DIV2;
  RCC_OscInitStruct.PLL.PLLQ = 7;
  if (HAL_RCC_OscConfig(&RCC_OscInitStruct) != HAL_OK)
  {
    Error_Handler();
  }

  /** Initializes the CPU, AHB and APB buses clocks
  */
  RCC_ClkInitStruct.ClockType = RCC_CLOCKTYPE_HCLK|RCC_CLOCKTYPE_SYSCLK
                              |RCC_CLOCKTYPE_PCLK1|RCC_CLOCKTYPE_PCLK2;
  RCC_ClkInitStruct.SYSCLKSource = RCC_SYSCLKSOURCE_PLLCLK;
  RCC_ClkInitStruct.AHBCLKDivider = RCC_SYSCLK_DIV1;
  RCC_ClkInitStruct.APB1CLKDivider = RCC_HCLK_DIV2;
  RCC_ClkInitStruct.APB2CLKDivider = RCC_HCLK_DIV1;

  if (HAL_RCC_ClockConfig(&RCC_ClkInitStruct, FLASH_LATENCY_3) != HAL_OK)
  {
    Error_Handler();
  }
}

/**
  * @brief ADC1 Initialization Function
  * @param None
  * @retval None
  */
static void MX_ADC1_Init(void)
{

  /* USER CODE BEGIN ADC1_Init 0 */

  /* USER CODE END ADC1_Init 0 */

  ADC_ChannelConfTypeDef sConfig = {0};

  /* USER CODE BEGIN ADC1_Init 1 */

  /* USER CODE END ADC1_Init 1 */

  /** Configure the global features of the ADC (Clock, Resolution, Data Alignment and number of conversion)
  */
  hadc1.Instance = ADC1;
  hadc1.Init.ClockPrescaler = ADC_CLOCK_SYNC_PCLK_DIV4;
  hadc1.Init.Resolution = ADC_RESOLUTION_10B;
  hadc1.Init.ScanConvMode = DISABLE;
  hadc1.Init.ContinuousConvMode = DISABLE;
  hadc1.Init.DiscontinuousConvMode = DISABLE;
  hadc1.Init.ExternalTrigConvEdge = ADC_EXTERNALTRIGCONVEDGE_NONE;
  hadc1.Init.ExternalTrigConv = ADC_SOFTWARE_START;
  hadc1.Init.DataAlign = ADC_DATAALIGN_RIGHT;
  hadc1.Init.NbrOfConversion = 1;
  hadc1.Init.DMAContinuousRequests = DISABLE;
  hadc1.Init.EOCSelection = ADC_EOC_SINGLE_CONV;
  if (HAL_ADC_Init(&hadc1) != HAL_OK)
  {
    Error_Handler();
  }

  /** Configure for the selected ADC regular channel its corresponding rank in the sequencer and its sample time.
  */
  sConfig.Channel = ADC_CHANNEL_0;
  sConfig.Rank = 1;
  sConfig.SamplingTime = ADC_SAMPLETIME_3CYCLES;
  if (HAL_ADC_ConfigChannel(&hadc1, &sConfig) != HAL_OK)
  {
    Error_Handler();
  }
  /* USER CODE BEGIN ADC1_Init 2 */

  /* USER CODE END ADC1_Init 2 */

}

/**
  * @brief CRC Initialization Function
  * @param None
  * @retval None
  */
static void MX_CRC_Init(void)
{

  /* USER CODE BEGIN CRC_Init 0 */

  /* USER CODE END CRC_Init 0 */

  /* USER CODE BEGIN CRC_Init 1 */

  /* USER CODE END CRC_Init 1 */
  hcrc.Instance = CRC;
  if (HAL_CRC_Init(&hcrc) != HAL_OK)
  {
    Error_Handler();
  }
  /* USER CODE BEGIN CRC_Init 2 */

  /* USER CODE END CRC_Init 2 */

}

/**
  * @brief USART2 Initialization Function
  * @param None
  * @retval None
  */
static void MX_USART2_UART_Init(void)
{

  /* USER CODE BEGIN USART2_Init 0 */

  /* USER CODE END USART2_Init 0 */

  /* USER CODE BEGIN USART2_Init 1 */

  /* USER CODE END USART2_Init 1 */
  huart2.Instance = USART2;
  huart2.Init.BaudRate = 115200;
  huart2.Init.WordLength = UART_WORDLENGTH_8B;
  huart2.Init.StopBits = UART_STOPBITS_1;
  huart2.Init.Parity = UART_PARITY_NONE;
  huart2.Init.Mode = UART_MODE_TX_RX;
  huart2.Init.HwFlowCtl = UART_HWCONTROL_NONE;
  huart2.Init.OverSampling = UART_OVERSAMPLING_16;
  if (HAL_UART_Init(&huart2) != HAL_OK)
  {
    Error_Handler();
  }
  /* USER CODE BEGIN USART2_Init 2 */

  /* USER CODE END USART2_Init 2 */

}

/**
  * @brief GPIO Initialization Function
  * @param None
  * @retval None
  */
static void MX_GPIO_Init(void)
{
  GPIO_InitTypeDef GPIO_InitStruct = {0};
/* USER CODE BEGIN MX_GPIO_Init_1 */
/* USER CODE END MX_GPIO_Init_1 */

  /* GPIO Ports Clock Enable */
  __HAL_RCC_GPIOC_CLK_ENABLE();
  __HAL_RCC_GPIOH_CLK_ENABLE();
  __HAL_RCC_GPIOA_CLK_ENABLE();
  __HAL_RCC_GPIOB_CLK_ENABLE();

  /*Configure GPIO pin Output Level */
  HAL_GPIO_WritePin(LD2_GPIO_Port, LD2_Pin, GPIO_PIN_RESET);

  /*Configure GPIO pin : B1_Pin */
  GPIO_InitStruct.Pin = B1_Pin;
  GPIO_InitStruct.Mode = GPIO_MODE_IT_FALLING;
  GPIO_InitStruct.Pull = GPIO_NOPULL;
  HAL_GPIO_Init(B1_GPIO_Port, &GPIO_InitStruct);

  /*Configure GPIO pin : LD2_Pin */
  GPIO_InitStruct.Pin = LD2_Pin;
  GPIO_InitStruct.Mode = GPIO_MODE_OUTPUT_PP;
  GPIO_InitStruct.Pull = GPIO_NOPULL;
  GPIO_InitStruct.Speed = GPIO_SPEED_FREQ_LOW;
  HAL_GPIO_Init(LD2_GPIO_Port, &GPIO_InitStruct);

/* USER CODE BEGIN MX_GPIO_Init_2 */
/* USER CODE END MX_GPIO_Init_2 */
}

/* USER CODE BEGIN 4 */
static void AI_Init(void)
{
  ai_error err;

  /* Create a local array with the addresses of the activations buffers */
  const ai_handle act_addr[] = { activations };
  /* Create an instance of the model */
  err = ai_network_create_and_init(&network, act_addr, NULL);
  if (err.type != AI_ERROR_NONE) {
    printf("ai_network_create error - type=%d code=%d\r\n", err.type, err.code);
    Error_Handler();
  }
  ai_input = ai_network_inputs_get(network, NULL);
  ai_output = ai_network_outputs_get(network, NULL);
}

static void AI_Run(float *pIn, float *pOut)
{
  ai_i32 batch;
  ai_error err;

  /* Update IO handlers with the data payload */
  ai_input[0].data = AI_HANDLE_PTR(pIn);
  ai_output[0].data = AI_HANDLE_PTR(pOut);

  batch = ai_network_run(network, ai_input, ai_output);
  if (batch != 1) {
    err = ai_network_get_error(network);
    printf("AI ai_network_run error - type=%d code=%d\r\n", err.type, err.code);
    Error_Handler();
  }
}

static uint32_t argmax(const float * values, uint32_t len)
{
  float max_value = values[0];
  uint32_t max_index = 0;
  for (uint32_t i = 1; i < len; i++) {
    if (values[i] > max_value) {
      max_value = values[i];
      max_index = i;
    }
  }
  return max_index;
}

//Now to compute the mfcc features
int __io_putchar(int ch) {
  HAL_UART_Transmit(&huart2, (uint8_t *)&ch, 1, HAL_MAX_DELAY);
  return ch;
}

// Pre-emphasis filter
void pre_emphasis(float32_t *input, float32_t *output, uint32_t length, float32_t alpha) {
// printf("\r\n Pre-Empahsis");
output[0] = input[0];
for (uint32_t i = 1; i < length; i++) {
output[i] = input[i] - alpha * input[i - 1];
// printf("After pre-emphasis[%d]: %0.2f \r\n",i,output[i]);
}
}

void apply_hamming_window(float32_t *frame, uint32_t length) {
  for (uint32_t i = 0; i < length; i++) {
    frame[i] *= 0.54 - 0.46 * cosf(2 * PI * i / (length - 1));
//    printf("After hamming window[%d]: %0.2f \r\n",i,frame[i]);
  }
}

void compute_mel_filter_banks(float32_t *power_spectrum, float32_t *mel_energies) {
  float32_t mel_min = 0;
  float32_t mel_max = 2595 * log10(1 + (SAMPLE_RATE / 2) / 700.0);
  float32_t mel_points[NUM_FILTER_BANKS + 2];
  float32_t hz_points[NUM_FILTER_BANKS + 2];
  uint32_t bin[NUM_FILTER_BANKS + 2];

  for (uint32_t i = 0; i < NUM_FILTER_BANKS + 2; i++) {
    mel_points[i] = mel_min + i * (mel_max - mel_min) / (NUM_FILTER_BANKS + 1);
    hz_points[i] = 700 * (pow(10, mel_points[i] / 2595) - 1);
    bin[i] = (uint32_t) floor((FRAME_SIZE + 1) * hz_points[i] / SAMPLE_RATE);
  }

  for (uint32_t i = 0; i < NUM_FILTER_BANKS; i++) {
    mel_energies[i] = 0;
    for (uint32_t j = bin[i]; j < bin[i + 1]; j++) {
      mel_energies[i] += (j - bin[i]) / (float32_t) (bin[i + 1] - bin[i]) * power_spectrum[j];
    }
    for (uint32_t j = bin[i + 1]; j < bin[i + 2]; j++) {
      mel_energies[i] += (bin[i + 2] - j) / (float32_t) (bin[i + 2] - bin[i + 1]) * power_spectrum[j];
    }
//    printf("After Mel-Energies[%d]: %0.2f \r\n",i,mel_energies[i]);

  }
}

void compute_dct(float32_t *input, float32_t *output, uint32_t length) {
  for (uint32_t k = 0; k < NUM_MFCC_COEFFS; k++) {
    output[k] = 0;
    for (uint32_t n = 0; n < length; n++) {
      output[k] += input[n] * cosf(PI * k * (2 * n + 1) / (2 * length));
    }
//    printf("After DCT[%d]: %0.2f \r\n",k,output[k]);
  }
}

void compute_mfcc(float32_t *input, uint32_t length, float32_t *mfcc) {
  float32_t pre_emphasized[length];
  pre_emphasis(input, pre_emphasized, length, 0.97);

  for (uint32_t i = 0; i < num_frames; i++) {
    memcpy(frame, &pre_emphasized[i * FRAME_SHIFT], FRAME_SIZE * sizeof(float32_t));
    apply_hamming_window(frame, FRAME_SIZE);

    arm_rfft_fast_f32(&fft_instance, frame, fft_output, 0);

    arm_cmplx_mag_squared_f32(fft_output, power_spectrum, FRAME_SIZE/2+1 );
    compute_mel_filter_banks(power_spectrum, mel_energies);

    for (uint32_t k = 0; k < NUM_FILTER_BANKS; k++) {
      log_mel_energies[k] = logf(mel_energies[k] + 1e-19);

      if(isinf(log_mel_energies[i])){
      	log_mel_energies[i] = -43.7f;
      }
    }
    compute_dct(log_mel_energies, dct_output, NUM_FILTER_BANKS);

    for (uint32_t m = 0; m < NUM_MFCC_COEFFS; m++) {
      mfcc_features[i * NUM_MFCC_COEFFS + m] = dct_output[m];
      printf("MFFC Features[%d]: %0.2f\r\n",m,mfcc_features[i * NUM_MFCC_COEFFS + m]);
    }
  }
}
/* USER CODE END 4 */

/**
  * @brief  This function is executed in case of error occurrence.
  * @retval None
  */
void Error_Handler(void)
{
  /* USER CODE BEGIN Error_Handler_Debug */
  /* User can add his own implementation to report the HAL error return state */
  __disable_irq();
  while (1)
  {
  }
  /* USER CODE END Error_Handler_Debug */
}

#ifdef  USE_FULL_ASSERT
/**
  * @brief  Reports the name of the source file and the source line number
  *         where the assert_param error has occurred.
  * @param  file: pointer to the source file name
  * @param  line: assert_param error line source number
  * @retval None
  */
void assert_failed(uint8_t *file, uint32_t line)
{
  /* USER CODE BEGIN 6 */
  /* User can add his own implementation to report the file name and line number,
     ex: printf("Wrong parameters value: file %s on line %d\r\n", file, line) */
  /* USER CODE END 6 */
}
#endif /* USE_FULL_ASSERT */