diff options
Diffstat (limited to 'noatun-plugins/tippercanoe/core.cpp')
| -rw-r--r-- | noatun-plugins/tippercanoe/core.cpp | 425 |
1 files changed, 425 insertions, 0 deletions
diff --git a/noatun-plugins/tippercanoe/core.cpp b/noatun-plugins/tippercanoe/core.cpp new file mode 100644 index 0000000..0fb6fa2 --- /dev/null +++ b/noatun-plugins/tippercanoe/core.cpp @@ -0,0 +1,425 @@ +/* Synaescope - a pretty noatun visualization (based on P. Harrison's Synaesthesia) + Copyright (C) 1997 Paul Francis Harrison <pfh@yoyo.cc.monash.edu.au> + 2001 Charles Samuels <charles@kde.org> + + This program is free software; you can redistribute it and/or modify it + under the terms of the GNU General Public License as published by the + Free Software Foundation; either version 2 of the License, or (at your + option) any later version. + + This program is distributed in the hope that it will be useful, but + WITHOUT ANY WARRANTY; without even the implied warranty of + MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU + General Public License for more details. + + You should have received a copy of the GNU General Public License along + with this program; if not, write to the Free Software Foundation, Inc., + 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA. +*/ + +#include <stdio.h> +#include <stdlib.h> +#include <math.h> +#include <string.h> +#include "syna.h" +#include <unistd.h> +Core *core; + +#define outputs unsigned char *Dlo=(unsigned char*)lastOutputBmp.data; \ +unsigned char *Dllo=(unsigned char*)lastLastOutputBmp.data; \ +unsigned char *Do=(unsigned char*)outputBmp.data; + +void Core::setupPalette(double) { interface->setupPalette(); } + +int Core::bitReverser(int i) +{ + int sum=0,j; + for(j=0;j<LogSize;j++) + { + sum = (i&1)+sum*2; + i >>= 1; + } + return sum; +} + +void Core::fft(double *x,double *y) +{ + int n2 = NumSamples, n1; + int twoToTheK; + for(twoToTheK=1;twoToTheK<NumSamples;twoToTheK*=2) + { + n1 = n2; + n2 /= 2; + for(int j=0;j<n2;j++) + { + double c = cosTable[j*twoToTheK&(NumSamples-1)], + s = negSinTable[j*twoToTheK&(NumSamples-1)]; + for(int i=j;i<NumSamples;i+=n1) + { + int l = i+n2; + double xt = x[i] - x[l]; + x[i] = (x[i] + x[l]); + double yt = y[i] - y[l]; + y[i] = (y[i] + y[l]); + x[l] = xt*c - yt*s; + y[l] = xt*s + yt*c; + } + } + } +} + +void Core::setStarSize(double size) +{ + //int factor = (fadeMode == Flame ? 100 : + // (fadeMode == Wave ? 150 : 200)); + + double fadeModeFudge = (fadeMode == Wave ? 0.4 : + (fadeMode == Flame ? 0.6 : 0.78)); + + int factor; + if (size > 0.0) + factor = int(exp(log(fadeModeFudge) / (size*8.0))*255); + else + factor = 0; + + if (factor > 255) factor = 255; + + for(int i=0;i<256;i++) + scaleDown[i] = i*factor>>8; + + maxStarRadius = 1; + for(int i=255;i;i = scaleDown[i]) + maxStarRadius++; +} + +inline void Core::addPixel(int x,int y,int br1,int br2) +{ + if (x < 0 || x >= outWidth || y < 0 || y >= outHeight) return; + + unsigned char *p = output()+x*2+y*outWidth*2; + if (p[0] < 255-br1) p[0] += br1; else p[0] = 255; + if (p[1] < 255-br2) p[1] += br2; else p[1] = 255; + //p += lastOutput()-output(); + //if (p[0] < 255-br1) p[0] += br1; else p[0] = 255; + //if (p[1] < 255-br2) p[1] += br2; else p[1] = 255; +} + +inline void Core::addPixelFast(unsigned char *p,int br1,int br2) +{ + if (p[0] < 255-br1) p[0] += br1; else p[0] = 255; + if (p[1] < 255-br2) p[1] += br2; else p[1] = 255; + //p += lastOutput()-output(); + //if (p[0] < 255-br1) p[0] += br1; else p[0] = 255; + //if (p[1] < 255-br2) p[1] += br2; else p[1] = 255; +} + +void Core::fadeFade() +{ + register unsigned long *ptr = (unsigned long*)output(); + int i = outWidth*outHeight*2/4; + do { + //Bytewize version was: *(ptr++) -= *ptr+(*ptr>>1)>>4; + if (*ptr) { + //if (*ptr & 0xf0f0f0f0ul) + *ptr -= ((*ptr & 0xf0f0f0f0ul) >> 4) + ((*ptr & 0xe0e0e0e0ul) >> 5); + ptr++; + } + + //else { + // *(ptr++) = (*ptr * 14 >> 4) & 0x0f0f0f0ful; + //} + else + ptr++; + } while(--i > 0); +} + +inline void Core::fadePixelWave(int x,int y,int where,int step) +{ +outputs + short j = + ( short(getPixel(x-1,y,where-2))+ + getPixel(x+1,y,where+2)+ + getPixel(x,y-1,where-step)+ + getPixel(x,y+1,where+step) + >> 2) + +Dlo[where]; + if (!j) { Do[where] = 0; return; } + j = j + -Dllo[where] + -1; + if (j < 0) Do[where] = 0; + else if (j & (255*256)) Do[where] = 255; + else Do[where] = j; +} + +void Core::fadeWave() +{ + unsigned short *t = lastLastOutputBmp.data; + lastLastOutputBmp.data = lastOutputBmp.data; + lastOutputBmp.data = outputBmp.data; + outputBmp.data = t; + + int x,y,i,j,start,end; + int step = outWidth*2; + for(x=0,i=0,j=outWidth*(outHeight-1)*2;x<outWidth;x++,i+=2,j+=2) { + fadePixelWave(x,0,i,step); + fadePixelWave(x,0,i+1,step); + fadePixelWave(x,outHeight-1,j,step); + fadePixelWave(x,outHeight-1,j+1,step); + } + + for(y=1,i=outWidth*2,j=outWidth*4-2;y<outHeight;y++,i+=step,j+=step) { + fadePixelWave(0,y,i,step); + fadePixelWave(0,y,i+1,step); + fadePixelWave(outWidth-1,y,j,step); + fadePixelWave(outWidth-1,y,j+1,step); + } +outputs + + for(y=1, + start=outWidth*2+2, + end=outWidth*4-2; y<outHeight-1; y++,start+=step,end+=step) { + int i = start; + do { + short j = + ( short(Dlo[i-2])+ + Dlo[i+2]+ + Dlo[i-step]+ + Dlo[i+step] + >> 2) + +Dlo[i]; + if (!j) { + Do[i] = 0; + } else { + j = j + -Dllo[i] + -1; + if (j < 0) Do[i] = 0; + else if (j & (255*256)) Do[i] = 255; + else Do[i] = j; + } + } while(++i < end); + } +} + +inline void Core::fadePixelHeat(int x,int y,int where,int step) +{ +outputs + short j = + ( short(getPixel(x-1,y,where-2))+ + getPixel(x+1,y,where+2)+ + getPixel(x,y-1,where-step)+ + getPixel(x,y+1,where+step) + >> 2) + +Dlo[where]; + if (!j) { Do[where] = 0; return; } + j = j + -Dllo[where] + -1; + if (j < 0) Do[where] = 0; + else if (j & (255*256)) Do[where] = 255; + else Do[where] = j; +} + +void Core::fadeHeat() +{ + unsigned short *t = lastLastOutputBmp.data; + lastLastOutputBmp.data = lastOutputBmp.data; + lastOutputBmp.data = outputBmp.data; + outputBmp.data = t; + + int x,y,i,j,start,end; + int step = outWidth*2; + for(x=0,i=0,j=outWidth*(outHeight-1)*2;x<outWidth;x++,i+=2,j+=2) + { + fadePixelHeat(x,0,i,step); + fadePixelHeat(x,0,i+1,step); + fadePixelHeat(x,outHeight-1,j,step); + fadePixelHeat(x,outHeight-1,j+1,step); + } + + for(y=1,i=outWidth*2,j=outWidth*4-2;y<outHeight;y++,i+=step,j+=step) + { + fadePixelHeat(0,y,i,step); + fadePixelHeat(0,y,i+1,step); + fadePixelHeat(outWidth-1,y,j,step); + fadePixelHeat(outWidth-1,y,j+1,step); + } + +outputs + for(y=1,start=outWidth*2+2, + end=outWidth*4-2; y<outHeight-1; y++,start+=step,end+=step) + { + int i = start; + do + { + short j = + ( short(Dlo[i-2])+ + Dlo[i+2]+ + +Dlo[i-step] + +Dlo[i+step] + >> 2) + +Dlo[i]; + if (!j) { + Do[i] = 0; + } + else + { + j = j + -Dllo[i] + +(Dllo[i] + -Dlo[i]>>2) + -1; + if (j < 0) Do[i] = 0; + else if (j & (255*256)) Do[i] = 255; + else Do[i] = j; + } + } while(++i < end); + } +} + +void Core::fade() +{ + switch(fadeMode) + { + case Stars : + fadeFade(); + break; + case Flame : + fadeHeat(); + break; + case Wave : + fadeWave(); + break; + default: + break; + } +} + +bool Core::calculate() +{ + double x[NumSamples], y[NumSamples]; + double a[NumSamples], b[NumSamples]; + int clarity[NumSamples]; //Surround sound + int i,j,k; +#ifndef LITTLEENDIAN + register sampleType temp; +#endif + + int brightFactor = int(Brightness * brightnessTwiddler /(starSize+0.01)); + + int num; + if ((num=read(0, data, NumSamples))==0) + return false; + + for(i=0;i<NumSamples;i++) + { +# ifdef LITTLEENDIAN + x[i] = data[i*2]; + y[i] = data[i*2+1]; +# else + // Need to convert to big-endian + temp = data[i*2]; + temp = (temp >> 8) | (temp << 8); + x[i] = temp; + temp = data[i*2+1]; + temp = (temp << 8) | (temp >> 8); + y[i] = temp; +# endif + } + + fft(x,y); + + for(i=0 +1;i<NumSamples;i++) + { + double x1 = x[bitReverse[i]], + y1 = y[bitReverse[i]], + x2 = x[bitReverse[NumSamples-i]], + y2 = y[bitReverse[NumSamples-i]], + aa,bb; + a[i] = sqrt(aa= (x1+x2)*(x1+x2) + (y1-y2)*(y1-y2) ); + b[i] = sqrt(bb= (x1-x2)*(x1-x2) + (y1+y2)*(y1+y2) ); + if (aa+bb != 0.0) + clarity[i] = (int)( + ( (x1+x2) * (x1-x2) + (y1+y2) * (y1-y2) )/(aa+bb) * 256 ); + else + clarity[i] = 0; + } + + int heightFactor = NumSamples/2 / outHeight + 1; + int actualHeight = NumSamples/2/heightFactor; + int heightAdd = outHeight + actualHeight >> 1; + + // Correct for window size + double brightFactor2 = (brightFactor/65536.0/NumSamples)* + sqrt(actualHeight*outWidth/(320.0*200.0)); + + for(i=1;i<NumSamples/2;i++) { + //int h = (int)( b[i]*280 / (a[i]+b[i]+0.0001)+20 ); + if (a[i] > 0 || b[i] > 0) { + int h = (int)( b[i]*outWidth / (a[i]+b[i]) ); + int br1, br2, br = (int)( + (a[i]+b[i])*i*brightFactor2 ); + br1 = br*(clarity[i]+128)>>8; + br2 = br*(128-clarity[i])>>8; + if (br1 < 0) br1 = 0; else if (br1 > 255) br1 = 255; + if (br2 < 0) br2 = 0; else if (br2 > 255) br2 = 255; + //unsigned char *p = output+ h*2+(164-((i<<8)>>m))*(outWidth*2); + int px = h, + py = heightAdd - i / heightFactor; + + if (pointsAreDiamonds) + { + addPixel(px,py,br1,br2); + br1=scaleDown[br1];br2=scaleDown[br2]; + + //TODO: Use addpixelfast + for(j=1;br1>0||br2>0;j++,br1=scaleDown[br1],br2=scaleDown[br2]) + { + for(k=0;k<j;k++) + { + addPixel(px-j+k,py-k,br1,br2); + addPixel(px+k,py-j+k,br1,br2); + addPixel(px+j-k,py+k,br1,br2); + addPixel(px-k,py+j-k,br1,br2); + } + } + } + else + { + if (px < maxStarRadius || py < maxStarRadius || + px > outWidth-maxStarRadius || py > outHeight-maxStarRadius) + { + addPixel(px,py,br1,br2); + for(j=1;br1>0||br2>0;j++,br1=scaleDown[br1],br2=scaleDown[br2]) + { + addPixel(px+j,py,br1,br2); + addPixel(px,py+j,br1,br2); + addPixel(px-j,py,br1,br2); + addPixel(px,py-j,br1,br2); + } + } + else + { + unsigned char *p = output()+px*2+py*outWidth*2, *p1=p, *p2=p, *p3=p, *p4=p; + addPixelFast(p,br1,br2); + for(;br1>0||br2>0;br1=scaleDown[br1],br2=scaleDown[br2]) + { + p1 += 2; + addPixelFast(p1,br1,br2); + p2 -= 2; + addPixelFast(p2,br1,br2); + p3 += outWidth*2; + addPixelFast(p3,br1,br2); + p4 -= outWidth*2; + addPixelFast(p4,br1,br2); + } + } + } + } + } + return true; +} + +#undef outputs + |