/[svn]/linuxsampler/trunk/src/engines/common/Resampler.h
ViewVC logotype

Diff of /linuxsampler/trunk/src/engines/common/Resampler.h

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 903 by persson, Sat Jul 22 14:22:53 2006 UTC revision 1485 by senoner, Thu Nov 15 23:35:45 2007 UTC
# Line 3  Line 3 
3   *   LinuxSampler - modular, streaming capable sampler                     *   *   LinuxSampler - modular, streaming capable sampler                     *
4   *                                                                         *   *                                                                         *
5   *   Copyright (C) 2003, 2004 by Benno Senoner and Christian Schoenebeck   *   *   Copyright (C) 2003, 2004 by Benno Senoner and Christian Schoenebeck   *
6   *   Copyright (C) 2005 Christian Schoenebeck                              *   *   Copyright (C) 2005 - 2007 Christian Schoenebeck                       *
7   *                                                                         *   *                                                                         *
8   *   This program is free software; you can redistribute it and/or modify  *   *   This program is free software; you can redistribute it and/or modify  *
9   *   it under the terms of the GNU General Public License as published by  *   *   it under the terms of the GNU General Public License as published by  *
# Line 21  Line 21 
21   *   MA  02111-1307  USA                                                   *   *   MA  02111-1307  USA                                                   *
22   ***************************************************************************/   ***************************************************************************/
23    
24    // Note: the assembly code is currently disabled, as it doesn't fit into
25    // the new synthesis core introduced by LS 0.4.0
26    
27  #ifndef __LS_RESAMPLER_H__  #ifndef __LS_RESAMPLER_H__
28  #define __LS_RESAMPLER_H__  #define __LS_RESAMPLER_H__
29    
30  #include "../../common/global.h"  #include "../../common/global_private.h"
31    
32  // TODO: cubic interpolation is not yet supported by the MMX/SSE(1) version though  // TODO: cubic interpolation is not yet supported by the MMX/SSE(1) version though
 // TODO: cubic interpolation is not supported for 24 bit samples  
