c@37: 
c@37: program test;
c@65: /*
c@37: af = load may.stream.audiofile;
c@37: plot = load may.plot;
c@37: cm = load may.matrix.complex;
c@37: mat = load may.matrix;
c@37: vec = load may.vector;
c@59: win = load may.signal.window;
c@59: mm = load may.mathmisc;
c@37: manipulate = load may.stream.manipulate;
c@37: syn = load may.stream.syntheticstream;
c@65: */
c@37: 
c@65: { runTests } = load may.test;
c@65: 
c@65: //{ makeKernel } = load cqtkernel;
c@65: //{ cqt } = load cqt;
c@37: 
c@60: // We want to test:
c@60: // 
c@60: // Kernel design -- check size (number of bins, number of atoms);
c@60: // check an example kernel against known data
c@60: //
c@60: // Time alignment -- feed a dirac train, check that peaks in all bins
c@60: // align
c@60: //
c@60: // Frequency discrimination -- feed a sinusoid, check peaks
c@60: //
c@60: // Latency compensation -- for dirac at 0, check peak can be found at
c@60: // 0 plus the declared latency
c@60: //
c@60: // Signal-noise ratio
c@60: //
c@60: // Specimen output for simple test case
c@60: 
c@65: tests = [
c@65: "kernel"    : load test_cqtkernel,
c@65: "frequency" : load test_frequency,
c@65: ];
c@59: 
c@65: bad = sum (mapHash do name testHash: runTests name testHash done tests);
c@59: 
c@65: if (bad > 0) then
c@65:     println "\n** \(bad) test(s) failed!";
c@65:     threadExit 1
c@65: else
c@65:     ()
c@65: fi
c@59: 
c@59: 
c@59: 
c@59: 
c@59: /*
c@37: //testStream = manipulate.withDuration 96000 (syn.sinusoid 48000 500);
c@37: //testStream = manipulate.withDuration 96000 (syn.pulseTrain 48000 4);
c@44: testStream = af.open "sweep-48000.wav";
c@37: //testStream = af.open "sweep.wav";
c@37: 
c@43: // So the stream is [ 0, 1, 0, -1, 0, 1, 0, -1, ... ] :
c@44: //testStream = manipulate.withDuration 64 (syn.sinusoid 8 2);
c@38: 
c@56: testStream = manipulate.withDuration 32 (syn.pulseTrain 8 0.001);
c@37: 
c@37: eprintln "have test stream";
c@37: 
c@44: cq = cqt { maxFreq = testStream.sampleRate/2, minFreq = 50, binsPerOctave = 24 } testStream;
c@37: 
c@40: eprintln "bin frequencies: \(cq.kernel.binFrequencies)";
c@40: 
c@40: bigM = mat.concatHorizontal (map cm.magnitudes cq.output);
c@37: 
c@38: eprintln "overall output size = \(mat.size bigM)";
c@38: 
c@39: mat.print bigM;
c@38: 
c@38: //\() (plot.plot [Contour bigM]);
c@56: \() (plot.plot [Grid bigM]);
c@59: */
c@65: //()
c@37: