You can not select more than 25 topics
Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.
426 lines
9.7 KiB
426 lines
9.7 KiB
/* 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()
|
|
{
|
|
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
|
|
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
|
|
|