33  #ifndef USE_LINEAR_INTERPOLATION  #ifndef USE_LINEAR_INTERPOLATION
34  # define USE_LINEAR_INTERPOLATION   1  ///< set to 0 if you prefer cubic interpolation (slower, better quality)  # define USE_LINEAR_INTERPOLATION   1  ///< set to 0 if you prefer cubic interpolation (slower, better quality)
35  #endif  #endif
# Line 75  namespace LinuxSampler { Line 77  namespace LinuxSampler {
77                  }                  }
78              }              }
79    
80  #if CONFIG_ASM && ARCH_X86  #if 0 // CONFIG_ASM && ARCH_X86
81              inline static void GetNext4SamplesMonoMMXSSE(sample_t* pSrc, void* Pos, float& Pitch) {              inline static void GetNext4SamplesMonoMMXSSE(sample_t* pSrc, void* Pos, float& Pitch) {
82                  if (INTERPOLATE) Interpolate4StepsMonoMMXSSE(pSrc, Pos, Pitch);                  if (INTERPOLATE) Interpolate4StepsMonoMMXSSE(pSrc, Pos, Pitch);
83                  else { // no pitch, so no interpolation necessary                  else { // no pitch, so no interpolation necessary
# Line 147  namespace LinuxSampler { Line 149  namespace LinuxSampler {
149    
150          protected:          protected:
151    
152              static int getSample(sample_t* src, int pos) {              inline static int32_t getSample(sample_t* src, int pos) {
153                  if (BITDEPTH24) {                  if (BITDEPTH24) {
154                      pos *= 3;                      pos *= 3;
155                        #if WORDS_BIGENDIAN
156                      unsigned char* p = (unsigned char*)src;                      unsigned char* p = (unsigned char*)src;
157                      return p[pos] << 8 | p[pos + 1] << 16 | p[pos + 2] << 24;                      return p[pos] << 8 | p[pos + 1] << 16 | p[pos + 2] << 24;
158                        #else
159                        // 24bit read optimization:
160                        // a misaligned 32bit read and subquent 8 bit shift is faster (on x86) than reading 3 single bytes and shifting them
161                        return (*((int32_t *)(&((char *)(src))[pos])))<<8;
162                        #endif
163                  } else {                  } else {
164                      return src[pos];                      return src[pos];
165                  }                  }
# Line 166  namespace LinuxSampler { Line 174  namespace LinuxSampler {
174                      int x2 = getSample(pSrc, pos_int + 1);                      int x2 = getSample(pSrc, pos_int + 1);
175                      float samplePoint  = (x1 + pos_fract * (x2 - x1));                      float samplePoint  = (x1 + pos_fract * (x2 - x1));
176                  #else // polynomial interpolation                  #else // polynomial interpolation
177                      float xm1 = pSrc[pos_int];                      float xm1 = getSample(pSrc, pos_int);
178                      float x0  = pSrc[pos_int+1];                      float x0  = getSample(pSrc, pos_int + 1);
179                      float x1  = pSrc[pos_int+2];                      float x1  = getSample(pSrc, pos_int + 2);
180                      float x2  = pSrc[pos_int+3];                      float x2  = getSample(pSrc, pos_int + 3);
181                      float a   = (3.0f * (x0 - x1) - xm1 + x2) * 0.5f;                      float a   = (3.0f * (x0 - x1) - xm1 + x2) * 0.5f;
182                      float b   = 2.0f * x1 + xm1 - (5.0f * x0 + x2) * 0.5f;                      float b   = 2.0f * x1 + xm1 - (5.0f * x0 + x2) * 0.5f;
183                      float c   = (x1 - xm1) * 0.5f;                      float c   = (x1 - xm1) * 0.5f;
# Line 198  namespace LinuxSampler { Line 206  namespace LinuxSampler {
206                      samplePoint.right = (x1 + pos_fract * (x2 - x1));                      samplePoint.right = (x1 + pos_fract * (x2 - x1));
207                  #else // polynomial interpolation                  #else // polynomial interpolation
208                      // calculate left channel                      // calculate left channel
209                      float xm1 = pSrc[pos_int];                      float xm1 = getSample(pSrc, pos_int);
210                      float x0  = pSrc[pos_int+2];                      float x0  = getSample(pSrc, pos_int + 2);
211                      float x1  = pSrc[pos_int+4];                      float x1  = getSample(pSrc, pos_int + 4);
212                      float x2  = pSrc[pos_int+6];                      float x2  = getSample(pSrc, pos_int + 6);
213                      float a   = (3.0f * (x0 - x1) - xm1 + x2) * 0.5f;                      float a   = (3.0f * (x0 - x1) - xm1 + x2) * 0.5f;
214                      float b   = 2.0f * x1 + xm1 - (5.0f * x0 + x2) * 0.5f;                      float b   = 2.0f * x1 + xm1 - (5.0f * x0 + x2) * 0.5f;
215                      float c   = (x1 - xm1) * 0.5f;                      float c   = (x1 - xm1) * 0.5f;
216                      samplePoint.left = (((a * pos_fract) + b) * pos_fract + c) * pos_fract + x0;                      samplePoint.left = (((a * pos_fract) + b) * pos_fract + c) * pos_fract + x0;
217    
218                      //calculate right channel                      //calculate right channel
219                      xm1 = pSrc[pos_int+1];                      xm1 = getSample(pSrc, pos_int + 1);
220                      x0  = pSrc[pos_int+3];                      x0  = getSample(pSrc, pos_int + 3);
221                      x1  = pSrc[pos_int+5];                      x1  = getSample(pSrc, pos_int + 5);
222                      x2  = pSrc[pos_int+7];                      x2  = getSample(pSrc, pos_int + 7);
223                      a   = (3.0f * (x0 - x1) - xm1 + x2) * 0.5f;                      a   = (3.0f * (x0 - x1) - xm1 + x2) * 0.5f;
224                      b   = 2.0f * x1 + xm1 - (5.0f * x0 + x2) * 0.5f;                      b   = 2.0f * x1 + xm1 - (5.0f * x0 + x2) * 0.5f;
225                      c   = (x1 - xm1) * 0.5f;                      c   = (x1 - xm1) * 0.5f;
# Line 222  namespace LinuxSampler { Line 230  namespace LinuxSampler {
230                  return samplePoint;                  return samplePoint;
231              }              }
232    
233  #if CONFIG_ASM && ARCH_X86  #if 0 // CONFIG_ASM && ARCH_X86
234              // TODO: no support for cubic interpolation yet              // TODO: no support for cubic interpolation yet
235              inline static void Interpolate4StepsMonoMMXSSE(sample_t* pSrc, void* Pos, float& Pitch) {              inline static void Interpolate4StepsMonoMMXSSE(sample_t* pSrc, void* Pos, float& Pitch) {
236                  /* calculate playback position of each of the 4 samples by adding the associated pitch */                  /* calculate playback position of each of the 4 samples by adding the associated pitch */
Legend:
Removed from v.903  
changed lines
  Added in v.1485
  ViewVC Help
Powered by ViewVC