/******************************************************************************

 *                                                                            *

 * Copyright (C) 2023 The Android Open Source Project

 *

 * Licensed under the Apache License, Version 2.0 (the "License");

 * you may not use this file except in compliance with the License.

 * You may obtain a copy of the License at:

 *

 * http://www.apache.org/licenses/LICENSE-2.0

 *

 * Unless required by applicable law or agreed to in writing, software

 * distributed under the License is distributed on an "AS IS" BASIS,

 * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.

 * See the License for the specific language governing permissions and

 * limitations under the License.

 *

 *****************************************************************************

 * Originally developed and contributed by Ittiam Systems Pvt. Ltd, Bangalore

 */

#include <math.h>

#include "ixheaac_type_def.h"

#include "ixheaace_mps_common_fix.h"

#include "ixheaace_mps_defines.h"

#include "ixheaace_mps_common_define.h"

FLOAT32

ixheaace_mps_212_sum_up_cplx_pow_2_dim_2(ixheaace_cmplx_str inp[MAX_ANA_TIME_SLOT][MAX_QMF_BANDS],

                                         const WORD32 start_dim_1, const WORD32 stop_dim_1,

                                         const WORD32 start_dim_2, const WORD32 stop_dim_2) {

  WORD32 idx_1, idx_2;

  FLOAT32 sum;

  sum = 0.0f;

  for (idx_1 = start_dim_1; idx_1 < stop_dim_1; idx_1++) {

    for (idx_2 = start_dim_2; idx_2 < stop_dim_2; idx_2++) {

      sum += inp[idx_1][idx_2].re * inp[idx_1][idx_2].re;

      sum += inp[idx_1][idx_2].im * inp[idx_1][idx_2].im;

    }

  }

  return (sum / 2);

}

VOID ixheaace_mps_212_cplx_scalar_product(

    ixheaace_cmplx_str *const out, ixheaace_cmplx_str inp_1[MAX_ANA_TIME_SLOT][MAX_QMF_BANDS],

    ixheaace_cmplx_str inp_2[MAX_ANA_TIME_SLOT][MAX_QMF_BANDS], const WORD32 start_dim_1,

    const WORD32 stop_dim_1, const WORD32 start_dim_2, const WORD32 stop_dim_2) {

  WORD32 idx_1, idx_2;

  FLOAT32 re_x, re_y, im_x, im_y, re, im;

  re = 0.0f;

  im = 0.0f;

  for (idx_1 = start_dim_1; idx_1 < stop_dim_1; idx_1++) {

    for (idx_2 = start_dim_2; idx_2 < stop_dim_2; idx_2++) {

      re_x = inp_1[idx_1][idx_2].re;

      im_x = inp_1[idx_1][idx_2].im;

      re_y = inp_2[idx_1][idx_2].re;

      im_y = inp_2[idx_1][idx_2].im;

      re += (re_x * re_y) + (im_x * im_y);

      im += (im_x * re_y) - (re_x * im_y);

    }

  }

  out->re = re / 2;

  out->im = im / 2;

}

FLOAT32 ixheaace_mps_212_sum_up_cplx_pow_2(const ixheaace_cmplx_str *const inp,

                                           const WORD32 len) {

  WORD32 idx;

  FLOAT32 sum;

  sum = 0.0f;

  for (idx = 0; idx < len; idx++) {

    sum += inp[idx].re * inp[idx].re;

    sum += inp[idx].im * inp[idx].im;

  }

  return (sum / 2);

}