dbf-halloween2015

annotate libs/vorbis/mapping0.c @ 3:c37fe5d8a4ed

find changesets by author, revision, files, or words in the commit message
windows port
author John Tsiombikas <nuclear@member.fsf.org>
date Sun, 01 Nov 2015 06:04:28 +0200
parents
children
rev   line source
nuclear@1 1 /********************************************************************
nuclear@1 2 * *
nuclear@1 3 * THIS FILE IS PART OF THE OggVorbis SOFTWARE CODEC SOURCE CODE. *
nuclear@1 4 * USE, DISTRIBUTION AND REPRODUCTION OF THIS LIBRARY SOURCE IS *
nuclear@1 5 * GOVERNED BY A BSD-STYLE SOURCE LICENSE INCLUDED WITH THIS SOURCE *
nuclear@1 6 * IN 'COPYING'. PLEASE READ THESE TERMS BEFORE DISTRIBUTING. *
nuclear@1 7 * *
nuclear@1 8 * THE OggVorbis SOURCE CODE IS (C) COPYRIGHT 1994-2010 *
nuclear@1 9 * by the Xiph.Org Foundation http://www.xiph.org/ *
nuclear@1 10 * *
nuclear@1 11 ********************************************************************
nuclear@1 12
nuclear@1 13 function: channel mapping 0 implementation
nuclear@1 14 last mod: $Id: mapping0.c 17022 2010-03-25 03:45:42Z xiphmont $
nuclear@1 15
nuclear@1 16 ********************************************************************/
nuclear@1 17
nuclear@1 18 #include <stdlib.h>
nuclear@1 19 #include <stdio.h>
nuclear@1 20 #include <string.h>
nuclear@1 21 #include <math.h>
nuclear@1 22 #include <ogg/ogg.h>
nuclear@1 23 #include "vorbis/codec.h"
nuclear@1 24 #include "codec_internal.h"
nuclear@1 25 #include "codebook.h"
nuclear@1 26 #include "window.h"
nuclear@1 27 #include "registry.h"
nuclear@1 28 #include "psy.h"
nuclear@1 29 #include "misc.h"
nuclear@1 30
nuclear@1 31 /* simplistic, wasteful way of doing this (unique lookup for each
nuclear@1 32 mode/submapping); there should be a central repository for
nuclear@1 33 identical lookups. That will require minor work, so I'm putting it
nuclear@1 34 off as low priority.
nuclear@1 35
nuclear@1 36 Why a lookup for each backend in a given mode? Because the
nuclear@1 37 blocksize is set by the mode, and low backend lookups may require
nuclear@1 38 parameters from other areas of the mode/mapping */
nuclear@1 39
nuclear@1 40 static void mapping0_free_info(vorbis_info_mapping *i){
nuclear@1 41 vorbis_info_mapping0 *info=(vorbis_info_mapping0 *)i;
nuclear@1 42 if(info){
nuclear@1 43 memset(info,0,sizeof(*info));
nuclear@1 44 _ogg_free(info);
nuclear@1 45 }
nuclear@1 46 }
nuclear@1 47
nuclear@1 48 static int ilog(unsigned int v){
nuclear@1 49 int ret=0;
nuclear@1 50 if(v)--v;
nuclear@1 51 while(v){
nuclear@1 52 ret++;
nuclear@1 53 v>>=1;
nuclear@1 54 }
nuclear@1 55 return(ret);
nuclear@1 56 }
nuclear@1 57
nuclear@1 58 static void mapping0_pack(vorbis_info *vi,vorbis_info_mapping *vm,
nuclear@1 59 oggpack_buffer *opb){
nuclear@1 60 int i;
nuclear@1 61 vorbis_info_mapping0 *info=(vorbis_info_mapping0 *)vm;
nuclear@1 62
nuclear@1 63 /* another 'we meant to do it this way' hack... up to beta 4, we
nuclear@1 64 packed 4 binary zeros here to signify one submapping in use. We
nuclear@1 65 now redefine that to mean four bitflags that indicate use of
nuclear@1 66 deeper features; bit0:submappings, bit1:coupling,
nuclear@1 67 bit2,3:reserved. This is backward compatable with all actual uses
nuclear@1 68 of the beta code. */
nuclear@1 69
nuclear@1 70 if(info->submaps>1){
nuclear@1 71 oggpack_write(opb,1,1);
nuclear@1 72 oggpack_write(opb,info->submaps-1,4);
nuclear@1 73 }else
nuclear@1 74 oggpack_write(opb,0,1);
nuclear@1 75
nuclear@1 76 if(info->coupling_steps>0){
nuclear@1 77 oggpack_write(opb,1,1);
nuclear@1 78 oggpack_write(opb,info->coupling_steps-1,8);
nuclear@1 79
nuclear@1 80 for(i=0;i<info->coupling_steps;i++){
nuclear@1 81 oggpack_write(opb,info->coupling_mag[i],ilog(vi->channels));
nuclear@1 82 oggpack_write(opb,info->coupling_ang[i],ilog(vi->channels));
nuclear@1 83 }
nuclear@1 84 }else
nuclear@1 85 oggpack_write(opb,0,1);
nuclear@1 86
nuclear@1 87 oggpack_write(opb,0,2); /* 2,3:reserved */
nuclear@1 88
nuclear@1 89 /* we don't write the channel submappings if we only have one... */
nuclear@1 90 if(info->submaps>1){
nuclear@1 91 for(i=0;i<vi->channels;i++)
nuclear@1 92 oggpack_write(opb,info->chmuxlist[i],4);
nuclear@1 93 }
nuclear@1 94 for(i=0;i<info->submaps;i++){
nuclear@1 95 oggpack_write(opb,0,8); /* time submap unused */
nuclear@1 96 oggpack_write(opb,info->floorsubmap[i],8);
nuclear@1 97 oggpack_write(opb,info->residuesubmap[i],8);
nuclear@1 98 }
nuclear@1 99 }
nuclear@1 100
nuclear@1 101 /* also responsible for range checking */
nuclear@1 102 static vorbis_info_mapping *mapping0_unpack(vorbis_info *vi,oggpack_buffer *opb){
nuclear@1 103 int i,b;
nuclear@1 104 vorbis_info_mapping0 *info=_ogg_calloc(1,sizeof(*info));
nuclear@1 105 codec_setup_info *ci=vi->codec_setup;
nuclear@1 106 memset(info,0,sizeof(*info));
nuclear@1 107
nuclear@1 108 b=oggpack_read(opb,1);
nuclear@1 109 if(b<0)goto err_out;
nuclear@1 110 if(b){
nuclear@1 111 info->submaps=oggpack_read(opb,4)+1;
nuclear@1 112 if(info->submaps<=0)goto err_out;
nuclear@1 113 }else
nuclear@1 114 info->submaps=1;
nuclear@1 115
nuclear@1 116 b=oggpack_read(opb,1);
nuclear@1 117 if(b<0)goto err_out;
nuclear@1 118 if(b){
nuclear@1 119 info->coupling_steps=oggpack_read(opb,8)+1;
nuclear@1 120 if(info->coupling_steps<=0)goto err_out;
nuclear@1 121 for(i=0;i<info->coupling_steps;i++){
nuclear@1 122 int testM=info->coupling_mag[i]=oggpack_read(opb,ilog(vi->channels));
nuclear@1 123 int testA=info->coupling_ang[i]=oggpack_read(opb,ilog(vi->channels));
nuclear@1 124
nuclear@1 125 if(testM<0 ||
nuclear@1 126 testA<0 ||
nuclear@1 127 testM==testA ||
nuclear@1 128 testM>=vi->channels ||
nuclear@1 129 testA>=vi->channels) goto err_out;
nuclear@1 130 }
nuclear@1 131
nuclear@1 132 }
nuclear@1 133
nuclear@1 134 if(oggpack_read(opb,2)!=0)goto err_out; /* 2,3:reserved */
nuclear@1 135
nuclear@1 136 if(info->submaps>1){
nuclear@1 137 for(i=0;i<vi->channels;i++){
nuclear@1 138 info->chmuxlist[i]=oggpack_read(opb,4);
nuclear@1 139 if(info->chmuxlist[i]>=info->submaps || info->chmuxlist[i]<0)goto err_out;
nuclear@1 140 }
nuclear@1 141 }
nuclear@1 142 for(i=0;i<info->submaps;i++){
nuclear@1 143 oggpack_read(opb,8); /* time submap unused */
nuclear@1 144 info->floorsubmap[i]=oggpack_read(opb,8);
nuclear@1 145 if(info->floorsubmap[i]>=ci->floors || info->floorsubmap[i]<0)goto err_out;
nuclear@1 146 info->residuesubmap[i]=oggpack_read(opb,8);
nuclear@1 147 if(info->residuesubmap[i]>=ci->residues || info->residuesubmap[i]<0)goto err_out;
nuclear@1 148 }
nuclear@1 149
nuclear@1 150 return info;
nuclear@1 151
nuclear@1 152 err_out:
nuclear@1 153 mapping0_free_info(info);
nuclear@1 154 return(NULL);
nuclear@1 155 }
nuclear@1 156
nuclear@1 157 #include "os.h"
nuclear@1 158 #include "lpc.h"
nuclear@1 159 #include "lsp.h"
nuclear@1 160 #include "envelope.h"
nuclear@1 161 #include "mdct.h"
nuclear@1 162 #include "psy.h"
nuclear@1 163 #include "scales.h"
nuclear@1 164
nuclear@1 165 #if 0
nuclear@1 166 static long seq=0;
nuclear@1 167 static ogg_int64_t total=0;
nuclear@1 168 static float FLOOR1_fromdB_LOOKUP[256]={
nuclear@1 169 1.0649863e-07F, 1.1341951e-07F, 1.2079015e-07F, 1.2863978e-07F,
nuclear@1 170 1.3699951e-07F, 1.4590251e-07F, 1.5538408e-07F, 1.6548181e-07F,
nuclear@1 171 1.7623575e-07F, 1.8768855e-07F, 1.9988561e-07F, 2.128753e-07F,
nuclear@1 172 2.2670913e-07F, 2.4144197e-07F, 2.5713223e-07F, 2.7384213e-07F,
nuclear@1 173 2.9163793e-07F, 3.1059021e-07F, 3.3077411e-07F, 3.5226968e-07F,
nuclear@1 174 3.7516214e-07F, 3.9954229e-07F, 4.2550680e-07F, 4.5315863e-07F,
nuclear@1 175 4.8260743e-07F, 5.1396998e-07F, 5.4737065e-07F, 5.8294187e-07F,
nuclear@1 176 6.2082472e-07F, 6.6116941e-07F, 7.0413592e-07F, 7.4989464e-07F,
nuclear@1 177 7.9862701e-07F, 8.5052630e-07F, 9.0579828e-07F, 9.6466216e-07F,
nuclear@1 178 1.0273513e-06F, 1.0941144e-06F, 1.1652161e-06F, 1.2409384e-06F,
nuclear@1 179 1.3215816e-06F, 1.4074654e-06F, 1.4989305e-06F, 1.5963394e-06F,
nuclear@1 180 1.7000785e-06F, 1.8105592e-06F, 1.9282195e-06F, 2.0535261e-06F,
nuclear@1 181 2.1869758e-06F, 2.3290978e-06F, 2.4804557e-06F, 2.6416497e-06F,
nuclear@1 182 2.8133190e-06F, 2.9961443e-06F, 3.1908506e-06F, 3.3982101e-06F,
nuclear@1 183 3.6190449e-06F, 3.8542308e-06F, 4.1047004e-06F, 4.3714470e-06F,
nuclear@1 184 4.6555282e-06F, 4.9580707e-06F, 5.2802740e-06F, 5.6234160e-06F,
nuclear@1 185 5.9888572e-06F, 6.3780469e-06F, 6.7925283e-06F, 7.2339451e-06F,
nuclear@1 186 7.7040476e-06F, 8.2047000e-06F, 8.7378876e-06F, 9.3057248e-06F,
nuclear@1 187 9.9104632e-06F, 1.0554501e-05F, 1.1240392e-05F, 1.1970856e-05F,
nuclear@1 188 1.2748789e-05F, 1.3577278e-05F, 1.4459606e-05F, 1.5399272e-05F,
nuclear@1 189 1.6400004e-05F, 1.7465768e-05F, 1.8600792e-05F, 1.9809576e-05F,
nuclear@1 190 2.1096914e-05F, 2.2467911e-05F, 2.3928002e-05F, 2.5482978e-05F,
nuclear@1 191 2.7139006e-05F, 2.8902651e-05F, 3.0780908e-05F, 3.2781225e-05F,
nuclear@1 192 3.4911534e-05F, 3.7180282e-05F, 3.9596466e-05F, 4.2169667e-05F,
nuclear@1 193 4.4910090e-05F, 4.7828601e-05F, 5.0936773e-05F, 5.4246931e-05F,
nuclear@1 194 5.7772202e-05F, 6.1526565e-05F, 6.5524908e-05F, 6.9783085e-05F,
nuclear@1 195 7.4317983e-05F, 7.9147585e-05F, 8.4291040e-05F, 8.9768747e-05F,
nuclear@1 196 9.5602426e-05F, 0.00010181521F, 0.00010843174F, 0.00011547824F,
nuclear@1 197 0.00012298267F, 0.00013097477F, 0.00013948625F, 0.00014855085F,
nuclear@1 198 0.00015820453F, 0.00016848555F, 0.00017943469F, 0.00019109536F,
nuclear@1 199 0.00020351382F, 0.00021673929F, 0.00023082423F, 0.00024582449F,
nuclear@1 200 0.00026179955F, 0.00027881276F, 0.00029693158F, 0.00031622787F,
nuclear@1 201 0.00033677814F, 0.00035866388F, 0.00038197188F, 0.00040679456F,
nuclear@1 202 0.00043323036F, 0.00046138411F, 0.00049136745F, 0.00052329927F,
nuclear@1 203 0.00055730621F, 0.00059352311F, 0.00063209358F, 0.00067317058F,
nuclear@1 204 0.00071691700F, 0.00076350630F, 0.00081312324F, 0.00086596457F,
nuclear@1 205 0.00092223983F, 0.00098217216F, 0.0010459992F, 0.0011139742F,
nuclear@1 206 0.0011863665F, 0.0012634633F, 0.0013455702F, 0.0014330129F,
nuclear@1 207 0.0015261382F, 0.0016253153F, 0.0017309374F, 0.0018434235F,
nuclear@1 208 0.0019632195F, 0.0020908006F, 0.0022266726F, 0.0023713743F,
nuclear@1 209 0.0025254795F, 0.0026895994F, 0.0028643847F, 0.0030505286F,
nuclear@1 210 0.0032487691F, 0.0034598925F, 0.0036847358F, 0.0039241906F,
nuclear@1 211 0.0041792066F, 0.0044507950F, 0.0047400328F, 0.0050480668F,
nuclear@1 212 0.0053761186F, 0.0057254891F, 0.0060975636F, 0.0064938176F,
nuclear@1 213 0.0069158225F, 0.0073652516F, 0.0078438871F, 0.0083536271F,
nuclear@1 214 0.0088964928F, 0.009474637F, 0.010090352F, 0.010746080F,
nuclear@1 215 0.011444421F, 0.012188144F, 0.012980198F, 0.013823725F,
nuclear@1 216 0.014722068F, 0.015678791F, 0.016697687F, 0.017782797F,
nuclear@1 217 0.018938423F, 0.020169149F, 0.021479854F, 0.022875735F,
nuclear@1 218 0.024362330F, 0.025945531F, 0.027631618F, 0.029427276F,
nuclear@1 219 0.031339626F, 0.033376252F, 0.035545228F, 0.037855157F,
nuclear@1 220 0.040315199F, 0.042935108F, 0.045725273F, 0.048696758F,
nuclear@1 221 0.051861348F, 0.055231591F, 0.058820850F, 0.062643361F,
nuclear@1 222 0.066714279F, 0.071049749F, 0.075666962F, 0.080584227F,
nuclear@1 223 0.085821044F, 0.091398179F, 0.097337747F, 0.10366330F,
nuclear@1 224 0.11039993F, 0.11757434F, 0.12521498F, 0.13335215F,
nuclear@1 225 0.14201813F, 0.15124727F, 0.16107617F, 0.17154380F,
nuclear@1 226 0.18269168F, 0.19456402F, 0.20720788F, 0.22067342F,
nuclear@1 227 0.23501402F, 0.25028656F, 0.26655159F, 0.28387361F,
nuclear@1 228 0.30232132F, 0.32196786F, 0.34289114F, 0.36517414F,
nuclear@1 229 0.38890521F, 0.41417847F, 0.44109412F, 0.46975890F,
nuclear@1 230 0.50028648F, 0.53279791F, 0.56742212F, 0.60429640F,
nuclear@1 231 0.64356699F, 0.68538959F, 0.72993007F, 0.77736504F,
nuclear@1 232 0.82788260F, 0.88168307F, 0.9389798F, 1.F,
nuclear@1 233 };
nuclear@1 234
nuclear@1 235 #endif
nuclear@1 236
nuclear@1 237
nuclear@1 238 static int mapping0_forward(vorbis_block *vb){
nuclear@1 239 vorbis_dsp_state *vd=vb->vd;
nuclear@1 240 vorbis_info *vi=vd->vi;
nuclear@1 241 codec_setup_info *ci=vi->codec_setup;
nuclear@1 242 private_state *b=vb->vd->backend_state;
nuclear@1 243 vorbis_block_internal *vbi=(vorbis_block_internal *)vb->internal;
nuclear@1 244 int n=vb->pcmend;
nuclear@1 245 int i,j,k;
nuclear@1 246
nuclear@1 247 int *nonzero = alloca(sizeof(*nonzero)*vi->channels);
nuclear@1 248 float **gmdct = _vorbis_block_alloc(vb,vi->channels*sizeof(*gmdct));
nuclear@1 249 int **iwork = _vorbis_block_alloc(vb,vi->channels*sizeof(*iwork));
nuclear@1 250 int ***floor_posts = _vorbis_block_alloc(vb,vi->channels*sizeof(*floor_posts));
nuclear@1 251
nuclear@1 252 float global_ampmax=vbi->ampmax;
nuclear@1 253 float *local_ampmax=alloca(sizeof(*local_ampmax)*vi->channels);
nuclear@1 254 int blocktype=vbi->blocktype;
nuclear@1 255
nuclear@1 256 int modenumber=vb->W;
nuclear@1 257 vorbis_info_mapping0 *info=ci->map_param[modenumber];
nuclear@1 258 vorbis_look_psy *psy_look=b->psy+blocktype+(vb->W?2:0);
nuclear@1 259
nuclear@1 260 vb->mode=modenumber;
nuclear@1 261
nuclear@1 262 for(i=0;i<vi->channels;i++){
nuclear@1 263 float scale=4.f/n;
nuclear@1 264 float scale_dB;
nuclear@1 265
nuclear@1 266 float *pcm =vb->pcm[i];
nuclear@1 267 float *logfft =pcm;
nuclear@1 268
nuclear@1 269 iwork[i]=_vorbis_block_alloc(vb,n/2*sizeof(**iwork));
nuclear@1 270 gmdct[i]=_vorbis_block_alloc(vb,n/2*sizeof(**gmdct));
nuclear@1 271
nuclear@1 272 scale_dB=todB(&scale) + .345; /* + .345 is a hack; the original
nuclear@1 273 todB estimation used on IEEE 754
nuclear@1 274 compliant machines had a bug that
nuclear@1 275 returned dB values about a third
nuclear@1 276 of a decibel too high. The bug
nuclear@1 277 was harmless because tunings
nuclear@1 278 implicitly took that into
nuclear@1 279 account. However, fixing the bug
nuclear@1 280 in the estimator requires
nuclear@1 281 changing all the tunings as well.
nuclear@1 282 For now, it's easier to sync
nuclear@1 283 things back up here, and
nuclear@1 284 recalibrate the tunings in the
nuclear@1 285 next major model upgrade. */
nuclear@1 286
nuclear@1 287 #if 0
nuclear@1 288 if(vi->channels==2){
nuclear@1 289 if(i==0)
nuclear@1 290 _analysis_output("pcmL",seq,pcm,n,0,0,total-n/2);
nuclear@1 291 else
nuclear@1 292 _analysis_output("pcmR",seq,pcm,n,0,0,total-n/2);
nuclear@1 293 }else{
nuclear@1 294 _analysis_output("pcm",seq,pcm,n,0,0,total-n/2);
nuclear@1 295 }
nuclear@1 296 #endif
nuclear@1 297
nuclear@1 298 /* window the PCM data */
nuclear@1 299 _vorbis_apply_window(pcm,b->window,ci->blocksizes,vb->lW,vb->W,vb->nW);
nuclear@1 300
nuclear@1 301 #if 0
nuclear@1 302 if(vi->channels==2){
nuclear@1 303 if(i==0)
nuclear@1 304 _analysis_output("windowedL",seq,pcm,n,0,0,total-n/2);
nuclear@1 305 else
nuclear@1 306 _analysis_output("windowedR",seq,pcm,n,0,0,total-n/2);
nuclear@1 307 }else{
nuclear@1 308 _analysis_output("windowed",seq,pcm,n,0,0,total-n/2);
nuclear@1 309 }
nuclear@1 310 #endif
nuclear@1 311
nuclear@1 312 /* transform the PCM data */
nuclear@1 313 /* only MDCT right now.... */
nuclear@1 314 mdct_forward(b->transform[vb->W][0],pcm,gmdct[i]);
nuclear@1 315
nuclear@1 316 /* FFT yields more accurate tonal estimation (not phase sensitive) */
nuclear@1 317 drft_forward(&b->fft_look[vb->W],pcm);
nuclear@1 318 logfft[0]=scale_dB+todB(pcm) + .345; /* + .345 is a hack; the
nuclear@1 319 original todB estimation used on
nuclear@1 320 IEEE 754 compliant machines had a
nuclear@1 321 bug that returned dB values about
nuclear@1 322 a third of a decibel too high.
nuclear@1 323 The bug was harmless because
nuclear@1 324 tunings implicitly took that into
nuclear@1 325 account. However, fixing the bug
nuclear@1 326 in the estimator requires
nuclear@1 327 changing all the tunings as well.
nuclear@1 328 For now, it's easier to sync
nuclear@1 329 things back up here, and
nuclear@1 330 recalibrate the tunings in the
nuclear@1 331 next major model upgrade. */
nuclear@1 332 local_ampmax[i]=logfft[0];
nuclear@1 333 for(j=1;j<n-1;j+=2){
nuclear@1 334 float temp=pcm[j]*pcm[j]+pcm[j+1]*pcm[j+1];
nuclear@1 335 temp=logfft[(j+1)>>1]=scale_dB+.5f*todB(&temp) + .345; /* +
nuclear@1 336 .345 is a hack; the original todB
nuclear@1 337 estimation used on IEEE 754
nuclear@1 338 compliant machines had a bug that
nuclear@1 339 returned dB values about a third
nuclear@1 340 of a decibel too high. The bug
nuclear@1 341 was harmless because tunings
nuclear@1 342 implicitly took that into
nuclear@1 343 account. However, fixing the bug
nuclear@1 344 in the estimator requires
nuclear@1 345 changing all the tunings as well.
nuclear@1 346 For now, it's easier to sync
nuclear@1 347 things back up here, and
nuclear@1 348 recalibrate the tunings in the
nuclear@1 349 next major model upgrade. */
nuclear@1 350 if(temp>local_ampmax[i])local_ampmax[i]=temp;
nuclear@1 351 }
nuclear@1 352
nuclear@1 353 if(local_ampmax[i]>0.f)local_ampmax[i]=0.f;
nuclear@1 354 if(local_ampmax[i]>global_ampmax)global_ampmax=local_ampmax[i];
nuclear@1 355
nuclear@1 356 #if 0
nuclear@1 357 if(vi->channels==2){
nuclear@1 358 if(i==0){
nuclear@1 359 _analysis_output("fftL",seq,logfft,n/2,1,0,0);
nuclear@1 360 }else{
nuclear@1 361 _analysis_output("fftR",seq,logfft,n/2,1,0,0);
nuclear@1 362 }
nuclear@1 363 }else{
nuclear@1 364 _analysis_output("fft",seq,logfft,n/2,1,0,0);
nuclear@1 365 }
nuclear@1 366 #endif
nuclear@1 367
nuclear@1 368 }
nuclear@1 369
nuclear@1 370 {
nuclear@1 371 float *noise = _vorbis_block_alloc(vb,n/2*sizeof(*noise));
nuclear@1 372 float *tone = _vorbis_block_alloc(vb,n/2*sizeof(*tone));
nuclear@1 373
nuclear@1 374 for(i=0;i<vi->channels;i++){
nuclear@1 375 /* the encoder setup assumes that all the modes used by any
nuclear@1 376 specific bitrate tweaking use the same floor */
nuclear@1 377
nuclear@1 378 int submap=info->chmuxlist[i];
nuclear@1 379
nuclear@1 380 /* the following makes things clearer to *me* anyway */
nuclear@1 381 float *mdct =gmdct[i];
nuclear@1 382 float *logfft =vb->pcm[i];
nuclear@1 383
nuclear@1 384 float *logmdct =logfft+n/2;
nuclear@1 385 float *logmask =logfft;
nuclear@1 386
nuclear@1 387 vb->mode=modenumber;
nuclear@1 388
nuclear@1 389 floor_posts[i]=_vorbis_block_alloc(vb,PACKETBLOBS*sizeof(**floor_posts));
nuclear@1 390 memset(floor_posts[i],0,sizeof(**floor_posts)*PACKETBLOBS);
nuclear@1 391
nuclear@1 392 for(j=0;j<n/2;j++)
nuclear@1 393 logmdct[j]=todB(mdct+j) + .345; /* + .345 is a hack; the original
nuclear@1 394 todB estimation used on IEEE 754
nuclear@1 395 compliant machines had a bug that
nuclear@1 396 returned dB values about a third
nuclear@1 397 of a decibel too high. The bug
nuclear@1 398 was harmless because tunings
nuclear@1 399 implicitly took that into
nuclear@1 400 account. However, fixing the bug
nuclear@1 401 in the estimator requires
nuclear@1 402 changing all the tunings as well.
nuclear@1 403 For now, it's easier to sync
nuclear@1 404 things back up here, and
nuclear@1 405 recalibrate the tunings in the
nuclear@1 406 next major model upgrade. */
nuclear@1 407
nuclear@1 408 #if 0
nuclear@1 409 if(vi->channels==2){
nuclear@1 410 if(i==0)
nuclear@1 411 _analysis_output("mdctL",seq,logmdct,n/2,1,0,0);
nuclear@1 412 else
nuclear@1 413 _analysis_output("mdctR",seq,logmdct,n/2,1,0,0);
nuclear@1 414 }else{
nuclear@1 415 _analysis_output("mdct",seq,logmdct,n/2,1,0,0);
nuclear@1 416 }
nuclear@1 417 #endif
nuclear@1 418
nuclear@1 419 /* first step; noise masking. Not only does 'noise masking'
nuclear@1 420 give us curves from which we can decide how much resolution
nuclear@1 421 to give noise parts of the spectrum, it also implicitly hands
nuclear@1 422 us a tonality estimate (the larger the value in the
nuclear@1 423 'noise_depth' vector, the more tonal that area is) */
nuclear@1 424
nuclear@1 425 _vp_noisemask(psy_look,
nuclear@1 426 logmdct,
nuclear@1 427 noise); /* noise does not have by-frequency offset
nuclear@1 428 bias applied yet */
nuclear@1 429 #if 0
nuclear@1 430 if(vi->channels==2){
nuclear@1 431 if(i==0)
nuclear@1 432 _analysis_output("noiseL",seq,noise,n/2,1,0,0);
nuclear@1 433 else
nuclear@1 434 _analysis_output("noiseR",seq,noise,n/2,1,0,0);
nuclear@1 435 }else{
nuclear@1 436 _analysis_output("noise",seq,noise,n/2,1,0,0);
nuclear@1 437 }
nuclear@1 438 #endif
nuclear@1 439
nuclear@1 440 /* second step: 'all the other crap'; all the stuff that isn't
nuclear@1 441 computed/fit for bitrate management goes in the second psy
nuclear@1 442 vector. This includes tone masking, peak limiting and ATH */
nuclear@1 443
nuclear@1 444 _vp_tonemask(psy_look,
nuclear@1 445 logfft,
nuclear@1 446 tone,
nuclear@1 447 global_ampmax,
nuclear@1 448 local_ampmax[i]);
nuclear@1 449
nuclear@1 450 #if 0
nuclear@1 451 if(vi->channels==2){
nuclear@1 452 if(i==0)
nuclear@1 453 _analysis_output("toneL",seq,tone,n/2,1,0,0);
nuclear@1 454 else
nuclear@1 455 _analysis_output("toneR",seq,tone,n/2,1,0,0);
nuclear@1 456 }else{
nuclear@1 457 _analysis_output("tone",seq,tone,n/2,1,0,0);
nuclear@1 458 }
nuclear@1 459 #endif
nuclear@1 460
nuclear@1 461 /* third step; we offset the noise vectors, overlay tone
nuclear@1 462 masking. We then do a floor1-specific line fit. If we're
nuclear@1 463 performing bitrate management, the line fit is performed
nuclear@1 464 multiple times for up/down tweakage on demand. */
nuclear@1 465
nuclear@1 466 #if 0
nuclear@1 467 {
nuclear@1 468 float aotuv[psy_look->n];
nuclear@1 469 #endif
nuclear@1 470
nuclear@1 471 _vp_offset_and_mix(psy_look,
nuclear@1 472 noise,
nuclear@1 473 tone,
nuclear@1 474 1,
nuclear@1 475 logmask,
nuclear@1 476 mdct,
nuclear@1 477 logmdct);
nuclear@1 478
nuclear@1 479 #if 0
nuclear@1 480 if(vi->channels==2){
nuclear@1 481 if(i==0)
nuclear@1 482 _analysis_output("aotuvM1_L",seq,aotuv,psy_look->n,1,1,0);
nuclear@1 483 else
nuclear@1 484 _analysis_output("aotuvM1_R",seq,aotuv,psy_look->n,1,1,0);
nuclear@1 485 }else{
nuclear@1 486 _analysis_output("aotuvM1",seq,aotuv,psy_look->n,1,1,0);
nuclear@1 487 }
nuclear@1 488 }
nuclear@1 489 #endif
nuclear@1 490
nuclear@1 491
nuclear@1 492 #if 0
nuclear@1 493 if(vi->channels==2){
nuclear@1 494 if(i==0)
nuclear@1 495 _analysis_output("mask1L",seq,logmask,n/2,1,0,0);
nuclear@1 496 else
nuclear@1 497 _analysis_output("mask1R",seq,logmask,n/2,1,0,0);
nuclear@1 498 }else{
nuclear@1 499 _analysis_output("mask1",seq,logmask,n/2,1,0,0);
nuclear@1 500 }
nuclear@1 501 #endif
nuclear@1 502
nuclear@1 503 /* this algorithm is hardwired to floor 1 for now; abort out if
nuclear@1 504 we're *not* floor1. This won't happen unless someone has
nuclear@1 505 broken the encode setup lib. Guard it anyway. */
nuclear@1 506 if(ci->floor_type[info->floorsubmap[submap]]!=1)return(-1);
nuclear@1 507
nuclear@1 508 floor_posts[i][PACKETBLOBS/2]=
nuclear@1 509 floor1_fit(vb,b->flr[info->floorsubmap[submap]],
nuclear@1 510 logmdct,
nuclear@1 511 logmask);
nuclear@1 512
nuclear@1 513 /* are we managing bitrate? If so, perform two more fits for
nuclear@1 514 later rate tweaking (fits represent hi/lo) */
nuclear@1 515 if(vorbis_bitrate_managed(vb) && floor_posts[i][PACKETBLOBS/2]){
nuclear@1 516 /* higher rate by way of lower noise curve */
nuclear@1 517
nuclear@1 518 _vp_offset_and_mix(psy_look,
nuclear@1 519 noise,
nuclear@1 520 tone,
nuclear@1 521 2,
nuclear@1 522 logmask,
nuclear@1 523 mdct,
nuclear@1 524 logmdct);
nuclear@1 525
nuclear@1 526 #if 0
nuclear@1 527 if(vi->channels==2){
nuclear@1 528 if(i==0)
nuclear@1 529 _analysis_output("mask2L",seq,logmask,n/2,1,0,0);
nuclear@1 530 else
nuclear@1 531 _analysis_output("mask2R",seq,logmask,n/2,1,0,0);
nuclear@1 532 }else{
nuclear@1 533 _analysis_output("mask2",seq,logmask,n/2,1,0,0);
nuclear@1 534 }
nuclear@1 535 #endif
nuclear@1 536
nuclear@1 537 floor_posts[i][PACKETBLOBS-1]=
nuclear@1 538 floor1_fit(vb,b->flr[info->floorsubmap[submap]],
nuclear@1 539 logmdct,
nuclear@1 540 logmask);
nuclear@1 541
nuclear@1 542 /* lower rate by way of higher noise curve */
nuclear@1 543 _vp_offset_and_mix(psy_look,
nuclear@1 544 noise,
nuclear@1 545 tone,
nuclear@1 546 0,
nuclear@1 547 logmask,
nuclear@1 548 mdct,
nuclear@1 549 logmdct);
nuclear@1 550
nuclear@1 551 #if 0
nuclear@1 552 if(vi->channels==2){
nuclear@1 553 if(i==0)
nuclear@1 554 _analysis_output("mask0L",seq,logmask,n/2,1,0,0);
nuclear@1 555 else
nuclear@1 556 _analysis_output("mask0R",seq,logmask,n/2,1,0,0);
nuclear@1 557 }else{
nuclear@1 558 _analysis_output("mask0",seq,logmask,n/2,1,0,0);
nuclear@1 559 }
nuclear@1 560 #endif
nuclear@1 561
nuclear@1 562 floor_posts[i][0]=
nuclear@1 563 floor1_fit(vb,b->flr[info->floorsubmap[submap]],
nuclear@1 564 logmdct,
nuclear@1 565 logmask);
nuclear@1 566
nuclear@1 567 /* we also interpolate a range of intermediate curves for
nuclear@1 568 intermediate rates */
nuclear@1 569 for(k=1;k<PACKETBLOBS/2;k++)
nuclear@1 570 floor_posts[i][k]=
nuclear@1 571 floor1_interpolate_fit(vb,b->flr[info->floorsubmap[submap]],
nuclear@1 572 floor_posts[i][0],
nuclear@1 573 floor_posts[i][PACKETBLOBS/2],
nuclear@1 574 k*65536/(PACKETBLOBS/2));
nuclear@1 575 for(k=PACKETBLOBS/2+1;k<PACKETBLOBS-1;k++)
nuclear@1 576 floor_posts[i][k]=
nuclear@1 577 floor1_interpolate_fit(vb,b->flr[info->floorsubmap[submap]],
nuclear@1 578 floor_posts[i][PACKETBLOBS/2],
nuclear@1 579 floor_posts[i][PACKETBLOBS-1],
nuclear@1 580 (k-PACKETBLOBS/2)*65536/(PACKETBLOBS/2));
nuclear@1 581 }
nuclear@1 582 }
nuclear@1 583 }
nuclear@1 584 vbi->ampmax=global_ampmax;
nuclear@1 585
nuclear@1 586 /*
nuclear@1 587 the next phases are performed once for vbr-only and PACKETBLOB
nuclear@1 588 times for bitrate managed modes.
nuclear@1 589
nuclear@1 590 1) encode actual mode being used
nuclear@1 591 2) encode the floor for each channel, compute coded mask curve/res
nuclear@1 592 3) normalize and couple.
nuclear@1 593 4) encode residue
nuclear@1 594 5) save packet bytes to the packetblob vector
nuclear@1 595
nuclear@1 596 */
nuclear@1 597
nuclear@1 598 /* iterate over the many masking curve fits we've created */
nuclear@1 599
nuclear@1 600 {
nuclear@1 601 int **couple_bundle=alloca(sizeof(*couple_bundle)*vi->channels);
nuclear@1 602 int *zerobundle=alloca(sizeof(*zerobundle)*vi->channels);
nuclear@1 603
nuclear@1 604 for(k=(vorbis_bitrate_managed(vb)?0:PACKETBLOBS/2);
nuclear@1 605 k<=(vorbis_bitrate_managed(vb)?PACKETBLOBS-1:PACKETBLOBS/2);
nuclear@1 606 k++){
nuclear@1 607 oggpack_buffer *opb=vbi->packetblob[k];
nuclear@1 608
nuclear@1 609 /* start out our new packet blob with packet type and mode */
nuclear@1 610 /* Encode the packet type */
nuclear@1 611 oggpack_write(opb,0,1);
nuclear@1 612 /* Encode the modenumber */
nuclear@1 613 /* Encode frame mode, pre,post windowsize, then dispatch */
nuclear@1 614 oggpack_write(opb,modenumber,b->modebits);
nuclear@1 615 if(vb->W){
nuclear@1 616 oggpack_write(opb,vb->lW,1);
nuclear@1 617 oggpack_write(opb,vb->nW,1);
nuclear@1 618 }
nuclear@1 619
nuclear@1 620 /* encode floor, compute masking curve, sep out residue */
nuclear@1 621 for(i=0;i<vi->channels;i++){
nuclear@1 622 int submap=info->chmuxlist[i];
nuclear@1 623 int *ilogmask=iwork[i];
nuclear@1 624
nuclear@1 625 nonzero[i]=floor1_encode(opb,vb,b->flr[info->floorsubmap[submap]],
nuclear@1 626 floor_posts[i][k],
nuclear@1 627 ilogmask);
nuclear@1 628 #if 0
nuclear@1 629 {
nuclear@1 630 char buf[80];
nuclear@1 631 sprintf(buf,"maskI%c%d",i?'R':'L',k);
nuclear@1 632 float work[n/2];
nuclear@1 633 for(j=0;j<n/2;j++)
nuclear@1 634 work[j]=FLOOR1_fromdB_LOOKUP[iwork[i][j]];
nuclear@1 635 _analysis_output(buf,seq,work,n/2,1,1,0);
nuclear@1 636 }
nuclear@1 637 #endif
nuclear@1 638 }
nuclear@1 639
nuclear@1 640 /* our iteration is now based on masking curve, not prequant and
nuclear@1 641 coupling. Only one prequant/coupling step */
nuclear@1 642
nuclear@1 643 /* quantize/couple */
nuclear@1 644 /* incomplete implementation that assumes the tree is all depth
nuclear@1 645 one, or no tree at all */
nuclear@1 646 _vp_couple_quantize_normalize(k,
nuclear@1 647 &ci->psy_g_param,
nuclear@1 648 psy_look,
nuclear@1 649 info,
nuclear@1 650 gmdct,
nuclear@1 651 iwork,
nuclear@1 652 nonzero,
nuclear@1 653 ci->psy_g_param.sliding_lowpass[vb->W][k],
nuclear@1 654 vi->channels);
nuclear@1 655
nuclear@1 656 #if 0
nuclear@1 657 for(i=0;i<vi->channels;i++){
nuclear@1 658 char buf[80];
nuclear@1 659 sprintf(buf,"res%c%d",i?'R':'L',k);
nuclear@1 660 float work[n/2];
nuclear@1 661 for(j=0;j<n/2;j++)
nuclear@1 662 work[j]=iwork[i][j];
nuclear@1 663 _analysis_output(buf,seq,work,n/2,1,0,0);
nuclear@1 664 }
nuclear@1 665 #endif
nuclear@1 666
nuclear@1 667 /* classify and encode by submap */
nuclear@1 668 for(i=0;i<info->submaps;i++){
nuclear@1 669 int ch_in_bundle=0;
nuclear@1 670 long **classifications;
nuclear@1 671 int resnum=info->residuesubmap[i];
nuclear@1 672
nuclear@1 673 for(j=0;j<vi->channels;j++){
nuclear@1 674 if(info->chmuxlist[j]==i){
nuclear@1 675 zerobundle[ch_in_bundle]=0;
nuclear@1 676 if(nonzero[j])zerobundle[ch_in_bundle]=1;
nuclear@1 677 couple_bundle[ch_in_bundle++]=iwork[j];
nuclear@1 678 }
nuclear@1 679 }
nuclear@1 680
nuclear@1 681 classifications=_residue_P[ci->residue_type[resnum]]->
nuclear@1 682 class(vb,b->residue[resnum],couple_bundle,zerobundle,ch_in_bundle);
nuclear@1 683
nuclear@1 684 ch_in_bundle=0;
nuclear@1 685 for(j=0;j<vi->channels;j++)
nuclear@1 686 if(info->chmuxlist[j]==i)
nuclear@1 687 couple_bundle[ch_in_bundle++]=iwork[j];
nuclear@1 688
nuclear@1 689 _residue_P[ci->residue_type[resnum]]->
nuclear@1 690 forward(opb,vb,b->residue[resnum],
nuclear@1 691 couple_bundle,zerobundle,ch_in_bundle,classifications,i);
nuclear@1 692 }
nuclear@1 693
nuclear@1 694 /* ok, done encoding. Next protopacket. */
nuclear@1 695 }
nuclear@1 696
nuclear@1 697 }
nuclear@1 698
nuclear@1 699 #if 0
nuclear@1 700 seq++;
nuclear@1 701 total+=ci->blocksizes[vb->W]/4+ci->blocksizes[vb->nW]/4;
nuclear@1 702 #endif
nuclear@1 703 return(0);
nuclear@1 704 }
nuclear@1 705
nuclear@1 706 static int mapping0_inverse(vorbis_block *vb,vorbis_info_mapping *l){
nuclear@1 707 vorbis_dsp_state *vd=vb->vd;
nuclear@1 708 vorbis_info *vi=vd->vi;
nuclear@1 709 codec_setup_info *ci=vi->codec_setup;
nuclear@1 710 private_state *b=vd->backend_state;
nuclear@1 711 vorbis_info_mapping0 *info=(vorbis_info_mapping0 *)l;
nuclear@1 712
nuclear@1 713 int i,j;
nuclear@1 714 long n=vb->pcmend=ci->blocksizes[vb->W];
nuclear@1 715
nuclear@1 716 float **pcmbundle=alloca(sizeof(*pcmbundle)*vi->channels);
nuclear@1 717 int *zerobundle=alloca(sizeof(*zerobundle)*vi->channels);
nuclear@1 718
nuclear@1 719 int *nonzero =alloca(sizeof(*nonzero)*vi->channels);
nuclear@1 720 void **floormemo=alloca(sizeof(*floormemo)*vi->channels);
nuclear@1 721
nuclear@1 722 /* recover the spectral envelope; store it in the PCM vector for now */
nuclear@1 723 for(i=0;i<vi->channels;i++){
nuclear@1 724 int submap=info->chmuxlist[i];
nuclear@1 725 floormemo[i]=_floor_P[ci->floor_type[info->floorsubmap[submap]]]->
nuclear@1 726 inverse1(vb,b->flr[info->floorsubmap[submap]]);
nuclear@1 727 if(floormemo[i])
nuclear@1 728 nonzero[i]=1;
nuclear@1 729 else
nuclear@1 730 nonzero[i]=0;
nuclear@1 731 memset(vb->pcm[i],0,sizeof(*vb->pcm[i])*n/2);
nuclear@1 732 }
nuclear@1 733
nuclear@1 734 /* channel coupling can 'dirty' the nonzero listing */
nuclear@1 735 for(i=0;i<info->coupling_steps;i++){
nuclear@1 736 if(nonzero[info->coupling_mag[i]] ||
nuclear@1 737 nonzero[info->coupling_ang[i]]){
nuclear@1 738 nonzero[info->coupling_mag[i]]=1;
nuclear@1 739 nonzero[info->coupling_ang[i]]=1;
nuclear@1 740 }
nuclear@1 741 }
nuclear@1 742
nuclear@1 743 /* recover the residue into our working vectors */
nuclear@1 744 for(i=0;i<info->submaps;i++){
nuclear@1 745 int ch_in_bundle=0;
nuclear@1 746 for(j=0;j<vi->channels;j++){
nuclear@1 747 if(info->chmuxlist[j]==i){
nuclear@1 748 if(nonzero[j])
nuclear@1 749 zerobundle[ch_in_bundle]=1;
nuclear@1 750 else
nuclear@1 751 zerobundle[ch_in_bundle]=0;
nuclear@1 752 pcmbundle[ch_in_bundle++]=vb->pcm[j];
nuclear@1 753 }
nuclear@1 754 }
nuclear@1 755
nuclear@1 756 _residue_P[ci->residue_type[info->residuesubmap[i]]]->
nuclear@1 757 inverse(vb,b->residue[info->residuesubmap[i]],
nuclear@1 758 pcmbundle,zerobundle,ch_in_bundle);
nuclear@1 759 }
nuclear@1 760
nuclear@1 761 /* channel coupling */
nuclear@1 762 for(i=info->coupling_steps-1;i>=0;i--){
nuclear@1 763 float *pcmM=vb->pcm[info->coupling_mag[i]];
nuclear@1 764 float *pcmA=vb->pcm[info->coupling_ang[i]];
nuclear@1 765
nuclear@1 766 for(j=0;j<n/2;j++){
nuclear@1 767 float mag=pcmM[j];
nuclear@1 768 float ang=pcmA[j];
nuclear@1 769
nuclear@1 770 if(mag>0)
nuclear@1 771 if(ang>0){
nuclear@1 772 pcmM[j]=mag;
nuclear@1 773 pcmA[j]=mag-ang;
nuclear@1 774 }else{
nuclear@1 775 pcmA[j]=mag;
nuclear@1 776 pcmM[j]=mag+ang;
nuclear@1 777 }
nuclear@1 778 else
nuclear@1 779 if(ang>0){
nuclear@1 780 pcmM[j]=mag;
nuclear@1 781 pcmA[j]=mag+ang;
nuclear@1 782 }else{
nuclear@1 783 pcmA[j]=mag;
nuclear@1 784 pcmM[j]=mag-ang;
nuclear@1 785 }
nuclear@1 786 }
nuclear@1 787 }
nuclear@1 788
nuclear@1 789 /* compute and apply spectral envelope */
nuclear@1 790 for(i=0;i<vi->channels;i++){
nuclear@1 791 float *pcm=vb->pcm[i];
nuclear@1 792 int submap=info->chmuxlist[i];
nuclear@1 793 _floor_P[ci->floor_type[info->floorsubmap[submap]]]->
nuclear@1 794 inverse2(vb,b->flr[info->floorsubmap[submap]],
nuclear@1 795 floormemo[i],pcm);
nuclear@1 796 }
nuclear@1 797
nuclear@1 798 /* transform the PCM data; takes PCM vector, vb; modifies PCM vector */
nuclear@1 799 /* only MDCT right now.... */
nuclear@1 800 for(i=0;i<vi->channels;i++){
nuclear@1 801 float *pcm=vb->pcm[i];
nuclear@1 802 mdct_backward(b->transform[vb->W][0],pcm,pcm);
nuclear@1 803 }
nuclear@1 804
nuclear@1 805 /* all done! */
nuclear@1 806 return(0);
nuclear@1 807 }
nuclear@1 808
nuclear@1 809 /* export hooks */
nuclear@1 810 const vorbis_func_mapping mapping0_exportbundle={
nuclear@1 811 &mapping0_pack,
nuclear@1 812 &mapping0_unpack,
nuclear@1 813 &mapping0_free_info,
nuclear@1 814 &mapping0_forward,
nuclear@1 815 &mapping0_inverse
nuclear@1 816 };