/********************************************************************

 *                                                                  *

 * THIS FILE IS PART OF THE OggVorbis SOFTWARE CODEC SOURCE CODE.   *

 * USE, DISTRIBUTION AND REPRODUCTION OF THIS LIBRARY SOURCE IS     *

 * GOVERNED BY A BSD-STYLE SOURCE LICENSE INCLUDED WITH THIS SOURCE *

 * IN 'COPYING'. PLEASE READ THESE TERMS BEFORE DISTRIBUTING.       *

 *                                                                  *

 * THE OggVorbis SOURCE CODE IS (C) COPYRIGHT 1994-2007             *

 * by the Xiph.Org Foundation http://www.xiph.org/                  *

 *                                                                  *

 ********************************************************************

 function: single-block PCM analysis mode dispatch

 last mod: $Id: analysis.c 16226 2009-07-08 06:43:49Z xiphmont $

 ********************************************************************/

#include <stdio.h>

#include <string.h>

#include <math.h>

#include <ogg/ogg.h>

#include "vorbis/codec.h"

#include "codec_internal.h"

#include "registry.h"

#include "scales.h"

#include "os.h"

#include "misc.h"

/* decides between modes, dispatches to the appropriate mapping. */

int vorbis_analysis(vorbis_block *vb, ogg_packet *op){

  int ret,i;

  vorbis_block_internal *vbi=vb->internal;

  vb->glue_bits=0;

  vb->time_bits=0;

  vb->floor_bits=0;

  vb->res_bits=0;

  /* first things first.  Make sure encode is ready */

  for(i=0;i

    oggpack_reset(vbi->packetblob[i]);

  /* we only have one mapping type (0), and we let the mapping code

     itself figure out what soft mode to use.  This allows easier

     bitrate management */

  if((ret=_mapping_P[0]->forward(vb)))

    return(ret);

  if(op){

    if(vorbis_bitrate_managed(vb))

      /* The app is using a bitmanaged mode... but not using the

         bitrate management interface. */

      return(OV_EINVAL);

    op->packet=oggpack_get_buffer(&vb->opb);

    op->bytes=oggpack_bytes(&vb->opb);

    op->b_o_s=0;

    op->e_o_s=vb->eofflag;

    op->granulepos=vb->granulepos;

    op->packetno=vb->sequence; /* for sake of completeness */

  }

  return(0);

}

#ifdef ANALYSIS

int analysis_noisy=1;

/* there was no great place to put this.... */

void _analysis_output_always(char *base,int i,float *v,int n,int bark,int dB,ogg_int64_t off){

  int j;

  FILE *of;

  char buffer[80];

  sprintf(buffer,"%s_%d.m",base,i);

  of=fopen(buffer,"w");

  if(!of)perror("failed to open data dump file");

  for(j=0;j

    if(bark){

      float b=toBARK((4000.f*j/n)+.25);

      fprintf(of,"%f ",b);

    }else

      if(off!=0)

        fprintf(of,"%f ",(double)(j+off)/8000.);

      else

        fprintf(of,"%f ",(double)j);

    if(dB){

      float val;

      if(v[j]==0.)

        val=-140.;

      else

        val=todB(v+j);

      fprintf(of,"%f\n",val);

    }else{

      fprintf(of,"%f\n",v[j]);

    }

  }

  fclose(of);

}

void _analysis_output(char *base,int i,float *v,int n,int bark,int dB,

                      ogg_int64_t off){

  if(analysis_noisy)_analysis_output_always(base,i,v,n,bark,dB,off);

}

#endif