void FileInMemory_showAsCode (FileInMemory me, const char32 *name, long numberOfBytesPerLine)
{
if (numberOfBytesPerLine <= 0) numberOfBytesPerLine = 20;
// autoNUMvector<unsigned char> data (0, my d_numberOfBytes); ????
MelderInfo_writeLine (U"\t\tstatic unsigned char ", name, U"_data[", my d_numberOfBytes+1, U"] = {");
for (long i = 0; i < my d_numberOfBytes; i++) {
unsigned char number = my d_data[i];
MelderInfo_write ((i % numberOfBytesPerLine == 0 ? U"\t\t\t" : U""), number, U",",
((i % numberOfBytesPerLine == (numberOfBytesPerLine - 1)) ? U"\n" : U" "));
}
MelderInfo_writeLine ((my d_numberOfBytes - 1) % numberOfBytesPerLine == (numberOfBytesPerLine - 1) ? U"\t\t\t0};" : U"0};");
MelderInfo_write (U"\t\tautoFileInMemory ", name, U" = FileInMemory_createWithData (");
MelderInfo_writeLine (my d_numberOfBytes, U", reinterpret_cast<const char *> (&", name, U"_data), \n\t\t\tU\"", my d_path, U"\", \n\t\t\tU\"", my d_id, U"\");");
}
void FileInMemory_showAsCode (FileInMemory me, const wchar_t *name, long numberOfBytesPerLine)
{
if (numberOfBytesPerLine <= 0) numberOfBytesPerLine = 20;
// autoNUMvector<unsigned char> data (0, my d_numberOfBytes); ????
MelderInfo_writeLine (L"\t\tstatic unsigned char ", name, L"_data[", Melder_integer (my d_numberOfBytes+1), L"] = {");
for (long i = 0; i < my d_numberOfBytes; i++) {
unsigned char number = my d_data[i];
MelderInfo_write ((i % numberOfBytesPerLine == 0 ? L"\t\t\t" : L""), Melder_integer (number), L",",
((i % numberOfBytesPerLine == (numberOfBytesPerLine - 1)) ? L"\n" : L" "));
}
MelderInfo_writeLine ((my d_numberOfBytes - 1) % numberOfBytesPerLine == (numberOfBytesPerLine - 1) ? L"\t\t\t0};" : L"0};");
MelderInfo_write (L"\t\tautoFileInMemory ", name, L" = FileInMemory_createWithData (");
MelderInfo_writeLine (Melder_integer (my d_numberOfBytes), L", reinterpret_cast<const char *> (&", name, L"_data), \n\t\t\tL\"", my d_path, L"\", \n\t\t\tL\"", my d_id, L"\");");
}
void structExcitation :: v_info () {
double *y = z [1];
long numberOfMaxima = 0;
structData :: v_info ();
MelderInfo_writeLine (L"Loudness: ", Melder_half (Excitation_getLoudness (this)), L" sones");
for (long i = 2; i < nx; i ++) if (y [i] > y [i - 1] && y [i] >= y [i + 1]) {
double i_real, formant_bark, strength;
if (++ numberOfMaxima > 15) break;
strength = NUMimproveMaximum (z [1], nx, i, NUM_PEAK_INTERPOLATE_SINC70, & i_real);
formant_bark = x1 + (i_real - 1) * dx;
MelderInfo_write (L"Peak at ", Melder_single (formant_bark), L" Bark");
MelderInfo_write (L", ", Melder_integer ((long) NUMbarkToHertz (formant_bark)), L" Hz");
MelderInfo_writeLine (L", ", Melder_half (strength), L" phon.");
}
}
void structLogisticRegression :: v_info () {
LogisticRegression_Parent :: v_info ();
MelderInfo_writeLine (U"Dependent 1: ", our dependent1);
MelderInfo_writeLine (U"Dependent 2: ", our dependent2);
MelderInfo_writeLine (U"Interpretation:");
MelderInfo_write (U" ln (P(", dependent2, U")/P(", dependent1, U")) " UNITEXT_ALMOST_EQUAL_TO U" ", Melder_fixed (intercept, 6));
for (long ivar = 1; ivar <= parameters.size; ivar ++) {
RegressionParameter parm = parameters.at [ivar];
MelderInfo_write (parm -> value < 0.0 ? U" - " : U" + ", Melder_fixed (fabs (parm -> value), 6), U" * ", parm -> label);
}
MelderInfo_writeLine (U"");
MelderInfo_writeLine (U"Log odds ratios:");
for (long ivar = 1; ivar <= parameters.size; ivar ++) {
RegressionParameter parm = parameters.at [ivar];
MelderInfo_writeLine (U" Log odds ratio of factor ", parm -> label, U": ", Melder_fixed ((parm -> maximum - parm -> minimum) * parm -> value, 6));
}
MelderInfo_writeLine (U"Odds ratios:");
for (long ivar = 1; ivar <= parameters.size; ivar ++) {
RegressionParameter parm = parameters.at [ivar];
MelderInfo_writeLine (U" Odds ratio of factor ", parm -> label, U": ", exp ((parm -> maximum - parm -> minimum) * parm -> value));
}
}
void structPitch :: v_info () {
long nVoiced;
autoNUMvector <double> frequencies (Sampled_getSortedValues (this, Pitch_LEVEL_FREQUENCY, kPitch_unit_HERTZ, & nVoiced), 1);
structDaata :: v_info ();
MelderInfo_writeLine (U"Time domain:");
MelderInfo_writeLine (U" Start time: ", xmin, U" seconds");
MelderInfo_writeLine (U" End time: ", xmax, U" seconds");
MelderInfo_writeLine (U" Total duration: ", xmax - xmin, U" seconds");
MelderInfo_writeLine (U"Time sampling:");
MelderInfo_writeLine (U" Number of frames: ", nx, U" (", nVoiced, U" voiced)");
MelderInfo_writeLine (U" Time step: ", dx, U" seconds");
MelderInfo_writeLine (U" First frame centred at: ", x1, U" seconds");
MelderInfo_writeLine (U"Ceiling at: ", ceiling, U" Hz");
if (nVoiced >= 1) { // quantiles
double quantile10, quantile16, quantile50, quantile84, quantile90;
quantile10 = NUMquantile (nVoiced, frequencies.peek(), 0.10);
quantile16 = NUMquantile (nVoiced, frequencies.peek(), 0.16);
quantile50 = NUMquantile (nVoiced, frequencies.peek(), 0.50); // median
quantile84 = NUMquantile (nVoiced, frequencies.peek(), 0.84);
quantile90 = NUMquantile (nVoiced, frequencies.peek(), 0.90);
MelderInfo_writeLine (U"\nEstimated quantiles:");
MelderInfo_write (U" 10% = ", Melder_single (quantile10), U" Hz = ", Melder_single (MEL (quantile10)), U" Mel = ");
MelderInfo_writeLine (Melder_single (SEMITONES (quantile10)), U" semitones above 100 Hz = ", Melder_single (ERB (quantile10)), U" ERB");
MelderInfo_write (U" 16% = ", Melder_single (quantile16), U" Hz = ", Melder_single (MEL (quantile16)), U" Mel = ");
MelderInfo_writeLine (Melder_single (SEMITONES (quantile16)), U" semitones above 100 Hz = ", Melder_single (ERB (quantile16)), U" ERB");
MelderInfo_write (U" 50% = ", Melder_single (quantile50), U" Hz = ", Melder_single (MEL (quantile50)), U" Mel = ");
MelderInfo_writeLine (Melder_single (SEMITONES (quantile50)), U" semitones above 100 Hz = ", Melder_single (ERB (quantile50)), U" ERB");
MelderInfo_write (U" 84% = ", Melder_single (quantile84), U" Hz = ", Melder_single (MEL (quantile84)), U" Mel = ");
MelderInfo_writeLine (Melder_single (SEMITONES (quantile84)), U" semitones above 100 Hz = ", Melder_single (ERB (quantile84)), U" ERB");
MelderInfo_write (U" 90% = ", Melder_single (quantile90), U" Hz = ", Melder_single (MEL (quantile90)), U" Mel = ");
MelderInfo_writeLine (Melder_single (SEMITONES (quantile90)), U" semitones above 100 Hz = ", Melder_single (ERB (quantile90)), U" ERB");
if (nVoiced > 1) {
double corr = sqrt (nVoiced / (nVoiced - 1.0));
MelderInfo_writeLine (U"\nEstimated spreading:");
MelderInfo_write (U" 84%-median = ", Melder_half ((quantile84 - quantile50) * corr), U" Hz = ", Melder_half ((MEL (quantile84) - MEL (quantile50)) * corr), U" Mel = ");
MelderInfo_writeLine (Melder_half ((SEMITONES (quantile84) - SEMITONES (quantile50)) * corr), U" semitones = ", Melder_half ((ERB (quantile84) - ERB (quantile50)) * corr), U" ERB");
MelderInfo_write (U" median-16% = ", Melder_half ((quantile50 - quantile16) * corr), U" Hz = ", Melder_half ((MEL (quantile50) - MEL (quantile16)) * corr), U" Mel = ");
MelderInfo_writeLine (Melder_half ((SEMITONES (quantile50) - SEMITONES (quantile16)) * corr), U" semitones = ", Melder_half ((ERB (quantile50) - ERB (quantile16)) * corr), U" ERB");
MelderInfo_write (U" 90%-10% = ", Melder_half ((quantile90 - quantile10) * corr), U" Hz = ", Melder_half ((MEL (quantile90) - MEL (quantile10)) * corr), U" Mel = ");
MelderInfo_writeLine (Melder_half ((SEMITONES (quantile90) - SEMITONES (quantile10)) * corr), U" semitones = ", Melder_half ((ERB (quantile90) - ERB (quantile10)) * corr), U" ERB");
}
}
if (nVoiced >= 1) { // extrema, range, mean and standard deviation
double minimum = Pitch_getMinimum (this, xmin, xmax, kPitch_unit_HERTZ, false);
double maximum = Pitch_getMaximum (this, xmin, xmax, kPitch_unit_HERTZ, false);
double meanHertz, meanMel, meanSemitones, meanErb;
MelderInfo_write (U"\nMinimum ", Melder_single (minimum), U" Hz = ", Melder_single (MEL (minimum)), U" Mel = ");
MelderInfo_writeLine (Melder_single (SEMITONES (minimum)), U" semitones above 100 Hz = ", Melder_single (ERB (minimum)), U" ERB");
MelderInfo_write (U"Maximum ", Melder_single (maximum), U" Hz = ", Melder_single (MEL (maximum)), U" Mel = ");
MelderInfo_writeLine (Melder_single (SEMITONES (maximum)), U" semitones above 100 Hz = ", Melder_single (ERB (maximum)), U" ERB");
MelderInfo_write (U"Range ", Melder_half (maximum - minimum), U" Hz = ", Melder_single (MEL (maximum) - MEL (minimum)), U" Mel = ");
MelderInfo_writeLine (Melder_half (SEMITONES (maximum) - SEMITONES (minimum)), U" semitones = ", Melder_half (ERB (maximum) - ERB (minimum)), U" ERB");
meanHertz = Pitch_getMean (this, 0, 0, kPitch_unit_HERTZ);
meanMel = Pitch_getMean (this, 0, 0, kPitch_unit_MEL);
meanSemitones = Pitch_getMean (this, 0, 0, kPitch_unit_SEMITONES_100);
meanErb = Pitch_getMean (this, 0, 0, kPitch_unit_ERB);
MelderInfo_write (U"Average: ", Melder_single (meanHertz), U" Hz = ", Melder_single (meanMel), U" Mel = ");
MelderInfo_writeLine (Melder_single (meanSemitones), U" semitones above 100 Hz = ", Melder_single (meanErb), U" ERB");
if (nVoiced >= 2) {
double stdevHertz = Pitch_getStandardDeviation (this, 0, 0, kPitch_unit_HERTZ);
double stdevMel = Pitch_getStandardDeviation (this, 0, 0, kPitch_unit_MEL);
double stdevSemitones = Pitch_getStandardDeviation (this, 0, 0, kPitch_unit_SEMITONES_100);
double stdevErb = Pitch_getStandardDeviation (this, 0, 0, kPitch_unit_ERB);
MelderInfo_write (U"Standard deviation: ", Melder_half (stdevHertz), U" Hz = ", Melder_half (stdevMel), U" Mel = ");
MelderInfo_writeLine (Melder_half (stdevSemitones), U" semitones = ", Melder_half (stdevErb), U" ERB");
}
}
if (nVoiced > 1) { // variability: mean absolute slope
double slopeHertz, slopeMel, slopeSemitones, slopeErb, slopeWithoutOctaveJumps;
Pitch_getMeanAbsoluteSlope (this, & slopeHertz, & slopeMel, & slopeSemitones, & slopeErb, & slopeWithoutOctaveJumps);
MelderInfo_write (U"\nMean absolute slope: ", Melder_half (slopeHertz), U" Hz/s = ", Melder_half (slopeMel), U" Mel/s = ");
MelderInfo_writeLine (Melder_half (slopeSemitones), U" semitones/s = ", Melder_half (slopeErb), U" ERB/s");
MelderInfo_writeLine (U"Mean absolute slope without octave jumps: ", Melder_half (slopeWithoutOctaveJumps), U" semitones/s");
}
}
void structSound :: v_info () {
structData :: v_info ();
const double rho_c = 400; /* rho = 1.14 kg m-3; c = 353 m s-1; [rho c] = kg m-2 s-1 */
double minimum = z [1] [1], maximum = minimum;
MelderInfo_writeLine (U"Number of channels: ", ny, ny == 1 ? U" (mono)" : ny == 2 ? U" (stereo)" : U"");
MelderInfo_writeLine (U"Time domain:");
MelderInfo_writeLine (U" Start time: ", xmin, U" seconds");
MelderInfo_writeLine (U" End time: ", xmax, U" seconds");
MelderInfo_writeLine (U" Total duration: ", xmax - xmin, U" seconds");
MelderInfo_writeLine (U"Time sampling:");
MelderInfo_writeLine (U" Number of samples: ", nx);
MelderInfo_writeLine (U" Sampling period: ", dx, U" seconds");
MelderInfo_writeLine (U" Sampling frequency: ", Melder_single (1.0 / dx), U" Hz");
MelderInfo_writeLine (U" First sample centred at: ", x1, U" seconds");
{// scope
double sum = 0.0, sumOfSquares = 0.0;
for (long channel = 1; channel <= ny; channel ++) {
double *amplitude = z [channel];
for (long i = 1; i <= nx; i ++) {
double value = amplitude [i];
sum += value;
sumOfSquares += value * value;
if (value < minimum) minimum = value;
if (value > maximum) maximum = value;
}
}
MelderInfo_writeLine (U"Amplitude:");
MelderInfo_writeLine (U" Minimum: ", Melder_single (minimum), U" Pascal");
MelderInfo_writeLine (U" Maximum: ", Melder_single (maximum), U" Pascal");
double mean = sum / (nx * ny);
MelderInfo_writeLine (U" Mean: ", Melder_single (mean), U" Pascal");
MelderInfo_writeLine (U" Root-mean-square: ", Melder_single (sqrt (sumOfSquares / (nx * ny))), U" Pascal");
double penergy = sumOfSquares * dx / ny; /* Pa2 s = kg2 m-2 s-3 */
MelderInfo_write (U"Total energy: ", Melder_single (penergy), U" Pascal\u00B2 sec");
double energy = penergy / rho_c; /* kg s-2 = Joule m-2 */
MelderInfo_writeLine (U" (energy in air: ", Melder_single (energy), U" Joule/m\u00B2)");
double power = energy / (dx * nx); /* kg s-3 = Watt/m2 */
MelderInfo_write (U"Mean power (intensity) in air: ", Melder_single (power), U" Watt/m\u00B2");
if (power != 0.0) {
MelderInfo_writeLine (U" = ", Melder_half (10 * log10 (power / 1e-12)), U" dB");
} else {
MelderInfo_writeLine (U"");
}
}
if (nx > 1) {
for (long channel = 1; channel <= ny; channel ++) {
double *amplitude = z [channel];
double sum = 0.0;
for (long i = 1; i <= nx; i ++) {
double value = amplitude [i];
sum += value;
}
double mean = sum / nx, stdev = 0.0;
for (long i = 1; i <= nx; i ++) {
double value = amplitude [i] - mean;
stdev += value * value;
}
stdev = sqrt (stdev / (nx - 1));
MelderInfo_writeLine (U"Standard deviation in channel ", channel, U": ", Melder_single (stdev), U" Pascal");
}
}
}
