Method writeWavetableToFile

Source/Utilities.cpp:971–1055

}

void Utilities::writeWavetableToFile(int index_wavetable) {

	float p_samplerate                                             = 44100.f;
	float wavetables[1][SUBTABLES_PER_WAVETABLE][WAVETABLE_LENGTH] = {0};

	float max = 0;

	// this flags hackishly if a waveform only uses sine components
	bool sine_only = m_fourier_coeffs[index_wavetable][0][0] > 0.5f ? true : false;

	double seed_freq = 27.5; // A0
	// float max = 0.f;

	// 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);

		// don't allow more than 800 harmonics (for big Samplerates this might
		// happen)
		number_of_harmonics = number_of_harmonics > NUMBER_OF_HARMONICS ? NUMBER_OF_HARMONICS : number_of_harmonics;

		for (int index_harmonics = 1; index_harmonics < number_of_harmonics; ++index_harmonics) {
			for (int index_position = 0; index_position < WAVETABLE_LENGTH; ++index_position) {

				// fill table with //sine harmonics
				wavetables[0][index_sub_table][index_position] +=
				    m_fourier_coeffs[index_wavetable][0][index_harmonics] *
				    sin(2.f * PI * index_position * index_harmonics / (float)WAVETABLE_LENGTH) *
				    m_fourier_coeffs[index_wavetable][1][0]; // last term is normalization
				if (!sine_only) {
					// cosine
					wavetables[0][index_sub_table][index_position] +=
					    m_fourier_coeffs[index_wavetable][1][index_harmonics] *
					    cos(2.f * PI * index_position * index_harmonics / (float)WAVETABLE_LENGTH) *
					    m_fourier_coeffs[index_wavetable][1][0]; // last term is normalization
				}
			}
		}

		for (int index_sub_table_arr = 0; index_sub_table_arr < SUBTABLES_PER_WAVETABLE; ++index_sub_table_arr) {
			for (int index_position = 0; index_position < WAVETABLE_LENGTH; ++index_position) {
				if (fabs(wavetables[0][index_sub_table_arr][index_position]) > max) {
					max = fabs(wavetables[0][index_sub_table_arr][index_position]);
				}
			}
		}

		// increment seed frequency by minor third = 2^(3/12)
		seed_freq *= 1.1892071150;
	}

	// avoid division by zero
	max = max < 0.0001 ? 1 : max;

	std::ofstream output_file;
	output_file.open("/home/frot/odin2/Source/audio/Oscillators/"