/ani/mrses

#ifndef _MRSES_IMPL_H

#define _MRSES_IMPL_H

#include <stdlib.h>

#include <math.h>

#include "mrses.h"

#include "hw_sched.h"

#include "tools.h"

#define HW_ALIGN 16		// For SPE support minimum 16 is required (should be multiple of MRSESDataType anyway)

#define SPE_BLOCK 16		// This indicates minimal block of processing by SPE code

#define SIMD_BLOCK 16		// In bytes

#define HW_ITERATE_BLOCKS 16

#define SPE_MAX_TRANSFER_SIZE 16384

//#define TRACE_TIMINGS 1

//#define USE_FAST_RANDOM 1		//* default random/rand somehow sequentalize threads

//#define FIX_RANDOM 1

//#define HW_USE_BLOCKED_MULTIPLY	//* Alternative outdated way of processing in SPU (don't use)

#define blas_asum cblas_sasum

#define blas_scal cblas_sscal

#define blas_axpy cblas_saxpy

#define blas_gemv cblas_sgemv

#define blas_gemm cblas_sgemm

#define blas_syrk cblas_ssyrk

#define blas_trsm cblas_strsm

#define blas_trsv cblas_strsv

#define blas_dot  cblas_sdot

#define lapack_potrf spotrf_

struct MRSESContextT {

    MRSESDistance dist;

    int palloc;			/**< Size of vectors in index array, iterate only */

    int properties;		/**< Number of properties */

    int alloc;			/**< Size of vectors in A and B matrixes for alignment */

    int nA, nB;			/**< Number of samples */

    MRSESDataType *A;		/**< Actually: (A - mean(A))/sqrt(nA) */

    MRSESDataType *B;		/**< Actually: (B - mean(B))/sqrt(nB) */

    MRSESDataType *mean;	/**< Actually: mean(A) - mean(B) */

    HWSched sched;

    int iterate_size;		/**< Number of blocks per block group */

    int max_block_size;

    int block_size;

    int cur_chunk;

    int width; 

    int iterations;   

    MRSESIntType *index;	/**< In compute: provided by client app, no allocation, in iteration: allocated, indicates that iterations were started before */

    MRSESIntType *ires;		/**< Array of resulting indeces */

    MRSESDataType *result;	/**< Provided by client app, no allocation */

};

typedef struct MRSESContextT MRSESContextS;

#endif /* _MRSES_IMPL_H */