public void run() {
runInit(); // superclass
// Haupt-Variablen fuer den Prozess
int ch, i, j, k, m;
double d1, d2;
SpectStreamSlot[] runInSlot = new SpectStreamSlot[2];
SpectStreamSlot runOutSlot;
SpectStream[] runInStream = new SpectStream[2];
SpectStream runOutStream;
SpectFrame[] runInFr = new SpectFrame[2];
SpectFrame runOutFr;
// Ziel-Frame Berechnung
int[] srcBands = new int[2];
int fftSize, fullFFTsize, complexFFTsize;
float[] convBuf1, convBuf2;
float[][] fftBuf;
float[] win;
int readDone, oldReadDone;
int[] phase = new int[2];
topLevel:
try {
// ------------------------------ Input-Slot ------------------------------
for (i = 0; i < 2; i++) {
runInSlot[i] = slots.elementAt(SLOT_INPUT1 + i);
if (runInSlot[i].getLinked() == null) {
runStop(); // threadDead = true -> folgendes for() wird uebersprungen
}
// diese while Schleife ist noetig, da beim initReader ein Pause eingelegt werden kann
// und die InterruptException ausgeloest wird; danach versuchen wir es erneut
for (boolean initDone = false; !initDone && !threadDead; ) {
try {
runInStream[i] = runInSlot[i].getDescr(); // throws InterruptedException
initDone = true;
srcBands[i] = runInStream[i].bands;
} catch (InterruptedException ignored) {
}
runCheckPause();
}
}
if (threadDead) break topLevel;
// ------------------------------ Output-Slot ------------------------------
runOutSlot = slots.elementAt(SLOT_OUTPUT);
runOutStream = new SpectStream(runInStream[0]);
runOutSlot.initWriter(runOutStream);
// ------------------------------ Vorberechnungen ------------------------------
fftSize = srcBands[0] - 1;
fullFFTsize = fftSize << 1;
complexFFTsize = fullFFTsize << 1;
fftBuf = new float[2][complexFFTsize];
win = new float[fullFFTsize];
d1 = 1.0 / (double) fullFFTsize * Math.PI;
for (i = 0; i < fullFFTsize; i++) {
d2 = Math.cos(i * d1);
win[i] = (float) (d2 * d2);
}
phase[0] = (int) (pr.para[PR_ANGLE].value / 100.0 * fullFFTsize + 0.5) % fullFFTsize;
phase[1] = (phase[0] + fftSize) % fullFFTsize;
// convBuf2 = fftBuf[0];
// Util.clear( convBuf2 );
// for( k = 0; k < 2; k++ ) {
// m = k * fftSize;
// for( i = 0; m < fullFFTsize; ) {
// convBuf2[ i++ ] += win[ m ];
// convBuf2[ i++ ] += win[ m++ ];
// }
// for( m = 0; i < complexFFTsize; ) {
// convBuf2[ i++ ] += win[ m ];
// convBuf2[ i++ ] += win[ m++ ];
// }
// }
// Debug.view( convBuf2, "win overlap" );
// ------------------------------ Hauptschleife ------------------------------
runSlotsReady();
mainLoop:
while (!threadDead) {
// ---------- Frame einlesen ----------
for (readDone = 0; (readDone < 2) && !threadDead; ) {
oldReadDone = readDone;
for (i = 0; i < 2; i++) {
try {
if (runInStream[i].framesReadable() > 0) {
runInFr[i] = runInSlot[i].readFrame();
readDone++;
}
} catch (InterruptedException ignored) {
} catch (EOFException e) {
break mainLoop;
}
runCheckPause();
}
if (oldReadDone == readDone) { // konnte nix gelesen werden
try {
Thread.sleep(500); // ...deshalb kurze Pause
} catch (InterruptedException ignored) {
} // mainLoop wird gleich automatisch verlassen
runCheckPause();
}
}
if (threadDead) break mainLoop;
runOutFr = runOutStream.allocFrame();
// ---------- Process: Ziel-Frame berechnen ----------
for (ch = 0; ch < runOutStream.chanNum; ch++) {
for (k = 0; k < 2; k++) {
convBuf1 = runInFr[k].data[ch];
convBuf2 = fftBuf[k];
// calculate complex log spectrum :
// Re( target ) = Log( Mag( source ))
// Im( target ) = Phase( source )
// convBuf1 is already in polar style
// -> fftBuf will be in rect style
for (i = 0; i <= fullFFTsize; ) {
convBuf2[i] = (float) Math.log(Math.max(1.0e-24, convBuf1[i])); // Re( target )
i++;
convBuf2[i] = convBuf1[i]; // Im( target )
i++;
}
// make full spectrum (negative frequencies = conjugate positive frequencies)
for (i = fullFFTsize + 2, j = fullFFTsize - 2; i < complexFFTsize; j -= 2) {
// bug but nice? - 1
convBuf2[i++] = convBuf2[j];
convBuf2[i++] = -convBuf2[j + 1];
}
// cepstrum domain
Fourier.complexTransform(convBuf2, fullFFTsize, Fourier.INVERSE);
// window
m = phase[k];
for (i = 0; m < fullFFTsize; ) {
convBuf2[i++] *= win[m];
convBuf2[i++] *= win[m++];
}
for (m = 0; i < complexFFTsize; ) {
convBuf2[i++] *= win[m];
convBuf2[i++] *= win[m++];
}
}
// mix cepstra
convBuf1 = fftBuf[0];
convBuf2 = runOutFr.data[ch];
Util.add(fftBuf[1], 0, convBuf1, 0, complexFFTsize);
// back to frequency domain
Fourier.complexTransform(convBuf1, fullFFTsize, Fourier.FORWARD);
// calculate real exponential spectrum :
// Mag( target ) = Exp( Re( source ))
// Phase( target ) = Im( source )
// ->convBuf2 shall be polar style, that makes things easy
for (i = 0; i <= fullFFTsize; ) {
convBuf2[i] = (float) Math.exp(convBuf1[i]);
i++;
convBuf2[i] = convBuf1[i];
i++;
}
}
// calculation done
runInSlot[0].freeFrame(runInFr[0]);
runInSlot[1].freeFrame(runInFr[1]);
for (boolean writeDone = false; (!writeDone) && !threadDead; ) {
try { // Unterbrechung
runOutSlot.writeFrame(runOutFr); // throws InterruptedException
writeDone = true;
runFrameDone(runOutSlot, runOutFr);
runOutStream.freeFrame(runOutFr);
} catch (InterruptedException ignored) {
} // mainLoop wird eh gleich verlassen
runCheckPause();
}
} // end of main loop
runInStream[0].closeReader();
runInStream[1].closeReader();
runOutStream.closeWriter();
} // break topLevel
catch (IOException e) {
runQuit(e);
return;
} catch (SlotAlreadyConnectedException e) {
runQuit(e);
return;
}
// catch( OutOfMemoryError e ) {
// abort( e );
// return;
// }
runQuit(null);
}
