/* abcio.cpp * * Copyright (C) 1992-2011,2015 Paul Boersma * * This program is free software; you can redistribute it and/or modify * it under the terms of the GNU General Public License as published by * the Free Software Foundation; either version 2 of the License, or (at * your option) any later version. * * This program is distributed in the hope that it will be useful, but * WITHOUT ANY WARRANTY; without even the implied warranty of * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. * See the GNU General Public License for more details. * * You should have received a copy of the GNU General Public License * along with this program; if not, write to the Free Software * Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA. */ /* * pb 2002/03/07 GPL * pb 2003/05/19 accept percent signs in getReal * pb 2004/10/01 Melder_double instead of %.17g * pb 2006/02/17 support for Intel-based Macs * pb 2006/02/20 corrected bingeti3, bingeti3LE, binputi3, binputi3LE * pb 2006/03/28 support for systems where a long is not 32 bits and a short is not 16 bits * pb 2007/07/21 MelderReadString * pb 2007/08/14 check for NULL pointer before Melder_isValidAscii * pb 2009/03/18 modern enums * fb 2010/02/26 UTF-16 via bin(get|put)utf16() * pb 2010/03/09 more support for Unicode values above 0xFFFF * pb 2010/12/23 corrected bingeti3 and bingeti3LE for 64-bit systems * pb 2011/03/30 C++ */ #include "melder.h" #include "NUM.h" #include #ifdef macintosh #include #endif #include "abcio.h" /********** text I/O **********/ static long getInteger (MelderReadText me) { char buffer [41]; char32 c; /* * Look for the first numeric character. */ for (c = MelderReadText_getChar (me); c != U'-' && ! isdigit ((int) c) && c != U'+'; c = MelderReadText_getChar (me)) { if (c == U'\0') Melder_throw (U"Early end of text detected while looking for an integer (line ", MelderReadText_getLineNumber (me), U")."); if (c == U'!') { // end-of-line comment? while ((c = MelderReadText_getChar (me)) != U'\n' && c != U'\r') { if (c == 0) Melder_throw (U"Early end of text detected in comment while looking for an integer (line ", MelderReadText_getLineNumber (me), U")."); } } if (c == U'\"') Melder_throw (U"Found a string while looking for an integer in text (line ", MelderReadText_getLineNumber (me), U")."); if (c == U'<') Melder_throw (U"Found an enumerated value while looking for an integer in text (line ", MelderReadText_getLineNumber (me), U")."); while (c != U' ' && c != U'\n' && c != U'\t' && c != U'\r') { if (c == U'\0') Melder_throw (U"Early end of text detected in comment (line ", MelderReadText_getLineNumber (me), U")."); c = MelderReadText_getChar (me); } } int i = 0; for (; i < 40; i ++) { if (c > 127) Melder_throw (U"Found strange text while looking for an integer in text (line ", MelderReadText_getLineNumber (me), U")."); buffer [i] = (char) (char8) c; // guarded conversion down c = MelderReadText_getChar (me); if (c == U'\0') { break; } // this may well be OK here if (c == U' ' || c == U'\n' || c == U'\t' || c == U'\r') break; } if (i >= 40) Melder_throw (U"Found long text while looking for an integer in text (line ", MelderReadText_getLineNumber (me), U")."); buffer [i + 1] = '\0'; return strtol (buffer, NULL, 10); } static unsigned long getUnsigned (MelderReadText me) { char buffer [41]; char32 c; for (c = MelderReadText_getChar (me); ! isdigit ((int) c) && c != U'+'; c = MelderReadText_getChar (me)) { if (c == U'\0') Melder_throw (U"Early end of text detected while looking for an unsigned integer (line ", MelderReadText_getLineNumber (me), U")."); if (c == U'!') { // end-of-line comment? while ((c = MelderReadText_getChar (me)) != '\n' && c != '\r') { if (c == U'\0') Melder_throw (U"Early end of text detected in comment while looking for an unsigned integer (line ", MelderReadText_getLineNumber (me), U")."); } } if (c == U'\"') Melder_throw (U"Found a string while looking for an unsigned integer in text (line ", MelderReadText_getLineNumber (me), U")."); if (c == U'<') Melder_throw (U"Found an enumerated value while looking for an unsigned integer in text (line ", MelderReadText_getLineNumber (me), U")."); if (c == U'-') Melder_throw (U"Found a negative value while looking for an unsigned integer in text (line ", MelderReadText_getLineNumber (me), U")."); while (c != U' ' && c != U'\n' && c != U'\t' && c != U'\r') { if (c == U'\0') Melder_throw (U"Early end of text detected in comment (line ", MelderReadText_getLineNumber (me), U")."); c = MelderReadText_getChar (me); } } int i = 0; for (i = 0; i < 40; i ++) { if (c > 127) Melder_throw (U"Found strange text while looking for an unsigned integer in text (line ", MelderReadText_getLineNumber (me), U")."); buffer [i] = (char) (char8) c; // guarded conversion down c = MelderReadText_getChar (me); if (c == U'\0') { break; } // this may well be OK here if (c == U' ' || c == U'\n' || c == U'\t' || c == U'\r') break; } if (i >= 40) Melder_throw (U"Found long text while searching for an unsigned integer in text (line ", MelderReadText_getLineNumber (me), U")."); buffer [i + 1] = '\0'; return strtoul (buffer, NULL, 10); } static double getReal (MelderReadText me) { int i; char buffer [41], *slash; char32 c; do { for (c = MelderReadText_getChar (me); c != U'-' && ! isdigit ((int) c) && c != U'+'; c = MelderReadText_getChar (me)) { if (c == U'\0') Melder_throw (U"Early end of text detected while looking for a real number (line ", MelderReadText_getLineNumber (me), U")."); if (c == U'!') { // end-of-line comment? while ((c = MelderReadText_getChar (me)) != U'\n' && c != U'\r') { if (c == U'\0') Melder_throw (U"Early end of text detected in comment while looking for a real number (line ", MelderReadText_getLineNumber (me), U")."); } } if (c == U'\"') Melder_throw (U"Found a string while looking for a real number in text (line ", MelderReadText_getLineNumber (me), U")."); if (c == U'<') Melder_throw (U"Found an enumerated value while looking for a real number in text (line ", MelderReadText_getLineNumber (me), U")."); while (c != U' ' && c != U'\n' && c != U'\t' && c != U'\r') { if (c == U'\0') Melder_throw (U"Early end of text detected in comment while looking for a real number (line ", MelderReadText_getLineNumber (me), U")."); c = MelderReadText_getChar (me); } } for (i = 0; i < 40; i ++) { if (c > 127) Melder_throw (U"Found strange text while looking for a real number in text (line ", MelderReadText_getLineNumber (me), U")."); buffer [i] = (char) (char8) c; // guarded conversion down c = MelderReadText_getChar (me); if (c == U'\0') { break; } // this may well be OK here if (c == U' ' || c == U'\n' || c == U'\t' || c == U'\r') break; } if (i >= 40) Melder_throw (U"Found long text while searching for a real number in text (line ", MelderReadText_getLineNumber (me), U")."); } while (i == 0 && buffer [0] == '+'); // guard against single '+' symbols, which occur in complex numbers buffer [i + 1] = '\0'; slash = strchr (buffer, '/'); if (slash) { double numerator, denominator; *slash = '\0'; numerator = Melder_a8tof (buffer), denominator = Melder_a8tof (slash + 1); if (numerator == HUGE_VAL || denominator == HUGE_VAL || denominator == 0.0) return HUGE_VAL; return numerator / denominator; } return Melder_a8tof (buffer); } static short getEnum (MelderReadText me, int (*getValue) (const char32 *)) { char32 buffer [41], c; for (c = MelderReadText_getChar (me); c != U'<'; c = MelderReadText_getChar (me)) { if (c == U'\0') Melder_throw (U"Early end of text detected while looking for an enumerated value (line ", MelderReadText_getLineNumber (me), U")."); if (c == U'!') { /* End-of-line comment? */ while ((c = MelderReadText_getChar (me)) != U'\n' && c != U'\r') { if (c == U'\0') Melder_throw (U"Early end of text detected in comment while looking for an enumerated value (line ", MelderReadText_getLineNumber (me), U")."); } } if (c == U'-' || isdigit ((int) c) || c == U'+') Melder_throw (U"Found a number while looking for an enumerated value in text (line ", MelderReadText_getLineNumber (me), U")."); if (c == U'\"') Melder_throw (U"Found a string while looking for an enumerated value in text (line ", MelderReadText_getLineNumber (me), U")."); while (c != U' ' && c != U'\n' && c != U'\t' && c != U'\r') { if (c == U'\0') Melder_throw (U"Early end of text detected in comment while looking for an enumerated value (line ", MelderReadText_getLineNumber (me), U")."); c = MelderReadText_getChar (me); } } int i = 0; for (; i < 40; i ++) { c = MelderReadText_getChar (me); // read past first '<' if (c == U'\0') Melder_throw (U"Early end of text detected while reading an enumerated value (line ", MelderReadText_getLineNumber (me), U")."); if (c == U' ' || c == U'\n' || c == U'\t' || c == U'\r') Melder_throw (U"No matching '>' while reading an enumerated value (line ", MelderReadText_getLineNumber (me), U")."); if (c == U'>') break; // the expected closing bracket; not added to the buffer buffer [i] = c; } if (i >= 40) Melder_throw (U"Found strange text while reading an enumerated value in text (line ", MelderReadText_getLineNumber (me), U")."); buffer [i] = U'\0'; int value = getValue (buffer); if (value < 0) Melder_throw (U"\"", buffer, U"\" is not a value of the enumerated type."); return value; } static char32 * getString (MelderReadText me) { static MelderString buffer { 0 }; MelderString_empty (& buffer); for (char32 c = MelderReadText_getChar (me); c != U'\"'; c = MelderReadText_getChar (me)) { if (c == U'\0') Melder_throw (U"Early end of text detected while looking for a string (line ", MelderReadText_getLineNumber (me), U")."); if (c == U'!') { // end-of-line comment? while ((c = MelderReadText_getChar (me)) != '\n' && c != '\r') { if (c == U'\0') Melder_throw (U"Early end of text detected in comment while looking for a string (line ", MelderReadText_getLineNumber (me), U")."); } } if (c == U'-' || isdigit ((int) c) || c == U'+') Melder_throw (U"Found a number while looking for a string in text (line ", MelderReadText_getLineNumber (me), U")."); if (c == U'<') Melder_throw (U"Found an enumerated value while looking for a string in text (line ", MelderReadText_getLineNumber (me), U")."); while (c != U' ' && c != U'\n' && c != U'\t' && c != U'\r') { if (c == U'\0') Melder_throw (U"Early end of text detected while looking for a string (line ", MelderReadText_getLineNumber (me), U")."); c = MelderReadText_getChar (me); } } for (int i = 0; 1; i ++) { char32 c = MelderReadText_getChar (me); // read past first '"' if (c == U'\0') Melder_throw (U"Early end of text detected while reading a string (line ", MelderReadText_getLineNumber (me), U")."); if (c == U'\"') { char32 next = MelderReadText_getChar (me); if (next == U'\0') { break; } // closing quote is last character in file: OK if (next != U'\"') { if (next == U' ' || next == U'\n' || next == U'\t' || next == U'\r') { // closing quote is followed by whitespace: it is OK to skip this whitespace (no need to "ungetChar") } else { char32 kar2 [2] = { next, U'\0' }; Melder_throw (U"Character ", kar2, U" following quote (line ", MelderReadText_getLineNumber (me), U"). End of string or undoubled quote?"); } break; // the expected closing double quote; not added to the buffer } // else: add only one of the two quotes to the buffer } MelderString_appendCharacter (& buffer, c); } return buffer. string; } #undef false #undef true #include "enums_getText.h" #include "abcio_enums.h" #include "enums_getValue.h" #include "abcio_enums.h" int texgeti1 (MelderReadText text) { try { long externalValue = getInteger (text); if (externalValue < -128 || externalValue > +127) Melder_throw (U"Value (", externalValue, U") out of range (-128 .. +127)."); return (int) externalValue; } catch (MelderError) { Melder_throw (U"Signed small integer not read from text file."); } } int texgeti2 (MelderReadText text) { try { long externalValue = getInteger (text); if (externalValue < -32768 || externalValue > +32767) Melder_throw (U"Value (", externalValue, U") out of range (-32768 .. +32767)."); return (int) externalValue; } catch (MelderError) { Melder_throw (U"Signed short integer not read from text file."); } } long texgeti4 (MelderReadText text) { try { long externalValue = getInteger (text); return externalValue; } catch (MelderError) { Melder_throw (U"Signed integer not read from text file."); } } unsigned int texgetu1 (MelderReadText text) { try { unsigned long externalValue = getUnsigned (text); if (externalValue > 255) Melder_throw (U"Value (", externalValue, U") out of range (0 .. 255)."); return (unsigned int) externalValue; } catch (MelderError) { Melder_throw (U"Unsigned small integer not read from text file."); } } unsigned int texgetu2 (MelderReadText text) { try { unsigned long externalValue = getUnsigned (text); if (externalValue > 65535) Melder_throw (U"Value (", externalValue, U") out of range (0 .. 65535)."); return (unsigned int) externalValue; } catch (MelderError) { Melder_throw (U"Unsigned short integer not read from text file."); } } unsigned long texgetu4 (MelderReadText text) { try { unsigned long externalValue = getUnsigned (text); return externalValue; } catch (MelderError) { Melder_throw (U"Unsigned integer not read from text file."); } } double texgetr4 (MelderReadText text) { return getReal (text); } double texgetr8 (MelderReadText text) { return getReal (text); } double texgetr10 (MelderReadText text) { return getReal (text); } fcomplex texgetc8 (MelderReadText text) { fcomplex z; z.re = getReal (text); z.im = getReal (text); return z; } dcomplex texgetc16 (MelderReadText text) { dcomplex z; z.re = getReal (text); z.im = getReal (text); return z; } short texgete1 (MelderReadText text, int (*getValue) (const char32 *)) { return getEnum (text, getValue); } short texgete2 (MelderReadText text, int (*getValue) (const char32 *)) { return getEnum (text, getValue); } short texgeteb (MelderReadText text) { return getEnum (text, kBoolean_getValue); } short texgeteq (MelderReadText text) { return getEnum (text, kQuestion_getValue); } short texgetex (MelderReadText text) { return getEnum (text, kExistence_getValue); } char *texgets2 (MelderReadText text) { return (char *) Melder_32to8 (getString (text)); } char *texgets4 (MelderReadText text) { return (char *) Melder_32to8 (getString (text)); } char32 *texgetw2 (MelderReadText text) { return Melder_dup (getString (text)); } char32 *texgetw4 (MelderReadText text) { return Melder_dup (getString (text)); } void texindent (MelderFile file) { file -> indent += 4; } void texexdent (MelderFile file) { file -> indent -= 4; } void texresetindent (MelderFile file) { file -> indent = 0; } void texputintro (MelderFile file, const char32 *s1, const char32 *s2, const char32 *s3, const char32 *s4, const char32 *s5, const char32 *s6) { if (file -> verbose) { MelderFile_write (file, U"\n"); for (int iindent = 1; iindent <= file -> indent; iindent ++) { MelderFile_write (file, U" "); } MelderFile_write (file, s1 && s1 [0] == U'd' && s1 [1] == U'_' ? & s1 [2] : & s1 [0], s2 && s2 [0] == U'd' && s2 [1] == U'_' ? & s2 [2] : & s2 [0], s3 && s3 [0] == U'd' && s3 [1] == U'_' ? & s3 [2] : & s3 [0], s4 && s4 [0] == U'd' && s4 [1] == U'_' ? & s4 [2] : & s4 [0], s5 && s5 [0] == U'd' && s5 [1] == U'_' ? & s5 [2] : & s5 [0], s6 && s6 [0] == U'd' && s6 [1] == U'_' ? & s6 [2] : & s6 [0]); } file -> indent += 4; } #define PUTLEADER \ MelderFile_write (file, U"\n"); \ if (file -> verbose) { \ for (int iindent = 1; iindent <= file -> indent; iindent ++) { \ MelderFile_write (file, U" "); \ } \ MelderFile_write (file, \ s1 && s1 [0] == U'd' && s1 [1] == U'_' ? & s1 [2] : & s1 [0], \ s2 && s2 [0] == U'd' && s2 [1] == U'_' ? & s2 [2] : & s2 [0], \ s3 && s3 [0] == U'd' && s3 [1] == U'_' ? & s3 [2] : & s3 [0], \ s4 && s4 [0] == U'd' && s4 [1] == U'_' ? & s4 [2] : & s4 [0], \ s5 && s5 [0] == U'd' && s5 [1] == U'_' ? & s5 [2] : & s5 [0], \ s6 && s6 [0] == U'd' && s6 [1] == U'_' ? & s6 [2] : & s6 [0]); \ } void texputi1 (MelderFile file, int i, const char32 *s1, const char32 *s2, const char32 *s3, const char32 *s4, const char32 *s5, const char32 *s6) { PUTLEADER MelderFile_write (file, file -> verbose ? U" = " : NULL, i, file -> verbose ? U" " : NULL); } void texputi2 (MelderFile file, int i, const char32 *s1, const char32 *s2, const char32 *s3, const char32 *s4, const char32 *s5, const char32 *s6) { PUTLEADER MelderFile_write (file, file -> verbose ? U" = " : NULL, i, file -> verbose ? U" " : NULL); } void texputi4 (MelderFile file, long i, const char32 *s1, const char32 *s2, const char32 *s3, const char32 *s4, const char32 *s5, const char32 *s6) { PUTLEADER MelderFile_write (file, file -> verbose ? U" = " : NULL, i, file -> verbose ? U" " : NULL); } void texputu1 (MelderFile file, unsigned int u, const char32 *s1, const char32 *s2, const char32 *s3, const char32 *s4, const char32 *s5, const char32 *s6) { PUTLEADER MelderFile_write (file, file -> verbose ? U" = " : NULL, u, file -> verbose ? U" " : NULL); } void texputu2 (MelderFile file, unsigned int u, const char32 *s1, const char32 *s2, const char32 *s3, const char32 *s4, const char32 *s5, const char32 *s6) { PUTLEADER MelderFile_write (file, file -> verbose ? U" = " : NULL, u, file -> verbose ? U" " : NULL); } void texputu4 (MelderFile file, unsigned long u, const char32 *s1, const char32 *s2, const char32 *s3, const char32 *s4, const char32 *s5, const char32 *s6) { PUTLEADER MelderFile_write (file, file -> verbose ? U" = " : NULL, u, file -> verbose ? U" " : NULL); } void texputr4 (MelderFile file, double x, const char32 *s1, const char32 *s2, const char32 *s3, const char32 *s4, const char32 *s5, const char32 *s6) { PUTLEADER MelderFile_write (file, file -> verbose ? U" = " : NULL, Melder_single (x), file -> verbose ? U" " : NULL); } void texputr8 (MelderFile file, double x, const char32 *s1, const char32 *s2, const char32 *s3, const char32 *s4, const char32 *s5, const char32 *s6) { PUTLEADER MelderFile_write (file, file -> verbose ? U" = " : NULL, x, file -> verbose ? U" " : NULL); } void texputc8 (MelderFile file, fcomplex z, const char32 *s1, const char32 *s2, const char32 *s3, const char32 *s4, const char32 *s5, const char32 *s6) { PUTLEADER MelderFile_write (file, file -> verbose ? U" = " : NULL, Melder_single (z.re), file -> verbose ? U" + " : U" ", Melder_single (z.im), file -> verbose ? U" i " : NULL); } void texputc16 (MelderFile file, dcomplex z, const char32 *s1, const char32 *s2, const char32 *s3, const char32 *s4, const char32 *s5, const char32 *s6) { PUTLEADER MelderFile_write (file, file -> verbose ? U" = " : NULL, z.re, file -> verbose ? U" + " : U" ", z.im, file -> verbose ? U" i " : NULL); } void texpute1 (MelderFile file, int i, const char32 * (*getText) (int), const char32 *s1, const char32 *s2, const char32 *s3, const char32 *s4, const char32 *s5, const char32 *s6) { PUTLEADER MelderFile_write (file, file -> verbose ? U" = <" : U"<", getText (i), file -> verbose ? U"> " : U">"); } void texpute2 (MelderFile file, int i, const char32 * (*getText) (int), const char32 *s1, const char32 *s2, const char32 *s3, const char32 *s4, const char32 *s5, const char32 *s6) { PUTLEADER MelderFile_write (file, file -> verbose ? U" = <" : U"<", getText (i), file -> verbose ? U"> " : U">"); } void texputeb (MelderFile file, bool i, const char32 *s1, const char32 *s2, const char32 *s3, const char32 *s4, const char32 *s5, const char32 *s6) { PUTLEADER MelderFile_write (file, file -> verbose ? U" = " : NULL, i ? U"" : U"", file -> verbose ? U" " : NULL); } void texputeq (MelderFile file, bool i, const char32 *s1, const char32 *s2, const char32 *s3, const char32 *s4, const char32 *s5, const char32 *s6) { PUTLEADER MelderFile_write (file, file -> verbose ? U"? " : NULL, i ? U"" : U"", file -> verbose ? U" " : NULL); } void texputex (MelderFile file, bool i, const char32 *s1, const char32 *s2, const char32 *s3, const char32 *s4, const char32 *s5, const char32 *s6) { PUTLEADER MelderFile_write (file, file -> verbose ? U"? " : NULL, i ? U"" : U"", file -> verbose ? U" " : NULL); } void texputs1 (MelderFile file, const char *s, const char32 *s1, const char32 *s2, const char32 *s3, const char32 *s4, const char32 *s5, const char32 *s6) { PUTLEADER MelderFile_write (file, file -> verbose ? U" = \"" : U"\""); if (s != NULL) { char c; while ((c = *s ++) != '\0') { MelderFile_writeCharacter (file, (char32) (char8) c); if (c == U'\"') MelderFile_writeCharacter (file, (char32) (char8) c); // double any internal quotes } } MelderFile_write (file, file -> verbose ? U"\" " : U"\""); } void texputs2 (MelderFile file, const char *s, const char32 *s1, const char32 *s2, const char32 *s3, const char32 *s4, const char32 *s5, const char32 *s6) { PUTLEADER MelderFile_write (file, file -> verbose ? U" = \"" : U"\""); if (s != NULL) { char c; while ((c = *s ++) != '\0') { MelderFile_writeCharacter (file, (char32) (char8) c); if (c == '\"') MelderFile_writeCharacter (file, (char32) (char8) c); // double any internal quotes } } MelderFile_write (file, file -> verbose ? U"\" " : U"\""); } void texputs4 (MelderFile file, const char *s, const char32 *s1, const char32 *s2, const char32 *s3, const char32 *s4, const char32 *s5, const char32 *s6) { PUTLEADER MelderFile_write (file, file -> verbose ? U" = \"" : U"\""); if (s != NULL) { char c; while ((c = *s ++) != '\0') { MelderFile_writeCharacter (file, (char32) (char8) c); if (c == '\"') MelderFile_writeCharacter (file, (char32) (char8) c); // double any internal quotes } } MelderFile_write (file, file -> verbose ? U"\" " : U"\""); } void texputw2 (MelderFile file, const char32 *s, const char32 *s1, const char32 *s2, const char32 *s3, const char32 *s4, const char32 *s5, const char32 *s6) { PUTLEADER MelderFile_write (file, file -> verbose ? U" = \"" : U"\""); if (s != NULL) { char32 c; while ((c = *s ++) != U'\0') { MelderFile_writeCharacter (file, c); if (c == U'\"') MelderFile_writeCharacter (file, c); // double any internal quotes } } MelderFile_write (file, file -> verbose ? U"\" " : U"\""); } void texputw4 (MelderFile file, const char32 *s, const char32 *s1, const char32 *s2, const char32 *s3, const char32 *s4, const char32 *s5, const char32 *s6) { PUTLEADER MelderFile_write (file, file -> verbose ? U" = \"" : U"\""); if (s != NULL) { char32 c; while ((c = *s ++) != '\0') { MelderFile_writeCharacter (file, c); if (c == '\"') MelderFile_writeCharacter (file, c); // double any internal quotes } } MelderFile_write (file, file -> verbose ? U"\" " : U"\""); } /********** machine-independent binary I/O **********/ /* Optimizations for machines for which some of the formats are native. */ /* On which machines is "short" a two's complement Big-Endian (MSB-first) 2-byte word? */ #if defined (macintosh) && TARGET_RT_BIG_ENDIAN == 1 #define binario_16bitBE 1 #define binario_16bitLE 0 #elif defined (_WIN32) || defined (macintosh) && TARGET_RT_LITTLE_ENDIAN == 1 #define binario_16bitBE 0 #define binario_16bitLE 1 #else #define binario_16bitBE 0 #define binario_16bitLE 0 #endif /* On which machines is "long" a two's complement Big-Endian (MSB-first) 4-byte word? */ #if defined (macintosh) && TARGET_RT_BIG_ENDIAN == 1 #define binario_32bitBE 1 #define binario_32bitLE 0 #elif defined (_WIN32) || defined (macintosh) && TARGET_RT_LITTLE_ENDIAN == 1 #define binario_32bitBE 0 #define binario_32bitLE 1 #else #define binario_32bitBE 0 #define binario_32bitLE 0 #endif /* On which machines is "float" IEEE, four bytes, Most Significant Bit first? */ #if defined (macintosh) && TARGET_RT_BIG_ENDIAN == 1 #define binario_floatIEEE4msb (sizeof (float) == 4) #define binario_floatIEEE4lsb 0 #elif defined (_WIN32) || defined (macintosh) && TARGET_RT_LITTLE_ENDIAN == 1 #define binario_floatIEEE4msb 0 #define binario_floatIEEE4lsb (sizeof (float) == 4) #else #define binario_floatIEEE4msb 0 #define binario_floatIEEE4lsb 0 #endif /* On which machines is "double" IEEE, eight bytes, Most Significant Bit first? */ #if defined (macintosh) && TARGET_RT_BIG_ENDIAN == 1 #define binario_doubleIEEE8msb (sizeof (double) == 8) #define binario_doubleIEEE8lsb 0 #elif defined (_WIN32) || defined (macintosh) && TARGET_RT_LITTLE_ENDIAN == 1 #define binario_doubleIEEE8msb 0 #define binario_doubleIEEE8lsb (sizeof (double) == 8) #else #define binario_doubleIEEE8msb 0 #define binario_doubleIEEE8lsb 0 #endif /* The routines bingetr4, bingetr8, binputr4, and binputr8, were implemented by Paul Boersma from the descriptions of the IEEE floating-point formats, as found in the MC68881/MC68882 User's Manual by Motorola (second edition, 1989). The following copyright notice only refers to the code of bingetr10 and binputr10. */ /* Copyright (C) 1988-1991 Apple Computer, Inc. * All rights reserved. * * Warranty Information * Even though Apple has reviewed this software, Apple makes no warranty * or representation, either express or implied, with respect to this * software, its quality, accuracy, merchantability, or fitness for a * particular purpose. As a result, this software is provided "as is," * and you, its user, are assuming the entire risk as to its quality * and accuracy. * * This code may be used and freely distributed as long as it includes * this copyright notice and the above warranty information. * * Machine-independent I/O routines for IEEE floating-point numbers. * * NaN's and infinities are converted to HUGE_VAL or HUGE, which * happens to be infinity on IEEE machines. Unfortunately, it is * impossible to preserve NaN's in a machine-independent way. * Infinities are, however, preserved on IEEE machines. * * These routines have been tested on the following machines: * Apple Macintosh, MPW 3.1 C compiler * Apple Macintosh, THINK C compiler * Silicon Graphics IRIS, MIPS compiler * Cray X/MP and Y/MP * Digital Equipment VAX * * * Implemented by Malcolm Slaney and Ken Turkowski. * * Malcolm Slaney contributions during 1988-1990 include big- and little- * endian file I/O, conversion to and from Motorola's extended 80-bit * floating-point format, and conversions to and from IEEE single- * precision floating-point format. * * In 1991, Ken Turkowski implemented the conversions to and from * IEEE double-precision format, added more precision to the extended * conversions, and accommodated conversions involving +/- infinity, * NaN's, and denormalized numbers. */ // QUESTION: do the following work correctly if a long is 64 bits? # define FloatToUnsigned(f) \ ((unsigned long)(((long)((f) - 2147483648.0)) + 2147483647L + 1)) # define UnsignedToFloat(u) \ (((double)((long)((u) - 2147483647L - 1))) + 2147483648.0) //#define FloatToUnsigned(f) (uint32_t) (f) //#define UnsignedToFloat(u) (double) (u) /**************************************************************** * Extended precision IEEE floating-point conversion routines. ****************************************************************/ /* * C O N V E R T T O I E E E E X T E N D E D */ /* * C O N V E R T F R O M I E E E E X T E N D E D */ /*************** End of Apple Computer intermezzo. ***************/ static void readError (FILE *f, const char32 *text) { Melder_throw (feof (f) ? U"Reached end of file" : U"Error in file", U" while trying to read ", text); } static void writeError (const char32 *text) { Melder_throw (U"Error in file while trying to write ", text); } unsigned int bingetu1 (FILE *f) { try { int externalValue = getc (f); // either -1 (EOF) or in the range 0..255 if (externalValue < 0) readError (f, U"a byte."); return (unsigned int) externalValue; } catch (MelderError) { Melder_throw (U"Unsigned integer not read from 1 byte in binary file."); } } void binputu1 (unsigned int u, FILE *f) { try { if (putc ((int) u, f) < 0) writeError (U"a byte."); } catch (MelderError) { Melder_throw (U"Unsigned integer not written to 1 byte in binary file."); } } int bingeti1 (FILE *f) { try { int externalValue = getc (f); if (externalValue < 0) readError (f, U"a byte."); return (signed char) externalValue; // this converts e.g. 200 to -56 } catch (MelderError) { Melder_throw (U"Signed integer not read from 1 byte in binary file."); } } void binputi1 (int u, FILE *f) { try { if (putc (u, f) < 0) writeError (U"a byte."); } catch (MelderError) { Melder_throw (U"Signed integer not written to 1 byte in binary file."); } } int bingete1 (FILE *f, int min, int max, const char32 *type) { try { int externalValue = getc (f); if (externalValue < 0) readError (f, U"a byte."); int result = (signed char) externalValue; // this converts e.g. 200 to -56, so the enumerated type is signed if (result < min || result > max) Melder_throw (result, U" is not a value of enumerated type <", type, U">."); return result; } catch (MelderError) { Melder_throw (U"Enumerated type not read from 1 byte in binary file."); } } void binpute1 (int value, FILE *f) { try { if (putc (value, f) < 0) writeError (U"a byte."); } catch (MelderError) { Melder_throw (U"Enumerated type not written to 1 byte in binary file."); } } static int bitsInReadBuffer = 0; static unsigned char readBuffer; #define macro_bingetb(nbits) \ unsigned int bingetb##nbits (FILE *f) { \ if (bitsInReadBuffer < nbits) { \ int externalValue = fgetc (f); \ if (externalValue < 0) readError (f, U"a bit."); \ readBuffer = (unsigned char) externalValue; \ bitsInReadBuffer = 8; \ } \ unsigned char result = (unsigned char) ((uint32_t) readBuffer << (8 - bitsInReadBuffer)); \ bitsInReadBuffer -= nbits; \ return result >> (8 - nbits); \ } macro_bingetb (1) macro_bingetb (2) macro_bingetb (3) macro_bingetb (4) macro_bingetb (5) macro_bingetb (6) macro_bingetb (7) void bingetb (FILE *f) { (void) f; bitsInReadBuffer = 0; } int16 bingeti2 (FILE *f) { try { if (binario_16bitBE && Melder_debug != 18) { int16 s; if (fread (& s, sizeof (int16), 1, f) != 1) readError (f, U"a signed 16-bit integer."); return s; } else { uint8 bytes [2]; if (fread (bytes, sizeof (uint8), 2, f) != 2) readError (f, U"two bytes."); return (int16) // reinterpret sign bit ((uint16) ((uint16) bytes [0] << 8) | (uint16) bytes [1]); } } catch (MelderError) { Melder_throw (U"Signed integer not read from 2 bytes in binary file."); } } int16 bingeti2LE (FILE *f) { try { if (binario_16bitLE && Melder_debug != 18) { int16 s; if (fread (& s, sizeof (int16), 1, f) != 1) readError (f, U"a signed 16-bit integer."); return s; } else { uint8 bytes [2]; if (fread (bytes, sizeof (uint8), 2, f) != 2) readError (f, U"two bytes."); return (int16) // reinterpret sign bit ((uint16) ((uint16) bytes [1] << 8) | (uint16) bytes [0]); } } catch (MelderError) { Melder_throw (U"Signed integer not read from 2 bytes in binary file."); } } uint16_t bingetu2 (FILE *f) { try { if (binario_16bitBE && Melder_debug != 18) { uint16 s; if (fread (& s, sizeof (uint16), 1, f) != 1) readError (f, U"an unsigned 16-bit integer."); return s; // without sign extension } else { uint8 bytes [2]; if (fread (bytes, sizeof (uint8), 2, f) != 2) readError (f, U"two bytes."); return (uint16) ((uint16) bytes [0] << 8) | (uint16) bytes [1]; } } catch (MelderError) { Melder_throw (U"Unsigned integer not read from 2 bytes in binary file."); } } uint16_t bingetu2LE (FILE *f) { try { if (binario_16bitLE && Melder_debug != 18) { uint16 s; if (fread (& s, sizeof (uint16), 1, f) != 1) readError (f, U"an unsigned 16-bit integer."); return s; // without sign extension } else { uint8 bytes [2]; if (fread (bytes, sizeof (uint8), 2, f) != 2) readError (f, U"two bytes."); return (uint16) ((uint16) bytes [1] << 8) | (uint16) bytes [0]; } } catch (MelderError) { Melder_throw (U"Unsigned integer not read from 2 bytes in binary file."); } } int bingete2 (FILE *f, int min, int max, const char32 *type) { try { int16 result; if (binario_16bitBE && Melder_debug != 18) { if (fread (& result, sizeof (int16), 1, f) != 1) readError (f, U"a signed 16-bit integer."); } else { uint8 bytes [2]; if (fread (bytes, sizeof (uint8), 2, f) != 2) readError (f, U"two bytes."); uint16 externalValue = (uint16) ((uint16) bytes [0] << 8) | (uint16) bytes [1]; result = (int16) externalValue; } if (result < min || result > max) Melder_throw (result, U" is not a value of enumerated type \"", type, U"\"."); return (int) result; } catch (MelderError) { Melder_throw (U"Enumerated value not read from 2 bytes in binary file."); } } int32 bingeti3 (FILE *f) { try { uint8 bytes [3]; if (fread (bytes, sizeof (uint8), 3, f) != 3) readError (f, U"three bytes."); uint32 externalValue = (uint32) ((uint32) bytes [0] << 16) | (uint32) ((uint32) bytes [1] << 8) | (uint32) bytes [2]; if ((bytes [0] & 128) != 0) // is the 24-bit sign bit on? externalValue |= 0xFF000000; // extend negative sign to 32 bits return (int32_t) externalValue; // reinterpret sign bit } catch (MelderError) { Melder_throw (U"Signed integer not read from 3 bytes in binary file."); } } int32 bingeti3LE (FILE *f) { try { uint8 bytes [3]; if (fread (bytes, sizeof (uint8), 3, f) != 3) readError (f, U"three bytes."); uint32 externalValue = (uint32) ((uint32) bytes [2] << 16) | (uint32) ((uint32) bytes [1] << 8) | (uint32) bytes [0]; if ((bytes [2] & 128) != 0) // is the 24-bit sign bit on? externalValue |= 0xFF000000; // extend negative sign to 32 bits return (int32_t) externalValue; // reinterpret sign bit } catch (MelderError) { Melder_throw (U"Signed integer not read from 3 bytes in binary file."); } } int32 bingeti4 (FILE *f) { try { if (binario_32bitBE && Melder_debug != 18) { int32 l; if (fread (& l, sizeof (int32), 1, f) != 1) readError (f, U"a signed 32-bit integer."); return l; } else { uint8 bytes [4]; if (fread (bytes, sizeof (uint8), 4, f) != 4) readError (f, U"four bytes."); return (int32) ((uint32) ((uint32) bytes [0] << 24) | (uint32) ((uint32) bytes [1] << 16) | (uint32) ((uint32) bytes [2] << 8) | (uint32) bytes [3]); } } catch (MelderError) { Melder_throw (U"Signed integer not read from 4 bytes in binary file."); } } int32 bingeti4LE (FILE *f) { try { if (binario_32bitLE && Melder_debug != 18) { int32 l; if (fread (& l, sizeof (int32), 1, f) != 1) readError (f, U"a signed 32-bit integer."); return l; } else { uint8 bytes [4]; if (fread (bytes, sizeof (uint8), 4, f) != 4) readError (f, U"four bytes."); return (int32) // reinterpret sign bit ((uint32) ((uint32) bytes [3] << 24) | (uint32) ((uint32) bytes [2] << 16) | (uint32) ((uint32) bytes [1] << 8) | (uint32) bytes [0]); } } catch (MelderError) { Melder_throw (U"Signed integer not read from 4 bytes in binary file."); } } uint32 bingetu4 (FILE *f) { try { if (binario_32bitBE && Melder_debug != 18) { uint32_t l; if (fread (& l, sizeof (uint32_t), 1, f) != 1) readError (f, U"an unsigned 32-bit integer."); return l; } else { uint8 bytes [4]; if (fread (bytes, sizeof (uint8), 4, f) != 4) readError (f, U"four bytes."); return (uint32) ((uint32) bytes [0] << 24) | (uint32) ((uint32) bytes [1] << 16) | (uint32) ((uint32) bytes [2] << 8) | (uint32) bytes [3]; } } catch (MelderError) { Melder_throw (U"Unsigned integer not read from 4 bytes in binary file."); } } uint32 bingetu4LE (FILE *f) { try { if (binario_32bitLE && Melder_debug != 18) { uint32 l; if (fread (& l, sizeof (uint32), 1, f) != 1) readError (f, U"an unsigned 32-bit integer."); return l; } else { uint8 bytes [4]; if (fread (bytes, sizeof (uint8), 4, f) != 4) readError (f, U"four bytes."); return (uint32) ((uint32) bytes [3] << 24) | (uint32) ((uint32) bytes [2] << 16) | (uint32) ((uint32) bytes [1] << 8) | (uint32) bytes [0]; } } catch (MelderError) { Melder_throw (U"Unsigned integer not read from 4 bytes in binary file."); } } double bingetr4 (FILE *f) { try { if (binario_floatIEEE4msb && Melder_debug != 18) { float x; if (fread (& x, sizeof (float), 1, f) != 1) readError (f, U"a 32-bit floating-point number."); return x; } else { uint8 bytes [4]; if (fread (bytes, sizeof (uint8), 4, f) != 4) readError (f, U"four bytes."); int32 exponent = (int32) ((uint32) ((uint32) ((uint32) bytes [0] & 0x0000007F) << 1) | (uint32) ((uint32) ((uint32) bytes [1] & 0x00000080) >> 7)); // between 0 and 255 (it's signed because we're going to subtract something) uint32 mantissa = (uint32) ((uint32) ((uint32) bytes [1] & 0x0000007F) << 16) | (uint32) ((uint32) bytes [2] << 8) | (uint32) bytes [3]; double x; if (exponent == 0) if (mantissa == 0) x = 0.0; else x = ldexp ((double) mantissa, exponent - 149); // denormalized else if (exponent == 0x000000FF) // Infinity or Not-a-Number x = HUGE_VAL; else // finite x = ldexp ((double) (mantissa | 0x00800000), exponent - 150); return bytes [0] & 0x80 ? - x : x; } } catch (MelderError) { Melder_throw (U"Floating-point number not read from 4 bytes in binary file."); } } double bingetr4LE (FILE *f) { try { if (binario_floatIEEE4lsb && Melder_debug != 18) { float x; if (fread (& x, sizeof (float), 1, f) != 1) readError (f, U"a 32-bit floating-point number."); return x; } else { uint8 bytes [4]; if (fread (bytes, sizeof (uint8), 4, f) != 4) readError (f, U"four bytes."); int32 exponent = (int32) ((uint32) ((uint32) ((uint32) bytes [3] & 0x0000007F) << 1) | (uint32) ((uint32) ((uint32) bytes [2] & 0x00000080) >> 7)); uint32 mantissa = (uint32) ((uint32) ((uint32) bytes [2] & 0x0000007F) << 16) | (uint32) ((uint32) bytes [1] << 8) | (uint32) bytes [0]; double x; if (exponent == 0) if (mantissa == 0) x = 0.0; else x = ldexp ((double) mantissa, exponent - 149); // denormalized else if (exponent == 0x000000FF) // Infinity or Not-a-Number x = HUGE_VAL; else // finite x = ldexp ((double) (mantissa | 0x00800000), exponent - 150); return bytes [3] & 0x80 ? - x : x; } } catch (MelderError) { Melder_throw (U"Floating-point number not read from 4 bytes in binary file."); } } double bingetr8 (FILE *f) { try { if (binario_doubleIEEE8msb && Melder_debug != 18) { double x; if (fread (& x, sizeof (double), 1, f) != 1) readError (f, U"a 64-bit floating-point number."); return x; } else { uint8 bytes [8]; if (fread (bytes, sizeof (uint8), 8, f) != 8) readError (f, U"eight bytes."); int32 exponent = (int32) ((uint32) ((uint32) ((uint32) bytes [0] & 0x0000007F) << 4) | (uint32) ((uint32) ((uint32) bytes [1] & 0x000000F0) >> 4)); uint32_t highMantissa = (uint32) ((uint32) ((uint32) bytes [1] & 0x0000000F) << 16) | (uint32) ((uint32) bytes [2] << 8) | (uint32) bytes [3]; uint32_t lowMantissa = (uint32) ((uint32) bytes [4] << 24) | (uint32) ((uint32) bytes [5] << 16) | (uint32) ((uint32) bytes [6] << 8) | (uint32) bytes [7]; double x; if (exponent == 0) if (highMantissa == 0 && lowMantissa == 0) x = 0.0; else x = ldexp ((double) highMantissa, exponent - 1042) + ldexp ((double) lowMantissa, exponent - 1074); // denormalized else if (exponent == 0x000007FF) // Infinity or Not-a-Number x = HUGE_VAL; else x = ldexp ((double) (highMantissa | 0x00100000), exponent - 1043) + ldexp ((double) lowMantissa, exponent - 1075); return bytes [0] & 0x80 ? - x : x; } } catch (MelderError) { Melder_throw (U"Floating-point number not read from 8 bytes in binary file."); } } double bingetr10 (FILE *f) { try { uint8 bytes [10]; if (fread (bytes, sizeof (uint8), 10, f) != 10) readError (f, U"ten bytes."); int32 exponent = (int32) ((uint32) ((uint32) ((uint32) bytes [0] & 0x0000007F) << 8) | (uint32) bytes [1]); // between 0 and 32767 uint32 highMantissa = (uint32) ((uint32) bytes [2] << 24) | (uint32) ((uint32) bytes [3] << 16) | (uint32) ((uint32) bytes [4] << 8) | (uint32) bytes [5]; uint32 lowMantissa = (uint32) ((uint32) bytes [6] << 24) | (uint32) ((uint32) bytes [7] << 16) | (uint32) ((uint32) bytes [8] << 8) | (uint32) bytes [9]; double x; if (exponent == 0 && highMantissa == 0 && lowMantissa == 0) x = 0; else if (exponent == 0x00007FFF) x = HUGE_VAL; // Infinity or NaN else { exponent -= 16383; // between -16382 and +16383 x = ldexp ((double) highMantissa, exponent - 31); x += ldexp ((double) lowMantissa, exponent - 63); } return bytes [0] & 0x80 ? - x : x; } catch (MelderError) { Melder_throw (U"Floating-point number not read from 10 bytes in binary file."); } } static int bitsInWriteBuffer = 0; static unsigned char writeBuffer = 0; #define macro_binputb(nbits) \ void binputb##nbits (unsigned int value, FILE *f) { \ if (bitsInWriteBuffer + nbits > 8) { \ if (fputc (writeBuffer, f) < 0) writeError (U"a bit."); \ bitsInWriteBuffer = 0; \ writeBuffer = 0; \ } \ writeBuffer |= (value << (8 - nbits)) >> bitsInWriteBuffer; \ bitsInWriteBuffer += nbits; \ } macro_binputb (1) macro_binputb (2) macro_binputb (3) macro_binputb (4) macro_binputb (5) macro_binputb (6) macro_binputb (7) void binputb (FILE *f) { if (bitsInWriteBuffer == 0) return; if (fputc (writeBuffer, f) < 0) writeError (U"a bit."); // flush bitsInWriteBuffer = 0; writeBuffer = 0; } void binputi2 (int16 i, FILE *f) { try { if (binario_16bitBE && Melder_debug != 18) { if (fwrite (& i, sizeof (short), 1, f) != 1) writeError (U"a signed 16-bit integer."); } else { char bytes [2]; bytes [0] = (char) (i >> 8); // truncate bytes [1] = (char) i; // truncate if (fwrite (bytes, sizeof (char), 2, f) != 2) writeError (U"two bytes."); } } catch (MelderError) { Melder_throw (U"Signed integer not written to 2 bytes in binary file."); } } void binputi2LE (int16 i, FILE *f) { try { if (binario_16bitLE && Melder_debug != 18) { if (fwrite (& i, sizeof (short), 1, f) != 1) writeError (U"a signed 16-bit integer."); } else { char bytes [2]; bytes [1] = (char) (i >> 8); // truncate bytes [0] = (char) i; // truncate if (fwrite (bytes, sizeof (char), 2, f) != 2) writeError (U"two bytes."); } } catch (MelderError) { Melder_throw (U"Signed integer not written to 2 bytes in binary file."); } } void binputu2 (uint16 u, FILE *f) { try { if (binario_16bitBE && Melder_debug != 18) { if (fwrite (& u, sizeof (uint16_t), 1, f) != 1) writeError (U"an unsigned 16-bit integer."); } else { char bytes [2]; bytes [0] = (char) (u >> 8); // truncate bytes [1] = (char) u; // truncate if (fwrite (bytes, sizeof (char), 2, f) != 2) writeError (U"two bytes."); } } catch (MelderError) { Melder_throw (U"Unsigned integer not written to 2 bytes in binary file."); } } void binputu2LE (uint16 u, FILE *f) { try { if (binario_16bitLE && Melder_debug != 18) { if (fwrite (& u, sizeof (uint16_t), 1, f) != 1) writeError (U"an unsigned 16-bit integer."); } else { char bytes [2]; bytes [1] = (char) (u >> 8); // truncate bytes [0] = (char) u; // truncate if (fwrite (bytes, sizeof (char), 2, f) != 2) writeError (U"two bytes."); } } catch (MelderError) { Melder_throw (U"Unsigned integer not written to 2 bytes in binary file."); } } void binpute2 (int value, FILE *f) { try { if (binario_16bitBE && Melder_debug != 18) { short s = value; if (fwrite (& s, sizeof (short), 1, f) != 1) writeError (U"a signed 16-bit integer"); } else { char bytes [2]; bytes [0] = (char) (value >> 8); // truncate bytes [1] = (char) value; // truncate if (fwrite (bytes, sizeof (char), 2, f) != 2) writeError (U"two bytes."); } } catch (MelderError) { Melder_throw (U"Enumerated value not written to 2 bytes in binary file."); } } void binputi3 (int32 i, FILE *f) { try { char bytes [3]; bytes [0] = (char) (i >> 16); // truncate bytes [1] = (char) (i >> 8); // truncate bytes [2] = (char) i; // truncate if (fwrite (bytes, sizeof (char), 3, f) != 3) writeError (U"three bytes"); } catch (MelderError) { Melder_throw (U"Signed integer not written to 3 bytes in binary file."); } } void binputi3LE (int32 i, FILE *f) { try { char bytes [3]; bytes [2] = (char) (i >> 16); // truncate bytes [1] = (char) (i >> 8); // truncate bytes [0] = (char) i; // truncate if (fwrite (bytes, sizeof (char), 3, f) != 3) writeError (U"three bytes"); } catch (MelderError) { Melder_throw (U"Signed integer not written to 3 bytes in binary file."); } } void binputi4 (int32 i, FILE *f) { try { if (binario_32bitBE && Melder_debug != 18) { if (fwrite (& i, sizeof (int32), 1, f) != 1) writeError (U"a signed 32-bit integer."); } else { char bytes [4]; bytes [0] = (char) (i >> 24); // truncate bytes [1] = (char) (i >> 16); // truncate bytes [2] = (char) (i >> 8); // truncate bytes [3] = (char) i; // truncate if (fwrite (bytes, sizeof (char), 4, f) != 4) writeError (U"four bytes."); } } catch (MelderError) { Melder_throw (U"Signed integer not written to 4 bytes in binary file."); } } void binputi4LE (int32 i, FILE *f) { try { if (binario_32bitLE && Melder_debug != 18) { if (fwrite (& i, sizeof (int32), 1, f) != 1) writeError (U"a signed 32-bit integer."); } else { char bytes [4]; bytes [3] = (char) (i >> 24); // truncate bytes [2] = (char) (i >> 16); // truncate bytes [1] = (char) (i >> 8); // truncate bytes [0] = (char) i; // truncate if (fwrite (bytes, sizeof (char), 4, f) != 4) writeError (U"four bytes."); } } catch (MelderError) { Melder_throw (U"Signed integer not written to 4 bytes in binary file."); } } void binputu4 (uint32 u, FILE *f) { try { if (binario_32bitBE && Melder_debug != 18) { if (fwrite (& u, sizeof (uint32), 1, f) != 1) writeError (U"an unsigned 32-bit integer."); } else { char bytes [4]; bytes [0] = (char) (u >> 24); // truncate bytes [1] = (char) (u >> 16); // truncate bytes [2] = (char) (u >> 8); // truncate bytes [3] = (char) u; // truncate if (fwrite (bytes, sizeof (char), 4, f) != 4) writeError (U"four bytes."); } } catch (MelderError) { Melder_throw (U"Unsigned integer not written to 4 bytes in binary file."); } } void binputu4LE (uint32 u, FILE *f) { try { if (binario_32bitLE && Melder_debug != 18) { if (fwrite (& u, sizeof (uint32), 1, f) != 1) writeError (U"an unsigned 32-bit integer."); } else { char bytes [4]; bytes [3] = (char) (u >> 24); // truncate bytes [2] = (char) (u >> 16); // truncate bytes [1] = (char) (u >> 8); // truncate bytes [0] = (char) u; // truncate if (fwrite (bytes, sizeof (char), 4, f) != 4) writeError (U"four bytes."); } } catch (MelderError) { Melder_throw (U"Unsigned integer not written to 4 bytes in binary file."); } } void binputr4 (double x, FILE *f) { try { if (binario_floatIEEE4msb && Melder_debug != 18) { float x4 = (float) x; // convert down, with loss of precision if (fwrite (& x4, sizeof (float), 1, f) != 1) writeError (U"a 32-bit floating-point number."); } else { uint8 bytes [4]; int sign, exponent; double fMantissa, fsMantissa; uint32 mantissa; if (x < 0.0) { sign = 0x0100; x *= -1; } else sign = 0; if (x == 0.0) { exponent = 0; mantissa = 0; } else { fMantissa = frexp (x, & exponent); if ((exponent > 128) || ! (fMantissa < 1)) // Infinity or Not-a-Number { exponent = sign | 0x00FF; mantissa = 0; } // Infinity else { // finite exponent += 126; // add bias if (exponent <= 0) { // denormalized fMantissa = ldexp (fMantissa, exponent - 1); exponent = 0; } exponent |= sign; fMantissa = ldexp (fMantissa, 24); fsMantissa = floor (fMantissa); mantissa = (uint32_t) fsMantissa & 0x007FFFFF; } } bytes [0] = (uint8) (exponent >> 1); // truncate: bits 2 through 9 (bit 9 is the sign bit) bytes [1] = (uint8) ((exponent << 7) | (mantissa >> 16)); // truncate bytes [2] = (uint8) (mantissa >> 8); // truncate bytes [3] = (uint8) mantissa; // truncate if (fwrite (bytes, sizeof (uint8), 4, f) != 4) writeError (U"four bytes."); } } catch (MelderError) { Melder_throw (U"Floating-point number not written to 4 bytes in binary file."); } } void binputr4LE (double x, FILE *f) { try { if (binario_floatIEEE4lsb && Melder_debug != 18) { float x4 = (float) x; // convert down, with loss of precision if (fwrite (& x4, sizeof (float), 1, f) != 1) writeError (U"a 32-bit floating-point number."); } else { uint8 bytes [4]; int sign, exponent; double fMantissa, fsMantissa; uint32 mantissa; if (x < 0.0) { sign = 0x0100; x *= -1; } else sign = 0; if (x == 0.0) { exponent = 0; mantissa = 0; } else { fMantissa = frexp (x, & exponent); if ((exponent > 128) || ! (fMantissa < 1)) // Infinity or Not-a-Number { exponent = sign | 0x00FF; mantissa = 0; } // Infinity else { // finite exponent += 126; // add bias if (exponent <= 0) { // denormalized fMantissa = ldexp (fMantissa, exponent - 1); exponent = 0; } exponent |= sign; fMantissa = ldexp (fMantissa, 24); fsMantissa = floor (fMantissa); mantissa = (uint32_t) fsMantissa & 0x007FFFFF; } } bytes [3] = (uint8) (exponent >> 1); bytes [2] = (uint8) ((exponent << 7) | (mantissa >> 16)); bytes [1] = (uint8) (mantissa >> 8); bytes [0] = (uint8) mantissa; if (fwrite (bytes, sizeof (uint8), 4, f) != 4) writeError (U"four bytes."); } } catch (MelderError) { Melder_throw (U"Floating-point number not written to 4 bytes in binary file."); } } void binputr8 (double x, FILE *f) { try { if (binario_doubleIEEE8msb && Melder_debug != 18) { if (fwrite (& x, sizeof (double), 1, f) != 1) writeError (U"a 64-bit floating-point number."); } else { uint8 bytes [8]; int sign, exponent; double fMantissa, fsMantissa; uint32 highMantissa, lowMantissa; if (x < 0.0) { sign = 0x0800; x *= -1; } else sign = 0; if (x == 0.0) { exponent = 0; highMantissa = 0; lowMantissa = 0; } else { fMantissa = frexp (x, & exponent); if ((exponent > 1024) || ! (fMantissa < 1)) // Infinity or Not-a-Number { exponent = sign | 0x07FF; highMantissa = 0; lowMantissa = 0; } // Infinity else { // finite exponent += 1022; // add bias if (exponent <= 0) { // denormalized fMantissa = ldexp (fMantissa, exponent - 1); exponent = 0; } exponent |= sign; fMantissa = ldexp (fMantissa, 21); fsMantissa = floor (fMantissa); highMantissa = (uint32_t) fsMantissa & 0x000FFFFF; fMantissa = ldexp (fMantissa - fsMantissa, 32); fsMantissa = floor (fMantissa); lowMantissa = (uint32_t) fsMantissa; } } bytes [0] = (uint8) (exponent >> 4); bytes [1] = (uint8) ((exponent << 4) | (highMantissa >> 16)); bytes [2] = (uint8) (highMantissa >> 8); bytes [3] = (uint8) highMantissa; bytes [4] = (uint8) (lowMantissa >> 24); bytes [5] = (uint8) (lowMantissa >> 16); bytes [6] = (uint8) (lowMantissa >> 8); bytes [7] = (uint8) lowMantissa; if (fwrite (bytes, sizeof (uint8), 8, f) != 8) writeError (U"eight bytes."); } } catch (MelderError) { Melder_throw (U"Floating-point number not written to 8 bytes in binary file."); } } void binputr10 (double x, FILE *f) { try { unsigned char bytes [10]; Melder_assert (sizeof (int) > 2); int sign, exponent; // these should be uint16_t, but frexp() expects an int double fMantissa, fsMantissa; uint32_t highMantissa, lowMantissa; if (x < 0.0) { sign = 0x8000; x *= -1; } else sign = 0; if (x == 0.0) { exponent = 0; highMantissa = 0; lowMantissa = 0; } else { fMantissa = frexp (x, & exponent); if ((exponent > 16384) || ! (fMantissa < 1)) // Infinity or Not-a-Number { exponent = sign | 0x7FFF; highMantissa = 0; lowMantissa = 0; } // Infinity else { // finite exponent += 16382; // add bias if (exponent < 0) { // denormalized fMantissa = ldexp (fMantissa, exponent); exponent = 0; } exponent |= sign; fMantissa = ldexp (fMantissa, 32); fsMantissa = floor (fMantissa); highMantissa = (uint32_t) fsMantissa; fMantissa = ldexp (fMantissa - fsMantissa, 32); fsMantissa = floor (fMantissa); lowMantissa = (uint32_t) fsMantissa; } } bytes [0] = (uint8) (exponent >> 8); bytes [1] = (uint8) exponent; bytes [2] = (uint8) (highMantissa >> 24); bytes [3] = (uint8) (highMantissa >> 16); bytes [4] = (uint8) (highMantissa >> 8); bytes [5] = (uint8) highMantissa; bytes [6] = (uint8) (lowMantissa >> 24); bytes [7] = (uint8) (lowMantissa >> 16); bytes [8] = (uint8) (lowMantissa >> 8); bytes [9] = (uint8) lowMantissa; if (fwrite (bytes, sizeof (uint8), 10, f) != 10) writeError (U"ten bytes."); } catch (MelderError) { Melder_throw (U"Floating-point number not written to 10 bytes in binary file."); } } fcomplex bingetc8 (FILE *f) { try { fcomplex result; result. re = bingetr4 (f); result. im = bingetr4 (f); return result; } catch (MelderError) { Melder_throw (U"Complex number not read from 8 bytes in binary file."); fcomplex result = { 0 }; return result; } } dcomplex bingetc16 (FILE *f) { try { dcomplex result; result. re = bingetr8 (f); result. im = bingetr8 (f); return result; } catch (MelderError) { Melder_throw (U"Complex number not read from 16 bytes in binary file."); dcomplex result = { 0 }; return result; } } void binputc8 (fcomplex z, FILE *f) { try { binputr4 (z. re, f); binputr4 (z. im, f); } catch (MelderError) { Melder_throw (U"Complex number not written to 8 bytes in binary file."); } } void binputc16 (dcomplex z, FILE *f) { try { binputr8 (z. re, f); binputr8 (z. im, f); } catch (MelderError) { Melder_throw (U"Complex number not written to 16 bytes in binary file."); } } char * bingets1 (FILE *f) { try { unsigned int length = bingetu1 (f); autostring8 result = Melder_malloc (char, length + 1); if (fread (result.peek(), sizeof (char), length, f) != length) Melder_throw (feof (f) ? U"Reached end of file" : U"Error in file", U" while trying to read ", length, U" one-byte characters."); result [length] = 0; // trailing null byte return result.transfer(); } catch (MelderError) { Melder_throw (U"Text not read from a binary file."); } } char * bingets2 (FILE *f) { try { uint16_t length = bingetu2 (f); autostring8 result = Melder_malloc (char, (int64_t) length + 1); if (fread (result.peek(), sizeof (char), length, f) != length) Melder_throw (feof (f) ? U"Reached end of file" : U"Error in file", U" while trying to read ", length, U" one-byte characters."); result [length] = 0; // trailing null byte return result.transfer(); } catch (MelderError) { Melder_throw (U"Text not read from a binary file."); } } char * bingets4 (FILE *f) { try { uint32_t length = bingetu4 (f); autostring8 result = Melder_malloc (char, (int64_t) length + 1); if (fread (result.peek(), sizeof (char), length, f) != length) Melder_throw (feof (f) ? U"Reached end of file" : U"Error in file", U" while trying to read ", length, U" one-byte characters."); result [length] = 0; // trailing null byte return result.transfer(); } catch (MelderError) { Melder_throw (U"Text not read from a binary file."); } } char32 * bingetw1 (FILE *f) { try { autostring32 result = NULL; unsigned short length = bingetu1 (f); if (length == 0xFF) { // an escape for encoding /* * UTF-16 */ length = bingetu1 (f); result.reset (Melder_malloc (char32, (int64_t) length + 1)); for (unsigned short i = 0; i < length; i ++) { uint16 kar = bingetu2 (f); if ((kar & 0xF800) == 0xD800) { if (kar > 0xDBFF) Melder_throw (U"Incorrect Unicode value (first surrogate member ", kar, U")."); uint16 kar2 = bingetu2 (f); if (kar2 < 0xDC00 || kar2 > 0xDFFF) Melder_throw (U"Incorrect Unicode value (second surrogate member ", kar2, U")."); result [i] = (((kar & 0x3FF) << 10) | (kar2 & 0x3FF)) + 0x10000; } else { result [i] = kar; } } } else { /* * ASCII */ result.reset (Melder_malloc (char32, (int64_t) length + 1)); for (unsigned short i = 0; i < length; i ++) { result [i] = bingetu1 (f); } } result [length] = U'\0'; return result.transfer(); } catch (MelderError) { Melder_throw (U"Text not read from a binary file."); } } char32 * bingetw2 (FILE *f) { try { autostring32 result = NULL; uint16 length = bingetu2 (f); if (length == 0xFFFF) { // an escape for encoding /* * UTF-16 */ length = bingetu2 (f); result.reset (Melder_malloc (char32, (int64_t) length + 1)); for (uint16 i = 0; i < length; i ++) { char32 kar = (char32) (char16) bingetu2 (f); if ((kar & 0x00F800) == 0x00D800) { if (kar > 0x00DBFF) Melder_throw (U"Incorrect Unicode value (first surrogate member ", kar, U")."); char32 kar2 = (char32) (char16) bingetu2 (f); if (kar2 < 0x00DC00 || kar2 > 0x00DFFF) Melder_throw (U"Incorrect Unicode value (second surrogate member ", kar2, U")."); result [i] = (((kar & 0x0003FF) << 10) | (kar2 & 0x0003FF)) + 0x010000; } else { result [i] = kar; } } } else { /* * ASCII */ result.reset (Melder_malloc (char32, length + 1)); for (unsigned short i = 0; i < length; i ++) { result [i] = (char32) (char8) bingetu1 (f); } } result [length] = U'\0'; return result.transfer(); } catch (MelderError) { Melder_throw (U"Text not read from a binary file."); } } char32 * bingetw4 (FILE *f) { try { autostring32 result = NULL; uint32_t length = bingetu4 (f); if (length == 0xFFFFFFFF) { // an escape for encoding /* * UTF-16 */ length = bingetu4 (f); result.reset (Melder_malloc (char32, (int64_t) length + 1)); for (uint32_t i = 0; i < length; i ++) { uint16_t kar = bingetu2 (f); if ((kar & 0xF800) == 0xD800) { if (kar > 0xDBFF) Melder_throw (U"Incorrect Unicode value (first surrogate member ", kar, U")."); uint16_t kar2 = bingetu2 (f); if (kar2 < 0xDC00 || kar2 > 0xDFFF) Melder_throw (U"Incorrect Unicode value (second surrogate member ", kar2, U")."); result [i] = (((kar & 0x3FF) << 10) | (kar2 & 0x3FF)) + 0x10000; } else { result [i] = kar; } } } else { /* * ASCII */ result.reset (Melder_malloc (char32, (int64_t) length + 1)); for (uint32_t i = 0; i < length; i ++) { result [i] = bingetu1 (f); } } result [length] = L'\0'; return result.transfer(); } catch (MelderError) { Melder_throw (U"Text not read from a binary file."); } } void binputs1 (const char *s, FILE *f) { try { if (s == NULL) { binputu1 (0, f); } else { size_t length = strlen (s); if (length > UINT8_MAX) { Melder_warning (U"Text of ", length, U" characters truncated to 255 characters."); length = UINT8_MAX; } binputu1 (length, f); if (fwrite (s, sizeof (char), length, f) != length) Melder_throw (U"Error in file while trying to write ", length, U" one-byte characters."); } } catch (MelderError) { Melder_throw (U"Text not written to a binary file."); } } void binputs2 (const char *s, FILE *f) { try { if (s == NULL) { binputu2 (0, f); } else { size_t length = strlen (s); if (length > UINT16_MAX) { Melder_warning (U"Text of ", length, U" characters truncated to 65535 characters."); length = UINT16_MAX; } binputu2 ((uint16_t) length, f); // safe conversion down if (fwrite (s, sizeof (char), length, f) != length) Melder_throw (U"Error in file while trying to write ", length, U" one-byte characters."); } } catch (MelderError) { Melder_throw (U"Text not written to a binary file."); } } void binputs4 (const char *s, FILE *f) { try { if (s == NULL) { binputu4 (0, f); } else { size_t length = strlen (s); if (length > UINT32_MAX) { Melder_warning (U"Text of ", length, U" characters truncated to 4,294,967,295 characters."); length = UINT32_MAX; } binputu4 (length, f); if (fwrite (s, sizeof (char), length, f) != length) Melder_throw (U"Error in file while trying to write ", length, U" one-byte characters."); } } catch (MelderError) { Melder_throw (U"Text not written to a binary file."); } } static inline void binpututf16 (char32 kar, FILE *f) { if (kar <= 0x00FFFF) { binputu2 ((uint16) kar, f); // truncate to lower 16 bits } else if (kar <= 0x10FFFF) { kar -= 0x010000; binputu2 ((uint16) (0x00D800 | (kar >> 10)), f); binputu2 ((uint16) (0x00DC00 | (kar & 0x0003FF)), f); } else { Melder_fatal (U"Impossible Unicode value."); } } void binputw1 (const char32 *s, FILE *f) { try { if (s == NULL) { binputu1 (0, f); } else { uint32 length = str32len (s); if (length > UINT8_MAX - 1) { Melder_warning (U"Text of ", length, U" characters truncated to 254 characters."); length = UINT8_MAX - 1; } if (Melder_isValidAscii (s)) { /* * ASCII */ binputu1 (length, f); for (size_t i = 0; i < length; i ++) { binputu1 ((unsigned int) (char) s [i], f); } } else { /* * UTF-16 */ binputu1 (0xFF, f); // an escape for multibyte encoding binputu1 (length, f); for (size_t i = 0; i < length; i ++) { binpututf16 (s [i], f); } } } } catch (MelderError) { Melder_throw (U"Text not written to a binary file."); } } void binputw2 (const char32 *s, FILE *f) { try { if (s == NULL) { binputu2 (0, f); } else { int64 length = str32len (s); if (length > UINT16_MAX - 1) { Melder_warning (U"Text of ", length, U" characters truncated to 65534 characters."); length = UINT16_MAX - 1; } if (Melder_isValidAscii (s)) { /* * ASCII */ binputu2 ((uint16) length, f); for (int64 i = 0; i < length; i ++) { binputu1 ((unsigned int) (char8) s [i], f); } } else { /* * UTF-16 */ binputu2 (0xFFFF, f); // an escape for multibyte encoding binputu2 ((uint16) length, f); for (int64 i = 0; i < length; i ++) { binpututf16 (s [i], f); } } } } catch (MelderError) { Melder_throw (U"Text not written to a binary file."); } } void binputw4 (const char32 *s, FILE *f) { try { if (s == NULL) { binputu4 (0, f); } else { int64 length = str32len (s); if (length > UINT32_MAX - 1) { Melder_warning (U"Text of ", length, U" characters truncated to 4,294,967,294 characters."); length = UINT32_MAX - 1; } if (Melder_isValidAscii (s)) { /* * ASCII */ binputu4 ((uint32) length, f); for (int64 i = 0; i < length; i ++) { binputu1 ((unsigned int) (char) s [i], f); } } else { /* * UTF-16 */ binputu4 (0xFFFFFFFF, f); // an escape for multibyte encoding binputu4 ((uint32) length, f); for (int64 i = 0; i < length; i ++) { binpututf16 (s [i], f); } } } } catch (MelderError) { Melder_throw (U"Text not written to a binary file."); } } /* End of file abcio.cpp */