/normxcorr/trunk

#include <stdio.h>

#include <stdlib.h>

 

#include <mex.h>

 

#include "local_sum.h"

 

#include "local_sum_kernel.cu"

 

#include <unistd.h>

 

int local_sum(TProcessingState *ps, 

 

int local_sum_validate(TProcessingState *ps, int icp, const mxArray *lsum, const mxArray *denom) {

    int res = 0;

#ifdef VALIDATE_LSUM

    int size = ps->fft_size;

    int size2 = size*size;

    int alloc_size = ps->fft_alloc_size;

 

    if (lsum&&denom) {

        float *tmp = (float*)malloc(size2*sizeof(float));

        

        cudaMemcpy(tmp, ps->cuda_lsum_cache + icp * alloc_size, size2*sizeof(float), cudaMemcpyDeviceToHost);

        float *real =  (float*)mxGetData(lsum);

        if (memcmp(tmp, real, size2*sizeof(float))) {

            printf("lsum fault: %i\n", 1);

            for (int i = 0; i < size2; i++) {

                if (tmp[i] != real[i]) {

                    res = 1;

                    printf("lsum fault: %i %i - %f %f\n", i / size, i % size, tmp[i], real[i]);

                    break;

                }

            }

        }

 

        cudaMemcpy(tmp, ps->cuda_denom_cache + icp * alloc_size, size2*sizeof(float), cudaMemcpyDeviceToHost);

        real =  (float*)mxGetData(denom);

        if (memcmp(tmp, real, size2*sizeof(float))) {

            for (int i = 0; i < size2; i++) {

                float diff = (tmp[i] == real[i])?0:2*fabs(tmp[i] - real[i])/(tmp[i] = real[i]);

                if (diff > 0.0001) {

                    res = 1;

                    printf("denom fault: %i %i - %f %f %f\n", i / size, i % size, tmp[i], real[i], diff);

                    break;

                }

            }

        }

 

        free(tmp);

    }

#endif /* VALIDATE_LSUM */

 

    return res;

}