1 /********************************************************************
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. *
8 * THE OggVorbis SOURCE CODE IS (C) COPYRIGHT 1994-2001 *
9 * by the XIPHOPHORUS Company http://www.xiph.org/ *
11 ********************************************************************
13 function: single-block PCM analysis mode dispatch
14 last mod: $Id: analysis.c,v 1.47 2001/12/12 09:45:24 xiphmont Exp $
16 ********************************************************************/
22 #include "vorbis/codec.h"
23 #include "codec_internal.h"
30 /* decides between modes, dispatches to the appropriate mapping. */
31 int vorbis_analysis(vorbis_block *vb, ogg_packet *op){
32 vorbis_dsp_state *vd=vb->vd;
33 backend_lookup_state *b=vd->backend_state;
34 vorbis_info *vi=vd->vi;
35 codec_setup_info *ci=vi->codec_setup;
44 /* first things first. Make sure encode is ready */
45 oggpack_reset(&vb->opb);
46 /* Encode the packet type */
47 oggpack_write(&vb->opb,0,1);
49 /* currently lazy. Short block dispatches to 0, long to 1. */
51 if(vb->W &&ci->modes>1)mode=1;
52 type=ci->map_type[ci->mode_param[mode]->mapping];
55 /* Encode frame mode, pre,post windowsize, then dispatch */
56 oggpack_write(&vb->opb,mode,b->modebits);
58 oggpack_write(&vb->opb,vb->lW,1);
59 oggpack_write(&vb->opb,vb->nW,1);
60 /*fprintf(stderr,"*");
62 fprintf(stderr,".");*/
65 if((ret=_mapping_P[type]->forward(vb,b->mode[mode])))
69 op->packet=oggpack_get_buffer(&vb->opb);
70 op->bytes=oggpack_bytes(&vb->opb);
72 op->e_o_s=vb->eofflag;
73 op->granulepos=vb->granulepos;
74 op->packetno=vb->sequence; /* for sake of completeness */
79 /* there was no great place to put this.... */
80 void _analysis_output_always(char *base,int i,float *v,int n,int bark,int dB){
86 sprintf(buffer,"%s_%d.m",base,i);
89 if(!of)perror("failed to open data dump file");
96 fprintf(of,"%g ",toBARK(22050.f*j/n));
98 fprintf(of,"%g ",(double)j);
101 fprintf(of,"%g\n",todB(v+j));
103 fprintf(of,"%g\n",v[j]);
111 void _analysis_output(char *base,int i,float *v,int n,int bark,int dB){
113 if(analysis_noisy)_analysis_output_always(base,i,v,n,bark,dB);