dsptest.c 4.32 KB
Newer Older
1 2
/*
 * MMX optimized DSP utils
3
 * Copyright (c) 2000, 2001, 2002 Fabrice Bellard.
4
 *
5 6 7 8
 * This library is free software; you can redistribute it and/or
 * modify it under the terms of the GNU Lesser General Public
 * License as published by the Free Software Foundation; either
 * version 2 of the License, or (at your option) any later version.
9
 *
10
 * This library is distributed in the hope that it will be useful,
11
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
12 13
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
 * Lesser General Public License for more details.
14
 *
15 16 17
 * You should have received a copy of the GNU Lesser General Public
 * License along with this library; if not, write to the Free Software
 * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
18 19
 */
#define TESTCPU_MAIN
20
#include "avcodec.h"
21
#include "dsputil.h"
22 23 24
#include "mpegvideo.h"
#include "mpeg12data.h"
#include "mpeg4data.h"
25 26
#include "../libavcodec/i386/cputest.c"
#include "../libavcodec/i386/dsputil_mmx.c"
27 28 29 30 31 32 33 34 35

#include "../libavcodec/i386/fdct_mmx.c"
#include "../libavcodec/i386/idct_mmx.c"
#include "../libavcodec/i386/motion_est_mmx.c"
#include "../libavcodec/i386/simple_idct_mmx.c"
#include "../libavcodec/dsputil.c"
#include "../libavcodec/simple_idct.c"
#include "../libavcodec/jfdctfst.c"

36 37 38 39 40 41 42 43 44
#undef TESTCPU_MAIN

#define PAD 0x10000
/*
 * for testing speed of various routine - should be probably extended
 * for a general purpose regression test later
 *
 * currently only for i386 - FIXME
 */
45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69

#define PIX_FUNC_C(a) \
    { #a "_c", a ## _c, 0 }, \
    { #a "_mmx", a ## _mmx, MM_MMX }, \
    { #a "_mmx2", a ## _mmx2, MM_MMXEXT | PAD }

#define PIX_FUNC(a) \
    { #a "_mmx", a ## _mmx, MM_MMX }, \
    { #a "_3dnow", a ## _3dnow, MM_3DNOW }, \
    { #a "_mmx2", a ## _mmx2, MM_MMXEXT | PAD }

#define PIX_FUNC_MMX(a) \
    { #a "_mmx", a ## _mmx, MM_MMX | PAD }

/*
    PIX_FUNC_C(pix_abs16x16),
    PIX_FUNC_C(pix_abs16x16_x2),
    PIX_FUNC_C(pix_abs16x16_y2),
    PIX_FUNC_C(pix_abs16x16_xy2),
    PIX_FUNC_C(pix_abs8x8),
    PIX_FUNC_C(pix_abs8x8_x2),
    PIX_FUNC_C(pix_abs8x8_y2),
    PIX_FUNC_C(pix_abs8x8_xy2),
*/

70 71 72 73 74 75
static const struct pix_func {
    char* name;
    op_pixels_func func;
    int mm_flags;
} pix_func[] = {

76
    PIX_FUNC_MMX(put_pixels),
77 78
    //PIX_FUNC_MMX(get_pixels),
    //PIX_FUNC_MMX(put_pixels_clamped),
Zdenek Kabelac's avatar
Zdenek Kabelac committed
79 80 81
#if 1
    PIX_FUNC(put_pixels_x2),
    PIX_FUNC(put_pixels_y2),
82 83 84 85 86 87 88 89 90 91
    PIX_FUNC_MMX(put_pixels_xy2),

    PIX_FUNC(put_no_rnd_pixels_x2),
    PIX_FUNC(put_no_rnd_pixels_y2),
    PIX_FUNC_MMX(put_no_rnd_pixels_xy2),

    PIX_FUNC(avg_pixels),
    PIX_FUNC(avg_pixels_x2),
    PIX_FUNC(avg_pixels_y2),
    PIX_FUNC(avg_pixels_xy2),
92 93 94 95 96

    PIX_FUNC_MMX(avg_no_rnd_pixels),
    PIX_FUNC_MMX(avg_no_rnd_pixels_x2),
    PIX_FUNC_MMX(avg_no_rnd_pixels_y2),
    PIX_FUNC_MMX(avg_no_rnd_pixels_xy2),
Zdenek Kabelac's avatar
Zdenek Kabelac committed
97
#endif
98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116
    { 0, 0 }
};

static inline long long rdtsc()
{
    long long l;
    asm volatile(   "rdtsc\n\t"
		    : "=A" (l)
		);
    return l;
}

static test_speed(int step)
{
    const struct pix_func* pix = pix_func;
    const int linesize = 720;
    char empty[32768];
    char* bu =(char*)(((long)empty + 32) & ~0xf);

117 118
    int sum = 0;

119 120 121 122 123 124 125
    while (pix->name)
    {
	int i;
        uint64_t te, ts;
        op_pixels_func func = pix->func;
	char* im = bu;

126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142
	if (pix->mm_flags & mm_flags)
	{
	    printf("%30s... ", pix->name);
	    fflush(stdout);
	    ts = rdtsc();
	    for(i=0; i<100000; i++){
		func(im, im + 1000, linesize, 16);
		im += step;
		if (im > bu + 20000)
		    im = bu;
	    }
	    te = rdtsc();
	    emms();
	    printf("% 9d\n", (int)(te - ts));
	    sum += (te - ts) / 100000;
	    if (pix->mm_flags & PAD)
		puts("");
143
	}
144
	pix++;
145
    }
146 147

    printf("Total sum: %d\n", sum);
148 149 150 151 152 153 154 155 156 157 158 159 160 161 162
}

int main(int argc, char* argv[])
{
    int step = 16;

    if (argc > 1)
    {
        // something simple for now
	if (argc > 2 && (strcmp("-s", argv[1]) == 0
			 || strcmp("-step", argv[1]) == 0))
            step = atoi(argv[2]);
    }

    mm_flags = mm_support();
Zdenek Kabelac's avatar
Zdenek Kabelac committed
163
    printf("%s: detected CPU flags:", argv[0]);
164 165 166 167 168 169 170 171 172 173 174 175 176 177 178
    if (mm_flags & MM_MMX)
        printf(" mmx");
    if (mm_flags & MM_MMXEXT)
        printf(" mmxext");
    if (mm_flags & MM_3DNOW)
        printf(" 3dnow");
    if (mm_flags & MM_SSE)
        printf(" sse");
    if (mm_flags & MM_SSE2)
        printf(" sse2");
    printf("\n");

    printf("Using step: %d\n", step);
    test_speed(step);
}