/*******************************************************************************
* Copyright 2015-2018 Intel Corporation.
*
* This software and the related documents are Intel copyrighted  materials,  and
* your use of  them is  governed by the  express license  under which  they were
* provided to you (License).  Unless the License provides otherwise, you may not
* use, modify, copy, publish, distribute,  disclose or transmit this software or
* the related documents without Intel's prior written permission.
*
* This software and the related documents  are provided as  is,  with no express
* or implied  warranties,  other  than those  that are  expressly stated  in the
* License.
*******************************************************************************/

#include <stdio.h>
#include "ipp.h"

/* Next two defines are created to simplify code reading and understanding */
#define EXIT_MAIN exitLine:                                  /* Label for Exit */
#define check_sts(st) if((st) != ippStsNoErr) goto exitLine; /* Go to Exit if Intel(R) Integrated Primitives (Intel(R) IPP) function returned status different from ippStsNoErr */

/* Results of ippMalloc() are not validated because Intel(R) IPP functions perform bad arguments check and will return an appropriate status  */

int main()
{
    int len = 512;
    Ipp64f pDst[512];
    Ipp64f magn = 4095;
    Ipp64f rLowFreq = 0.2;
    Ipp64f rHighFreq = 0.3;
    int tapslen = 27;
    int numIters = 512;
    IppsFIRSpec_64f* pSpec = NULL;
    IppStatus status = ippStsNoErr;
    int i = 0, bufSize = 0, specSize = 0;
    Ipp8u* pBuffer = NULL;
    IppAlgType algType = ippAlgDirect;
    Ipp64f*  pDlySrc = NULL;
    Ipp64f*  pDlyDst = NULL;

    Ipp64f* FIRDst = ippsMalloc_64f(numIters * sizeof(Ipp64f));
    Ipp64f* taps = ippsMalloc_64f(tapslen*sizeof(Ipp64f));

    pDlySrc = ippsMalloc_64f(tapslen*sizeof(Ipp64f));
    check_sts( status = ippsZero_64f(pDlySrc, tapslen) )

    /* generate source vector */
    check_sts( status = ippsVectorJaehne_64f(pDst, len, magn) ) /* create a Jaehne vector */
    printf("\nSource vector\n");
    for (i = 0; i < 32; i++) printf("%f,", pDst[i]);

    /* computes tapsLen coefficients for lowpass FIR filter */
    check_sts( status = ippsFIRGenGetBufferSize(tapslen, &bufSize) )

    pBuffer     = ippsMalloc_8u(bufSize);

    check_sts( status = ippsFIRGenBandpass_64f(rLowFreq, rHighFreq, taps, tapslen, ippWinBartlett, ippTrue, pBuffer) )

    printf("\nGenerated taps:\n");
    for (i = 0; i < tapslen; i++) printf("%f,", taps[i]);

    check_sts( status = ippsFIRSRGetSize(tapslen, ipp64f, &specSize, &bufSize) )
    pSpec       = (IppsFIRSpec_64f*)ippsMalloc_8u(specSize);
    pBuffer     =                     ippsMalloc_8u(bufSize);

    check_sts( status = ippsFIRSRInit_64f(taps, tapslen, algType, pSpec) )

    /*filter an input vector*/
    check_sts( status = ippsFIRSR_64f(pDst, FIRDst, numIters, pSpec, pDlySrc, pDlyDst, pBuffer) )

    printf("\nAfter LowPass filter\n");
    for (i = 0; i < 32; i++) printf("%f,", FIRDst[i]);
    printf("\n\n");

EXIT_MAIN
    ippsFree(FIRDst);
    ippsFree(taps);
    printf("Exit status %d (%s)\n", (int)status, ippGetStatusString(status));
    return 0;
}