{

public:

virtual void SetUp() {

subMat0.push_back(af_make_seq(0, 4, 1));

subMat0.push_back(af_make_seq(0, 2, 1));

subMat0.push_back(af_make_seq(0, 1, 1));

subMat1.push_back(af_make_seq(1, 4, 1));

subMat1.push_back(af_make_seq(0, 2, 1));

subMat1.push_back(af_make_seq(0, 1, 1));

subMat2.push_back(af_make_seq(1, 4, 1));

subMat2.push_back(af_make_seq(0, 3, 1));

subMat2.push_back(af_make_seq(0, 2, 1));

}

vector<af_seq> subMat0;

vector<af_seq> subMat1;

vector<af_seq> subMat2;

{

if (noDoubleTests<T>()) return;

vector<af::dim4> numDims;

vector<vector<T> > in;

vector<vector<T> > tests;

readTests<T,T,int>(pTestFile,numDims,in,tests);

af::dim4 dims = numDims[0];

T *outData;

af_array inArray = 0;

af_array outArray = 0;

af_array tempArray = 0;

// Get input array

if (isSubRef) {

ASSERT_EQ(AF_SUCCESS, af_create_array(&tempArray, &(in[0].front()), dims.ndims(), dims.get(), (af_dtype) af::dtype_traits<T>::af_type));

ASSERT_EQ(AF_SUCCESS, af_index(&inArray, tempArray, seqv->size(), &seqv->front()));

} else {

ASSERT_EQ(AF_SUCCESS, af_create_array(&inArray, &(in[0].front()), dims.ndims(), dims.get(), (af_dtype) af::dtype_traits<T>::af_type));

}

// Run diff2

ASSERT_EQ(AF_SUCCESS, af_diff2(&outArray, inArray, dim));

// Get result

outData = new T[dims.elements()];

ASSERT_EQ(AF_SUCCESS, af_get_data_ptr((void*)outData, outArray));

// Compare result

for (size_t testIter = 0; testIter < tests.size(); ++testIter) {

vector<T> currGoldBar = tests[testIter];

size_t nElems = currGoldBar.size();

for (size_t elIter = 0; elIter < nElems; ++elIter) {

ASSERT_EQ(currGoldBar[elIter], outData[elIter]) << "at: " << elIter << std::endl;

}

}

// Delete

delete[] outData;

if(inArray != 0) af_release_array(inArray);

if(outArray != 0) af_release_array(outArray);

if(tempArray != 0) af_release_array(tempArray);

{

if (noDoubleTests<T>()) return;

vector<af::dim4> numDims;

vector<vector<T> > in;

vector<vector<T> > tests;

readTests<T,T,int>(pTestFile,numDims,in,tests);

af::dim4 dims = numDims[0];

af_array inArray = 0;

af_array outArray = 0;

ASSERT_EQ(AF_SUCCESS, af_create_array(&inArray, &(in[0].front()), dims.ndims(), dims.get(), (af_dtype) af::dtype_traits<T>::af_type));

ASSERT_EQ(AF_ERR_ARG, af_diff2(&outArray, inArray, -1));

ASSERT_EQ(AF_ERR_ARG, af_diff2(&outArray, inArray, 5));

if(inArray != 0) af_release_array(inArray);

if(outArray != 0) af_release_array(outArray);

{

if (noDoubleTests<float>()) return;

const unsigned dim = 1;

vector<af::dim4> numDims;

vector<vector<float> > in;

vector<vector<float> > tests;

readTests<float,float,int>(string(TEST_DIR"/diff2/matrix1.test"),numDims,in,tests);

af::dim4 dims = numDims[0];

af::array input(dims, &(in[0].front()));

af::array output = af::diff2(input, dim);

float *outData = new float[dims.elements()];

output.host((void*)outData);

// Compare result

for (size_t testIter = 0; testIter < tests.size(); ++testIter) {

vector<float> currGoldBar = tests[testIter];

size_t nElems = currGoldBar.size();

for (size_t elIter = 0; elIter < nElems; ++elIter) {

ASSERT_EQ(currGoldBar[elIter], outData[elIter]) << "at: " << elIter << std::endl;

}

}

// Delete

delete[] outData;