2023.11.26

VERSION

@@ -1,2 +1,2 @@
-November 13 2023
+November 26 2023
 

doc/midistats.1

@@ -1,4 +1,4 @@
-.TH MIDISTATS 1 "1 November 2023"
+.TH MIDISTATS 1 "17 November 2023"
 .SH NAME
 \fBmidistats\fP \- program to summarize the statistical properties of a midi file
 .SH SYNOPSIS

doc/readme.txt

@@ -6,7 +6,7 @@ abc2abc  version 2.20 February 07 2023
 yaps     version 1.92 January 06 2023
 abcmatch version 1.82 June 14 2022
 midicopy version 1.39 November 08 2022
-midistats version 0.80 November 13 2023
+midistats version 0.81 November 26 2023
 
 24th January 2002
 Copyright James Allwright

midistats.c

@@ -18,7 +18,7 @@
  * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301, USA
  */
  
-#define VERSION "0.80 November 13 2023 midistats"
+#define VERSION "0.81 November 26 2023 midistats"
 
 #include <limits.h>
 /* Microsoft Visual C++ Version 6.0 or higher */
@@ -49,7 +49,7 @@ extern char* strchr();
 #include "midifile.h"
 void initfuncs();
 void stats_finish();
-float histogram_entropy (int *histogram, int size); 
+float histogram_perplexity (int *histogram, int size); 
 void stats_noteoff(int chan,int pitch,int vol);
 void stats_eot ();
 #define max(a,b)  (( a > b ? a : b))
@@ -67,6 +67,7 @@ int lasttrack = 0; /* lasttrack */
 int division;    /* pulses per quarter note defined in MIDI header    */
 int quietLimit;  /* minimum number of pulses with no activity */
 long tempo = 500000; /* the default tempo is 120 quarter notes/minute */
+int bpm = 120; /*default tempo */
 long laston = 0; /* length of MIDI track in pulses or ticks           */
 int key[12];
 int sharps;
