/* USER CODE BEGIN Header */
/**
  ******************************************************************************
  * @file           : main.c
  * @brief          : Main program body
  ******************************************************************************
  * @attention
  *
  * Copyright (c) 2024 STecgroelectronics.
  * 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 "arm_math.h"
#include "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 2048

#define AXIS_NUMBER 1

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

#define FRAME_SIZE 1024

#define FRAME_SHIFT 512

#define NUM_MFCC_COEFFS 13

#define NUM_FILTER_BANKS 26
/* USER CODE END PM */

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

UART_HandleTypeDef huart2;

/* USER CODE BEGIN PV */
int n = 0;
uint16_t ecg_x = 0.0;
float ecg_buffer[DATA_USER * AXIS_NUMBER] = {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 * ((FRAME_SIZE - 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];
/* 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);
/* 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);
/* USER CODE END PFP */

/* Private user code ---------------------------------------------------------*/
/* USER CODE BEGIN 0 */
//void HAL_ADC_ConvCpltCallback(ADC_HandleTypeDef *hadc){
//	  ecg_x = HAL_ADC_GetValue(&hadc1);
//	  printf("The converted value is:");
//	  printf("%0.2f", ecg_x);
//	  printf("\r\n");
//	  HAL_ADC_Start(&hadc1);
//}

//void HAL_TIM_PeriodElapsedCallback(TIM_HandleTypeDef *htim){
////	HAL_ADC_Start_IT(&hadc1);
//	HAL_TIM_Base_Start(&htim2);
//}
/* USER CODE END 0 */

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

  /* USER CODE BEGIN 1 */

  /* 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();
  /* USER CODE BEGIN 2 */
  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);
//  arm_dct4_init_f32(&dct_instance, &fft_instance, &S_CFFT, DATA_USER, DATA_USER / 2, 0.5f);
  /* USER CODE END 2 */

  /* Infinite loop */
  /* USER CODE BEGIN WHILE */
  while (1)
  {
	  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;
//	        printf("\r\n");
	        compute_mfcc(ecg_buffer, DATA_USER, mfcc_features);

//
//	        pre_emphasis(ecg_buffer, pre_emphasized, DATA_USER, 0.97);
//
//	        for (uint32_t i = 0; i < num_frames; i++) {
//	        	printf("%d times\r\n",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 );
//	            for(int i =0;i<FRAME_SIZE;i++){
//	            	printf("After Power Spectrum[%d]: %0.2f\r\n",i,power_spectrum[i]);
//	            }
//	            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;
//	              }
//
//	              printf("After Log Mel-Energies[%d]: %0.9f \r\n",k,log_mel_energies[k]);
//	            }
//
//	            compute_dct(log_mel_energies, dct_output, NUM_FILTER_BANKS);

	          }
//	      HAL_Delay(1000);
  }
    /* USER CODE END WHILE */

    /* USER CODE BEGIN 3 */
//	compute_mfcc(ecg_buffer, DATA_USER, mfcc_features);
//	// Transmit or process the MFCC features as needed
//	for (int i = 0; i < NUM_MFCC_COEFFS; i++) {
//		printf("MFCC[%d]: %0.2f\r\n", i, mfcc_features[i]);
//	}
//	HAL_Delay(1000);

  /* 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 = 16;
  RCC_OscInitStruct.PLL.PLLN = 336;
  RCC_OscInitStruct.PLL.PLLP = RCC_PLLP_DIV4;
  RCC_OscInitStruct.PLL.PLLQ = 4;
  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_2) != 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 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 */
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_dct(float32_t *input, float32_t *output, uint32_t length) {
//    // Initialize the twiddle factors and normalization factor
//    float32_t pState[length];
//    float32_t pInlineBuffer[2 * length];
//
//    // Copy input to the buffer
//    memcpy(pInlineBuffer, input, length * sizeof(float32_t));
//
//    // Perform DCT
//    arm_dct4_f32(&dct_instance, pState, pInlineBuffer);
//
//    // Copy the result to output
//    memcpy(output, pInlineBuffer, length * sizeof(float32_t));
//
//    printf("%d the length of the buffer",sizeof(output));
//}



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++) {
	printf("times: %d\r\n",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]);
    }
  }
  printf("\r\n NEW \r\n");
}
/* 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 */