| @@ -21,8 +21,7 @@ · 20 unmodified lines |
| 21 | // The C-embedded Audio Library. | 21 | // The C-embedded Audio Library. |
| 22 | #define TWO_TO_16 65536.f | 22 | #define TWO_TO_16 65536.f |
| 23 | | 23 | |
| 24 | #define EXPONENTIAL_TABLE_SIZE 65536 | 24 | #define EXPONENTIAL_TABLE_SIZE 65536 |
| 25 | | | |
| 26 | | 25 | |
| 27 | float interpolate3max(float *buf, const int peakindex) | 26 | float interpolate3max(float *buf, const int peakindex) |
| 28 | { | 27 | { |
| @@ -71,7 +70,7 @@ float fastabs(float f) · 42 unmodified lines |
| 71 | | 70 | |
| 72 | alias.f = f; | 71 | alias.f = f; |
| 73 | alias.ui &= 0x7fffffff; | 72 | alias.ui &= 0x7fffffff; |
| 74 | return alias.f; | 73 | return alias.f; |
| 75 | } | 74 | } |
| 76 | | 75 | |
| 77 | // dope af | 76 | // dope af |
| @@ -207,8 +206,7 @@ int LEAF_clipInt(int min, int val, int max) · 129 unmodified lines |
| 207 | } else { | 206 | } else { |
| 208 | return val; | 207 | return val; |
| 209 | } | 208 | } |
| 210 | } | 209 | } |
| 211 | | | |
| 212 | | 210 | |
| 213 | oBool LEAF_isPrime(uint64_t number ) | 211 | oBool LEAF_isPrime(uint64_t number ) |
| 214 | { | 212 | { |
| @@ -232,91 +230,91 @@ float LEAF_tanh(float x) · 17 unmodified lines |
| 232 | else | 230 | else |
| 233 | return x * ( 27 + x * x ) / ( 27 + 9 * x * x ); | 231 | return x * ( 27 + x * x ) / ( 27 + 9 * x * x ); |
| 234 | } | 232 | } |
| 235 | | 233 | |
| 236 | | 234 | |
| 237 | void LEAF_generate_sine(float* buffer, int size) | 235 | void LEAF_generate_sine(float* buffer, int size) |
| 238 | { | 236 | { |
| 239 | float phase; | 237 | float phase; |
| 240 | for (int i = 0; i < size; i++) | 238 | for (int i = 0; i < size; i++) |
| 241 | { | 239 | { |
| 242 | phase = (float) i / (float) size; | 240 | phase = (float) i / (float) size; |
| 243 | buffer[i] = sinf(phase * TWO_PI); | 241 | buffer[i] = sinf(phase * TWO_PI); |
| 244 | } | 242 | } |
| 245 | } | 243 | } |
| 246 | | 244 | |
| 247 | void LEAF_generate_sawtooth(float* buffer, float basefreq, int size) | 245 | void LEAF_generate_sawtooth(float* buffer, float basefreq, int size) |
| 248 | { | 246 | { |
| 249 | int harmonic = 1; | 247 | int harmonic = 1; |
| 250 | float phase = 0.0f; | 248 | float phase = 0.0f; |
| 251 | float freq = harmonic * basefreq; | 249 | float freq = harmonic * basefreq; |
| 252 | float amp; | 250 | float amp; |
| 253 | | 251 | |
| 254 | while (freq < (leaf.sampleRate * 0.5)) | 252 | while (freq < (leaf.sampleRate * 0.5)) |
| 255 | { | 253 | { |
| 256 | amp = 1.0f / harmonic; | 254 | amp = 1.0f / harmonic; |
| 257 | for (int i = 0; i < size; i++) | 255 | for (int i = 0; i < size; i++) |
| 258 | { | 256 | { |
| 259 | phase = (float) i / (float) size; | 257 | phase = (float) i / (float) size; |
| 260 | buffer[i] += (amp * sinf(harmonic * phase * TWO_PI)); | 258 | buffer[i] += (amp * sinf(harmonic * phase * TWO_PI)); |
| 261 | } | 259 | } |
| 262 | | 260 | |
| 263 | harmonic++; | 261 | harmonic++; |
| 264 | freq = harmonic * basefreq; | 262 | freq = harmonic * basefreq; |
| 265 | } | 263 | } |
| 266 | } | 264 | } |
| 267 | | 265 | |
| 268 | | 266 | |
| 269 | void LEAF_generate_triangle(float* buffer, float basefreq, int size) | 267 | void LEAF_generate_triangle(float* buffer, float basefreq, int size) |
| 270 | { | 268 | { |
| 271 | int harmonic = 1; | 269 | int harmonic = 1; |
| 272 | float phase = 0.0f; | 270 | float phase = 0.0f; |
| 273 | float freq = harmonic * basefreq; | 271 | float freq = harmonic * basefreq; |
| 274 | float amp = 1.0f; | 272 | float amp = 1.0f; |
| 275 | | 273 | |
| 276 | int count = 0; | 274 | int count = 0; |
| 277 | float mult = 1.0f; | 275 | float mult = 1.0f; |
| 278 | | 276 | |
| 279 | while (freq < (leaf.sampleRate * 0.5)) | 277 | while (freq < (leaf.sampleRate * 0.5)) |
| 280 | { | 278 | { |
| 281 | amp = 1.0f / (float)(harmonic * harmonic); | 279 | amp = 1.0f / (float)(harmonic * harmonic); |
| 282 | | 280 | |
| 283 | if (count % 2) mult = -1.0f; | 281 | if (count % 2) mult = -1.0f; |
| 284 | else mult = 1.0f; | 282 | else mult = 1.0f; |
| 285 | | 283 | |
| 286 | for (int i = 0; i < size; i++) | 284 | for (int i = 0; i < size; i++) |
| 287 | { | 285 | { |
| 288 | phase = (float) i / (float) size; | 286 | phase = (float) i / (float) size; |
| 289 | buffer[i] += (mult * amp * sinf(harmonic * phase * TWO_PI)); | 287 | buffer[i] += (mult * amp * sinf(harmonic * phase * TWO_PI)); |
| 290 | } | 288 | } |
| 291 | | 289 | |
| 292 | count++; | 290 | count++; |
| 293 | harmonic += 2; | 291 | harmonic += 2; |
| 294 | freq = harmonic * basefreq; | 292 | freq = harmonic * basefreq; |
| 295 | } | 293 | } |
| 296 | } | 294 | } |
| 297 | | 295 | |
| 298 | void LEAF_generate_square(float* buffer, float basefreq, int size) | 296 | void LEAF_generate_square(float* buffer, float basefreq, int size) |
| 299 | { | 297 | { |
| 300 | int harmonic = 1; | 298 | int harmonic = 1; |
| 301 | float phase = 0.0f; | 299 | float phase = 0.0f; |
| 302 | float freq = harmonic * basefreq; | 300 | float freq = harmonic * basefreq; |
| 303 | float amp = 1.0f; | 301 | float amp = 1.0f; |
| 304 | | 302 | |
| 305 | while (freq < (leaf.sampleRate * 0.5)) | 303 | while (freq < (leaf.sampleRate * 0.5)) |
| 306 | { | 304 | { |
| 307 | amp = 1.0f / (float)(harmonic); | 305 | amp = 1.0f / (float)(harmonic); |
| 308 | | 306 | |
| 309 | for (int i = 0; i < size; i++) | 307 | for (int i = 0; i < size; i++) |
| 310 | { | 308 | { |
| 311 | phase = (float) i / (float) size; | 309 | phase = (float) i / (float) size; |
| 312 | buffer[i] += (amp * sinf(harmonic * phase * TWO_PI)); | 310 | buffer[i] += (amp * sinf(harmonic * phase * TWO_PI)); |
| 313 | } | 311 | } |
| 314 | | 312 | |
| 315 | harmonic += 2; | 313 | harmonic += 2; |
| 316 | freq = harmonic * basefreq; | 314 | freq = harmonic * basefreq; |
| 317 | } | 315 | } |
| 318 | } | 316 | } |
| 319 | | 317 | |
| 320 | | 318 | |
| 321 | //----------------------------------------------------------------------------- | 319 | //----------------------------------------------------------------------------- |
| 322 | // name: mtof() | 320 | // name: mtof() |
| @@ -326,35 +324,34 @@ float LEAF_midiToFrequency(float f) · 3 unmodified lines |
| 326 | { | 324 | { |
| 327 | if( f <= -1500.0f ) return (0); | 325 | if( f <= -1500.0f ) return (0); |
| 328 | else if( f > 1499.0f ) return (LEAF_midiToFrequency(1499.0f)); | 326 | else if( f > 1499.0f ) return (LEAF_midiToFrequency(1499.0f)); |
| 329 | else return ( powf(2.0f, (f - 69.0f) / 12.0f) * 440.0f ); | 327 | else return ( powf(2.0f, (f - 69.0f) * 0.083333333333333f) * 440.0f ); |
| 330 | } | 328 | } |
| 331 | | 329 | |
| 332 | | 330 | |
| 333 | // alpha, [0.0, 1.0] | 331 | // alpha, [0.0, 1.0] |
| 334 | float LEAF_interpolate_hermite (float A, float B, float C, float D, float alpha) | 332 | float LEAF_interpolate_hermite (float A, float B, float C, float D, float alpha) |
| 335 | { | 333 | { |
| 336 | alpha = LEAF_clip(0.0f, alpha, 1.0f); | 334 | alpha = LEAF_clip(0.0f, alpha, 1.0f); |
| 337 | | 335 | |
| 338 | float a = -A/2.0f + (3.0f*B)/2.0f - (3.0f*C)/2.0f + D/2.0f; | 336 | float a = -A/2.0f + (3.0f*B)/2.0f - (3.0f*C)/2.0f + D/2.0f; |
| 339 | float b = A - (5.0f*B)/2.0f + 2.0f*C - D / 2.0f; | 337 | float b = A - (5.0f*B)/2.0f + 2.0f*C - D / 2.0f; |
| 340 | float c = -A/2.0f + C/2.0f; | 338 | float c = -A/2.0f + C/2.0f; |
| 341 | float d = B; | 339 | float d = B; |
| 342 | | 340 | |
| 343 | return a*alpha*alpha*alpha + b*alpha*alpha + c*alpha + d; | 341 | return a*alpha*alpha*alpha + b*alpha*alpha + c*alpha + d; |
| 344 | } | 342 | } |
| 345 | | 343 | |
| 346 | // alpha, [0.0, 1.0] | 344 | // alpha, [0.0, 1.0] |
| 347 | float LEAF_interpolation_linear (float A, float B, float alpha) | 345 | float LEAF_interpolation_linear (float A, float B, float alpha) |
| 348 | { | 346 | { |
| 349 | alpha = LEAF_clip(0.0f, alpha, 1.0f); | 347 | alpha = LEAF_clip(0.0f, alpha, 1.0f); |
| 350 | | 348 | |
| 351 | float omAlpha = 1.0f - alpha; | 349 | float omAlpha = 1.0f - alpha; |
| 352 | | 350 | |
| 353 | // First 1/2 of interpolation | 351 | // First 1/2 of interpolation |
| 354 | float out = A * omAlpha; | 352 | float out = A * omAlpha; |
| 355 | | | |
| 356 | out += B * alpha; | | |
| 357 | | | |
| 358 | return out; | | |
| 359 | } | | |
| 360 | | 353 | |
| | 354 | out += B * alpha; |
| | 355 | |
| | 356 | return out; |
| | 357 | } |