static void gui_button_cb_insert (I, GuiButtonEvent event) {
(void) event;
iam (StringsEditor);
Strings strings = (Strings) my data;
/*
* Find the first selected item.
*/
long numberOfSelected, *selected = GuiList_getSelectedPositions (my list, & numberOfSelected);
long position = selected == NULL ? strings -> numberOfStrings + 1 : selected [1];
NUMvector_free (selected, 1);
wchar_t *text = GuiText_getString (my text);
/*
* Change the data.
*/
Strings_insert (strings, position, text);
/*
* Change the list.
*/
GuiList_insertItem (my list, text, position);
GuiList_deselectAllItems (my list);
GuiList_selectItem (my list, position);
/*
* Clean up.
*/
Melder_free (text);
my broadcastDataChanged ();
}
void TableOfReal_insertRow (TableOfReal me, long rowNumber) {
try {
if (rowNumber < 1 || rowNumber > my numberOfRows + 1)
Melder_throw ("Cannot create row ", rowNumber, ".");
autoNUMmatrix <double> data (1, my numberOfRows + 1, 1, my numberOfColumns);
autoNUMvector <wchar_t *> rowLabels (1, my numberOfRows + 1); // not autostringvector...
for (long irow = 1; irow < rowNumber; irow ++) {
rowLabels [irow] = my rowLabels [irow]; // ...because this is a dangling copy
for (long icol = 1; icol <= my numberOfColumns; icol ++)
data [irow] [icol] = my data [irow] [icol];
}
for (long irow = my numberOfRows + 1; irow > rowNumber; irow --) {
rowLabels [irow] = my rowLabels [irow - 1];
for (long icol = 1; icol <= my numberOfColumns; icol ++)
data [irow] [icol] = my data [irow - 1] [icol];
}
/*
* Change without error.
*/
NUMvector_free (my rowLabels, 1);
my rowLabels = rowLabels.transfer();
NUMmatrix_free (my data, 1, 1);
my data = data.transfer();
my numberOfRows ++;
} catch (MelderError) {
Melder_throw (me, ": row at position ", rowNumber, " not inserted.");
}
}
void PointProcess_addPoint (PointProcess me, double t) {
try {
if (t == NUMundefined)
Melder_throw (U"Cannot add a point at an undefined time.");
if (my nt >= my maxnt) {
/*
* Create without change.
*/
autoNUMvector <double> dum (1, 2 * my maxnt);
NUMvector_copyElements (my t, dum.peek(), 1, my nt);
/*
* Change without error.
*/
NUMvector_free (my t, 1);
my t = dum.transfer();
my maxnt *= 2;
}
if (my nt == 0 || t >= my t [my nt]) { // special case that often occurs in practice
my t [++ my nt] = t;
} else {
long left = PointProcess_getLowIndex (me, t);
if (left == 0 || my t [left] != t) {
for (long i = my nt; i > left; i --) my t [i + 1] = my t [i];
my nt ++;
my t [left + 1] = t;
}
}
} catch (MelderError) {
Melder_throw (me, U": point not added.");
}
}
static void gui_button_cb_insert (StringsEditor me, GuiButtonEvent /* event */) {
Strings strings = (Strings) my data;
/*
* Find the first selected item.
*/
long numberOfSelected, *selected = GuiList_getSelectedPositions (my list, & numberOfSelected);
long position = selected ? selected [1] : strings -> numberOfStrings + 1;
NUMvector_free (selected, 1);
char32 *text = GuiText_getString (my text);
/*
* Change the data.
*/
Strings_insert (strings, position, text);
/*
* Change the list.
*/
GuiList_insertItem (my list, text, position);
GuiList_deselectAllItems (my list);
GuiList_selectItem (my list, position);
/*
* Clean up.
*/
Melder_free (text);
Editor_broadcastDataChanged (me);
}
void TableOfReal_insertColumn (TableOfReal me, long columnNumber) {
try {
if (columnNumber < 1 || columnNumber > my numberOfColumns + 1)
Melder_throw ("Cannot create column ", columnNumber, ".");
autoNUMmatrix <double> data (1, my numberOfRows, 1, my numberOfColumns + 1);
autoNUMvector <wchar_t*> columnLabels (1, my numberOfColumns + 1); // not autostringvector...
for (long j = 1; j < columnNumber; j ++) {
columnLabels [j] = my columnLabels [j]; // ...because this is a dangling copy
for (long i = 1; i <= my numberOfRows; i ++) data [i] [j] = my data [i] [j];
}
for (long j = my numberOfColumns + 1; j > columnNumber; j --) {
columnLabels [j] = my columnLabels [j - 1];
for (long i = 1; i <= my numberOfRows; i ++) data [i] [j] = my data [i] [j - 1];
}
/*
* Change without error.
*/
NUMvector_free (my columnLabels, 1);
my columnLabels = columnLabels.transfer();
NUMmatrix_free (my data, 1, 1);
my data = data.transfer();
my numberOfColumns ++;
} catch (MelderError) {
Melder_throw (me, ": column at position ", columnNumber, " not inserted.");
}
}
static int tryAdoption (int **distance, int *path, int numberOfCities, long *totalDistance)
{
int *help = NUMvector <int> (0, numberOfCities);
int i, maximumGainLeft, result = 0;
/* Compute maximum distance between two successive cities. */
int city1 = path [0], city2 = path [1];
int maximumDistance = distance [city1] [city2];
for (i = 2; i <= numberOfCities; i ++) {
city1 = city2;
city2 = path [i];
if (distance [city1] [city2] > maximumDistance)
maximumDistance = distance [city1] [city2];
}
maximumGainLeft = maximumDistance;
for (i = 1; i <= numberOfCities; i ++) {
int cont = 1, b1, b2, distance_b1_b2, d1nr = 3, cc, e1nrMax = 6;
int numberOfCitiesMinus1 = numberOfCities - 1, j;
for (j = 0; j <= numberOfCitiesMinus1; j ++) path [j] = path [j + 1];
path [numberOfCities] = path [0];
b1 = path [0];
b2 = path [1];
distance_b1_b2 = distance [b1] [b2];
cc = path [2];
while (d1nr < numberOfCitiesMinus1 && cont) {
int d1 = path [d1nr];
int gain1 = distance_b1_b2 + distance [d1] [cc] - distance [d1] [b2];
if (gain1 + maximumGainLeft > 0) {
int e1nr = d1nr + 1;
int dn = path [d1nr];
if (e1nrMax > numberOfCitiesMinus1) e1nrMax = numberOfCitiesMinus1;
while (e1nr < e1nrMax && cont) {
int e1 = path [e1nr];
int gain = gain1 + distance [dn] [e1] - distance [dn] [b1] - distance [cc] [e1];
if (gain > 0) {
int nAdoption = e1nr - d1nr;
int dnnr = e1nr - 1;
cont = 0;
*totalDistance -= gain;
for (j = 0; j <= dnnr - 1; j ++) help [j] = path [j + 1];
for (j = 1; j <= nAdoption; j ++) path [j] = help [dnnr - j];
for (j = 0; j <= d1nr - 2; j ++) path [nAdoption + j + 1] = help [j];
}
dn = e1;
e1nr ++;
}
}
e1nrMax ++;
cc = d1;
d1nr ++;
}
result |= ! cont;
}
NUMvector_free (help, 0);
return result;
}
static void gui_button_cb_remove (StringsEditor me, GuiButtonEvent /* event */) {
long numberOfSelected, *selected = GuiList_getSelectedPositions (my list, & numberOfSelected);
for (long iselected = numberOfSelected; iselected >= 1; iselected --) {
Strings_remove ((Strings) my data, selected [iselected]);
}
NUMvector_free (selected, 1);
updateList (me);
Editor_broadcastDataChanged (me);
}
void Sound_playPart (Sound me, double tmin, double tmax,
int (*callback) (void *closure, int phase, double tmin, double tmax, double t), void *closure)
{
try {
long ifsamp = lround (1.0 / my dx), bestSampleRate = MelderAudio_getOutputBestSampleRate (ifsamp);
if (ifsamp == bestSampleRate) {
struct SoundPlay *thee = (struct SoundPlay *) & thePlayingSound;
double *fromLeft = my z [1], *fromRight = my ny > 1 ? my z [2] : NULL;
MelderAudio_stopPlaying (MelderAudio_IMPLICIT);
long i1, i2;
if ((thy numberOfSamples = Matrix_getWindowSamplesX (me, tmin, tmax, & i1, & i2)) < 1) return;
thy tmin = tmin;
thy tmax = tmax;
thy dt = my dx;
thy t1 = my x1;
thy callback = callback;
thy closure = closure;
thy silenceBefore = (long) (ifsamp * MelderAudio_getOutputSilenceBefore ());
thy silenceAfter = (long) (ifsamp * MelderAudio_getOutputSilenceAfter ());
int numberOfChannels = my ny;
NUMvector_free (thy buffer, 1); // just in case
thy buffer = NUMvector <short> (1, (i2 - i1 + 1 + thy silenceBefore + thy silenceAfter) * numberOfChannels);
thy i1 = i1;
thy i2 = i2;
short *to = thy buffer + thy silenceBefore * numberOfChannels;
if (numberOfChannels > 2) {
for (long i = i1; i <= i2; i ++) {
for (long chan = 1; chan <= my ny; chan ++) {
long value = (long) round (my z [chan] [i] * 32768.0);
* ++ to = value < -32768 ? -32768 : value > 32767 ? 32767 : value;
}
}
} else if (numberOfChannels == 2) {
for (long i = i1; i <= i2; i ++) {
long valueLeft = (long) round (fromLeft [i] * 32768.0);
* ++ to = valueLeft < -32768 ? -32768 : valueLeft > 32767 ? 32767 : valueLeft;
long valueRight = (long) round (fromRight [i] * 32768.0);
* ++ to = valueRight < -32768 ? -32768 : valueRight > 32767 ? 32767 : valueRight;
}
} else {
for (long i = i1; i <= i2; i ++) {
long value = (long) round (fromLeft [i] * 32768.0);
* ++ to = value < -32768 ? -32768 : value > 32767 ? 32767 : value;
}
}
if (thy callback) thy callback (thy closure, 1, tmin, tmax, tmin);
MelderAudio_play16 (thy buffer + 1, ifsamp,
thy silenceBefore + thy numberOfSamples + thy silenceAfter, numberOfChannels, melderPlayCallback, thee);
} else {
autoSound resampled = Sound_resample (me, bestSampleRate, 1);
Sound_playPart (resampled.peek(), tmin, tmax, callback, closure); // recursively
}
} catch (MelderError) {
Melder_throw (me, U": not played.");
}
}
void Configuration_rotateToPrincipalDirections (Configuration me) {
try {
autoNUMmatrix<double> m (NUMmatrix_copy (my data, 1, my numberOfRows, 1, my numberOfColumns), 1, 1);
NUMdmatrix_into_principalComponents (my data, my numberOfRows, my numberOfColumns, my numberOfColumns, m.peek());
NUMvector_free (my data, 1);
my data = m.transfer();
} catch (MelderError) {
Melder_throw (me, U": not rotated to principal directions.");
}
}
void Pitch_Frame_init (Pitch_Frame me, int nCandidates) {
/*
* Create without change.
*/
autoNUMvector <structPitch_Candidate> candidate (1, nCandidates);
/*
* Change without error.
*/
NUMvector_free (my candidate, 1);
my candidate = candidate.transfer();
my nCandidates = nCandidates;
}
void structVDSmagtMinimizer :: v_destroy () {
NUMvector_free (dp, 1);
NUMvector_free (pc, 1);
NUMvector_free (gc, 1);
NUMvector_free (g0, 1);
NUMvector_free (s, 1);
NUMvector_free (srst, 1);
NUMvector_free (grst, 1);
VDSmagtMinimizer_Parent :: v_destroy ();
}
static bool melderPlayCallback (void *closure, long samplesPlayed) {
struct SoundPlay *me = (struct SoundPlay *) closure;
int phase = 2;
double t = samplesPlayed <= my silenceBefore ? my tmin :
samplesPlayed >= my silenceBefore + my numberOfSamples ? my tmax :
my t1 + (my i1 - 1.5 + samplesPlayed - my silenceBefore) * my dt;
if (! MelderAudio_isPlaying) {
NUMvector_free (my buffer, 1), my buffer = NULL;
phase = 3;
}
if (my callback)
return my callback (my closure, phase, my tmin, my tmax, t);
return true;
}
static long Pitch_getMeanAbsoluteSlope (Pitch me,
double *out_hertz, double *out_mel, double *out_semitones, double *out_erb, double *out_withoutOctaveJumps)
{
long firstVoicedFrame = 0, lastVoicedFrame = 0, nVoiced = 0, i;
double *frequencies = NUMvector <double> (1, my nx);
for (i = 1; i <= my nx; i ++) {
double frequency = my frame [i]. candidate [1]. frequency;
frequencies [i] = frequency > 0.0 && frequency < my ceiling ? frequency : 0.0;
if (frequencies [i]) nVoiced ++;
}
for (i = 1; i <= my nx; i ++) /* Look for first voiced frame. */
if (frequencies [i] != 0.0) { firstVoicedFrame = i; break; }
for (i = my nx; i >= 1; i --) /* Look for last voiced frame. */
if (frequencies [i] != 0.0) { lastVoicedFrame = i; break; }
if (nVoiced > 1) {
int ilast = firstVoicedFrame;
double span = (lastVoicedFrame - firstVoicedFrame) * my dx, flast = frequencies [ilast];
double slopeHz = 0, slopeMel = 0, slopeSemitones = 0, slopeErb = 0, slopeRobust = 0;
for (i = firstVoicedFrame + 1; i <= lastVoicedFrame; i ++) if (frequencies [i] != 0.0) {
double localStepSemitones = fabs (SEMITONES (frequencies [i]) - SEMITONES (flast));
slopeHz += fabs (frequencies [i] - flast);
slopeMel += fabs (MEL (frequencies [i]) - MEL (flast));
slopeSemitones += localStepSemitones;
slopeErb += fabs (ERB (frequencies [i]) - ERB (flast));
while (localStepSemitones >= 12.0) localStepSemitones -= 12.0; /* Kill octave jumps. */
if (localStepSemitones > 6.0) localStepSemitones = 12.0 - localStepSemitones;
slopeRobust += localStepSemitones;
ilast = i;
flast = frequencies [i];
}
if (out_hertz) *out_hertz = slopeHz / span;
if (out_mel) *out_mel = slopeMel / span;
if (out_semitones) *out_semitones = slopeSemitones / span;
if (out_erb) *out_erb = slopeErb / span;
if (out_withoutOctaveJumps) *out_withoutOctaveJumps = slopeRobust / span;
} else {
if (out_hertz) *out_hertz = NUMundefined;
if (out_mel) *out_mel = NUMundefined;
if (out_semitones) *out_semitones = NUMundefined;
if (out_erb) *out_erb = NUMundefined;
if (out_withoutOctaveJumps) *out_withoutOctaveJumps = NUMundefined;
}
NUMvector_free (frequencies, 1);
return nVoiced;
}
void NUMvector_insert (long elementSize, void **v, long lo, long *hi, long position) {
try {
char *result;
if (! *v) {
result = reinterpret_cast <char *> (NUMvector (elementSize, lo, lo));
*hi = lo;
Melder_assert (position == lo);
} else {
result = reinterpret_cast <char *> (NUMvector (elementSize, lo, *hi + 1));
Melder_assert (position >= lo && position <= *hi + 1);
NUMvector_copyElements (elementSize, *v, result, lo, position - 1);
NUMvector_copyElements (elementSize, *v, result + elementSize, position, *hi);
NUMvector_free (elementSize, *v, lo);
(*hi) ++;
}
*v = result;
} catch (MelderError) {
Melder_throw (U"Vector: element not inserted.");
}
}
void ExperimentMFC_start (ExperimentMFC me) {
try {
long maximumStimulusPlaySamples, maximumResponsePlaySamples, maximumPlaySamples;
long stimulusCarrierBeforeSamples = 0, stimulusCarrierAfterSamples = 0, maximumStimulusSamples = 0;
long responseCarrierBeforeSamples = 0, responseCarrierAfterSamples = 0, maximumResponseSamples = 0;
Melder_warningOff ();
my trial = 0;
NUMvector_free <long> (my stimuli, 1);
NUMvector_free <long> (my responses, 1);
NUMvector_free <double> (my goodnesses, 1);
NUMvector_free <double> (my reactionTimes, 1);
my playBuffer.reset();
my pausing = false;
my numberOfTrials = my numberOfDifferentStimuli * my numberOfReplicationsPerStimulus;
my stimuli = NUMvector <long> (1, my numberOfTrials);
my responses = NUMvector <long> (1, my numberOfTrials);
my goodnesses = NUMvector <double> (1, my numberOfTrials);
my reactionTimes = NUMvector <double> (1, my numberOfTrials);
/*
* Read all the sounds. They must all have the same sampling frequency and number of channels.
*/
my samplePeriod = 0.0;
my numberOfChannels = 0;
if (my stimuliAreSounds) {
if (my stimulusCarrierBefore. name && my stimulusCarrierBefore. name [0]) {
readSound (me, my stimulusFileNameHead, my stimulusFileNameTail, my stimulusMedialSilenceDuration,
& my stimulusCarrierBefore. name, & my stimulusCarrierBefore. sound);
stimulusCarrierBeforeSamples = my stimulusCarrierBefore. sound -> nx;
}
if (my stimulusCarrierAfter. name && my stimulusCarrierAfter. name [0]) {
readSound (me, my stimulusFileNameHead, my stimulusFileNameTail, my stimulusMedialSilenceDuration,
& my stimulusCarrierAfter. name, & my stimulusCarrierAfter. sound);
stimulusCarrierAfterSamples = my stimulusCarrierAfter. sound -> nx;
}
for (long istim = 1; istim <= my numberOfDifferentStimuli; istim ++) {
readSound (me, my stimulusFileNameHead, my stimulusFileNameTail, my stimulusMedialSilenceDuration,
& my stimulus [istim]. name, & my stimulus [istim]. sound);
if (my stimulus [istim]. sound -> nx > maximumStimulusSamples)
maximumStimulusSamples = my stimulus [istim]. sound -> nx;
}
}
if (my responsesAreSounds) {
if (my responseCarrierBefore. name && my responseCarrierBefore. name [0]) {
readSound (me, my responseFileNameHead, my responseFileNameTail, my responseMedialSilenceDuration,
& my responseCarrierBefore. name, & my responseCarrierBefore. sound);
responseCarrierBeforeSamples = my responseCarrierBefore. sound -> nx;
}
if (my responseCarrierAfter. name && my responseCarrierAfter. name [0]) {
readSound (me, my responseFileNameHead, my responseFileNameTail, my responseMedialSilenceDuration,
& my responseCarrierAfter. name, & my responseCarrierAfter. sound);
responseCarrierAfterSamples = my responseCarrierAfter. sound -> nx;
}
for (long iresp = 1; iresp <= my numberOfDifferentResponses; iresp ++) {
readSound (me, my responseFileNameHead, my responseFileNameTail, my responseMedialSilenceDuration,
& my response [iresp]. name, & my response [iresp]. sound);
if (my response [iresp]. sound -> nx > maximumResponseSamples)
maximumResponseSamples = my response [iresp]. sound -> nx;
}
}
/*
* Create the play buffer.
*/
maximumStimulusPlaySamples =
lround (my stimulusInitialSilenceDuration / my samplePeriod)
+ lround (my stimulusFinalSilenceDuration / my samplePeriod)
+ stimulusCarrierBeforeSamples + maximumStimulusSamples + stimulusCarrierAfterSamples + 2;
maximumResponsePlaySamples =
lround (my responseInitialSilenceDuration / my samplePeriod)
+ lround (my responseFinalSilenceDuration / my samplePeriod)
+ responseCarrierBeforeSamples + maximumResponseSamples + responseCarrierAfterSamples + 2;
maximumPlaySamples = maximumStimulusPlaySamples > maximumResponsePlaySamples ? maximumStimulusPlaySamples : maximumResponsePlaySamples;
my playBuffer = Sound_create (my numberOfChannels, 0.0, maximumPlaySamples * my samplePeriod,
maximumPlaySamples, my samplePeriod, 0.5 * my samplePeriod);
/*
* Determine the order in which the stimuli will be presented to the subject.
*/
if (my randomize == kExperiment_randomize_CYCLIC_NON_RANDOM) {
for (long itrial = 1; itrial <= my numberOfTrials; itrial ++)
my stimuli [itrial] = (itrial - 1) % my numberOfDifferentStimuli + 1;
} else if (my randomize == kExperiment_randomize_PERMUTE_ALL) {
for (long itrial = 1; itrial <= my numberOfTrials; itrial ++)
my stimuli [itrial] = (itrial - 1) % my numberOfDifferentStimuli + 1;
permuteRandomly (me, 1, my numberOfTrials);
} else if (my randomize == kExperiment_randomize_PERMUTE_BALANCED) {
for (long ireplica = 1; ireplica <= my numberOfReplicationsPerStimulus; ireplica ++) {
long offset = (ireplica - 1) * my numberOfDifferentStimuli;
for (long istim = 1; istim <= my numberOfDifferentStimuli; istim ++)
my stimuli [offset + istim] = istim;
permuteRandomly (me, offset + 1, offset + my numberOfDifferentStimuli);
}
} else if (my randomize == kExperiment_randomize_PERMUTE_BALANCED_NO_DOUBLETS) {
for (long ireplica = 1; ireplica <= my numberOfReplicationsPerStimulus; ireplica ++) {
long offset = (ireplica - 1) * my numberOfDifferentStimuli;
for (long istim = 1; istim <= my numberOfDifferentStimuli; istim ++)
my stimuli [offset + istim] = istim;
do {
permuteRandomly (me, offset + 1, offset + my numberOfDifferentStimuli);
} while (ireplica != 1 && my stimuli [offset + 1] == my stimuli [offset] && my numberOfDifferentStimuli > 1);
}
} else if (my randomize == kExperiment_randomize_WITH_REPLACEMENT) {
for (long itrial = 1; itrial <= my numberOfTrials; itrial ++)
my stimuli [itrial] = NUMrandomInteger (1, my numberOfDifferentStimuli);
}
Melder_warningOn ();
} catch (MelderError) {
Melder_warningOn ();
my numberOfTrials = 0;
NUMvector_free (my stimuli, 1); my stimuli = nullptr;
Melder_throw (me, U": not started.");
}
}
void structCategoriesEditorCommand :: v_destroy () {
NUMvector_free (selection, 1);
CategoriesEditorCommand_Parent :: v_destroy ();
}
int Praat_tests (int itest, char32 *arg1, char32 *arg2, char32 *arg3, char32 *arg4) {
int64 n = Melder_atoi (arg1);
double t = 0.0;
(void) arg1;
(void) arg2;
(void) arg3;
(void) arg4;
Melder_clearInfo ();
Melder_stopwatch ();
switch (itest) {
case kPraatTests_TIME_RANDOM_FRACTION: {
for (int64 i = 1; i <= n; i ++)
(void) NUMrandomFraction ();
t = Melder_stopwatch ();
} break;
case kPraatTests_TIME_RANDOM_GAUSS: {
for (int64 i = 1; i <= n; i ++)
(void) NUMrandomGauss (0.0, 1.0);
t = Melder_stopwatch ();
} break;
case kPraatTests_TIME_SORT: {
long m = Melder_atoi (arg2);
long *array = NUMvector <long> (1, m);
for (int64 i = 1; i <= m; i ++)
array [i] = NUMrandomInteger (1, 100);
Melder_stopwatch ();
for (int64 i = 1; i <= n; i ++)
NUMsort_l (m, array);
t = Melder_stopwatch ();
NUMvector_free (array, 1);
} break;
case kPraatTests_TIME_INTEGER: {
int64 sum = 0;
for (int64 i = 1; i <= n; i ++)
sum += i * (i - 1) * (i - 2);
t = Melder_stopwatch ();
MelderInfo_writeLine (sum);
} break;
case kPraatTests_TIME_FLOAT: {
double sum = 0.0, fn = n;
for (double fi = 1.0; fi <= fn; fi = fi + 1.0)
sum += fi * (fi - 1.0) * (fi - 2.0);
t = Melder_stopwatch ();
MelderInfo_writeLine (sum);
} break;
case kPraatTests_TIME_FLOAT_TO_UNSIGNED_BUILTIN: {
uint64_t sum = 0;
double fn = n;
for (double fi = 1.0; fi <= fn; fi = fi + 1.0)
sum += (uint32) fi;
t = Melder_stopwatch (); // 2.59 // 1.60
MelderInfo_writeLine (sum);
} break;
case kPraatTests_TIME_FLOAT_TO_UNSIGNED_EXTERN: {
uint64_t sum = 0;
double fn = n;
for (double fi = 1.0; fi <= fn; fi = fi + 1.0)
sum += (uint32) ((int32) (fi - 2147483648.0) + 2147483647L + 1);
t = Melder_stopwatch (); // 1.60
MelderInfo_writeLine (sum);
} break;
case kPraatTests_TIME_UNSIGNED_TO_FLOAT_BUILTIN: {
double sum = 0.0;
uint32 nu = (uint32) n;
for (uint32 iu = 1; iu <= nu; iu ++)
sum += (double) iu;
t = Melder_stopwatch (); // 1.35
MelderInfo_writeLine (sum);
} break;
case kPraatTests_TIME_UNSIGNED_TO_FLOAT_EXTERN: {
double sum = 0.0;
uint32 nu = (uint32) n;
for (uint32 iu = 1; iu <= nu; iu ++)
sum += (double) (int32) (iu - 2147483647L - 1) + 2147483648.0;
t = Melder_stopwatch (); // 0.96
MelderInfo_writeLine (sum);
} break;
case kPraatTests_TIME_STRING_MELDER_32: {
autoMelderString string;
char32 word [] { U"abc" };
word [2] = NUMrandomInteger ('a', 'z');
for (int64 i = 1; i <= n; i ++) {
MelderString_copy (& string, word);
for (int j = 1; j <= 30; j ++)
MelderString_append (& string, word);
}
t = Melder_stopwatch ();
} break;
case kPraatTests_TIME_STRING_MELDER_32_ALLOC: {
char32 word [] { U"abc" };
word [2] = NUMrandomInteger ('a', 'z');
for (int64 i = 1; i <= n; i ++) {
autoMelderString string;
MelderString_copy (& string, word);
for (int j = 1; j <= 30; j ++)
MelderString_append (& string, word);
}
t = Melder_stopwatch ();
} break;
case kPraatTests_TIME_STRING_CPP_S: {
std::string s = "";
char word [] { "abc" };
word [2] = (char) NUMrandomInteger ('a', 'z');
for (int64 i = 1; i <= n; i ++) {
s = word;
for (int j = 1; j <= 30; j ++)
s += word;
}
t = Melder_stopwatch ();
} break;
case kPraatTests_TIME_STRING_CPP_C: {
std::basic_string<char> s = "";
char word [] { "abc" };
word [2] = (char) NUMrandomInteger ('a', 'z');
for (int64 i = 1; i <= n; i ++) {
s = word;
for (int j = 1; j <= 30; j ++)
s += word;
}
t = Melder_stopwatch ();
} break;
case kPraatTests_TIME_STRING_CPP_WS: {
std::wstring s = L"";
wchar_t word [] { L"abc" };
word [2] = NUMrandomInteger ('a', 'z');
for (int64 i = 1; i <= n; i ++) {
s = word;
for (int j = 1; j <= 30; j ++)
s += word;
}
t = Melder_stopwatch ();
} break;
case kPraatTests_TIME_STRING_CPP_WC: {
std::basic_string<wchar_t> s = L"";
wchar_t word [] { L"abc" };
word [2] = NUMrandomInteger ('a', 'z');
for (int64 i = 1; i <= n; i ++) {
s = word;
for (int j = 1; j <= 30; j ++)
s += word;
}
t = Melder_stopwatch ();
} break;
case kPraatTests_TIME_STRING_CPP_32: {
std::basic_string<char32_t> s = U"";
char32 word [] { U"abc" };
word [2] = NUMrandomInteger ('a', 'z');
for (int64 i = 1; i <= n; i ++) {
s = word;
for (int j = 1; j <= 30; j ++)
s += word;
}
t = Melder_stopwatch ();
} break;
case kPraatTests_TIME_STRING_CPP_U32STRING: {
std::u32string s = U"";
char32 word [] { U"abc" };
word [2] = NUMrandomInteger ('a', 'z');
for (int64 i = 1; i <= n; i ++) {
s = word;
for (int j = 1; j <= 30; j ++)
s += word;
}
t = Melder_stopwatch ();
} break;
case kPraatTests_TIME_STRCPY: {
char buffer [100];
char word [] { "abc" };
word [2] = (char) NUMrandomInteger ('a', 'z');
for (int64 i = 1; i <= n; i ++) {
strcpy (buffer, word);
for (int j = 1; j <= 30; j ++)
strcpy (buffer + strlen (buffer), word);
}
t = Melder_stopwatch ();
MelderInfo_writeLine (Melder_peek8to32 (buffer));
} break;
case kPraatTests_TIME_WCSCPY: {
wchar_t buffer [100];
wchar_t word [] { L"abc" };
word [2] = NUMrandomInteger ('a', 'z');
for (int64 i = 1; i <= n; i ++) {
wcscpy (buffer, word);
for (int j = 1; j <= 30; j ++)
wcscpy (buffer + wcslen (buffer), word);
}
t = Melder_stopwatch ();
} break;
case kPraatTests_TIME_STR32CPY: {
char32 buffer [100];
char32 word [] { U"abc" };
word [2] = NUMrandomInteger ('a', 'z');
for (int64 i = 1; i <= n; i ++) {
str32cpy (buffer, word);
for (int j = 1; j <= 30; j ++)
str32cpy (buffer + str32len (buffer), word);
}
t = Melder_stopwatch ();
MelderInfo_writeLine (buffer);
} break;
case kPraatTests_TIME_GRAPHICS_TEXT_TOP: {
autoPraatPicture picture;
for (int64 i = 1; i <= n; i ++) {
Graphics_textTop (GRAPHICS, false, U"hello world");
}
t = Melder_stopwatch ();
} break;
case kPraatTests_THING_AUTO: {
int numberOfThingsBefore = theTotalNumberOfThings;
{
Melder_casual (U"1\n");
autoDaata data = Thing_new (Daata);
Thing_setName (data.get(), U"hello");
Melder_casual (U"2\n");
testData (data.get());
testAutoData (data.move());
autoDaata data18 = Thing_new (Daata);
testAutoData (data18.move());
fprintf (stderr, "3\n");
autoDaata data2 = newAutoData ();
fprintf (stderr, "4\n");
autoDaata data3 = newAutoData ();
fprintf (stderr, "5\n");
//data2 = data; // disabled l-value copy assignment from same class
fprintf (stderr, "6\n");
autoOrdered ordered = Thing_new (Ordered);
fprintf (stderr, "7\n");
//data = ordered; // disabled l-value copy assignment from subclass
data = ordered.move();
//ordered = data; // disabled l-value copy assignment from superclass
//ordered = data.move(); // assignment from superclass to subclass is rightfully refused by compiler
fprintf (stderr, "8\n");
data2 = newAutoData ();
fprintf (stderr, "8a\n");
autoDaata data5 = newAutoData ();
fprintf (stderr, "8b\n");
data2 = data5.move();
fprintf (stderr, "9\n");
//ordered = data; // rightfully refused by compiler
fprintf (stderr, "10\n");
//autoOrdered ordered2 = Thing_new (Daata); // rightfully refused by compiler
fprintf (stderr, "11\n");
autoDaata data4 = Thing_new (Ordered); // constructor
fprintf (stderr, "12\n");
//autoDaata data6 = data4; // disabled l-value copy constructor from same class
fprintf (stderr, "13\n");
autoDaata data7 = data4.move();
fprintf (stderr, "14\n");
autoOrdered ordered3 = Thing_new (Ordered);
autoDaata data8 = ordered3.move();
fprintf (stderr, "15\n");
//autoDaata data9 = ordered; // disabled l-value copy constructor from subclass
fprintf (stderr, "16\n");
autoDaata data10 = data7.move();
fprintf (stderr, "17\n");
autoDaata data11 = Thing_new (Daata); // constructor, move assignment, null destructor
fprintf (stderr, "18\n");
data11 = Thing_new (Ordered);
fprintf (stderr, "19\n");
testAutoDataRef (data11);
fprintf (stderr, "20\n");
//data11 = nullptr; // disabled implicit assignment of pointer to autopointer
fprintf (stderr, "21\n");
}
int numberOfThingsAfter = theTotalNumberOfThings;
fprintf (stderr, "Number of things: before %d, after %d\n", numberOfThingsBefore, numberOfThingsAfter);
#if 1
MelderCallback<void,structDaata>::FunctionType f;
typedef void (*DataFunc) (Daata);
typedef void (*OrderedFunc) (Ordered);
DataFunc dataFun;
OrderedFunc orderedFun;
MelderCallback<void,structDaata> dataFun2 (dataFun);
MelderCallback<void,structOrdered> orderedFun2 (orderedFun);
MelderCallback<void,structDaata> dataFun3 (orderedFun);
//MelderCallback<void,structOrdered> orderedFun3 (dataFun); // rightfully refused by compiler
autoDaata data = Thing_new (Daata);
dataFun3 (data.get());
#endif
} break;
}
MelderInfo_writeLine (Melder_single (t / n * 1e9), U" nanoseconds");
MelderInfo_close ();
return 1;
}
void structLineMinimizer :: v_destroy () {
NUMvector_free (ptry, 1);
NUMvector_free (direction, 1);
LineMinimizer_Parent :: v_destroy ();
}
void structMinimizer :: v_destroy () {
NUMvector_free (p, 1);
NUMvector_free (history, 1);
Minimizer_Parent :: v_destroy ();
}
void structDelta :: v_destroy () {
NUMvector_free (this -> tube, 1);
Delta_Parent :: v_destroy ();
}
void structPitch :: v_info () {
long nVoiced;
double *frequencies = Sampled_getSortedValues (this, Pitch_LEVEL_FREQUENCY, kPitch_unit_HERTZ, & nVoiced);
structData :: v_info ();
MelderInfo_writeLine (L"Time domain:");
MelderInfo_writeLine (L" Start time: ", Melder_double (xmin), L" seconds");
MelderInfo_writeLine (L" End time: ", Melder_double (xmax), L" seconds");
MelderInfo_writeLine (L" Total duration: ", Melder_double (xmax - xmin), L" seconds");
MelderInfo_writeLine (L"Time sampling:");
MelderInfo_writeLine (L" Number of frames: ", Melder_integer (nx), L" (", Melder_integer (nVoiced), L" voiced)");
MelderInfo_writeLine (L" Time step: ", Melder_double (dx), L" seconds");
MelderInfo_writeLine (L" First frame centred at: ", Melder_double (x1), L" seconds");
MelderInfo_writeLine (L"Ceiling at: ", Melder_double (ceiling), L" Hz");
if (nVoiced >= 1) { /* Quantiles. */
double quantile10, quantile16, quantile50, quantile84, quantile90;
quantile10 = NUMquantile (nVoiced, frequencies, 0.10);
quantile16 = NUMquantile (nVoiced, frequencies, 0.16);
quantile50 = NUMquantile (nVoiced, frequencies, 0.50); /* Median. */
quantile84 = NUMquantile (nVoiced, frequencies, 0.84);
quantile90 = NUMquantile (nVoiced, frequencies, 0.90);
MelderInfo_writeLine (L"\nEstimated quantiles:");
MelderInfo_write (L" 10% = ", Melder_single (quantile10), L" Hz = ", Melder_single (MEL (quantile10)), L" Mel = ");
MelderInfo_writeLine (Melder_single (SEMITONES (quantile10)), L" semitones above 100 Hz = ", Melder_single (ERB (quantile10)), L" ERB");
MelderInfo_write (L" 16% = ", Melder_single (quantile16), L" Hz = ", Melder_single (MEL (quantile16)), L" Mel = ");
MelderInfo_writeLine (Melder_single (SEMITONES (quantile16)), L" semitones above 100 Hz = ", Melder_single (ERB (quantile16)), L" ERB");
MelderInfo_write (L" 50% = ", Melder_single (quantile50), L" Hz = ", Melder_single (MEL (quantile50)), L" Mel = ");
MelderInfo_writeLine (Melder_single (SEMITONES (quantile50)), L" semitones above 100 Hz = ", Melder_single (ERB (quantile50)), L" ERB");
MelderInfo_write (L" 84% = ", Melder_single (quantile84), L" Hz = ", Melder_single (MEL (quantile84)), L" Mel = ");
MelderInfo_writeLine (Melder_single (SEMITONES (quantile84)), L" semitones above 100 Hz = ", Melder_single (ERB (quantile84)), L" ERB");
MelderInfo_write (L" 90% = ", Melder_single (quantile90), L" Hz = ", Melder_single (MEL (quantile90)), L" Mel = ");
MelderInfo_writeLine (Melder_single (SEMITONES (quantile90)), L" semitones above 100 Hz = ", Melder_single (ERB (quantile90)), L" ERB");
if (nVoiced > 1) {
double corr = sqrt (nVoiced / (nVoiced - 1.0));
MelderInfo_writeLine (L"\nEstimated spreading:");
MelderInfo_write (L" 84%-median = ", Melder_half ((quantile84 - quantile50) * corr), L" Hz = ", Melder_half ((MEL (quantile84) - MEL (quantile50)) * corr), L" Mel = ");
MelderInfo_writeLine (Melder_half ((SEMITONES (quantile84) - SEMITONES (quantile50)) * corr), L" semitones = ", Melder_half ((ERB (quantile84) - ERB (quantile50)) * corr), L" ERB");
MelderInfo_write (L" median-16% = ", Melder_half ((quantile50 - quantile16) * corr), L" Hz = ", Melder_half ((MEL (quantile50) - MEL (quantile16)) * corr), L" Mel = ");
MelderInfo_writeLine (Melder_half ((SEMITONES (quantile50) - SEMITONES (quantile16)) * corr), L" semitones = ", Melder_half ((ERB (quantile50) - ERB (quantile16)) * corr), L" ERB");
MelderInfo_write (L" 90%-10% = ", Melder_half ((quantile90 - quantile10) * corr), L" Hz = ", Melder_half ((MEL (quantile90) - MEL (quantile10)) * corr), L" Mel = ");
MelderInfo_writeLine (Melder_half ((SEMITONES (quantile90) - SEMITONES (quantile10)) * corr), L" semitones = ", Melder_half ((ERB (quantile90) - ERB (quantile10)) * corr), L" 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 (L"\nMinimum ", Melder_single (minimum), L" Hz = ", Melder_single (MEL (minimum)), L" Mel = ");
MelderInfo_writeLine (Melder_single (SEMITONES (minimum)), L" semitones above 100 Hz = ", Melder_single (ERB (minimum)), L" ERB");
MelderInfo_write (L"Maximum ", Melder_single (maximum), L" Hz = ", Melder_single (MEL (maximum)), L" Mel = ");
MelderInfo_writeLine (Melder_single (SEMITONES (maximum)), L" semitones above 100 Hz = ", Melder_single (ERB (maximum)), L" ERB");
MelderInfo_write (L"Range ", Melder_half (maximum - minimum), L" Hz = ", Melder_single (MEL (maximum) - MEL (minimum)), L" Mel = ");
MelderInfo_writeLine (Melder_half (SEMITONES (maximum) - SEMITONES (minimum)), L" semitones = ", Melder_half (ERB (maximum) - ERB (minimum)), L" 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 (L"Average: ", Melder_single (meanHertz), L" Hz = ", Melder_single (meanMel), L" Mel = ");
MelderInfo_writeLine (Melder_single (meanSemitones), L" semitones above 100 Hz = ", Melder_single (meanErb), L" 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 (L"Standard deviation: ", Melder_half (stdevHertz), L" Hz = ", Melder_half (stdevMel), L" Mel = ");
MelderInfo_writeLine (Melder_half (stdevSemitones), L" semitones = ", Melder_half (stdevErb), L" ERB");
}
}
NUMvector_free (frequencies, 1);
if (nVoiced > 1) { /* Variability: mean absolute slope. */
double slopeHertz, slopeMel, slopeSemitones, slopeErb, slopeWithoutOctaveJumps;
Pitch_getMeanAbsoluteSlope (this, & slopeHertz, & slopeMel, & slopeSemitones, & slopeErb, & slopeWithoutOctaveJumps);
MelderInfo_write (L"\nMean absolute slope: ", Melder_half (slopeHertz), L" Hz/s = ", Melder_half (slopeMel), L" Mel/s = ");
MelderInfo_writeLine (Melder_half (slopeSemitones), L" semitones/s = ", Melder_half (slopeErb), L" ERB/s");
MelderInfo_writeLine (L"Mean absolute slope without octave jumps: ", Melder_half (slopeWithoutOctaveJumps), L" semitones/s");
}
}
int Praat_tests (int itest, char32 *arg1, char32 *arg2, char32 *arg3, char32 *arg4) {
int64 n = Melder_atoi (arg1);
double t;
(void) arg1;
(void) arg2;
(void) arg3;
(void) arg4;
Melder_clearInfo ();
Melder_stopwatch ();
switch (itest) {
case kPraatTests_TIME_RANDOM_FRACTION: {
for (int64 i = 1; i <= n; i ++)
(void) NUMrandomFraction ();
t = Melder_stopwatch ();
} break;
case kPraatTests_TIME_RANDOM_GAUSS: {
for (int64 i = 1; i <= n; i ++)
(void) NUMrandomGauss (0.0, 1.0);
t = Melder_stopwatch ();
} break;
case kPraatTests_TIME_SORT: {
long m = Melder_atoi (arg2);
long *array = NUMvector <long> (1, m);
for (int64 i = 1; i <= m; i ++)
array [i] = NUMrandomInteger (1, 100);
Melder_stopwatch ();
for (int64 i = 1; i <= n; i ++)
NUMsort_l (m, array);
t = Melder_stopwatch ();
NUMvector_free (array, 1);
} break;
case kPraatTests_TIME_INTEGER: {
double sum = 0;
for (int64 i = 1; i <= n; i ++)
sum += i * (i - 1) * (i - 2);
t = Melder_stopwatch ();
MelderInfo_writeLine (sum);
} break;
case kPraatTests_TIME_FLOAT: {
double sum = 0.0, fn = n;
for (double fi = 1.0; fi <= fn; fi = fi + 1.0)
sum += fi * (fi - 1.0) * (fi - 2.0);
t = Melder_stopwatch ();
MelderInfo_writeLine (sum);
} break;
case kPraatTests_TIME_FLOAT_TO_UNSIGNED_BUILTIN: {
uint64_t sum = 0;
double fn = n;
for (double fi = 1.0; fi <= fn; fi = fi + 1.0)
sum += (uint32_t) fi;
t = Melder_stopwatch (); // 2.59 // 1.60
MelderInfo_writeLine (sum);
} break;
case kPraatTests_TIME_FLOAT_TO_UNSIGNED_EXTERN: {
uint64_t sum = 0;
double fn = n;
for (double fi = 1.0; fi <= fn; fi = fi + 1.0)
sum += (uint32_t) ((int32_t) (fi - 2147483648.0) + 2147483647L + 1);
t = Melder_stopwatch (); // 1.60
MelderInfo_writeLine (sum);
} break;
case kPraatTests_TIME_UNSIGNED_TO_FLOAT_BUILTIN: {
double sum = 0.0;
uint32_t nu = (uint32_t) n;
for (uint32_t iu = 1; iu <= nu; iu ++)
sum += (double) iu;
t = Melder_stopwatch (); // 1.35
MelderInfo_writeLine (sum);
} break;
case kPraatTests_TIME_UNSIGNED_TO_FLOAT_EXTERN: {
double sum = 0.0;
uint32_t nu = (uint32_t) n;
for (uint32_t iu = 1; iu <= nu; iu ++)
sum += (double) (int32_t) (iu - 2147483647L - 1) + 2147483648.0;
t = Melder_stopwatch (); // 0.96
MelderInfo_writeLine (sum);
} break;
case kPraatTests_TIME_STRING_MELDER_32: {
autoMelderString string;
char32 word [] { U"abc" };
word [2] = NUMrandomInteger ('a', 'z');
for (int64 i = 1; i <= n; i ++) {
MelderString_copy (& string, word);
for (int j = 1; j <= 30; j ++)
MelderString_append (& string, word);
}
t = Melder_stopwatch ();
} break;
case kPraatTests_TIME_STRING_CPP_S: {
std::string s = "";
char word [] { "abc" };
word [2] = (char) NUMrandomInteger ('a', 'z');
for (int64 i = 1; i <= n; i ++) {
s = word;
for (int j = 1; j <= 30; j ++)
s += word;
}
t = Melder_stopwatch ();
} break;
case kPraatTests_TIME_STRING_CPP_C: {
std::basic_string<char> s = "";
char word [] { "abc" };
word [2] = (char) NUMrandomInteger ('a', 'z');
for (int64 i = 1; i <= n; i ++) {
s = word;
for (int j = 1; j <= 30; j ++)
s += word;
}
t = Melder_stopwatch ();
} break;
case kPraatTests_TIME_STRING_CPP_WS: {
std::wstring s = L"";
wchar_t word [] { L"abc" };
word [2] = NUMrandomInteger ('a', 'z');
for (int64 i = 1; i <= n; i ++) {
s = word;
for (int j = 1; j <= 30; j ++)
s += word;
}
t = Melder_stopwatch ();
} break;
case kPraatTests_TIME_STRING_CPP_WC: {
std::basic_string<wchar_t> s = L"";
wchar_t word [] { L"abc" };
word [2] = NUMrandomInteger ('a', 'z');
for (int64 i = 1; i <= n; i ++) {
s = word;
for (int j = 1; j <= 30; j ++)
s += word;
}
t = Melder_stopwatch ();
} break;
case kPraatTests_TIME_STRING_CPP_32: {
std::basic_string<char32_t> s = U"";
char32 word [] { U"abc" };
word [2] = NUMrandomInteger ('a', 'z');
for (int64 i = 1; i <= n; i ++) {
s = word;
for (int j = 1; j <= 30; j ++)
s += word;
}
t = Melder_stopwatch ();
} break;
case kPraatTests_TIME_STRING_CPP_U32STRING: {
#if ! defined (macintosh) || ! useCarbon
std::u32string s = U"";
char32 word [] { U"abc" };
word [2] = NUMrandomInteger ('a', 'z');
for (int64 i = 1; i <= n; i ++) {
s = word;
for (int j = 1; j <= 30; j ++)
s += word;
}
#endif
t = Melder_stopwatch ();
} break;
case kPraatTests_TIME_STRCPY: {
char buffer [100];
char word [] { "abc" };
word [2] = (char) NUMrandomInteger ('a', 'z');
for (int64 i = 1; i <= n; i ++) {
strcpy (buffer, word);
for (int j = 1; j <= 30; j ++)
strcpy (buffer + strlen (buffer), word);
}
t = Melder_stopwatch ();
MelderInfo_writeLine (Melder_peek8to32 (buffer));
} break;
case kPraatTests_TIME_WCSCPY: {
wchar_t buffer [100];
wchar_t word [] { L"abc" };
word [2] = NUMrandomInteger ('a', 'z');
for (int64 i = 1; i <= n; i ++) {
wcscpy (buffer, word);
for (int j = 1; j <= 30; j ++)
wcscpy (buffer + wcslen (buffer), word);
}
t = Melder_stopwatch ();
} break;
case kPraatTests_TIME_STR32CPY: {
char32 buffer [100];
char32 word [] { U"abc" };
word [2] = NUMrandomInteger ('a', 'z');
for (int64 i = 1; i <= n; i ++) {
str32cpy (buffer, word);
for (int j = 1; j <= 30; j ++)
str32cpy (buffer + str32len (buffer), word);
}
t = Melder_stopwatch ();
MelderInfo_writeLine (buffer);
} break;
case kPraatTests_TIME_GRAPHICS_TEXT_TOP: {
autoPraatPicture picture;
for (int64 i = 1; i <= n; i ++) {
Graphics_textTop (GRAPHICS, false, U"hello world");
}
t = Melder_stopwatch ();
} break;
}
MelderInfo_writeLine (Melder_single (t / n * 1e9), U" nanoseconds");
MelderInfo_close ();
return 1;
}
void Melder_freeTokens (char32 ***tokens) {
NUMvector_free (*tokens, 1);
*tokens = nullptr;
}
