lib/analysis.c

   1 /********************************************************************
   2  *                                                                  *
   3  * THIS FILE IS PART OF THE OggVorbis SOFTWARE CODEC SOURCE CODE.   *
   4  * USE, DISTRIBUTION AND REPRODUCTION OF THIS LIBRARY SOURCE IS     *
   5  * GOVERNED BY A BSD-STYLE SOURCE LICENSE INCLUDED WITH THIS SOURCE *
   6  * IN 'COPYING'. PLEASE READ THESE TERMS BEFORE DISTRIBUTING.       *
   7  *                                                                  *
   8  * THE OggVorbis SOURCE CODE IS (C) COPYRIGHT 1994-2002             *
   9  * by the XIPHOPHORUS Company http://www.xiph.org/                  *
  10  *                                                                  *
  11  ********************************************************************
  12
  13  function: single-block PCM analysis mode dispatch
  14  last mod: $Id: analysis.c,v 1.48 2002/01/22 08:06:06 xiphmont Exp $
  15
  16  ********************************************************************/
  17
  18 #include <stdio.h>
  19 #include <string.h>
  20 #include <math.h>
  21 #include <ogg/ogg.h>
  22 #include "vorbis/codec.h"
  23 #include "codec_internal.h"
  24 #include "registry.h"
  25 #include "scales.h"
  26 #include "os.h"
  27 #include "misc.h"
  28
  29 int analysis_noisy=1;
  30
  31 /* decides between modes, dispatches to the appropriate mapping. */
  32 int vorbis_analysis(vorbis_block *vb, ogg_packet *op){
  33   vorbis_dsp_state     *vd=vb->vd;
  34   backend_lookup_state *b=vd->backend_state;
  35   vorbis_info          *vi=vd->vi;
  36   codec_setup_info     *ci=vi->codec_setup;
  37   int                   type,ret;
  38   int                   mode=0;
  39
  40   vb->glue_bits=0;
  41   vb->time_bits=0;
  42   vb->floor_bits=0;
  43   vb->res_bits=0;
  44
  45   /* first things first.  Make sure encode is ready */
  46   oggpack_reset(&vb->opb);
  47   /* Encode the packet type */
  48   oggpack_write(&vb->opb,0,1);
  49
  50   /* currently lazy.  Short block dispatches to 0, long to 1. */
  51
  52   if(vb->W &&ci->modes>1)mode=1;
  53   type=ci->map_type[ci->mode_param[mode]->mapping];
  54   vb->mode=mode;
  55
  56   /* Encode frame mode, pre,post windowsize, then dispatch */
  57   oggpack_write(&vb->opb,mode,b->modebits);
  58   if(vb->W){
  59     oggpack_write(&vb->opb,vb->lW,1);
  60     oggpack_write(&vb->opb,vb->nW,1);
  61     /*fprintf(stderr,"*");
  62   }else{
  63   fprintf(stderr,".");*/
  64   }
  65
  66   if((ret=_mapping_P[type]->forward(vb,b->mode[mode])))
  67     return(ret);
  68
  69   if(op){
  70     op->packet=oggpack_get_buffer(&vb->opb);
  71     op->bytes=oggpack_bytes(&vb->opb);
  72     op->b_o_s=0;
  73     op->e_o_s=vb->eofflag;
  74     op->granulepos=vb->granulepos;
  75     op->packetno=vb->sequence; /* for sake of completeness */
  76   }
  77   return(0);
  78 }
  79
  80 /* there was no great place to put this.... */
  81 void _analysis_output_always(char *base,int i,float *v,int n,int bark,int dB){
  82   int j;
  83   FILE *of;
  84   char buffer[80];
  85
  86   /*  if(i==5870){*/
  87     sprintf(buffer,"%s_%d.m",base,i);
  88     of=fopen(buffer,"w");
  89
  90     if(!of)perror("failed to open data dump file");
  91
  92     for(j=0;j<n;j++){
  93       if(dB && v[j]==0)
  94         fprintf(of,"\n\n");
  95       else{
  96         if(bark)
  97           fprintf(of,"%g ",toBARK(22050.f*j/n));
  98         else
  99           fprintf(of,"%g ",(double)j);
 100
 101         if(dB){
 102           fprintf(of,"%g\n",todB(v+j));
 103         }else{
 104           fprintf(of,"%g\n",v[j]);
 105         }
 106       }
 107     }
 108     fclose(of);
 109     /*  } */
 110 }
 111
 112 void _analysis_output(char *base,int i,float *v,int n,int bark,int dB){
 113 #ifdef ANALYSIS
 114   if(analysis_noisy)_analysis_output_always(base,i,v,n,bark,dB);
 115 #endif
 116 }
 117