BLIT schiffler_binarize(void *src, void *dst, int bytes, LinkList<ParameterInstance> ¶ms) {
unsigned char v;
{
v = (unsigned int) *(float*)(params.front()->get()); // only one value
}
SDL_imageFilterBinarizeUsingThreshold
((unsigned char*)src, (unsigned char*)dst, bytes, v);
}
BLIT blue_mask(void *src, void *dst, int bytes, LinkList<ParameterInstance> ¶ms) {
register int c;
register uint32_t *s = (uint32_t*)src;
register uint32_t *d = (uint32_t*)dst;
for(c = bytes >> 2; c > 0; c--, s++, d++)
*d |= *s & blue_bitmask;
unsigned char v;
{
v = (unsigned int) *(float*)(params.front()->get()); // only one value
}
SDL_imageFilterBinarizeUsingThreshold((unsigned char*)dst, (unsigned char*)dst, bytes, v);
}
void Draw (SDL_Surface *screen, char *bmpfile)
{
char messageText[128];
Uint32 rmask, gmask, bmask, amask;
SDL_Surface *picture;
SDL_Surface *mapped_picture;
SDL_Surface *rotozoom_picture;
SDL_PixelFormat *pixelFormat;
Uint8 *grayscale, *map, *curmap;
double *unrelaxed, *relaxed, *currelax;
int mapsize, relaxsize;
int rowskip;
char *pixel;
Uint32 p;
int x, y;
Uint8 r, g, b, a;
double dy;
double t, s;
int i;
double c1, c2, c3, c4, ca;
Uint8 lookupTable[256];
double zw, zh, zf;
/* SDL interprets each pixel as a 32-bit number, so our masks must depend
on the endianness (byte order) of the machine */
#if SDL_BYTEORDER == SDL_BIG_ENDIAN
rmask = 0xff000000;
gmask = 0x00ff0000;
bmask = 0x0000ff00;
amask = 0x000000ff;
#else
rmask = 0x000000ff;
gmask = 0x0000ff00;
bmask = 0x00ff0000;
amask = 0xff000000;
#endif
/* Load the image into a surface */
printf("Loading picture: %s\n", bmpfile);
picture = SDL_LoadBMP(bmpfile);
if ( picture == NULL ) {
fprintf(stderr, "Couldn't load %s: %s\n", bmpfile, SDL_GetError());
return;
}
/* Convert image to a brightness map */
printf("Allocating workbuffers\n");
mapsize = picture->w * picture->h;
relaxsize = mapsize * sizeof(double);
grayscale = (Uint8 *)malloc(mapsize);
map = (Uint8 *)malloc(mapsize);
unrelaxed = (double *)malloc(relaxsize);
relaxed = (double *)malloc(relaxsize);
if ((grayscale == NULL) || (map == NULL) || (unrelaxed == NULL) || (relaxed == NULL))
{
fprintf(stderr, "Out of memory\n");
return;
}
memset(grayscale, 0, mapsize);
memset(map, 0, mapsize);
memset(unrelaxed, 0, relaxsize);
memset(relaxed, 0, relaxsize);
printf("Converting image to brightness map\n");
pixel = picture->pixels;
pixelFormat = picture->format;
rowskip = (picture->pitch - picture->w * picture->format->BytesPerPixel);
curmap = grayscale;
for (y=0; y < picture->h; y++) {
for (x=0; x < picture->w; x++) {
// Get RGB
p = getPixel(picture, x, y);
SDL_GetRGBA(p, pixelFormat, &r, &g, &b, &a);
// Calculate luminance (Y = 0.3R + 0.59G + 0.11B;)
dy = (0.30 * r) + (0.59 * g) + (0.11 * b);
if (dy<0.0) {
dy=0.0;
} else if (dy>255.0) {
dy=255.0;
}
*curmap = (Uint8)dy;
// Next pixel
pixel += picture->format->BytesPerPixel;
curmap++;
}
// Next row
pixel += rowskip;
}
/* --- Prepare relaxation loop --- */
/* Binarize luminance map */
SDL_imageFilterBinarizeUsingThreshold(grayscale, map, mapsize, threshold);
/* Create cos^5 based lookup table */
t = 0.0;
for (i = 0; i < 256; i++)
{
s = 1.0 - 0.5 * (1.0 + cos(t));
s = 255.0 * s * s * s * s * s;
lookupTable[i] = (Uint8)s;
t += ((double)contours*2.0*3.141592653589/128.0);
}
/* Create new surface for relaxed image */
mapped_picture = SDL_CreateRGBSurface(SDL_SWSURFACE, picture->w, picture->h, 32,
rmask, gmask, bmask, amask);
if (mapped_picture == NULL) {
fprintf(stderr, "CreateRGBSurface failed: %s\n", SDL_GetError());
return;
}
/* Apply relaxation algorithm */
printf("Applying Laplacian relaxation: %i iterations\n", iterations);
// Iterate to relax
for (i = 0; i <= iterations; i++)
{
// Backup original
memcpy(unrelaxed, relaxed, relaxsize);
// Average into relaxed
for (x=0; x < picture->w; x++) {
for (y=0; y < picture->h; y++) {
// up
c1 = GetPotential(map, grayscale, unrelaxed, x, y-1, picture->w, picture->h);
// down
c2 = GetPotential(map, grayscale, unrelaxed, x, y+1, picture->w, picture->h);
// left
c3 = GetPotential(map, grayscale, unrelaxed, x-1, y, picture->w, picture->h);
// right
c4 = GetPotential(map, grayscale, unrelaxed, x+1, y, picture->w, picture->h);
// average and store
ca = ((c1 + c2 + c3 + c4)/4.0);
relaxed[x+y*picture->w] = ca;
}
}
// Draw only sometimes
if (((i % 10)==0) || (i==iterations)) {
/* --- Create image with contour map --- */
/* Fill output surface via colormap */
currelax = relaxed;
for (y=0; y<mapped_picture->h; y++) {
for (x=0; x<mapped_picture->w; x++) {
if (map[x+y*picture->w]!=0) {
r = g = b = grayscale[x+y*picture->w];
} else {
r = g = b = lookupTable[(Uint8)*currelax];
}
pixelRGBA(mapped_picture, x, y, r, g, b, 255);
currelax++;
}
}
/* --- Scale and draw to screen --- */
/* Scale to screen size */
zw = (double)screen->w/(double)mapped_picture->w;
zh = (double)screen->h/(double)mapped_picture->h;
zf = (zw < zh) ? zw : zh;
if ((rotozoom_picture=zoomSurface(mapped_picture, zf, zf, 1))==NULL) {
fprintf(stderr, "Rotozoom failed: %s\n", SDL_GetError());
return;
}
/* Draw surface to screen */
if ( SDL_BlitSurface(rotozoom_picture, NULL, screen, NULL) < 0 ) {
fprintf(stderr, "Blit failed: %s\n", SDL_GetError());
return;
}
SDL_FreeSurface(rotozoom_picture);
/* Info */
if (i != iterations) {
sprintf(messageText,"%i", i);
stringRGBA(screen, 8, 8, messageText, 255, 255, 255, 255);
}
/* Display by flipping screens */
SDL_Flip(screen);
}
/* Maybe quit */
HandleEvent();
}
/* Save final picture */
if (SDL_SaveBMP(mapped_picture, "result.bmp") <0) {
fprintf(stderr, "Save BMP failed: %s\n", SDL_GetError());
return;
}
free(map);
free(grayscale);
free(unrelaxed);
free(relaxed);
SDL_FreeSurface(picture);
SDL_FreeSurface(mapped_picture);
return;
}
int main(int argc, char *argv[])
{
unsigned char src1[15],src2[15],dstm[15],dstc[15];
/* SDL_imageFilter Test */
printf ("TestImageFilter\n\n");
printf ("Testing an array of 15 bytes - first 8 bytes should be processed\n");
printf ("by MMX or C code, the last 7 bytes only by C code.\n\n");
print_line();
SDL_imageFilterMMXon();
setup_src(src1, src2);
SDL_imageFilterBitAnd ((unsigned char *)src1,(unsigned char *)src2,(unsigned char *)dstm,15);
print_result ("MMX BitAnd", src1, src2, dstm);
SDL_imageFilterMMXoff();
setup_src(src1, src2);
SDL_imageFilterBitAnd ((unsigned char *)src1,(unsigned char *)src2,(unsigned char *)dstc,15);
print_result (" C BitAnd", src1, src2, dstc);
print_compare(dstm,dstc);
print_line();
SDL_imageFilterMMXon();
setup_src(src1, src2);
SDL_imageFilterBitOr ((unsigned char *)src1,(unsigned char *)src2,(unsigned char *)dstm,15);
print_result ("MMX BitOr", src1, src2, dstm);
SDL_imageFilterMMXoff();
setup_src(src1, src2);
SDL_imageFilterBitOr ((unsigned char *)src1,(unsigned char *)src2,(unsigned char *)dstc,15);
print_result (" C BitOr", src1, src2, dstc);
print_compare(dstm,dstc);
print_line();
SDL_imageFilterMMXon();
setup_src(src1, src2);
SDL_imageFilterAdd ((unsigned char *)src1,(unsigned char *)src2,(unsigned char *)dstm,15);
print_result ("MMX Add", src1, src2, dstm);
SDL_imageFilterMMXoff();
setup_src(src1, src2);
SDL_imageFilterAdd ((unsigned char *)src1,(unsigned char *)src2,(unsigned char *)dstc,15);
print_result (" C Add", src1, src2, dstc);
print_compare(dstm,dstc);
print_line();
SDL_imageFilterMMXon();
setup_src(src1, src2);
SDL_imageFilterAbsDiff ((unsigned char *)src1,(unsigned char *)src2,(unsigned char *)dstm,15);
print_result ("MMX AbsDiff", src1, src2, dstm);
SDL_imageFilterMMXoff();
setup_src(src1, src2);
SDL_imageFilterAbsDiff ((unsigned char *)src1,(unsigned char *)src2,(unsigned char *)dstc,15);
print_result (" C AbsDiff", src1, src2, dstc);
print_compare(dstm,dstc);
print_line();
SDL_imageFilterMMXon();
setup_src(src1, src2);
SDL_imageFilterMean ((unsigned char *)src1,(unsigned char *)src2,(unsigned char *)dstm,15);
print_result ("MMX Mean", src1, src2, dstm);
SDL_imageFilterMMXoff();
setup_src(src1, src2);
SDL_imageFilterMean ((unsigned char *)src1,(unsigned char *)src2,(unsigned char *)dstc,15);
print_result (" C Mean", src1, src2, dstc);
print_compare(dstm,dstc);
print_line();
SDL_imageFilterMMXon();
setup_src(src1, src2);
SDL_imageFilterSub ((unsigned char *)src1,(unsigned char *)src2,(unsigned char *)dstm,15);
print_result ("MMX Sub", src1, src2, dstm);
SDL_imageFilterMMXoff();
setup_src(src1, src2);
SDL_imageFilterSub ((unsigned char *)src1,(unsigned char *)src2,(unsigned char *)dstc,15);
print_result (" C Sub", src1, src2, dstc);
print_compare(dstm,dstc);
print_line();
SDL_imageFilterMMXon();
setup_src(src1, src2);
SDL_imageFilterMult ((unsigned char *)src1,(unsigned char *)src2,(unsigned char *)dstm,15);
print_result ("MMX Mult", src1, src2, dstm);
SDL_imageFilterMMXoff();
setup_src(src1, src2);
SDL_imageFilterMult ((unsigned char *)src1,(unsigned char *)src2,(unsigned char *)dstc,15);
print_result (" C Mult", src1, src2, dstc);
print_compare(dstm,dstc);
print_line();
SDL_imageFilterMMXon();
setup_src(src1, src2);
SDL_imageFilterMultNor ((unsigned char *)src1,(unsigned char *)src2,(unsigned char *)dstm,15);
print_result ("ASM MultNor", src1, src2, dstm);
SDL_imageFilterMMXoff();
setup_src(src1, src2);
SDL_imageFilterMultNor ((unsigned char *)src1,(unsigned char *)src2,(unsigned char *)dstc,15);
print_result (" C MultNor", src1, src2, dstc);
print_compare(dstm,dstc);
print_line();
SDL_imageFilterMMXon();
setup_src(src1, src2);
SDL_imageFilterMultDivby2 ((unsigned char *)src1,(unsigned char *)src2,(unsigned char *)dstm,15);
print_result ("MMX MultDivby2", src1, src2, dstm);
SDL_imageFilterMMXoff();
setup_src(src1, src2);
SDL_imageFilterMultDivby2 ((unsigned char *)src1,(unsigned char *)src2,(unsigned char *)dstc,15);
print_result (" C MultDivby2", src1, src2, dstc);
print_compare(dstm,dstc);
print_line();
SDL_imageFilterMMXon();
setup_src(src1, src2);
SDL_imageFilterMultDivby4 ((unsigned char *)src1,(unsigned char *)src2,(unsigned char *)dstm,15);
print_result ("MMX MultDivby4", src1, src2, dstm);
SDL_imageFilterMMXoff();
setup_src(src1, src2);
SDL_imageFilterMultDivby4 ((unsigned char *)src1,(unsigned char *)src2,(unsigned char *)dstc,15);
print_result (" C MultDivby4", src1, src2, dstc);
print_compare(dstm,dstc);
print_line();
SDL_imageFilterMMXon();
setup_src(src1, src2);
SDL_imageFilterDiv ((unsigned char *)src1,(unsigned char *)src2,(unsigned char *)dstm,15);
print_result ("ASM Div", src1, src2, dstm);
SDL_imageFilterMMXoff();
setup_src(src1, src2);
SDL_imageFilterDiv ((unsigned char *)src1,(unsigned char *)src2,(unsigned char *)dstc,15);
print_result (" C Div", src1, src2, dstc);
print_compare(dstm,dstc);
print_line();
SDL_imageFilterMMXon();
setup_src(src1, src2);
SDL_imageFilterBitNegation ((unsigned char *)src1,(unsigned char *)dstm,15);
print_result ("MMX BitNegation", src1, NULL, dstm);
SDL_imageFilterMMXoff();
setup_src(src1, src2);
SDL_imageFilterBitNegation ((unsigned char *)src1,(unsigned char *)dstc,15);
print_result (" C BitNegation", src1, NULL, dstc);
print_compare(dstm,dstc);
print_line();
SDL_imageFilterMMXon();
setup_src(src1, src2);
SDL_imageFilterAddByte ((unsigned char *)src1,(unsigned char *)dstm,15, 3);
print_result ("MMX AddByte(3)", src1, NULL, dstm);
SDL_imageFilterMMXoff();
setup_src(src1, src2);
SDL_imageFilterAddByte ((unsigned char *)src1,(unsigned char *)dstc,15, 3);
print_result (" C AddByte(3)", src1, NULL, dstc);
print_compare(dstm,dstc);
print_line();
SDL_imageFilterMMXon();
setup_src(src1, src2);
SDL_imageFilterAddByteToHalf ((unsigned char *)src1,(unsigned char *)dstm,15, 3);
print_result ("MMX AddByteToHalf(3)", src1, NULL, dstm);
SDL_imageFilterMMXoff();
setup_src(src1, src2);
SDL_imageFilterAddByteToHalf ((unsigned char *)src1,(unsigned char *)dstc,15, 3);
print_result (" C AddByteToHalf(3)", src1, NULL, dstc);
print_compare(dstm,dstc);
print_line();
SDL_imageFilterMMXon();
setup_src(src1, src2);
SDL_imageFilterSubByte ((unsigned char *)src1,(unsigned char *)dstm,15, 3);
print_result ("MMX SubByte(3)", src1, NULL, dstm);
SDL_imageFilterMMXoff();
setup_src(src1, src2);
SDL_imageFilterSubByte ((unsigned char *)src1,(unsigned char *)dstc,15, 3);
print_result (" C SubByte(3)", src1, NULL, dstc);
print_compare(dstm,dstc);
print_line();
SDL_imageFilterMMXon();
setup_src(src1, src2);
SDL_imageFilterShiftRight ((unsigned char *)src1,(unsigned char *)dstm,15, 1);
print_result ("MMX ShiftRight(1)", src1, NULL, dstm);
SDL_imageFilterMMXoff();
setup_src(src1, src2);
SDL_imageFilterShiftRight ((unsigned char *)src1,(unsigned char *)dstc,15, 1);
print_result (" C ShiftRight(1)", src1, NULL, dstc);
print_compare(dstm,dstc);
print_line();
SDL_imageFilterMMXon();
setup_src(src1, src2);
SDL_imageFilterMultByByte ((unsigned char *)src1,(unsigned char *)dstm,15, 3);
print_result ("MMX MultByByte(3)", src1, NULL, dstm);
SDL_imageFilterMMXoff();
setup_src(src1, src2);
SDL_imageFilterMultByByte ((unsigned char *)src1,(unsigned char *)dstc,15, 3);
print_result (" C MultByByte(3)", src1, NULL, dstc);
print_compare(dstm,dstc);
print_line();
SDL_imageFilterMMXon();
setup_src(src1, src2);
SDL_imageFilterShiftRightAndMultByByte ((unsigned char *)src1,(unsigned char *)dstm,15, 1, 3);
print_result ("MMX ShiftRightAndMultByByte(1,3)", src1, NULL, dstm);
SDL_imageFilterMMXoff();
setup_src(src1, src2);
SDL_imageFilterShiftRightAndMultByByte ((unsigned char *)src1,(unsigned char *)dstc,15, 1, 3);
print_result (" C ShuftRightAndMultByByte(1,3)", src1, NULL, dstc);
print_compare(dstm,dstc);
print_line();
SDL_imageFilterMMXon();
setup_src(src1, src2);
SDL_imageFilterShiftLeftByte ((unsigned char *)src1,(unsigned char *)dstm,15, 3);
print_result ("MMX ShiftLeftByte(3)", src1, NULL, dstm);
SDL_imageFilterMMXoff();
setup_src(src1, src2);
SDL_imageFilterShiftLeftByte ((unsigned char *)src1,(unsigned char *)dstc,15, 3);
print_result (" C ShiftLeftByte(3)", src1, NULL, dstc);
print_compare(dstm,dstc);
print_line();
SDL_imageFilterMMXon();
setup_src(src1, src2);
SDL_imageFilterShiftLeft ((unsigned char *)src1,(unsigned char *)dstm,15, 3);
print_result ("MMX ShiftLeft(3)", src1, NULL, dstm);
SDL_imageFilterMMXoff();
setup_src(src1, src2);
SDL_imageFilterShiftLeft ((unsigned char *)src1,(unsigned char *)dstc,15, 3);
print_result (" C ShiftLeft(3)", src1, NULL, dstc);
print_compare(dstm,dstc);
print_line();
SDL_imageFilterMMXon();
setup_src(src1, src2);
SDL_imageFilterBinarizeUsingThreshold ((unsigned char *)src1,(unsigned char *)dstm,15, 2);
print_result ("MMX BinarizeUsingThreshold(2)", src1, NULL, dstm);
SDL_imageFilterMMXoff();
setup_src(src1, src2);
SDL_imageFilterBinarizeUsingThreshold ((unsigned char *)src1,(unsigned char *)dstc,15, 2);
print_result (" C BinarizeUsingThreshold(2)", src1, NULL, dstc);
print_compare(dstm,dstc);
print_line();
SDL_imageFilterMMXon();
setup_src(src1, src2);
SDL_imageFilterClipToRange ((unsigned char *)src1,(unsigned char *)dstm,15, 1,7);
print_result ("MMX ClipToRange(1,7)", src1, NULL, dstm);
SDL_imageFilterMMXoff();
setup_src(src1, src2);
SDL_imageFilterClipToRange ((unsigned char *)src1,(unsigned char *)dstc,15, 1,7);
print_result (" C ClipToRange(1,7)", src1, NULL, dstc);
print_compare(dstm,dstc);
print_line();
SDL_imageFilterMMXon();
setup_src(src1, src2);
SDL_imageFilterNormalizeLinear ((unsigned char *)src1,(unsigned char *)dstm,15, 0,33,0,255);
print_result ("MMX NormalizeLinear(0,33,0,255)", src1, NULL, dstm);
SDL_imageFilterMMXoff();
setup_src(src1, src2);
SDL_imageFilterNormalizeLinear ((unsigned char *)src1,(unsigned char *)dstc,15, 0,33,0,255);
print_result (" C NormalizeLinear(0,33,0,255)", src1, NULL, dstc);
print_compare(dstm,dstc);
print_line();
/* Uint functions */
/* Disabled, since broken */
/*
SDL_imageFilterMMXon();
setup_src(src1, src2);
SDL_imageFilterAddUint ((unsigned char *)src1,(unsigned char *)dstm,15, 0x01020304);
print_result ("MMX AddUint(0x01020304)", src1, NULL, dstm);
SDL_imageFilterMMXoff();
setup_src(src1, src2);
SDL_imageFilterAddUint ((unsigned char *)src1,(unsigned char *)dstc,15, 0x01020304);
print_result (" C AddUint(0x01020304)", src1, NULL, dstc);
print_compare(dstm,dstc);
print_line();
SDL_imageFilterMMXon();
setup_src(src1, src2);
SDL_imageFilterSubUint ((unsigned char *)src1,(unsigned char *)dstm,15, 0x01020304);
print_result ("MMX SubUint(0x01020304)", src1, NULL, dstm);
SDL_imageFilterMMXoff();
setup_src(src1, src2);
SDL_imageFilterSubUint ((unsigned char *)src1,(unsigned char *)dstc,15, 0x01020304);
print_result (" C SubUint(0x01020304)", src1, NULL, dstc);
print_compare(dstm,dstc);
print_line();
SDL_imageFilterMMXon();
setup_src(src1, src2);
SDL_imageFilterShiftRightUint ((unsigned char *)src1,(unsigned char *)dstm,15, 4);
print_result ("MMX ShiftRightUint(4)", src1, NULL, dstm);
SDL_imageFilterMMXoff();
setup_src(src1, src2);
SDL_imageFilterShiftRightUint ((unsigned char *)src1,(unsigned char *)dstc,15, 4);
print_result (" C ShiftRightUint(4)", src1, NULL, dstc);
print_compare(dstm,dstc);
print_line();
SDL_imageFilterMMXon();
setup_src(src1, src2);
SDL_imageFilterShiftLeftUint ((unsigned char *)src1,(unsigned char *)dstm,15, 4);
print_result ("MMX ShiftLeftUint(4)", src1, NULL, dstm);
SDL_imageFilterMMXoff();
setup_src(src1, src2);
SDL_imageFilterShiftLeftUint ((unsigned char *)src1,(unsigned char *)dstc,15, 4);
print_result (" C ShiftLeftUint(4)", src1, NULL, dstc);
print_compare(dstm,dstc);
print_line();
*/
#ifdef USE_MMX
SDL_imageFilterMMXon();
if (SDL_imageFilterMMXdetect())
{
printf("MMX was detected\n\n");
}
else
{
printf("MMX was NOT detected\n\n");
}
#else
printf("MMX support disabled in SDL_gfx.\n\n");
#endif
printf ("Result: %i of %i passed OK.\n", ok_count, total_count);
#ifdef WIN32
printf("Press Enter to continue ...");
pause();
#endif
exit(0);
}
