~ubuntu-branches/ubuntu/utopic/ffmpeg-debian/utopic

#undef PMINUB

#undef PMAXUB

 

#define REAL_PAVGB(a,b) "pavgb " #a ", " #b " \n\t"

#define REAL_PAVGB(a,b) "pavgusb " #a ", " #b " \n\t"

#endif

#define PAVGB(a,b)  REAL_PAVGB(a,b)

 

#define PMINUB(a,b,t) "pminub " #a ", " #b " \n\t"

#define PMINUB(b,a,t) \

    "movq " #a ", " #t " \n\t"\

    "psubusb " #b ", " #t " \n\t"\

    "psubb " #t ", " #a " \n\t"

#endif

 

#define PMAXUB(a,b) "pmaxub " #a ", " #b " \n\t"

#define PMAXUB(a,b) \

    "psubusb " #a ", " #b " \n\t"\

    "paddb " #a ", " #b " \n\t"

#endif

 

//FIXME? |255-0| = 1 (should not be a problem ...)

/**

 * Check if the middle 8x8 Block in the given 8x16 block is flat

 */

        "psubusb %%mm3, %%mm4                   \n\t"

 

        "                                       \n\t"

        "pxor %%mm7, %%mm7                      \n\t"

        "psadbw %%mm7, %%mm0                    \n\t"

#else

 * Do a vertical low pass filter on the 8x16 block (only write to the 8x8 block in the middle)

 * using the 9-Tap Filter (1,1,2,2,4,2,2,1,1)/16

 */

static inline void RENAME(doVertLowPass)(uint8_t *src, int stride, PPContext *c)

{

    src+= stride*3;

    __asm__ volatile(        //"movv %0 %1 %2\n\t"

        "movq %2, %%mm0                         \n\t"  // QP,..., QP

        : "r" (src), "r" ((x86_reg)stride), "m" (c->pQPb)

        : "%"REG_a, "%"REG_c

    );

    const int l1= stride;

    const int l2= stride + l1;

    const int l3= stride + l2;

 

        src++;

    }

}

#endif //HAVE_ALTIVEC

 

 */

static inline void RENAME(vertRK1Filter)(uint8_t *src, int stride, int QP)

{

    src+= stride*3;

// FIXME rounding

    __asm__ volatile(

        : "r" (src), "r" ((x86_reg)stride)

        : "%"REG_a, "%"REG_c

    );

    const int l1= stride;

    const int l2= stride + l1;

    const int l3= stride + l2;

        }

    }

 

}

#endif //0

 

 */

static inline void RENAME(vertX1Filter)(uint8_t *src, int stride, PPContext *co)

{

    src+= stride*3;

 

    __asm__ volatile(

        : "r" (src), "r" ((x86_reg)stride), "m" (co->pQPb)

        : "%"REG_a, "%"REG_c

    );

 

    const int l1= stride;

    const int l2= stride + l1;

        }

        src++;

    }

}

 

static inline void RENAME(doVertDefFilter)(uint8_t src[], int stride, PPContext *c)

{

/*

    uint8_t tmp[16];

    const int l1= stride;

    }

}

*/

    src+= stride*4;

    __asm__ volatile(

        "pxor %%mm7, %%mm7                      \n\t"

        "movq (%%"REG_c"), %%mm2                \n\t" // 2L0 - 5L1 + 5L2 - 2L3

        "movq 8(%%"REG_c"), %%mm3               \n\t" // 2H0 - 5H1 + 5H2 - 2H3

 

        "movq %%mm7, %%mm6                      \n\t" // 0

        "psubw %%mm0, %%mm6                     \n\t"

        "pmaxsw %%mm6, %%mm0                    \n\t" // |2L4 - 5L5 + 5L6 - 2L7|

        "psubw %%mm6, %%mm3                     \n\t" // |2H0 - 5H1 + 5H2 - 2H3|

#endif

 

        "pminsw %%mm2, %%mm0                    \n\t"

        "pminsw %%mm3, %%mm1                    \n\t"

#else

        "pand %%mm2, %%mm4                      \n\t"

        "pand %%mm3, %%mm5                      \n\t"

 

        "pminsw %%mm0, %%mm4                    \n\t"

        "pminsw %%mm1, %%mm5                    \n\t"

#else

        : "r" ((x86_reg)stride), "m" (c->pQPb)

        : "%"REG_a, "%"REG_c

    );

    const int l1= stride;

    const int l2= stride + l1;

    const int l3= stride + l2;

        }

        src++;

    }

}

