/** Build with g++ ipp_dtest.cpp -lippcore -lippi -lipps -lippcv -o ipp_dtest **/

#include <iostream>
#include <stdio.h>
#include <stdlib.h>
#include "ippcore.h"
#include "ippcv.h"
#include "ipps.h"
#include "ippi.h"

using namespace std;

int main()
{
int WIDTH = 10; /* image width */
int HEIGHT = 10; /* image height */
IppStatus status = ippStsNoErr;
Ipp8u* pSrc=NULL, *pDst = NULL; /* Pointers to source/destination images */
int srcStep = 0, dstStep = 0; /* Steps, in bytes, through the source/destination images */
IppiSize roiSize = { WIDTH, HEIGHT }; /* Size of source/destination ROI in pixels */
Ipp32u kernelSize = 3;
Ipp32f sigma = 0.56f;
Ipp8u *pBuffer = NULL; /* Pointer to the work buffer */
IppFilterGaussianSpec* pSpec = NULL; /* context structure */
int iTmpBufSize = 0, iSpecSize = 0; /* Common work buffer size */
IppiBorderType borderType = ippBorderRepl;
//IppiBorderType borderType = ippBorderMirrorR;
Ipp8u borderValue = 0;
int numChannels = 1;

Ipp8u src_pixels[HEIGHT][WIDTH] = {{5,7,9,11,13,15,17,19,21,23},
	{9,11,13,15,17,19,21,23,25,27},
	{18,20,22,24,26,28,30,32,34,36},
	{34,36,38,40,42,44,46,48,50,52},
	{59,61,63,65,67,69,71,73,75,77},
	{95,97,99,101,103,105,107,109,111,113},
	{144,146,148,150,152,154,156,158,160,162},
	{208,210,212,214,216,218,220,222,224,226},
	{59,61,63,65,67,69,71,73,75,77},
	{9,11,13,15,17,19,21,23,25,27}};

pSrc=ippiMalloc_8u_C1(roiSize.width, roiSize.height, &srcStep);
pDst = ippiMalloc_8u_C1(roiSize.width, roiSize.height, &dstStep);

/* fill source image with some values*/
for(int i=0;i<HEIGHT; i++)
for(int j=0; j<WIDTH; j++)
*(pSrc + (i * srcStep/sizeof(char)) + j ) = src_pixels[i][j];

status = ippiFilterGaussianGetBufferSize(roiSize, kernelSize, ipp8u,
numChannels, &iSpecSize, &iTmpBufSize);
if(status != ippStsNoErr)
return 0;

pSpec = (IppFilterGaussianSpec *)ippsMalloc_8u(iSpecSize);
pBuffer = ippsMalloc_8u(iTmpBufSize);

status = ippiFilterGaussianInit(roiSize, kernelSize, sigma,
borderType, ipp8u, numChannels, pSpec, pBuffer);
if(status != ippStsNoErr)
return 0;

status = ippiFilterGaussianBorder_8u_C1R(pSrc, srcStep, pDst, dstStep,
roiSize, borderValue, pSpec, pBuffer);
if(status != ippStsNoErr)
return 0;

// Print output
cout<<"\nsrc"<<endl;
cout<<"------"<<endl;
for(int i=0;i<HEIGHT; i++)
{
for(int j=0; j<WIDTH; j++)
{
cout<<" "<<int(*(pSrc + (i * srcStep/sizeof(char)) + j ));
}
cout<<endl;
}
cout<<"\ndst"<<endl;
cout<<"------"<<endl;
for(int i=0;i<HEIGHT; i++)
{
for(int j=0; j<WIDTH; j++)
{
cout<<" "<<int(*(pDst + (i * srcStep/sizeof(char)) + j ));
}
cout<<endl;
}

ippsFree(pBuffer);
ippsFree(pSpec);
ippiFree(pSrc);
ippiFree(pDst);

return (int)status;
}