/********************************************************************
* *
- * THIS FILE IS PART OF THE Ogg Vorbis SOFTWARE CODEC SOURCE CODE. *
- * USE, DISTRIBUTION AND REPRODUCTION OF THIS SOURCE IS GOVERNED BY *
- * THE GNU PUBLIC LICENSE 2, WHICH IS INCLUDED WITH THIS SOURCE. *
- * PLEASE READ THESE TERMS DISTRIBUTING. *
+ * 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 OggSQUISH SOURCE CODE IS (C) COPYRIGHT 1994-1999 *
- * by 1999 Monty <monty@xiph.org> and The XIPHOPHORUS Company *
- * http://www.xiph.org/ *
+ * THE OggVorbis SOURCE CODE IS (C) COPYRIGHT 1994-2007 *
+ * by the Xiph.Org Foundation http://www.xiph.org/ *
* *
********************************************************************
- function: single-block PCM analysis
- author: Monty <xiphmont@mit.edu>
- modifications by: Monty
- last modification date: Oct 15 1999
+ function: single-block PCM analysis mode dispatch
+ last mod: $Id$
********************************************************************/
#include <stdio.h>
#include <string.h>
#include <math.h>
-#include "lpc.h"
-#include "lsp.h"
-#include "codec.h"
-#include "envelope.h"
-#include "mdct.h"
-#include "psy.h"
-#include "bitwise.h"
-#include "spectrum.h"
-
-/* this code is still seriously abbreviated. I'm filling in pieces as
- we go... --Monty 19991004 */
-
-int vorbis_analysis(vorbis_block *vb,ogg_packet *op){
- int i;
- double *window=vb->vd->window[vb->W][vb->lW][vb->nW];
- psy_lookup *vp=&vb->vd->vp[vb->W];
- lpc_lookup *vl=&vb->vd->vl[vb->W];
- vorbis_dsp_state *vd=vb->vd;
- vorbis_info *vi=vd->vi;
- oggpack_buffer *opb=&vb->opb;
-
- int n=vb->pcmend;
- int spectral_order=vi->floororder[vb->W];
-
- vb->gluebits=0;
- vb->time_envelope_bits=0;
- vb->spectral_envelope_bits=0;
- vb->spectral_residue_bits=0;
-
- /*lpc_lookup *vbal=&vb->vd->vbal[vb->W];
- double balance_v[vbal->m];
- double balance_amp;*/
-
- /* first things first. Make sure encode is ready*/
- _oggpack_reset(opb);
- /* Encode the packet type */
- _oggpack_write(opb,0,1);
- /* Encode the block size */
- _oggpack_write(opb,vb->W,1);
-
- /* we have the preecho metrics; decide what to do with them */
- _ve_envelope_sparsify(vb);
- _ve_envelope_apply(vb,0);
-
- /* Encode the envelope */
- if(_ve_envelope_encode(vb))return(-1);
-
- /* time domain PCM -> MDCT domain */
- for(i=0;i<vi->channels;i++)
- mdct_forward(&vd->vm[vb->W],vb->pcm[i],vb->pcm[i],window);
-
- /* no balance yet */
-
- /* extract the spectral envelope and residue */
- /* just do by channel. No coupling yet */
- {
- for(i=0;i<vi->channels;i++){
- static int frameno=0;
- int j;
- double floor[n/2];
- double curve[n/2];
- double *lpc=vb->lpc[i];
- double *lsp=vb->lsp[i];
-
- memset(floor,0,sizeof(double)*n/2);
-
- _vp_noise_floor(vp,vb->pcm[i],floor);
+#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<PACKETBLOBS;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
- {
- FILE *out;
- char buffer[80];
-
- sprintf(buffer,"spectrum.m");
- out=fopen(buffer,"w+");
- for(j=0;j<n/2;j++)
- fprintf(out,"%g\n",vb->pcm[i][j]);
- fclose(out);
-
- sprintf(buffer,"noise.m");
- out=fopen(buffer,"w+");
- for(j=0;j<n/2;j++)
- fprintf(out,"%g\n",floor[j]);
- fclose(out);
- }
-#endif
+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<n;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);
+}
- _vp_mask_floor(vp,vb->pcm[i],floor);
+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);
+}
-#ifdef ANALYSIS
- {
- FILE *out;
- char buffer[80];
-
- sprintf(buffer,"premask.m");
- out=fopen(buffer,"w+");
- for(j=0;j<n/2;j++)
- fprintf(out,"%g\n",floor[j]);
- fclose(out);
- }
#endif
- /* Convert our floor to a set of lpc coefficients */
- vb->amp[i]=sqrt(vorbis_curve_to_lpc(floor,lpc,vl));
- /* LSP <-> LPC is orthogonal and LSP quantizes more stably */
- vorbis_lpc_to_lsp(lpc,lsp,vl->m);
- /* code the spectral envelope; mutates the lsp coeffs to reflect
- what was actually encoded */
- _vs_spectrum_encode(vb,vb->amp[i],lsp);
- /* Generate residue from the decoded envelope, which will be
- slightly different to the pre-encoding floor due to
- quantization. Slow, yes, but perhaps more accurate */
- vorbis_lsp_to_lpc(lsp,lpc,vl->m);
- vorbis_lpc_to_curve(curve,lpc,vb->amp[i],vl);
-
- /* this may do various interesting massaging too...*/
- _vs_residue_quantize(vb->pcm[i],curve,vi,n/2);
-#ifdef ANALYSIS
- {
- FILE *out;
- char buffer[80];
-
- sprintf(buffer,"mask.m");
- out=fopen(buffer,"w+");
- for(j=0;j<n/2;j++)
- fprintf(out,"%g\n",curve[j]);
- fclose(out);
-
- sprintf(buffer,"res.m");
- out=fopen(buffer,"w+");
- for(j=0;j<n/2;j++)
- fprintf(out,"%g\n",vb->pcm[i][j]);
- fclose(out);
- }
-#endif
- /* encode the residue */
- _vs_residue_encode(vb,vb->pcm[i]);
-
- }
- }
-
- /* set up the packet wrapper */
-
- op->packet=opb->buffer;
- op->bytes=_oggpack_bytes(opb);
- op->b_o_s=0;
- op->e_o_s=vb->eofflag;
- op->frameno=vb->frameno;
-
- return(0);
-}
-/* commented out, relocated balance stuff */
- /*{
- double *C=vb->pcm[0];
- double *D=vb->pcm[1];
-
- balance_amp=_vp_balance_compute(D,C,balance_v,vbal);
-
- {
- FILE *out;
- char buffer[80];
-
- sprintf(buffer,"com%d.m",frameno);
- out=fopen(buffer,"w+");
- for(i=0;i<n/2;i++){
- fprintf(out," 0. 0.\n");
- fprintf(out,"%g %g\n",C[i],D[i]);
- fprintf(out,"\n");
- }
- fclose(out);
-
- sprintf(buffer,"L%d.m",frameno);
- out=fopen(buffer,"w+");
- for(i=0;i<n/2;i++){
- fprintf(out,"%g\n",C[i]);
- }
- fclose(out);
- sprintf(buffer,"R%d.m",frameno);
- out=fopen(buffer,"w+");
- for(i=0;i<n/2;i++){
- fprintf(out,"%g\n",D[i]);
- }
- fclose(out);
-
- }
-
- _vp_balance_apply(D,C,balance_v,balance_amp,vbal,1);
-
- {
- FILE *out;
- char buffer[80];
-
- sprintf(buffer,"bal%d.m",frameno);
- out=fopen(buffer,"w+");
- for(i=0;i<n/2;i++){
- fprintf(out," 0. 0.\n");
- fprintf(out,"%g %g\n",C[i],D[i]);
- fprintf(out,"\n");
- }
- fclose(out);
- sprintf(buffer,"C%d.m",frameno);
- out=fopen(buffer,"w+");
- for(i=0;i<n/2;i++){
- fprintf(out,"%g\n",C[i]);
- }
- fclose(out);
- sprintf(buffer,"D%d.m",frameno);
- out=fopen(buffer,"w+");
- for(i=0;i<n/2;i++){
- fprintf(out,"%g\n",D[i]);
- }
- fclose(out);
-
- }
- }*/