#endif //HAVE_ALTIVEC

 

static inline void RENAME(dering)(uint8_t src[], int stride, PPContext *c)

{

    __asm__ volatile(

        "pxor %%mm6, %%mm6                      \n\t"

        "pcmpeqb %%mm7, %%mm7                   \n\t"

//        %0        eax        eax+%1        eax+2%1        %0+4%1        edx        edx+%1        edx+2%1        %0+8%1        edx+4%1

 

#undef FIND_MIN_MAX

#define REAL_FIND_MIN_MAX(addr)\

        "movq " #addr ", %%mm0                  \n\t"\

        "pminub %%mm0, %%mm7                    \n\t"\

 

        "movq %%mm7, %%mm4                      \n\t"

        "psrlq $8, %%mm7                        \n\t"

        "pminub %%mm4, %%mm7                    \n\t" // min of pixels

        "pshufw $0xF9, %%mm7, %%mm4             \n\t"

        "pminub %%mm4, %%mm7                    \n\t" // min of pixels

 

        "movq %%mm6, %%mm4                      \n\t"

        "psrlq $8, %%mm6                        \n\t"

        "pmaxub %%mm4, %%mm6                    \n\t" // max of pixels

        "pshufw $0xF9, %%mm6, %%mm4             \n\t"

        "pmaxub %%mm4, %%mm6                    \n\t"

        : : "r" (src), "r" ((x86_reg)stride), "m" (c->pQPb), "m"(c->pQPb2)

        : "%"REG_a, "%"REG_d, "%"REG_c

    );

    int y;

    int min=255;

    int max=0;

//        src[0] = src[7]=src[stride*7]=src[stride*7 + 7]=255;

    }

#endif

}

#endif //HAVE_ALTIVEC

 

 */

static inline void RENAME(deInterlaceInterpolateLinear)(uint8_t src[], int stride)

{

    src+= 4*stride;

    __asm__ volatile(

        "lea (%0, %1), %%"REG_a"                \n\t"

 */

static inline void RENAME(deInterlaceInterpolateCubic)(uint8_t src[], int stride)

{

    src+= stride*3;

    __asm__ volatile(

        "lea (%0, %1), %%"REG_a"                \n\t"

        : : "r" (src), "r" ((x86_reg)stride)

        : "%"REG_a, "%"REG_d, "%"REG_c

    );

    int x;

    src+= stride*3;

    for(x=0; x<8; x++){

        src[stride*9] = CLIP((-src[stride*6] + 9*src[stride*8] + 9*src[stride*10] - src[stride*12])>>4);

        src++;

    }

}

 

/**

 */

static inline void RENAME(deInterlaceFF)(uint8_t src[], int stride, uint8_t *tmp)

{

    src+= stride*4;

    __asm__ volatile(

        "lea (%0, %1), %%"REG_a"                \n\t"

        : : "r" (src), "r" ((x86_reg)stride), "r"(tmp)

        : "%"REG_a, "%"REG_d

    );

    int x;

    src+= stride*4;

    for(x=0; x<8; x++){

 

        src++;

    }

}

 

/**

 */

static inline void RENAME(deInterlaceL5)(uint8_t src[], int stride, uint8_t *tmp, uint8_t *tmp2)

{

    src+= stride*4;

    __asm__ volatile(

        "lea (%0, %1), %%"REG_a"                \n\t"

        : : "r" (src), "r" ((x86_reg)stride), "r"(tmp), "r"(tmp2)

        : "%"REG_a, "%"REG_d

    );

    int x;

    src+= stride*4;

    for(x=0; x<8; x++){

 

        src++;

    }

}

 

/**

 */

static inline void RENAME(deInterlaceBlendLinear)(uint8_t src[], int stride, uint8_t *tmp)

{

    src+= 4*stride;

    __asm__ volatile(

        "lea (%0, %1), %%"REG_a"                \n\t"

        : : "r" (src), "r" ((x86_reg)stride), "r" (tmp)

        : "%"REG_a, "%"REG_d

    );

    int a, b, c, x;

    src+= 4*stride;

 

        src += 4;

        tmp += 4;

    }

}

 

/**

 */

static inline void RENAME(deInterlaceMedian)(uint8_t src[], int stride)

{

    src+= 4*stride;

    __asm__ volatile(

        "lea (%0, %1), %%"REG_a"                \n\t"

        "lea (%%"REG_a", %1, 4), %%"REG_d"      \n\t"

#endif //HAVE_MMX

}

 

/**

 * transposes and shift the given 8x8 Block into dst1 and dst2

 */

#endif //HAVE_MMX

//static long test=0;

 

static inline void RENAME(tempNoiseReducer)(uint8_t *src, int stride,

                                    uint8_t *tempBlurred, uint32_t *tempBlurredPast, int *maxNoise)

{

 

#define FAST_L2_DIFF

//#define L1_DIFF //u should change the thresholds too if u try that one

    __asm__ volatile(

        "lea (%2, %2, 2), %%"REG_a"             \n\t" // 3*stride

        "lea (%2, %2, 4), %%"REG_d"             \n\t" // 5*stride

        :: "r" (src), "r" (tempBlurred), "r"((x86_reg)stride), "m" (tempBlurredPast)

        : "%"REG_a, "%"REG_d, "%"REG_c, "memory"

    );

{

    int y;

    int d=0;

        }

    }

}

}

