/* test_streams - Simple test pattern generator * Copyright (C) 2000,2001,2002,2003,2004,2005,2006,2007 Josh Coalson * * This program is free software; you can redistribute it and/or * modify it under the terms of the GNU General Public License * as published by the Free Software Foundation; either version 2 * of the License, or (at your option) any later version. * * This program is distributed in the hope that it will be useful, * but WITHOUT ANY WARRANTY; without even the implied warranty of * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the * GNU General Public License for more details. * * You should have received a copy of the GNU General Public License * along with this program; if not, write to the Free Software * Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA. */ #if HAVE_CONFIG_H # include #endif #include #include #include #if defined _MSC_VER || defined __MINGW32__ #include #else #include #endif #include "FLAC/assert.h" #include "FLAC/ordinals.h" #ifndef M_PI /* math.h in VC++ doesn't seem to have this (how Microsoft is that?) */ #define M_PI 3.14159265358979323846 #endif #if defined _WIN32 || defined __EMX__ static const char *mode = "wb"; #else static const char *mode = "w"; #endif #if !defined _MSC_VER && !defined __MINGW32__ #define GET_RANDOM_BYTE (((unsigned)random()) & 0xff) #else #define GET_RANDOM_BYTE (((unsigned)rand()) & 0xff) #endif static FLAC__bool is_big_endian_host; static FLAC__bool write_little_endian_uint16(FILE *f, FLAC__uint16 x) { return fputc(x, f) != EOF && fputc(x >> 8, f) != EOF ; } static FLAC__bool write_little_endian_int16(FILE *f, FLAC__int16 x) { return write_little_endian_uint16(f, (FLAC__uint16)x); } static FLAC__bool write_little_endian_uint24(FILE *f, FLAC__uint32 x) { return fputc(x, f) != EOF && fputc(x >> 8, f) != EOF && fputc(x >> 16, f) != EOF ; } static FLAC__bool write_little_endian_int24(FILE *f, FLAC__int32 x) { return write_little_endian_uint24(f, (FLAC__uint32)x); } static FLAC__bool write_little_endian_uint32(FILE *f, FLAC__uint32 x) { return fputc(x, f) != EOF && fputc(x >> 8, f) != EOF && fputc(x >> 16, f) != EOF && fputc(x >> 24, f) != EOF ; } #if 0 /* @@@ not used (yet) */ static FLAC__bool write_little_endian_int32(FILE *f, FLAC__int32 x) { return write_little_endian_uint32(f, (FLAC__uint32)x); } #endif static FLAC__bool write_big_endian_uint16(FILE *f, FLAC__uint16 x) { return fputc(x >> 8, f) != EOF && fputc(x, f) != EOF ; } #if 0 /* @@@ not used (yet) */ static FLAC__bool write_big_endian_int16(FILE *f, FLAC__int16 x) { return write_big_endian_uint16(f, (FLAC__uint16)x); } #endif #if 0 /* @@@ not used (yet) */ static FLAC__bool write_big_endian_uint24(FILE *f, FLAC__uint32 x) { return fputc(x >> 16, f) != EOF && fputc(x >> 8, f) != EOF && fputc(x, f) != EOF ; } #endif #if 0 /* @@@ not used (yet) */ static FLAC__bool write_big_endian_int24(FILE *f, FLAC__int32 x) { return write_big_endian_uint24(f, (FLAC__uint32)x); } #endif static FLAC__bool write_big_endian_uint32(FILE *f, FLAC__uint32 x) { return fputc(x >> 24, f) != EOF && fputc(x >> 16, f) != EOF && fputc(x >> 8, f) != EOF && fputc(x, f) != EOF ; } #if 0 /* @@@ not used (yet) */ static FLAC__bool write_big_endian_int32(FILE *f, FLAC__int32 x) { return write_big_endian_uint32(f, (FLAC__uint32)x); } #endif static FLAC__bool write_sane_extended(FILE *f, unsigned val) /* Write to 'f' a SANE extended representation of 'val'. Return false if * the write succeeds; return true otherwise. * * SANE extended is an 80-bit IEEE-754 representation with sign bit, 15 bits * of exponent, and 64 bits of significand (mantissa). Unlike most IEEE-754 * representations, it does not imply a 1 above the MSB of the significand. * * Preconditions: * val!=0U */ { unsigned int shift, exponent; FLAC__ASSERT(val!=0U); /* handling 0 would require a special case */ for(shift= 0U; (val>>(31-shift))==0U; ++shift) ; val<<= shift; exponent= 63U-(shift+32U); /* add 32 for unused second word */ if(!write_big_endian_uint16(f, (FLAC__uint16)(exponent+0x3FFF))) return false; if(!write_big_endian_uint32(f, val)) return false; if(!write_big_endian_uint32(f, 0)) /* unused second word */ return false; return true; } /* a mono one-sample 16bps stream */ static FLAC__bool generate_01(void) { FILE *f; FLAC__int16 x = -32768; if(0 == (f = fopen("test01.raw", mode))) return false; if(!write_little_endian_int16(f, x)) goto foo; fclose(f); return true; foo: fclose(f); return false; } /* a stereo one-sample 16bps stream */ static FLAC__bool generate_02(void) { FILE *f; FLAC__int16 xl = -32768, xr = 32767; if(0 == (f = fopen("test02.raw", mode))) return false; if(!write_little_endian_int16(f, xl)) goto foo; if(!write_little_endian_int16(f, xr)) goto foo; fclose(f); return true; foo: fclose(f); return false; } /* a mono five-sample 16bps stream */ static FLAC__bool generate_03(void) { FILE *f; FLAC__int16 x[] = { -25, 0, 25, 50, 100 }; unsigned i; if(0 == (f = fopen("test03.raw", mode))) return false; for(i = 0; i < 5; i++) if(!write_little_endian_int16(f, x[i])) goto foo; fclose(f); return true; foo: fclose(f); return false; } /* a stereo five-sample 16bps stream */ static FLAC__bool generate_04(void) { FILE *f; FLAC__int16 x[] = { -25, 500, 0, 400, 25, 300, 50, 200, 100, 100 }; unsigned i; if(0 == (f = fopen("test04.raw", mode))) return false; for(i = 0; i < 10; i++) if(!write_little_endian_int16(f, x[i])) goto foo; fclose(f); return true; foo: fclose(f); return false; } /* a mono full-scale deflection 8bps stream */ static FLAC__bool generate_fsd8(const char *fn, const int pattern[], unsigned reps) { FILE *f; unsigned rep, p; FLAC__ASSERT(pattern != 0); if(0 == (f = fopen(fn, mode))) return false; for(rep = 0; rep < reps; rep++) { for(p = 0; pattern[p]; p++) { signed char x = pattern[p] > 0? 127 : -128; if(fwrite(&x, sizeof(x), 1, f) < 1) goto foo; } } fclose(f); return true; foo: fclose(f); return false; } /* a mono full-scale deflection 16bps stream */ static FLAC__bool generate_fsd16(const char *fn, const int pattern[], unsigned reps) { FILE *f; unsigned rep, p; FLAC__ASSERT(pattern != 0); if(0 == (f = fopen(fn, mode))) return false; for(rep = 0; rep < reps; rep++) { for(p = 0; pattern[p]; p++) { FLAC__int16 x = pattern[p] > 0? 32767 : -32768; if(!write_little_endian_int16(f, x)) goto foo; } } fclose(f); return true; foo: fclose(f); return false; } /* a stereo wasted-bits-per-sample 16bps stream */ static FLAC__bool generate_wbps16(const char *fn, unsigned samples) { FILE *f; unsigned sample; if(0 == (f = fopen(fn, mode))) return false; for(sample = 0; sample < samples; sample++) { FLAC__int16 l = (sample % 2000) << 2; FLAC__int16 r = (sample % 1000) << 3; if(!write_little_endian_int16(f, l)) goto foo; if(!write_little_endian_int16(f, r)) goto foo; } fclose(f); return true; foo: fclose(f); return false; } /* a mono full-scale deflection 24bps stream */ static FLAC__bool generate_fsd24(const char *fn, const int pattern[], unsigned reps) { FILE *f; unsigned rep, p; FLAC__ASSERT(pattern != 0); if(0 == (f = fopen(fn, mode))) return false; for(rep = 0; rep < reps; rep++) { for(p = 0; pattern[p]; p++) { FLAC__int32 x = pattern[p] > 0? 8388607 : -8388608; if(!write_little_endian_int24(f, x)) goto foo; } } fclose(f); return true; foo: fclose(f); return false; } /* a mono sine-wave 8bps stream */ static FLAC__bool generate_sine8_1(const char *fn, const double sample_rate, const unsigned samples, const double f1, const double a1, const double f2, const double a2) { const FLAC__int8 full_scale = 127; const double delta1 = 2.0 * M_PI / ( sample_rate / f1); const double delta2 = 2.0 * M_PI / ( sample_rate / f2); FILE *f; double theta1, theta2; unsigned i; if(0 == (f = fopen(fn, mode))) return false; for(i = 0, theta1 = theta2 = 0.0; i < samples; i++, theta1 += delta1, theta2 += delta2) { double val = (a1*sin(theta1) + a2*sin(theta2))*(double)full_scale; FLAC__int8 v = (FLAC__int8)(val + 0.5); if(fwrite(&v, sizeof(v), 1, f) < 1) goto foo; } fclose(f); return true; foo: fclose(f); return false; } /* a stereo sine-wave 8bps stream */ static FLAC__bool generate_sine8_2(const char *fn, const double sample_rate, const unsigned samples, const double f1, const double a1, const double f2, const double a2, double fmult) { const FLAC__int8 full_scale = 127; const double delta1 = 2.0 * M_PI / ( sample_rate / f1); const double delta2 = 2.0 * M_PI / ( sample_rate / f2); FILE *f; double theta1, theta2; unsigned i; if(0 == (f = fopen(fn, mode))) return false; for(i = 0, theta1 = theta2 = 0.0; i < samples; i++, theta1 += delta1, theta2 += delta2) { double val = (a1*sin(theta1) + a2*sin(theta2))*(double)full_scale; FLAC__int8 v = (FLAC__int8)(val + 0.5); if(fwrite(&v, sizeof(v), 1, f) < 1) goto foo; val = -(a1*sin(theta1*fmult) + a2*sin(theta2*fmult))*(double)full_scale; v = (FLAC__int8)(val + 0.5); if(fwrite(&v, sizeof(v), 1, f) < 1) goto foo; } fclose(f); return true; foo: fclose(f); return false; } /* a mono sine-wave 16bps stream */ static FLAC__bool generate_sine16_1(const char *fn, const double sample_rate, const unsigned samples, const double f1, const double a1, const double f2, const double a2) { const FLAC__int16 full_scale = 32767; const double delta1 = 2.0 * M_PI / ( sample_rate / f1); const double delta2 = 2.0 * M_PI / ( sample_rate / f2); FILE *f; double theta1, theta2; unsigned i; if(0 == (f = fopen(fn, mode))) return false; for(i = 0, theta1 = theta2 = 0.0; i < samples; i++, theta1 += delta1, theta2 += delta2) { double val = (a1*sin(theta1) + a2*sin(theta2))*(double)full_scale; FLAC__int16 v = (FLAC__int16)(val + 0.5); if(!write_little_endian_int16(f, v)) goto foo; } fclose(f); return true; foo: fclose(f); return false; } /* a stereo sine-wave 16bps stream */ static FLAC__bool generate_sine16_2(const char *fn, const double sample_rate, const unsigned samples, const double f1, const double a1, const double f2, const double a2, double fmult) { const FLAC__int16 full_scale = 32767; const double delta1 = 2.0 * M_PI / ( sample_rate / f1); const double delta2 = 2.0 * M_PI / ( sample_rate / f2); FILE *f; double theta1, theta2; unsigned i; if(0 == (f = fopen(fn, mode))) return false; for(i = 0, theta1 = theta2 = 0.0; i < samples; i++, theta1 += delta1, theta2 += delta2) { double val = (a1*sin(theta1) + a2*sin(theta2))*(double)full_scale; FLAC__int16 v = (FLAC__int16)(val + 0.5); if(!write_little_endian_int16(f, v)) goto foo; val = -(a1*sin(theta1*fmult) + a2*sin(theta2*fmult))*(double)full_scale; v = (FLAC__int16)(val + 0.5); if(!write_little_endian_int16(f, v)) goto foo; } fclose(f); return true; foo: fclose(f); return false; } /* a mono sine-wave 24bps stream */ static FLAC__bool generate_sine24_1(const char *fn, const double sample_rate, const unsigned samples, const double f1, const double a1, const double f2, const double a2) { const FLAC__int32 full_scale = 0x7fffff; const double delta1 = 2.0 * M_PI / ( sample_rate / f1); const double delta2 = 2.0 * M_PI / ( sample_rate / f2); FILE *f; double theta1, theta2; unsigned i; if(0 == (f = fopen(fn, mode))) return false; for(i = 0, theta1 = theta2 = 0.0; i < samples; i++, theta1 += delta1, theta2 += delta2) { double val = (a1*sin(theta1) + a2*sin(theta2))*(double)full_scale; FLAC__int32 v = (FLAC__int32)(val + 0.5); if(!write_little_endian_int24(f, v)) goto foo; } fclose(f); return true; foo: fclose(f); return false; } /* a stereo sine-wave 24bps stream */ static FLAC__bool generate_sine24_2(const char *fn, const double sample_rate, const unsigned samples, const double f1, const double a1, const double f2, const double a2, double fmult) { const FLAC__int32 full_scale = 0x7fffff; const double delta1 = 2.0 * M_PI / ( sample_rate / f1); const double delta2 = 2.0 * M_PI / ( sample_rate / f2); FILE *f; double theta1, theta2; unsigned i; if(0 == (f = fopen(fn, mode))) return false; for(i = 0, theta1 = theta2 = 0.0; i < samples; i++, theta1 += delta1, theta2 += delta2) { double val = (a1*sin(theta1) + a2*sin(theta2))*(double)full_scale; FLAC__int32 v = (FLAC__int32)(val + 0.5); if(!write_little_endian_int24(f, v)) goto foo; val = -(a1*sin(theta1*fmult) + a2*sin(theta2*fmult))*(double)full_scale; v = (FLAC__int32)(val + 0.5); if(!write_little_endian_int24(f, v)) goto foo; } fclose(f); return true; foo: fclose(f); return false; } static FLAC__bool generate_noise(const char *fn, unsigned bytes) { FILE *f; unsigned b; if(0 == (f = fopen(fn, mode))) return false; for(b = 0; b < bytes; b++) { #if !defined _MSC_VER && !defined __MINGW32__ FLAC__byte x = (FLAC__byte)(((unsigned)random()) & 0xff); #else FLAC__byte x = (FLAC__byte)(((unsigned)rand()) & 0xff); #endif if(fwrite(&x, sizeof(x), 1, f) < 1) goto foo; } fclose(f); return true; foo: fclose(f); return false; } static FLAC__bool generate_aiff(const char *filename, unsigned sample_rate, unsigned channels, unsigned bps, unsigned samples) { const unsigned true_size = channels * ((bps+7)/8) * samples; const unsigned padded_size = (true_size + 1) & (~1u); const unsigned shift = 8 - (bps%8); FILE *f; unsigned i; if(0 == (f = fopen(filename, mode))) return false; if(fwrite("FORM", 1, 4, f) < 4) goto foo; if(!write_big_endian_uint32(f, padded_size + 46)) goto foo; if(fwrite("AIFFCOMM\000\000\000\022", 1, 12, f) < 12) goto foo; if(!write_big_endian_uint16(f, (FLAC__uint16)channels)) goto foo; if(!write_big_endian_uint32(f, samples)) goto foo; if(!write_big_endian_uint16(f, (FLAC__uint16)bps)) goto foo; if(!write_sane_extended(f, sample_rate)) goto foo; if(fwrite("SSND", 1, 4, f) < 4) goto foo; if(!write_big_endian_uint32(f, true_size + 8)) goto foo; if(fwrite("\000\000\000\000\000\000\000\000", 1, 8, f) < 8) goto foo; for(i = 0; i < true_size; i++) if(fputc(GET_RANDOM_BYTE< 2 || (bps%8)); /* ^^^^^^^ * (bps%8) allows 24 bps which is technically supposed to be WAVEFORMATEXTENSIBLE but we * write 24bps as WAVEFORMATEX since it's unambiguous and matches how flac writes it */ const unsigned true_size = channels * ((bps+7)/8) * samples; const unsigned padded_size = (true_size + 1) & (~1u); const unsigned shift = (bps%8)? 8 - (bps%8) : 0; FILE *f; unsigned i; if(0 == (f = fopen(filename, mode))) return false; if(fwrite("RIFF", 1, 4, f) < 4) goto foo; if(!write_little_endian_uint32(f, padded_size + (waveformatextensible?60:36))) goto foo; if(fwrite("WAVEfmt ", 1, 8, f) < 8) goto foo; if(!write_little_endian_uint32(f, waveformatextensible?40:16)) goto foo; if(!write_little_endian_uint16(f, (FLAC__uint16)(waveformatextensible?65534:1))) goto foo; if(!write_little_endian_uint16(f, (FLAC__uint16)channels)) goto foo; if(!write_little_endian_uint32(f, sample_rate)) goto foo; if(!write_little_endian_uint32(f, sample_rate * channels * ((bps+7)/8))) goto foo; if(!write_little_endian_uint16(f, (FLAC__uint16)(channels * ((bps+7)/8)))) /* block align */ goto foo; if(!write_little_endian_uint16(f, (FLAC__uint16)(bps+shift))) goto foo; if(waveformatextensible) { if(!write_little_endian_uint16(f, (FLAC__uint16)22)) /* cbSize */ goto foo; if(!write_little_endian_uint16(f, (FLAC__uint16)bps)) /* validBitsPerSample */ goto foo; if(!write_little_endian_uint32(f, 0)) /* channelMask */ goto foo; /* GUID = {0x00000001, 0x0000, 0x0010, {0x80, 0x00, 0x00, 0xaa, 0x00, 0x38, 0x9b, 0x71}} */ if(fwrite("\x01\x00\x00\x00\x00\x00\x10\x00\x80\x00\x00\xaa\x00\x38\x9b\x71", 1, 16, f) != 16) goto foo; } if(fwrite("data", 1, 4, f) < 4) goto foo; if(!write_little_endian_uint32(f, true_size)) goto foo; for(i = 0; i < true_size; i++) if(fputc(GET_RANDOM_BYTE<