@@ -475,15 +476,17 @@ int i;
 /* [SS] 2023-10-30 */
 void stats_interpret_pulseCounter () {
 int i,j;
-int maxcount,ncounts;
+int maxcount,ncounts,npeaks,npositives,peaklimit;
 int maxloc;
 float threshold,peak;
 int decimate;
 float tripletsCriterion8,tripletsCriterion4;
 int resolution = 12;
+int nzeros;
 threshold = 10.0/(float) division;
 maxcount = 0;
 ncounts = 0;
+npeaks = 0;
 for (i=0;i<division;i++) {
   ncounts = ncounts + pulseCounter[i];
   if (pulseCounter[i] > maxcount) {
@@ -491,6 +494,10 @@ for (i=0;i<division;i++) {
        maxcount = pulseCounter[i];
        } 
   } 
+peaklimit = (int) (ncounts * 0.020);
+for (i=0;i<division;i++) {
+  if (pulseCounter[i] > peaklimit) npeaks++;
+  }
 for (i = 0; i < resolution; i++) pulseDistribution[i] = 0;
 decimate = division/resolution;
 for (i = 0; i < division; i++) {
@@ -498,14 +505,28 @@ for (i = 0; i < division; i++) {
    pulseDistribution[j] += pulseCounter[i];
   }
 
+/* count zeros */
+nzeros = 0;
+for (i=0;i<resolution;i++) if((float) pulseDistribution[i]/(float) ncounts < 0.015 ) nzeros++;
+npositives = resolution - nzeros;
+if (nzeros > 3 && (float) pulseDistribution[resolution-1]/(float) ncounts < 0.1) {printf("clean_quantization\n");
+ } else if ((float) pulseDistribution[resolution-1]/(float) ncounts > 0.09 ||
+            npeaks > npositives) {printf("dithered_quantization\n");
+ } else {
 peak = (float) maxcount/ (float) ncounts;
-/*printf("maxcount = %d ncounts = %d peak = %f threshold = %f\n",maxcount,ncounts,peak,threshold); */
 if (peak < threshold)  printf("unquantized\n");
+}
+
 tripletsCriterion8 = (float) pulseDistribution[8]/ (float) ncounts;
 tripletsCriterion4 = (float) pulseDistribution[4]/ (float) ncounts;
-/*printf("tripletsCriterion = %f\n",tripletsCriterion);*/
 if (tripletsCriterion8 > 0.10 || tripletsCriterion4 > 0.10) printf("triplets\n");
 if (pulseDistribution[0]/(float) ncounts > 0.95) printf("qnotes");
+/*
+printf("pulseDistribution:");
+for (i=0;i<resolution;i++) printf("%6.3f",(float) pulseDistribution[i]/(float) ncounts);
+printf("\n");
+printf("nzeros = %d npeaks = %d \n",nzeros,npeaks);
+*/
 }
 
 void stats_finish()
@@ -574,7 +595,7 @@ else
 printf("\ntrkact ");
   lasttrack++;
   for (i=0;i<lasttrack;i++) printf("% 5d",trkactivity[i]);
-printf("\npitchentropy %f\n",histogram_entropy(pitchclass_activity,12));
+printf("\npitchperplexity %f\n",histogram_perplexity(pitchclass_activity,12));
 printf("totalrhythmpatterns =%d\n",nrpatterns);
 printf("collisions = %d\n",ncollisions);
 if (hasLyrics) printf("Lyrics\n");
@@ -584,8 +605,9 @@ printf("\n");
 
 
 
-float histogram_entropy (int *histogram, int size) 
+float histogram_perplexity (int *histogram, int size) 
   {
+  /* The perplexity is 2 to the power of the entropy */
   int i;
   int total;
   float entropy;
@@ -604,7 +626,7 @@ float histogram_entropy (int *histogram, int size)
     entropy = entropy + e;
     } 
   //printf("\n");
-  return -entropy/log(2.0);
+  return pow(2.0,-entropy/log(2.0));
   } 
 
 
@@ -693,7 +715,7 @@ void stats_trackend()
    if (channel_used_in_track[chan] > 0) trkdata.quietTime[chan] += (trkdata.npulses[0] - trkdata.lastNoteOff[chan]); 
  for (chan=0;chan<16;chan++) {  /* 2023-09-13 */
    if (chan == 9 || channel_used_in_track[chan+1] == 0) continue;
-   trkdata.pitchEntropy[chan+1] = histogram_entropy(chanpitchhistogram +chan*12,11);
+   trkdata.pitchEntropy[chan+1] = histogram_perplexity(chanpitchhistogram +chan*12,11);
    }
  output_track_summary(); 
 }
@@ -938,6 +960,12 @@ void record_noteoff(int chan,int pitch,int vol)
 void record_trackend()
 {
 }
+void record_tempo(long ltempo)
+{ 
+tempo = ltempo;
+if (bpm == 120) bpm = 60000000.0/tempo;
+tempocount++;
+}
 
 int int_compare_events(const void *a, const void *b) {
    struct eventstruc *ia  = (struct eventstruc *)a;
@@ -1006,7 +1034,7 @@ void initfunc_for_loadNoteEvents()
     Mf_eot = no_op0;
     Mf_timesig = no_op4;
     Mf_smpte = no_op5;
-    Mf_tempo = no_op1;
+    Mf_tempo = record_tempo;
     Mf_keysig = no_op2;
     Mf_seqspecific = no_op3;
     Mf_text = no_op3;
@@ -1231,7 +1259,8 @@ int nchannels;
 nchannels = 0;
 for (i=1;i<17;i++) 
   if (channel_active[i] > 0) nchannels++;
-printf("%d\t%d\t%d\t%d\t%d\n",lasttrack,nchannels, division,lastEvent,lastBeat);
+printf("%d\t%d\t%d\t%d\t%d\t%d\n",lasttrack,nchannels, division,bpm,lastEvent,lastBeat);
+/*printf("%d\n",tempocount);*/
 }