ADX (ファイル形式)

 #define M_PI acos(-1.0) 
double a , b , c ; a = sqrt ( 2.0 ) - cos ( 2.0 * M_PI * (( double ) adx_header -> highpass_frequency / adx_header -> sample_rate )); b = sqrt ( 2.0 ) - 1.0 ; c = ( a - sqrt (( a + b ) * ( a - b ))) / b ; //(a+b)*(ab) = a*ab*b、ただし、より単純な式では浮動小数点の精度が失われます// double factor[2]; factor [ 0 ] = c * 2.0 ; factor [ 1 ] = - ( c * c );    
           
     
              
 
 
     
     

 static int32_t * past_samples ; // 各チャネルから以前にデコードされたサンプル。開始時にゼロに設定されます (size = 2*channel_count) static uint_fast32_t sample_index = 0 ; // sample_index は、次にデコードする必要があるサンプル セットのインデックスですstatic ADX_header * adx_header ; // buffer は、デコードされたサンプルが配置される場所です// samples_needed は、バッファーを埋めるためにデコードする必要があるサンプル 'セット' (各チャネルから 1 つのサンプル) の数を示します// looping_enabled は、組み込みループの使用を制御するブール フラグです// バッファー内で埋めることができなかったサンプル 'セット' の数 (EOS) を返しますunsigned decode_adx_standard ( int16_t * buffer , unsigned samples_needed , bool looping_enabled ) { unsigned const samples_per_block = ( adx_header -> block_size - 2 ) * 8 / adx_header -> sample_bitdepth ; int16_t scale [ adx_header -> channel_count ]; if ( looping_enabled && ! adx_header -> loop_enabled ) looping_enabled = false ; // 要求された数のサンプルがデコードされるか、ファイルの終わりに達するまでループしますwhile ( samples_needed > 0 && sample_index < adx_header -> total_samples ) { // 現在のブロックでデコードされていないサンプルの数を計算しますunsigned sample_offset = sample_index % samples_per_block ; unsigned samples_can_get = samples_per_block - sample_offset ;        
      
     
 
 
 
 
 
         
 
            
     
 
     
       
 
  
         
  
     
          
          

     // サンプルがバッファーに収まらない場合は、この実行中に取得できるサンプル数をクランプします。
if ( samples_can_get > samples_needed ) samples_can_get = samples_needed ; // ストリームの長さが十分でないか、ループ トリガーが近くにある場合は、取得するサンプル数をクランプします。if ( looping_enabled && sample_index + samples_can_get > adx_header - > loop_end_index ) samples_can_get = adx_header - > loop_end_index - sample_index ; else if ( sample_index + samples_can_get > adx_header - > total_samples ) samples_can_get = adx_header - > total_samples - sample_index ; // sample_index が存在するフレームの開始のビット アドレスを計算し、その位置を記録しますunsigned long started_at = ( adx_header -> copyright_offset + 4 + \
 sample_index / samples_per_block * adx_header -> block_size * adx_header -> channel_count ) * 8 ; // このフレームの各ブロックの先頭からスケール値を読み取りますfor ( unsigned i = 0 ; i < adx_header -> channel_count ; ++ i ) { bitstream_seek ( started_at + adx_header -> block_size * i * 8 ); scale [ i ] = ntohs ( bitstream_read ( 16 ) ); } // sample_offset の停止値を事前計算しますunsigned sample_endoffset = sample_offset + samples_can_get ; // フレームの最初のブロックのスケール直後の最初のサンプルのビットストリームアドレスを保存します。started_at += 16 ; while ( sample_offset        
          
 
     
            
            
           
            
 
     
                                          
 
     
               
     
                
              
     
 
     
          
 
     
       
        < sample_endoffset ) { for ( unsigned i = 0 ; i < adx_header -> channel_count ; ++ i ) { // 次のサンプルを予測しますdouble sample_prediction = factor [ 0 ] * past_samples [ i * 2 + 0 ] + factor [ 1 ] * past_samples [ i * 2 + 1 ]; // サンプル オフセットまでシークし、読み取って 32 ビット整数に符号拡張します// 符号拡張の実装は読者の演習として残されています// 8 ビットを超える場合は、符号拡張にエンディアン調整も含める必要がありますbitstream_seek ( started_at + adx_header -> sample_bitdepth * sample_offset + \
 adx_header -> block_size * 8 * i ); int_fast32_t sample_error = bitstream_read ( adx_header -> sample_bitdepth ); sample_error = sign_extend ( sample_error , adx_header -> sample_bitdepth ); // エラー訂正値をスケールするsample_error *= scale [ i ]; // 予測とエラー訂正を組み合わせてサンプルを計算するint_fast32_t sample = sample_error + ( int_fast32_t ) sample_prediction ; // 過去のサンプルを新しいサンプルで更新するpast_samples [ i * 2 + 1 ] = past_samples [ i * 2 + 0 ]; past_samples [ i * 2 + 0 ] = sample ; // デコードされたサンプルを 16 ビット整数の有効範囲にクランプするif ( sample > 32767 ) sample =  
     
                  
        
           
                        
 
           
           
           
                                                  
                
                
 
           
             
 
           
                
 
           
                 
               
 
           
              
                32767 ; 
else if ( sample < -32768 ) sample = -32768 ; // サンプルをバッファーに保存し、1 つ進めます* buffer ++ = sample ; } ++ sample_offset ; // すべてのブロックから 1 つのサンプルをデコードしたので、ブロック オフセットを 1 つ進めます++ sample_index ; // これは、ストリームに 1 サンプル進んでいることも意味します-- samples_needed ; // したがって、デコードする必要があるサンプル セットが 1 つ少なくなります} // ループ終了マーカーに達したかどうかを確認します。達した場合は、ループの開始にジャンプする必要がありますif ( looping_enabled && sample_index == adx_header -> loop_end_index ) sample_index = adx_header -> loop_start_index ; } return samples_needed ; }