static void classMinimizer_afterHook (Minimizer me, Thing /* boss */) {
if (my success || ! my gmonitor) {
return;
}
if (my start == 1) {
Minimizer_drawHistory (me, my gmonitor, 0, my maxNumOfIterations, 0.0, 1.1 * my history[1], 1);
Graphics_textTop (my gmonitor, false, Melder_cat (U"Dimension of search space: ", my nParameters));
}
Graphics_setInner (my gmonitor);
Graphics_line (my gmonitor, my iteration, my history[my iteration], my iteration, my history[my iteration]);
Graphics_unsetInner (my gmonitor);
Melder_monitor ((double) (my iteration) / my maxNumOfIterations, U"Iterations: ", my iteration,
U", Function calls: ", my funcCalls, U", Cost: ", my minimum);
}
void ERP_drawChannel_number (ERP me, Graphics graphics, long channelNumber, double tmin, double tmax, double vmin, double vmax, bool garnish) {
if (channelNumber < 1 || channelNumber > my ny) return;
/*
* Automatic domain.
*/
if (tmin == tmax) {
tmin = my xmin;
tmax = my xmax;
}
/*
* Domain expressed in sample numbers.
*/
long ixmin, ixmax;
Matrix_getWindowSamplesX (me, tmin, tmax, & ixmin, & ixmax);
/*
* Automatic vertical range.
*/
if (vmin == vmax) {
Matrix_getWindowExtrema (me, ixmin, ixmax, channelNumber, channelNumber, & vmin, & vmax);
if (vmin == vmax) {
vmin -= 1.0;
vmax += 1.0;
}
}
/*
* Set coordinates for drawing.
*/
Graphics_setInner (graphics);
Graphics_setWindow (graphics, tmin, tmax, vmin, vmax);
Graphics_function (graphics, my z [channelNumber], ixmin, ixmax, Matrix_columnToX (me, ixmin), Matrix_columnToX (me, ixmax));
Graphics_unsetInner (graphics);
if (garnish) {
Graphics_drawInnerBox (graphics);
Graphics_textTop (graphics, true, Melder_wcscat (L"Channel ", my channelNames [channelNumber]));
Graphics_textBottom (graphics, true, L"Time (s)");
Graphics_marksBottom (graphics, 2, true, true, false);
if (0.0 > tmin && 0.0 < tmax)
Graphics_markBottom (graphics, 0.0, true, true, true, NULL);
Graphics_markLeft (graphics, vmin, true, true, false, NULL);
Graphics_markLeft (graphics, vmax, true, true, false, NULL);
Graphics_markBottom (graphics, 0.0, true, true, true, NULL);
if (vmin != 0.0 && vmax != 0.0 && (vmin > 0.0) != (vmax > 0.0)) {
Graphics_markLeft (graphics, 0.0, true, true, true, NULL);
}
}
}
void SPINET_spectralRepresentation (SPINET me, Graphics g, double fromTime, double toTime, double fromErb, double toErb, double minimum, double maximum, int enhanced, int garnish) {
double **z = enhanced ? my s : my y;
autoMatrix thee = Matrix_create (my xmin, my xmax, my nx, my dx, my x1, my ymin, my ymax, my ny, my dy, my y1);
for (long j = 1; j <= my ny; j++) {
for (long i = 1; i <= my nx; i++) {
thy z[j][i] = z[j][i];
}
}
Matrix_paintCells (thee.peek(), g, fromTime, toTime, fromErb, toErb, minimum, maximum);
if (garnish) {
Graphics_drawInnerBox (g);
Graphics_textBottom (g, true, U"Time (s)");
Graphics_marksBottom (g, 2, true, true, false);
Graphics_textLeft (g, true, U"Frequency (ERB)");
Graphics_marksLeft (g, 2, true, true, false);
Graphics_textTop (g, false, enhanced ? U"Cooperative interaction output" : U"Gammatone filterbank output");
}
}
void Artword_draw (Artword me, Graphics g, int feature, bool garnish) {
long numberOfTargets = my data [feature]. numberOfTargets;
if (numberOfTargets > 0) {
autoNUMvector <double> x (1, numberOfTargets);
autoNUMvector <double> y (1, numberOfTargets);
Graphics_setInner (g);
Graphics_setWindow (g, 0, my totalTime, -1.0, 1.0);
for (int i = 1; i <= numberOfTargets; i ++) {
x [i] = my data [feature]. times [i];
y [i] = my data [feature]. targets [i];
}
Graphics_polyline (g, numberOfTargets, & x [1], & y [1]);
Graphics_unsetInner (g);
}
if (garnish) {
Graphics_drawInnerBox (g);
Graphics_marksBottom (g, 2, true, true, false);
Graphics_marksLeft (g, 3, true, true, true);
Graphics_textTop (g, false, kArt_muscle_getText (feature));
Graphics_textBottom (g, true, U"Time (s)");
}
}
static void classMinimizer_after (I, Any aclosure) {
iam (Minimizer);
(void) aclosure;
if (my success || ! my gmonitor) {
return;
}
if (my start == 1) {
wchar_t s[35];
Minimizer_drawHistory (me, my gmonitor, 0, my maxNumOfIterations, 0, 1.1 * my history[1], 1);
swprintf (s, 35, L"Dimension of search space: %6ld", my nParameters);
Graphics_textTop (my gmonitor, 0, s);
}
Graphics_setInner (my gmonitor);
Graphics_line (my gmonitor, my iteration, my history[my iteration],
my iteration, my history[my iteration]);
Graphics_unsetInner (my gmonitor);
Melder_monitor ( (double) (my iteration) / my maxNumOfIterations,
L"Iterations: ", Melder_integer (my iteration),
L", Function calls: ", Melder_integer (my funcCalls),
L", Cost: ", Melder_double (my minimum));
}
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;
}
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;
}