#define SIZE 256

/*
 * Feautrier schedule :
 * S1[i, j] -> [1 + j, i]
 * S0[i] -> [0, i]
 */
void BlockBasedFIR(double y[SIZE],double x[SIZE],
        double c[SIZE],int BLOCK, int NTAP){
    int i, j;
#pragma scop BlockBasedFIR, dims = (s,j,i:PAR)
    for (i=0; i<BLOCK; i++) {
#pragma  scop_schedule_statement (0,0,i)
        y[i] = 0;
        for (j=0; j<NTAP; j++) {
            if(i>=j)
#pragma  scop_schedule_statement (1,j,i)
                y[i] = y[i] + x[i-j]*c[j];
        }
    }
}