#endif //HAVE_ALTIVEC

 

/**

 * accurate deblock filter

 */

            "movq (%%"REG_c"), %%mm2                \n\t" // 2L0 - 5L1 + 5L2 - 2L3

            "movq 8(%%"REG_c"), %%mm3               \n\t" // 2H0 - 5H1 + 5H2 - 2H3

 

            "movq %%mm7, %%mm6                      \n\t" // 0

            "psubw %%mm0, %%mm6                     \n\t"

            "pmaxsw %%mm6, %%mm0                    \n\t" // |2L4 - 5L5 + 5L6 - 2L7|

            "psubw %%mm6, %%mm3                     \n\t" // |2H0 - 5H1 + 5H2 - 2H3|

#endif

 

            "pminsw %%mm2, %%mm0                    \n\t"

            "pminsw %%mm3, %%mm1                    \n\t"

#else

            "pand %%mm2, %%mm4                      \n\t"

            "pand %%mm3, %%mm5                      \n\t"

 

            "pminsw %%mm0, %%mm4                    \n\t"

            "pminsw %%mm1, %%mm5                    \n\t"

#else

static inline void RENAME(blockCopy)(uint8_t dst[], int dstStride, const uint8_t src[], int srcStride,

                                     int levelFix, int64_t *packedOffsetAndScale)

