void TableOfReal_writeToHeaderlessSpreadsheetFile (TableOfReal me, MelderFile file) {
try {
autoMelderString buffer;
MelderString_copy (& buffer, L"rowLabel");
for (long icol = 1; icol <= my numberOfColumns; icol ++) {
MelderString_appendCharacter (& buffer, '\t');
wchar_t *s = my columnLabels [icol];
MelderString_append (& buffer, s != NULL && s [0] != '\0' ? s : L"?");
}
MelderString_appendCharacter (& buffer, '\n');
for (long irow = 1; irow <= my numberOfRows; irow ++) {
wchar_t *s = my rowLabels [irow];
MelderString_append (& buffer, s != NULL && s [0] != '\0' ? s : L"?");
for (long icol = 1; icol <= my numberOfColumns; icol ++) {
MelderString_appendCharacter (& buffer, '\t');
double x = my data [irow] [icol];
MelderString_append (& buffer, Melder_double (x));
}
MelderString_appendCharacter (& buffer, '\n');
}
MelderFile_writeText (file, buffer.string, Melder_getOutputEncoding ());
} catch (MelderError) {
Melder_throw (me, ": not saved to tab-separated file.");
}
}
const char32 * Melder_pad (const char32 *string, int64 width) {
if (++ iPadBuffer == NUMBER_OF_BUFFERS) iPadBuffer = 0;
int64 length = str32len (string);
int64 tooShort = width - length;
if (tooShort <= 0) return string;
MelderString_copy (& thePadBuffers [iPadBuffer], string);
for (int64 i = 0; i < tooShort; i ++)
MelderString_appendCharacter (& thePadBuffers [iPadBuffer], U' ');
return thePadBuffers [iPadBuffer]. string;
}
autostring32 FFNet_createNameFromTopology (FFNet me) {
autoMelderString name;
MelderString_copy (&name, my nUnitsInLayer[0]);
for (long i = 1; i <= my nLayers; i++) {
MelderString_appendCharacter (&name, U'-');
MelderString_append (&name, my nUnitsInLayer[i]);
}
autostring32 naam = Melder_dup (name.string);
return naam.transfer();
}
const char32 * Melder_padOrTruncate (const char32 *string, int64 width) {
if (++ iPadBuffer == NUMBER_OF_BUFFERS) iPadBuffer = 0;
int64 length = str32len (string);
int64 tooLong = length - width;
if (tooLong == 0) return string;
if (tooLong < 0) {
int64 tooShort = - tooLong;
MelderString_copy (& thePadBuffers [iPadBuffer], string);
for (int64 i = 0; i < tooShort; i ++)
MelderString_appendCharacter (& thePadBuffers [iPadBuffer], U' ');
} else {
MelderString_ncopy (& thePadBuffers [iPadBuffer], string, width);
}
return thePadBuffers [iPadBuffer]. string;
}
static Strings Strings_createAsFileOrDirectoryList (const wchar_t *path, int type) {
#if USE_STAT
/*
* Initialize.
*/
DIR *d = NULL;
try {
autoMelderString filePath, searchDirectory, left, right;
/*
* Parse the path.
* Example: in /Users/paul/sounds/h*.wav",
* the search directory is "/Users/paul/sounds",
* the left environment is "h", and the right environment is ".wav".
*/
MelderString_copy (& searchDirectory, path);
wchar_t *asterisk = wcsrchr (searchDirectory. string, '*');
if (asterisk != NULL) {
*asterisk = '\0';
searchDirectory. length = asterisk - searchDirectory. string; // probably superfluous, but correct
wchar_t *lastSlash = wcsrchr (searchDirectory. string, Melder_DIRECTORY_SEPARATOR);
if (lastSlash != NULL) {
*lastSlash = '\0'; // This fixes searchDirectory.
searchDirectory. length = lastSlash - searchDirectory. string; // probably superfluous, but correct
MelderString_copy (& left, lastSlash + 1);
} else {
MelderString_copy (& left, searchDirectory. string); // quickly save...
MelderString_empty (& searchDirectory); // ...before destruction
}
MelderString_copy (& right, asterisk + 1);
}
char buffer8 [1+kMelder_MAXPATH];
Melder_wcsTo8bitFileRepresentation_inline (searchDirectory. string, buffer8);
d = opendir (buffer8 [0] ? buffer8 : ".");
if (d == NULL)
Melder_throw ("Cannot open directory ", searchDirectory. string, ".");
//Melder_casual ("opened");
autoStrings me = Thing_new (Strings);
my strings = NUMvector <wchar_t *> (1, 1000000);
struct dirent *entry;
while ((entry = readdir (d)) != NULL) {
MelderString_copy (& filePath, searchDirectory. string [0] ? searchDirectory. string : L".");
MelderString_appendCharacter (& filePath, Melder_DIRECTORY_SEPARATOR);
wchar_t bufferW [1+kMelder_MAXPATH];
Melder_8bitFileRepresentationToWcs_inline (entry -> d_name, bufferW);
MelderString_append (& filePath, bufferW);
//Melder_casual ("read %s", filePath. string);
Melder_wcsTo8bitFileRepresentation_inline (filePath. string, buffer8);
struct stat stats;
if (stat (buffer8, & stats) != 0) {
Melder_throw ("Cannot look at file ", filePath. string, ".");
//stats. st_mode = -1L;
}
//Melder_casual ("statted %s", filePath. string);
//Melder_casual ("file %s mode %s", filePath. string, Melder_integer (stats. st_mode / 4096));
if ((type == Strings_createAsFileOrDirectoryList_TYPE_FILE && S_ISREG (stats. st_mode)) ||
(type == Strings_createAsFileOrDirectoryList_TYPE_DIRECTORY && S_ISDIR (stats. st_mode)))
{
Melder_8bitFileRepresentationToWcs_inline (entry -> d_name, bufferW);
unsigned long length = wcslen (bufferW);
if (bufferW [0] != '.' &&
(left. length == 0 || wcsnequ (bufferW, left. string, left. length)) &&
(right. length == 0 || (length >= right. length && wcsequ (bufferW + (length - right. length), right. string))))
{
my strings [++ my numberOfStrings] = Melder_wcsdup (bufferW);
}
}
}
closedir (d);
Strings_sort (me.peek());
return me.transfer();
} catch (MelderError) {
if (d) closedir (d); // "finally"
throw;
}
#elif defined (_WIN32)
try {
wchar_t searchPath [1+kMelder_MAXPATH];
int len = wcslen (path), hasAsterisk = wcschr (path, '*') != NULL, endsInSeparator = len != 0 && path [len - 1] == '\\';
autoStrings me = Thing_new (Strings);
my strings = NUMvector <wchar_t *> (1, 1000000);
swprintf (searchPath, 1+kMelder_MAXPATH, L"%ls%ls%ls", path, hasAsterisk || endsInSeparator ? L"" : L"\\", hasAsterisk ? L"" : L"*");
WIN32_FIND_DATAW findData;
HANDLE searchHandle = FindFirstFileW (searchPath, & findData);
if (searchHandle != INVALID_HANDLE_VALUE) {
do {
if ((type == Strings_createAsFileOrDirectoryList_TYPE_FILE &&
(findData. dwFileAttributes & FILE_ATTRIBUTE_DIRECTORY) == 0)
|| (type == Strings_createAsFileOrDirectoryList_TYPE_DIRECTORY &&
(findData. dwFileAttributes & FILE_ATTRIBUTE_DIRECTORY) != 0))
{
if (findData. cFileName [0] != '.') {
my strings [++ my numberOfStrings] = Melder_wcsdup (findData. cFileName);
}
}
} while (FindNextFileW (searchHandle, & findData));
FindClose (searchHandle);
}
Strings_sort (me.peek());
return me.transfer();
} catch (MelderError) {
throw;
}
#endif
}
autoEEG EEG_readFromBdfFile (MelderFile file) {
try {
autofile f = Melder_fopen (file, "rb");
char buffer [81];
fread (buffer, 1, 8, f); buffer [8] = '\0';
bool is24bit = buffer [0] == (char) 255;
fread (buffer, 1, 80, f); buffer [80] = '\0';
trace (U"Local subject identification: \"", Melder_peek8to32 (buffer), U"\"");
fread (buffer, 1, 80, f); buffer [80] = '\0';
trace (U"Local recording identification: \"", Melder_peek8to32 (buffer), U"\"");
fread (buffer, 1, 8, f); buffer [8] = '\0';
trace (U"Start date of recording: \"", Melder_peek8to32 (buffer), U"\"");
fread (buffer, 1, 8, f); buffer [8] = '\0';
trace (U"Start time of recording: \"", Melder_peek8to32 (buffer), U"\"");
fread (buffer, 1, 8, f); buffer [8] = '\0';
long numberOfBytesInHeaderRecord = atol (buffer);
trace (U"Number of bytes in header record: ", numberOfBytesInHeaderRecord);
fread (buffer, 1, 44, f); buffer [44] = '\0';
trace (U"Version of data format: \"", Melder_peek8to32 (buffer), U"\"");
fread (buffer, 1, 8, f); buffer [8] = '\0';
long numberOfDataRecords = strtol (buffer, nullptr, 10);
trace (U"Number of data records: ", numberOfDataRecords);
fread (buffer, 1, 8, f); buffer [8] = '\0';
double durationOfDataRecord = atof (buffer);
trace (U"Duration of a data record: ", durationOfDataRecord);
fread (buffer, 1, 4, f); buffer [4] = '\0';
long numberOfChannels = atol (buffer);
trace (U"Number of channels in data record: ", numberOfChannels);
if (numberOfBytesInHeaderRecord != (numberOfChannels + 1) * 256)
Melder_throw (U"Number of bytes in header record (", numberOfBytesInHeaderRecord,
U") doesn't match number of channels (", numberOfChannels, U").");
autostring32vector channelNames (1, numberOfChannels);
for (long ichannel = 1; ichannel <= numberOfChannels; ichannel ++) {
fread (buffer, 1, 16, f); buffer [16] = '\0'; // labels of the channels
/*
* Strip all final spaces.
*/
for (int i = 15; i >= 0; i --) {
if (buffer [i] == ' ') {
buffer [i] = '\0';
} else {
break;
}
}
channelNames [ichannel] = Melder_8to32 (buffer);
trace (U"Channel <<", channelNames [ichannel], U">>");
}
bool hasLetters = str32equ (channelNames [numberOfChannels], U"EDF Annotations");
double samplingFrequency = NUMundefined;
for (long channel = 1; channel <= numberOfChannels; channel ++) {
fread (buffer, 1, 80, f); buffer [80] = '\0'; // transducer type
}
for (long channel = 1; channel <= numberOfChannels; channel ++) {
fread (buffer, 1, 8, f); buffer [8] = '\0'; // physical dimension of channels
}
autoNUMvector <double> physicalMinimum (1, numberOfChannels);
for (long ichannel = 1; ichannel <= numberOfChannels; ichannel ++) {
fread (buffer, 1, 8, f); buffer [8] = '\0';
physicalMinimum [ichannel] = atof (buffer);
}
autoNUMvector <double> physicalMaximum (1, numberOfChannels);
for (long ichannel = 1; ichannel <= numberOfChannels; ichannel ++) {
fread (buffer, 1, 8, f); buffer [8] = '\0';
physicalMaximum [ichannel] = atof (buffer);
}
autoNUMvector <double> digitalMinimum (1, numberOfChannels);
for (long ichannel = 1; ichannel <= numberOfChannels; ichannel ++) {
fread (buffer, 1, 8, f); buffer [8] = '\0';
digitalMinimum [ichannel] = atof (buffer);
}
autoNUMvector <double> digitalMaximum (1, numberOfChannels);
for (long ichannel = 1; ichannel <= numberOfChannels; ichannel ++) {
fread (buffer, 1, 8, f); buffer [8] = '\0';
digitalMaximum [ichannel] = atof (buffer);
}
for (long channel = 1; channel <= numberOfChannels; channel ++) {
fread (buffer, 1, 80, f); buffer [80] = '\0'; // prefiltering
}
long numberOfSamplesPerDataRecord = 0;
for (long channel = 1; channel <= numberOfChannels; channel ++) {
fread (buffer, 1, 8, f); buffer [8] = '\0'; // number of samples in each data record
long numberOfSamplesInThisDataRecord = atol (buffer);
if (samplingFrequency == NUMundefined) {
numberOfSamplesPerDataRecord = numberOfSamplesInThisDataRecord;
samplingFrequency = numberOfSamplesInThisDataRecord / durationOfDataRecord;
}
if (numberOfSamplesInThisDataRecord / durationOfDataRecord != samplingFrequency)
Melder_throw (U"Number of samples per data record in channel ", channel,
U" (", numberOfSamplesInThisDataRecord,
U") doesn't match sampling frequency of channel 1 (", samplingFrequency, U").");
}
for (long channel = 1; channel <= numberOfChannels; channel ++) {
fread (buffer, 1, 32, f); buffer [32] = '\0'; // reserved
}
double duration = numberOfDataRecords * durationOfDataRecord;
autoEEG him = EEG_create (0, duration);
his numberOfChannels = numberOfChannels;
autoSound me = Sound_createSimple (numberOfChannels, duration, samplingFrequency);
Melder_assert (my nx == numberOfSamplesPerDataRecord * numberOfDataRecords);
autoNUMvector <unsigned char> dataBuffer (0L, 3 * numberOfSamplesPerDataRecord - 1);
for (long record = 1; record <= numberOfDataRecords; record ++) {
for (long channel = 1; channel <= numberOfChannels; channel ++) {
double factor = channel == numberOfChannels ? 1.0 : physicalMinimum [channel] / digitalMinimum [channel];
if (channel < numberOfChannels - EEG_getNumberOfExtraSensors (him.peek())) factor /= 1000000.0;
if (is24bit) {
fread (& dataBuffer [0], 3, numberOfSamplesPerDataRecord, f);
unsigned char *p = & dataBuffer [0];
for (long i = 1; i <= numberOfSamplesPerDataRecord; i ++) {
long sample = i + (record - 1) * numberOfSamplesPerDataRecord;
Melder_assert (sample <= my nx);
uint8_t lowByte = *p ++, midByte = *p ++, highByte = *p ++;
uint32_t externalValue = ((uint32_t) highByte << 16) | ((uint32_t) midByte << 8) | (uint32_t) lowByte;
if ((highByte & 128) != 0) // is the 24-bit sign bit on?
externalValue |= 0xFF000000; // extend negative sign to 32 bits
my z [channel] [sample] = (int32_t) externalValue * factor;
}
} else {
fread (& dataBuffer [0], 2, numberOfSamplesPerDataRecord, f);
unsigned char *p = & dataBuffer [0];
for (long i = 1; i <= numberOfSamplesPerDataRecord; i ++) {
long sample = i + (record - 1) * numberOfSamplesPerDataRecord;
Melder_assert (sample <= my nx);
uint8 lowByte = *p ++, highByte = *p ++;
uint16 externalValue = (uint16) ((uint16) highByte << 8) | (uint16) lowByte;
my z [channel] [sample] = (int16) externalValue * factor;
}
}
}
}
int numberOfStatusBits = 8;
for (long i = 1; i <= my nx; i ++) {
unsigned long value = (long) my z [numberOfChannels] [i];
if (value & 0x0000FF00) {
numberOfStatusBits = 16;
}
}
autoTextGrid thee;
if (hasLetters) {
thee = TextGrid_create (0, duration, U"Mark Trigger", U"Mark Trigger");
autoMelderString letters;
double time = NUMundefined;
for (long i = 1; i <= my nx; i ++) {
unsigned long value = (long) my z [numberOfChannels] [i];
for (int byte = 1; byte <= numberOfStatusBits / 8; byte ++) {
unsigned long mask = byte == 1 ? 0x000000ff : 0x0000ff00;
char32 kar = byte == 1 ? (value & mask) : (value & mask) >> 8;
if (kar != U'\0' && kar != 20) {
MelderString_appendCharacter (& letters, kar);
} else if (letters. string [0] != U'\0') {
if (letters. string [0] == U'+') {
if (NUMdefined (time)) {
try {
TextGrid_insertPoint (thee.peek(), 1, time, U"");
} catch (MelderError) {
Melder_throw (U"Did not insert empty mark (", letters. string, U") on Mark tier.");
}
time = NUMundefined; // defensive
}
time = Melder_atof (& letters. string [1]);
MelderString_empty (& letters);
} else {
if (! NUMdefined (time)) {
Melder_throw (U"Undefined time for label at sample ", i, U".");
}
try {
if (Melder_nequ (letters. string, U"Trigger-", 8)) {
try {
TextGrid_insertPoint (thee.peek(), 2, time, & letters. string [8]);
} catch (MelderError) {
Melder_clearError ();
trace (U"Duplicate trigger at ", time, U" seconds: ", & letters. string [8]);
}
} else {
TextGrid_insertPoint (thee.peek(), 1, time, & letters. string [0]);
}
} catch (MelderError) {
Melder_throw (U"Did not insert mark (", letters. string, U") on Trigger tier.");
}
time = NUMundefined; // crucial
MelderString_empty (& letters);
}
}
}
}
if (NUMdefined (time)) {
TextGrid_insertPoint (thee.peek(), 1, time, U"");
time = NUMundefined; // defensive
}
} else {
thee = TextGrid_create (0, duration,
numberOfStatusBits == 8 ? U"S1 S2 S3 S4 S5 S6 S7 S8" : U"S1 S2 S3 S4 S5 S6 S7 S8 S9 S10 S11 S12 S13 S14 S15 S16", U"");
for (int bit = 1; bit <= numberOfStatusBits; bit ++) {
unsigned long bitValue = 1 << (bit - 1);
IntervalTier tier = (IntervalTier) thy tiers -> item [bit];
for (long i = 1; i <= my nx; i ++) {
unsigned long previousValue = i == 1 ? 0 : (long) my z [numberOfChannels] [i - 1];
unsigned long thisValue = (long) my z [numberOfChannels] [i];
if ((thisValue & bitValue) != (previousValue & bitValue)) {
double time = i == 1 ? 0.0 : my x1 + (i - 1.5) * my dx;
if (time != 0.0)
TextGrid_insertBoundary (thee.peek(), bit, time);
if ((thisValue & bitValue) != 0)
TextGrid_setIntervalText (thee.peek(), bit, tier -> intervals -> size, U"1");
}
}
}
}
f.close (file);
his channelNames = channelNames.transfer();
his sound = me.move();
his textgrid = thee.move();
if (EEG_getNumberOfCapElectrodes (him.peek()) == 32) {
EEG_setChannelName (him.peek(), 1, U"Fp1");
EEG_setChannelName (him.peek(), 2, U"AF3");
EEG_setChannelName (him.peek(), 3, U"F7");
EEG_setChannelName (him.peek(), 4, U"F3");
EEG_setChannelName (him.peek(), 5, U"FC1");
EEG_setChannelName (him.peek(), 6, U"FC5");
EEG_setChannelName (him.peek(), 7, U"T7");
EEG_setChannelName (him.peek(), 8, U"C3");
EEG_setChannelName (him.peek(), 9, U"CP1");
EEG_setChannelName (him.peek(), 10, U"CP5");
EEG_setChannelName (him.peek(), 11, U"P7");
EEG_setChannelName (him.peek(), 12, U"P3");
EEG_setChannelName (him.peek(), 13, U"Pz");
EEG_setChannelName (him.peek(), 14, U"PO3");
EEG_setChannelName (him.peek(), 15, U"O1");
EEG_setChannelName (him.peek(), 16, U"Oz");
EEG_setChannelName (him.peek(), 17, U"O2");
EEG_setChannelName (him.peek(), 18, U"PO4");
EEG_setChannelName (him.peek(), 19, U"P4");
EEG_setChannelName (him.peek(), 20, U"P8");
EEG_setChannelName (him.peek(), 21, U"CP6");
EEG_setChannelName (him.peek(), 22, U"CP2");
EEG_setChannelName (him.peek(), 23, U"C4");
EEG_setChannelName (him.peek(), 24, U"T8");
EEG_setChannelName (him.peek(), 25, U"FC6");
EEG_setChannelName (him.peek(), 26, U"FC2");
EEG_setChannelName (him.peek(), 27, U"F4");
EEG_setChannelName (him.peek(), 28, U"F8");
EEG_setChannelName (him.peek(), 29, U"AF4");
EEG_setChannelName (him.peek(), 30, U"Fp2");
EEG_setChannelName (him.peek(), 31, U"Fz");
EEG_setChannelName (him.peek(), 32, U"Cz");
} else if (EEG_getNumberOfCapElectrodes (him.peek()) == 64) {
EEG_setChannelName (him.peek(), 1, U"Fp1");
EEG_setChannelName (him.peek(), 2, U"AF7");
EEG_setChannelName (him.peek(), 3, U"AF3");
EEG_setChannelName (him.peek(), 4, U"F1");
EEG_setChannelName (him.peek(), 5, U"F3");
EEG_setChannelName (him.peek(), 6, U"F5");
EEG_setChannelName (him.peek(), 7, U"F7");
EEG_setChannelName (him.peek(), 8, U"FT7");
EEG_setChannelName (him.peek(), 9, U"FC5");
EEG_setChannelName (him.peek(), 10, U"FC3");
EEG_setChannelName (him.peek(), 11, U"FC1");
EEG_setChannelName (him.peek(), 12, U"C1");
EEG_setChannelName (him.peek(), 13, U"C3");
EEG_setChannelName (him.peek(), 14, U"C5");
EEG_setChannelName (him.peek(), 15, U"T7");
EEG_setChannelName (him.peek(), 16, U"TP7");
EEG_setChannelName (him.peek(), 17, U"CP5");
EEG_setChannelName (him.peek(), 18, U"CP3");
EEG_setChannelName (him.peek(), 19, U"CP1");
EEG_setChannelName (him.peek(), 20, U"P1");
EEG_setChannelName (him.peek(), 21, U"P3");
EEG_setChannelName (him.peek(), 22, U"P5");
EEG_setChannelName (him.peek(), 23, U"P7");
EEG_setChannelName (him.peek(), 24, U"P9");
EEG_setChannelName (him.peek(), 25, U"PO7");
EEG_setChannelName (him.peek(), 26, U"PO3");
EEG_setChannelName (him.peek(), 27, U"O1");
EEG_setChannelName (him.peek(), 28, U"Iz");
EEG_setChannelName (him.peek(), 29, U"Oz");
EEG_setChannelName (him.peek(), 30, U"POz");
EEG_setChannelName (him.peek(), 31, U"Pz");
EEG_setChannelName (him.peek(), 32, U"CPz");
EEG_setChannelName (him.peek(), 33, U"Fpz");
EEG_setChannelName (him.peek(), 34, U"Fp2");
EEG_setChannelName (him.peek(), 35, U"AF8");
EEG_setChannelName (him.peek(), 36, U"AF4");
EEG_setChannelName (him.peek(), 37, U"AFz");
EEG_setChannelName (him.peek(), 38, U"Fz");
EEG_setChannelName (him.peek(), 39, U"F2");
EEG_setChannelName (him.peek(), 40, U"F4");
EEG_setChannelName (him.peek(), 41, U"F6");
EEG_setChannelName (him.peek(), 42, U"F8");
EEG_setChannelName (him.peek(), 43, U"FT8");
EEG_setChannelName (him.peek(), 44, U"FC6");
EEG_setChannelName (him.peek(), 45, U"FC4");
EEG_setChannelName (him.peek(), 46, U"FC2");
EEG_setChannelName (him.peek(), 47, U"FCz");
EEG_setChannelName (him.peek(), 48, U"Cz");
EEG_setChannelName (him.peek(), 49, U"C2");
EEG_setChannelName (him.peek(), 50, U"C4");
EEG_setChannelName (him.peek(), 51, U"C6");
EEG_setChannelName (him.peek(), 52, U"T8");
EEG_setChannelName (him.peek(), 53, U"TP8");
EEG_setChannelName (him.peek(), 54, U"CP6");
EEG_setChannelName (him.peek(), 55, U"CP4");
EEG_setChannelName (him.peek(), 56, U"CP2");
EEG_setChannelName (him.peek(), 57, U"P2");
EEG_setChannelName (him.peek(), 58, U"P4");
EEG_setChannelName (him.peek(), 59, U"P6");
EEG_setChannelName (him.peek(), 60, U"P8");
EEG_setChannelName (him.peek(), 61, U"P10");
EEG_setChannelName (him.peek(), 62, U"PO8");
EEG_setChannelName (him.peek(), 63, U"PO4");
EEG_setChannelName (him.peek(), 64, U"O2");
}
return him;
} catch (MelderError) {
static autoStrings Strings_createAsFileOrDirectoryList (const char32 *path /* cattable */, int type) {
#if USE_STAT
/*
* Initialize.
*/
DIR *d = nullptr;
try {
autoMelderString filePath, searchDirectory, left, right;
/*
* Parse the path.
* Example: in /Users/paul/sounds/h*.wav",
* the search directory is "/Users/paul/sounds",
* the left environment is "h", and the right environment is ".wav".
*/
MelderString_copy (& searchDirectory, path);
char32 *asterisk = str32rchr (searchDirectory. string, '*');
if (asterisk) {
*asterisk = '\0';
searchDirectory. length = asterisk - searchDirectory. string; // probably superfluous, but correct
char32 *lastSlash = str32rchr (searchDirectory. string, Melder_DIRECTORY_SEPARATOR);
if (lastSlash) {
*lastSlash = '\0'; // This fixes searchDirectory.
searchDirectory. length = lastSlash - searchDirectory. string; // probably superfluous, but correct
MelderString_copy (& left, lastSlash + 1);
} else {
MelderString_copy (& left, searchDirectory. string); // quickly save...
MelderString_empty (& searchDirectory); // ...before destruction
}
MelderString_copy (& right, asterisk + 1);
}
char buffer8 [kMelder_MAXPATH+1];
Melder_str32To8bitFileRepresentation_inline (searchDirectory. string, buffer8);
d = opendir (buffer8 [0] ? buffer8 : ".");
if (! d)
Melder_throw (U"Cannot open directory ", searchDirectory. string, U".");
//Melder_casual (U"opened");
autoStrings me = Thing_new (Strings);
my strings = NUMvector <char32 *> (1, 1000000);
struct dirent *entry;
while (!! (entry = readdir (d))) {
MelderString_copy (& filePath, searchDirectory. string [0] ? searchDirectory. string : U".");
MelderString_appendCharacter (& filePath, Melder_DIRECTORY_SEPARATOR);
char32 buffer32 [kMelder_MAXPATH+1];
Melder_8bitFileRepresentationToStr32_inline (entry -> d_name, buffer32);
MelderString_append (& filePath, buffer32);
//Melder_casual (U"read ", filePath. string);
Melder_str32To8bitFileRepresentation_inline (filePath. string, buffer8);
struct stat stats;
if (stat (buffer8, & stats) != 0) {
//Melder_throw (U"Cannot look at file ", filePath. string, U".");
//stats. st_mode = -1L;
}
//Melder_casual (U"statted ", filePath. string);
//Melder_casual (U"file ", filePath. string, U" mode ", stats. st_mode / 4096);
if ((type == Strings_createAsFileOrDirectoryList_TYPE_FILE && S_ISREG (stats. st_mode)) ||
(type == Strings_createAsFileOrDirectoryList_TYPE_DIRECTORY && S_ISDIR (stats. st_mode)))
{
Melder_8bitFileRepresentationToStr32_inline (entry -> d_name, buffer32);
int64 length = str32len (buffer32);
if (buffer32 [0] != U'.' &&
(left. length == 0 || str32nequ (buffer32, left. string, left. length)) &&
(right. length == 0 || (length >= right. length && str32equ (buffer32 + (length - right. length), right. string))))
{
my strings [++ my numberOfStrings] = Melder_dup (buffer32);
}
}
}
closedir (d);
Strings_sort (me.get());
return me;
} catch (MelderError) {
if (d) closedir (d); // "finally"
throw;
}
#elif defined (_WIN32)
try {
char32 searchPath [kMelder_MAXPATH+1];
int len = str32len (path);
bool hasAsterisk = !! str32chr (path, U'*');
bool endsInSeparator = ( len != 0 && path [len - 1] == U'\\' );
autoStrings me = Thing_new (Strings);
my strings = NUMvector <char32 *> (1, 1000000);
Melder_sprint (searchPath, kMelder_MAXPATH+1, path, hasAsterisk || endsInSeparator ? U"" : U"\\", hasAsterisk ? U"" : U"*");
WIN32_FIND_DATAW findData;
HANDLE searchHandle = FindFirstFileW (Melder_peek32toW (searchPath), & findData);
if (searchHandle != INVALID_HANDLE_VALUE) {
do {
if ((type == Strings_createAsFileOrDirectoryList_TYPE_FILE &&
(findData. dwFileAttributes & FILE_ATTRIBUTE_DIRECTORY) == 0)
|| (type == Strings_createAsFileOrDirectoryList_TYPE_DIRECTORY &&
(findData. dwFileAttributes & FILE_ATTRIBUTE_DIRECTORY) != 0))
{
if (findData. cFileName [0] != L'.') {
my strings [++ my numberOfStrings] = Melder_dup (Melder_peekWto32 (findData. cFileName));
}
}
} while (FindNextFileW (searchHandle, & findData));
FindClose (searchHandle);
}
Strings_sort (me.get());
return me;
} catch (MelderError) {
throw;
}
#endif
}
TableOfReal TableOfReal_readFromHeaderlessSpreadsheetFile (MelderFile file) {
try {
autostring string = MelderFile_readText (file);
long nrow, ncol, nelements;
/*
* Count columns.
*/
ncol = 0;
wchar_t *p = & string [0];
for (;;) {
wchar_t kar = *p++;
if (kar == '\n' || kar == '\0') break;
if (kar == ' ' || kar == '\t') continue;
ncol ++;
do { kar = *p++; } while (kar != ' ' && kar != '\t' && kar != '\n' && kar != '\0');
if (kar == '\n' || kar == '\0') break;
}
ncol --;
if (ncol < 1) Melder_throw ("No columns.");
/*
* Count elements.
*/
p = & string [0];
nelements = 0;
for (;;) {
wchar_t kar = *p++;
if (kar == '\0') break;
if (kar == ' ' || kar == '\t' || kar == '\n') continue;
nelements ++;
do { kar = *p++; } while (kar != ' ' && kar != '\t' && kar != '\n' && kar != '\0');
if (kar == '\0') break;
}
/*
* Check if all columns are complete.
*/
if (nelements == 0 || nelements % (ncol + 1) != 0)
Melder_throw ("The number of elements (", nelements, ") is not a multiple of the number of columns plus 1 (", ncol + 1, ").");
/*
* Create empty table.
*/
nrow = nelements / (ncol + 1) - 1;
autoTableOfReal me = TableOfReal_create (nrow, ncol);
/*
* Read elements.
*/
p = & string [0];
while (*p == ' ' || *p == '\t') { Melder_assert (*p != '\0'); p ++; }
while (*p != ' ' && *p != '\t') { Melder_assert (*p != '\0'); p ++; } // ignore the header of the zeroth column ("rowLabel" perhaps)
for (long icol = 1; icol <= ncol; icol ++) {
while (*p == ' ' || *p == '\t') { Melder_assert (*p != '\0'); p ++; }
static MelderString buffer = { 0 };
MelderString_empty (& buffer);
while (*p != ' ' && *p != '\t' && *p != '\n') {
MelderString_appendCharacter (& buffer, *p);
p ++;
}
TableOfReal_setColumnLabel (me.peek(), icol, buffer.string);
MelderString_empty (& buffer);
}
for (long irow = 1; irow <= nrow; irow ++) {
while (*p == ' ' || *p == '\t' || *p == '\n') { Melder_assert (*p != '\0'); p ++; }
static MelderString buffer = { 0 };
MelderString_empty (& buffer);
while (*p != ' ' && *p != '\t') {
MelderString_appendCharacter (& buffer, *p);
p ++;
}
TableOfReal_setRowLabel (me.peek(), irow, buffer.string);
MelderString_empty (& buffer);
for (long icol = 1; icol <= ncol; icol ++) {
while (*p == ' ' || *p == '\t' || *p == '\n') { Melder_assert (*p != '\0'); p ++; }
MelderString_empty (& buffer);
while (*p != ' ' && *p != '\t' && *p != '\n' && *p != '\0') {
MelderString_appendCharacter (& buffer, *p);
p ++;
}
my data [irow] [icol] = Melder_atof (buffer.string); /* If cell contains a string, this will be 0. */
MelderString_empty (& buffer);
}
}
return me.transfer();
} catch (MelderError) {
Melder_throw ("TableOfReal: tab-separated file ", file, " not read.");
}
}
