Method writeSampleTableToFile

Source/Utilities.cpp:1209–1339

}

void Utilities::writeSampleTableToFile(const std::string &p_filename) {

	// load in data
	std::string filename = "/home/frot/odin2/samples/" + p_filename;
	File file(filename);
	WavAudioFormat wavAudioFormat;
	AudioFormatReader *audioFormatReader = wavAudioFormat.createReaderFor(file.createInputStream().get(), true);
	float *wavData[2]; // new int *[3];
	wavData[0] = new float[audioFormatReader->lengthInSamples];
	wavData[1] = new float[audioFormatReader->lengthInSamples];
	audioFormatReader->read(wavData, 2, 0, audioFormatReader->lengthInSamples);
	int n_samples = audioFormatReader->lengthInSamples;
#define WAVE wavData[0]

	// remove offsets from beginning start by adding linear * cos function
	float difference = WAVE[n_samples - 1] - WAVE[0];

	for (int i = 0; i < n_samples; ++i) {
		WAVE[i] -= difference * (cos((2 * M_PI * i) / (float)n_samples)) * (float)i / (float)(n_samples);
	}

	// CREATE WAVEFILE

	float wavedraw_coefficients[SIN_AND_COS][NUMBER_OF_HARMONICS];

	float step_width = 2 * PI / (float)n_samples;

	for (int harmonic = 1; harmonic < NUMBER_OF_HARMONICS; ++harmonic) {

		float coeff_sine   = 0.f;
		float coeff_cosine = 0.f;

		for (int segment = 0; segment < n_samples; ++segment) {
			// use either const sections or linear sections

			// wrap function value at end to start
			float segment_end_value = (segment == n_samples - 1) ? WAVE[0] : WAVE[segment + 1];

			coeff_sine += lin_segment_one_overtone_sine(
			    segment * step_width, (segment + 1) * step_width, WAVE[segment], segment_end_value, harmonic);
			coeff_cosine += lin_segment_one_overtone_cosine(
			    segment * step_width, (segment + 1) * step_width, WAVE[segment], segment_end_value, harmonic);
		}
		wavedraw_coefficients[0][harmonic] = coeff_sine;
		wavedraw_coefficients[1][harmonic] = coeff_cosine;
	}

	// now create the wavetable from the fourier coefficients
	double seed_freq                                                 = 27.5; // A0
	float max                                                        = 0.f;
	float wavedraw_tables[SUBTABLES_PER_WAVETABLE][WAVETABLE_LENGTH] = {0};
	float p_samplerate                                               = 44100;

	// loop over subtables
	for (int index_sub_table = 0; index_sub_table < SUBTABLES_PER_WAVETABLE; ++index_sub_table) {

		// how many harmonics are needed for this subtable
		int number_of_harmonics = (int)((p_samplerate * 0.5f / seed_freq) - 1);