+ return pulses;
+}
+
+void output_cycle_plot(const std::vector<pulse> &pulses, double calibration_factor)
+{
+ FILE *fp = fopen("cycles.plot", "w");
+ for (unsigned i = 0; i < pulses.size(); ++i) {
+ double cycles = pulses[i].len * calibration_factor * C64_FREQUENCY;
+ fprintf(fp, "%f %f\n", pulses[i].time, cycles);
+ }
+ fclose(fp);
+}
+
+int main(int argc, char **argv)
+{
+ parse_options(argc, argv);
+
+ make_lanczos_weight_table();
+ std::vector<float> pcm;
+ int sample_rate;
+ if (!read_audio_file(argv[optind], &pcm, &sample_rate)) {
+ exit(1);
+ }
+
+ if (do_crop) {
+ pcm = crop(pcm, crop_start, crop_end, sample_rate);
+ }
+
+ if (use_filter) {
+ pcm = do_filter(pcm, filter_coeff);
+ }
+
+#if 0
+ for (int i = 0; i < LEN; ++i) {
+ in[i] += rand() % 10000;
+ }
+#endif
+
+#if 0
+ for (int i = 0; i < LEN; ++i) {
+ printf("%d\n", in[i]);
+ }
+#endif
+
+ std::vector<pulse> pulses = detect_pulses(pcm, sample_rate);