{

    int i;

#endif

    if(levelFix){

    __asm__ volatile(

        "movq (%%"REG_a"), %%mm2        \n\t" // packedYOffset

        "movq 8(%%"REG_a"), %%mm3       \n\t" // packedYScale

        "lea (%2,%4), %%"REG_a"         \n\t"

        "lea (%3,%5), %%"REG_d"         \n\t"

        "pxor %%mm4, %%mm4              \n\t"

#define REAL_SCALED_CPY(src1, src2, dst1, dst2)                                                \

        "movq " #src1 ", %%mm0          \n\t"\

        "movq " #src1 ", %%mm5          \n\t"\

                &(src[srcStride*i]), BLOCK_SIZE);

#endif //HAVE_MMX

    }else{

    __asm__ volatile(

        "lea (%0,%2), %%"REG_a"                 \n\t"

        "lea (%1,%3), %%"REG_d"                 \n\t"

 */

static inline void RENAME(duplicate)(uint8_t src[], int stride)

{

    __asm__ volatile(

        "movq (%0), %%mm0               \n\t"

        "add %1, %0                     \n\t"

    int QPCorrecture= 256*256;

 

    int copyAhead;

    int i;

#endif

 

    uint8_t * const tempDst= dstStride > 0 ? c.tempDst : c.tempDst - 23*dstStride;

    //const int mbWidth= isColor ? (width+7)>>3 : (width+15)>>4;

 

    for(i=0; i<57; i++){

        int offset= ((i*c.ppMode.baseDcDiff)>>8) + 1;

        int threshold= offset*2 + 1;

 

        scale= (double)(c.ppMode.maxAllowedY - c.ppMode.minAllowedY) / (double)(white-black);

 

        c.packedYScale= (uint16_t)(scale*256.0 + 0.5);

        c.packedYOffset= (((black*c.packedYScale)>>8) - c.ppMode.minAllowedY) & 0xFFFF;

#else

        // with the L1 Cache of the P4 ... or only a few blocks at a time or soemthing

        for(x=0; x<width; x+=BLOCK_SIZE){

 

/*

            prefetchnta(srcBlock + (((x>>2)&6) + 5)*srcStride + 32);

            prefetchnta(srcBlock + (((x>>2)&6) + 6)*srcStride + 32);

                : "%"REG_a, "%"REG_d

            );

 

//FIXME check if this is faster on an 3dnow chip or if it is faster without the prefetch or ...

/*          prefetch(srcBlock + (((x>>3)&3) + 5)*srcStride + 32);

            prefetch(srcBlock + (((x>>3)&3) + 9)*srcStride + 32);

        //1% speedup if these are here instead of the inner loop

        const uint8_t *srcBlock= &(src[y*srcStride]);

        uint8_t *dstBlock= &(dst[y*dstStride]);

        uint8_t *tempBlock1= c.tempBlocks;

        uint8_t *tempBlock2= c.tempBlocks + 8;

#endif

        // with the L1 Cache of the P4 ... or only a few blocks at a time or soemthing

        for(x=0; x<width; x+=BLOCK_SIZE){

            const int stride= dstStride;

            uint8_t *tmpXchg;

#endif

            if(isColor){

                yHistogram[ srcBlock[srcStride*12 + 4] ]++;

            }

            c.QP= QP;

            __asm__ volatile(

                "movd %1, %%mm7         \n\t"

                "packuswb %%mm7, %%mm7  \n\t" // 0, 0, 0, QP, 0, 0, 0, QP

#endif

 

 

/*

            prefetchnta(srcBlock + (((x>>2)&6) + 5)*srcStride + 32);

            prefetchnta(srcBlock + (((x>>2)&6) + 6)*srcStride + 32);

                : "%"REG_a, "%"REG_d

            );

 

//FIXME check if this is faster on an 3dnow chip or if it is faster without the prefetch or ...

/*          prefetch(srcBlock + (((x>>3)&3) + 5)*srcStride + 32);

            prefetch(srcBlock + (((x>>3)&3) + 9)*srcStride + 32);

                }

            }

 

            RENAME(transpose1)(tempBlock1, tempBlock2, dstBlock, dstStride);

#endif

            /* check if we have a previous block to deblock it with dstBlock */

            if(x - 8 >= 0){

                if(mode & H_X1_FILTER)

                        RENAME(vertX1Filter)(tempBlock1, 16, &c);

                else if(mode & H_DEBLOCK){

                if(mode & H_X1_FILTER)

                    horizX1Filter(dstBlock-4, stride, QP);

                else if(mode & H_DEBLOCK){

                    DECLARE_ALIGNED(16, unsigned char, tempBlock[272]);

                    transpose_16x8_char_toPackedAlign_altivec(tempBlock, dstBlock - (4 + 1), stride);

 

            dstBlock+=8;

            srcBlock+=8;

 

            tmpXchg= tempBlock1;

            tempBlock1= tempBlock2;

            tempBlock2 = tmpXchg;

                + dstBlock[x +14*dstStride] + dstBlock[x +15*dstStride];

        }*/

    }

    __asm__ volatile("femms");

    __asm__ volatile("emms");

#endif