/*******************************************************************************
* 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;
}