libs/libsndfile/src/paf.c

   1 /*
   2 ** Copyright (C) 1999-2005 Erik de Castro Lopo <erikd@mega-nerd.com>
   3 **
   4 ** This program is free software; you can redistribute it and/or modify
   5 ** it under the terms of the GNU Lesser General Public License as published by
   6 ** the Free Software Foundation; either version 2.1 of the License, or
   7 ** (at your option) any later version.
   8 **
   9 ** This program is distributed in the hope that it will be useful,
  10 ** but WITHOUT ANY WARRANTY; without even the implied warranty of
  11 ** MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  12 ** GNU Lesser General Public License for more details.
  13 **
  14 ** You should have received a copy of the GNU Lesser General Public License
  15 ** along with this program; if not, write to the Free Software
  16 ** Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA.
  17 */
  18
  19 #include "sfconfig.h"
  20
  21 #include <stdio.h>
  22 #include <stdlib.h>
  23 #include <fcntl.h>
  24 #include <string.h>
  25 #include <ctype.h>
  26
  27 #include "sndfile.h"
  28 #include "sfendian.h"
  29 #include "float_cast.h"
  30 #include "common.h"
  31
  32 /*------------------------------------------------------------------------------
  33 ** Macros to handle big/little endian issues.
  34 */
  35
  36 #define FAP_MARKER      (MAKE_MARKER ('f', 'a', 'p', ' '))
  37 #define PAF_MARKER      (MAKE_MARKER (' ', 'p', 'a', 'f'))
  38
  39 /*------------------------------------------------------------------------------
  40 ** Other defines.
  41 */
  42
  43 #define PAF_HEADER_LENGTH                       2048
  44
  45 #define PAF24_SAMPLES_PER_BLOCK         10
  46 #define PAF24_BLOCK_SIZE                        32
  47
  48 /*------------------------------------------------------------------------------
  49 ** Typedefs.
  50 */
  51
  52 typedef struct
  53 {       int     version ;
  54         int     endianness ;
  55     int samplerate ;
  56     int format ;
  57         int     channels ;
  58         int     source ;
  59 } PAF_FMT ;
  60
  61 typedef struct
  62 {       int                             max_blocks, channels, samplesperblock, blocksize ;
  63         int                             read_block, write_block, read_count, write_count ;
  64         sf_count_t              sample_count ;
  65         int                             *samples ;
  66         unsigned char   *block ;
  67 #if HAVE_FLEXIBLE_ARRAY
  68         int                             data [] ; /* ISO C99 struct flexible array. */
  69 #else
  70         int                             data [1] ; /* This is a hack and may not work. */
  71 #endif
  72 } PAF24_PRIVATE ;
  73
  74 /*------------------------------------------------------------------------------
  75 ** Private static functions.
  76 */
  77
  78 static int paf24_init (SF_PRIVATE *psf) ;
  79
  80 static int      paf_read_header (SF_PRIVATE *psf) ;
  81 static int      paf_write_header (SF_PRIVATE *psf, int calc_length) ;
  82
  83 static sf_count_t paf24_read_s (SF_PRIVATE *psf, short *ptr, sf_count_t len) ;
  84 static sf_count_t paf24_read_i (SF_PRIVATE *psf, int *ptr, sf_count_t len) ;
  85 static sf_count_t paf24_read_f (SF_PRIVATE *psf, float *ptr, sf_count_t len) ;
  86 static sf_count_t paf24_read_d (SF_PRIVATE *psf, double *ptr, sf_count_t len) ;
  87
  88 static sf_count_t paf24_write_s (SF_PRIVATE *psf, const short *ptr, sf_count_t len) ;
  89 static sf_count_t paf24_write_i (SF_PRIVATE *psf, const int *ptr, sf_count_t len) ;
  90 static sf_count_t paf24_write_f (SF_PRIVATE *psf, const float *ptr, sf_count_t len) ;
  91 static sf_count_t paf24_write_d (SF_PRIVATE *psf, const double *ptr, sf_count_t len) ;
  92
  93 static sf_count_t paf24_seek (SF_PRIVATE *psf, int mode, sf_count_t offset) ;
  94
  95 enum
  96 {       PAF_PCM_16 = 0,
  97         PAF_PCM_24 = 1,
  98         PAF_PCM_S8 = 2
  99 } ;
 100
 101 /*------------------------------------------------------------------------------
 102 ** Public function.
 103 */
 104
 105 int
 106 paf_open        (SF_PRIVATE *psf)
 107 {       int             subformat, error, endian ;
 108
 109         psf->dataoffset = PAF_HEADER_LENGTH ;
 110
 111         if (psf->mode == SFM_READ || (psf->mode == SFM_RDWR && psf->filelength > 0))
 112         {       if ((error = paf_read_header (psf)))
 113                         return error ;
 114                 } ;
 115
 116         subformat = psf->sf.format & SF_FORMAT_SUBMASK ;
 117
 118         if (psf->mode == SFM_WRITE || psf->mode == SFM_RDWR)
 119         {       if ((psf->sf.format & SF_FORMAT_TYPEMASK) != SF_FORMAT_PAF)
 120                         return  SFE_BAD_OPEN_FORMAT ;
 121
 122                 endian = psf->sf.format & SF_FORMAT_ENDMASK ;
 123
 124                 /* PAF is by default big endian. */
 125                 psf->endian = SF_ENDIAN_BIG ;
 126
 127                 if (endian == SF_ENDIAN_LITTLE || (CPU_IS_LITTLE_ENDIAN && (endian == SF_ENDIAN_CPU)))
 128                         psf->endian = SF_ENDIAN_LITTLE ;
 129
 130                 if ((error = paf_write_header (psf, SF_FALSE)))
 131                         return error ;
 132
 133                 psf->write_header = paf_write_header ;
 134                 } ;
 135
 136         switch (subformat)
 137         {       case SF_FORMAT_PCM_S8 :
 138                                         psf->bytewidth = 1 ;
 139                                         error = pcm_init (psf) ;
 140                                         break ;
 141
 142                 case SF_FORMAT_PCM_16 :
 143                                         psf->bytewidth = 2 ;
 144                                         error = pcm_init (psf) ;
 145                                         break ;
 146
 147                 case SF_FORMAT_PCM_24 :
 148                                         /* No bytewidth because of whacky 24 bit encoding. */
 149                                         error = paf24_init (psf) ;
 150                                         break ;
 151
 152                 default : return SFE_PAF_UNKNOWN_FORMAT ;
 153                 } ;
 154
 155         return error ;
 156 } /* paf_open */
 157
 158 /*------------------------------------------------------------------------------
 159 */
 160
 161 static int
 162 paf_read_header (SF_PRIVATE *psf)
 163 {       PAF_FMT         paf_fmt ;
 164         int                     marker ;
 165
 166         memset (&paf_fmt, 0, sizeof (paf_fmt)) ;
 167         psf_binheader_readf (psf, "pm", 0, &marker) ;
 168
 169         psf_log_printf (psf, "Signature   : '%M'\n", marker) ;
 170
 171         if (marker == PAF_MARKER)
 172         {       psf_binheader_readf (psf, "E444444", &(paf_fmt.version), &(paf_fmt.endianness),
 173                         &(paf_fmt.samplerate), &(paf_fmt.format), &(paf_fmt.channels), &(paf_fmt.source)) ;
 174                 }
 175         else if (marker == FAP_MARKER)
 176         {       psf_binheader_readf (psf, "e444444", &(paf_fmt.version), &(paf_fmt.endianness),
 177                         &(paf_fmt.samplerate), &(paf_fmt.format), &(paf_fmt.channels), &(paf_fmt.source)) ;
 178                 }
 179         else
 180                 return SFE_PAF_NO_MARKER ;
 181
 182         psf_log_printf (psf, "Version     : %d\n", paf_fmt.version) ;
 183
 184         if (paf_fmt.version != 0)
 185         {       psf_log_printf (psf, "*** Bad version number. should be zero.\n") ;
 186                 return SFE_PAF_VERSION ;
 187                 } ;
 188
 189         psf_log_printf (psf, "Sample Rate : %d\n", paf_fmt.samplerate) ;
 190         psf_log_printf (psf, "Channels    : %d\n", paf_fmt.channels) ;
 191
 192         psf_log_printf (psf, "Endianness  : %d => ", paf_fmt.endianness) ;
 193         if (paf_fmt.endianness)
 194         {       psf_log_printf (psf, "Little\n", paf_fmt.endianness) ;
 195                 psf->endian = SF_ENDIAN_LITTLE ;
 196                 }
 197         else
 198         {       psf_log_printf (psf, "Big\n", paf_fmt.endianness) ;
 199                 psf->endian = SF_ENDIAN_BIG ;
 200                 } ;
 201
 202         if (psf->filelength < PAF_HEADER_LENGTH)
 203                 return SFE_PAF_SHORT_HEADER ;
 204
 205         psf->datalength = psf->filelength - psf->dataoffset ;
 206
 207         psf_binheader_readf (psf, "p", (int) psf->dataoffset) ;
 208
 209         psf->sf.samplerate      = paf_fmt.samplerate ;
 210         psf->sf.channels        = paf_fmt.channels ;
 211
 212         /* Only fill in type major. */
 213         psf->sf.format = SF_FORMAT_PAF ;
 214
 215         psf_log_printf (psf, "Format      : %d => ", paf_fmt.format) ;
 216
 217         /* PAF is by default big endian. */
 218         psf->sf.format |= paf_fmt.endianness ? SF_ENDIAN_LITTLE : SF_ENDIAN_BIG ;
 219
 220         switch (paf_fmt.format)
 221         {       case PAF_PCM_S8 :
 222                                         psf_log_printf (psf, "8 bit linear PCM\n") ;
 223                                         psf->bytewidth = 1 ;
 224
 225                                         psf->sf.format |= SF_FORMAT_PCM_S8 ;
 226
 227                                         psf->blockwidth = psf->bytewidth * psf->sf.channels ;
 228                                         psf->sf.frames = psf->datalength / psf->blockwidth ;
 229                                         break ;
 230
 231                 case PAF_PCM_16 :
 232                                         psf_log_printf (psf, "16 bit linear PCM\n") ;
 233                                         psf->bytewidth = 2 ;
 234
 235                                         psf->sf.format |= SF_FORMAT_PCM_16 ;
 236
 237                                         psf->blockwidth = psf->bytewidth * psf->sf.channels ;
 238                                         psf->sf.frames = psf->datalength / psf->blockwidth ;
 239                                         break ;
 240
 241                 case PAF_PCM_24 :
 242                                         psf_log_printf (psf, "24 bit linear PCM\n") ;
 243                                         psf->bytewidth = 3 ;
 244
 245                                         psf->sf.format |= SF_FORMAT_PCM_24 ;
 246
 247                                         psf->blockwidth = 0 ;
 248                                         psf->sf.frames = PAF24_SAMPLES_PER_BLOCK * psf->datalength /
 249                                                                                         (PAF24_BLOCK_SIZE * psf->sf.channels) ;
 250                                         break ;
 251
 252                 default :       psf_log_printf (psf, "Unknown\n") ;
 253                                         return SFE_PAF_UNKNOWN_FORMAT ;
 254                                         break ;
 255                 } ;
 256
 257         psf_log_printf (psf, "Source      : %d => ", paf_fmt.source) ;
 258
 259         switch (paf_fmt.source)
 260         {       case 1 : psf_log_printf (psf, "Analog Recording\n") ;
 261                                         break ;
 262                 case 2 : psf_log_printf (psf, "Digital Transfer\n") ;
 263                                         break ;
 264                 case 3 : psf_log_printf (psf, "Multi-track Mixdown\n") ;
 265                                         break ;
 266                 case 5 : psf_log_printf (psf, "Audio Resulting From DSP Processing\n") ;
 267                                         break ;
 268                 default : psf_log_printf (psf, "Unknown\n") ;
 269                                         break ;
 270                 } ;
 271
 272         return 0 ;
 273 } /* paf_read_header */
 274
 275 static int
 276 paf_write_header (SF_PRIVATE *psf, int calc_length)
 277 {       int                     paf_format ;
 278
 279         /* PAF header already written so no need to re-write. */
 280         if (psf_ftell (psf) >= PAF_HEADER_LENGTH)
 281                 return 0 ;
 282
 283         psf->dataoffset = PAF_HEADER_LENGTH ;
 284
 285         psf->dataoffset = PAF_HEADER_LENGTH ;
 286
 287         /* Prevent compiler warning. */
 288         calc_length = calc_length ;
 289
 290         switch (psf->sf.format & SF_FORMAT_SUBMASK)
 291         {       case SF_FORMAT_PCM_S8 :
 292                                         paf_format = PAF_PCM_S8 ;
 293                                         break ;
 294
 295                 case SF_FORMAT_PCM_16 :
 296                                         paf_format = PAF_PCM_16 ;
 297                                         break ;
 298
 299                 case SF_FORMAT_PCM_24 :
 300                                         paf_format = PAF_PCM_24 ;
 301                                         break ;
 302
 303                 default : return SFE_PAF_UNKNOWN_FORMAT ;
 304                 } ;
 305
 306         /* Reset the current header length to zero. */
 307         psf->header [0] = 0 ;
 308         psf->headindex = 0 ;
 309
 310         if (psf->endian == SF_ENDIAN_BIG)
 311         {       /* Marker, version, endianness, samplerate */
 312                 psf_binheader_writef (psf, "Em444", PAF_MARKER, 0, 0, psf->sf.samplerate) ;
 313                 /* format, channels, source */
 314                 psf_binheader_writef (psf, "E444", paf_format, psf->sf.channels, 0) ;
 315                 }
 316         else if (psf->endian == SF_ENDIAN_LITTLE)
 317         {       /* Marker, version, endianness, samplerate */
 318                 psf_binheader_writef (psf, "em444", FAP_MARKER, 0, 1, psf->sf.samplerate) ;
 319                 /* format, channels, source */
 320                 psf_binheader_writef (psf, "e444", paf_format, psf->sf.channels, 0) ;
 321                 } ;
 322
 323         /* Zero fill to dataoffset. */
 324         psf_binheader_writef (psf, "z", (size_t) (psf->dataoffset - psf->headindex)) ;
 325
 326         psf_fwrite (psf->header, psf->headindex, 1, psf) ;
 327
 328         return psf->error ;
 329 } /* paf_write_header */
 330
 331 /*===============================================================================
 332 **      24 bit PAF files have a really weird encoding.
 333 **  For a mono file, 10 samples (each being 3 bytes) are packed into a 32 byte
 334 **      block. The 8 ints in this 32 byte block are then endian swapped (as ints)
 335 **      if necessary before being written to disk.
 336 **  For a stereo file, blocks of 10 samples from the same channel are encoded
 337 **  into 32 bytes as for the mono case. The 32 byte blocks are then interleaved
 338 **      on disk.
 339 **      Reading has to reverse the above process :-).
 340 **      Weird!!!
 341 **
 342 **      The code below attempts to gain efficiency while maintaining readability.
 343 */
 344
 345 static int paf24_read_block (SF_PRIVATE *psf, PAF24_PRIVATE *ppaf24) ;
 346 static int paf24_write_block (SF_PRIVATE *psf, PAF24_PRIVATE *ppaf24) ;
 347 static int paf24_close (SF_PRIVATE *psf) ;
 348
 349
 350 static int
 351 paf24_init (SF_PRIVATE *psf)
 352 {       PAF24_PRIVATE   *ppaf24 ;
 353         int     paf24size ;
 354
 355         paf24size = sizeof (PAF24_PRIVATE) + psf->sf.channels *
 356                                         (PAF24_BLOCK_SIZE + PAF24_SAMPLES_PER_BLOCK * sizeof (int)) ;
 357
 358         /*
 359         **      Not exatly sure why this needs to be here but the tests
 360         **      fail without it.
 361         */
 362         psf->last_op = 0 ;
 363
 364         if (! (psf->fdata = malloc (paf24size)))
 365                 return SFE_MALLOC_FAILED ;
 366
 367         ppaf24 = (PAF24_PRIVATE*) psf->fdata ;
 368         memset (ppaf24, 0, paf24size) ;
 369
 370         ppaf24->channels        = psf->sf.channels ;
 371         ppaf24->samples         = ppaf24->data ;
 372         ppaf24->block           = (unsigned char*) (ppaf24->data + PAF24_SAMPLES_PER_BLOCK * ppaf24->channels) ;
 373
 374         ppaf24->blocksize = PAF24_BLOCK_SIZE * ppaf24->channels ;
 375         ppaf24->samplesperblock = PAF24_SAMPLES_PER_BLOCK ;
 376
 377         if (psf->mode == SFM_READ || psf->mode == SFM_RDWR)
 378         {       paf24_read_block (psf, ppaf24) ;        /* Read first block. */
 379
 380                 psf->read_short         = paf24_read_s ;
 381                 psf->read_int           = paf24_read_i ;
 382                 psf->read_float         = paf24_read_f ;
 383                 psf->read_double        = paf24_read_d ;
 384                 } ;
 385
 386         if (psf->mode == SFM_WRITE || psf->mode == SFM_RDWR)
 387         {       psf->write_short        = paf24_write_s ;
 388                 psf->write_int          = paf24_write_i ;
 389                 psf->write_float        = paf24_write_f ;
 390                 psf->write_double       = paf24_write_d ;
 391                 } ;
 392
 393         psf->seek       = paf24_seek ;
 394         psf->container_close    = paf24_close ;
 395
 396         psf->filelength = psf_get_filelen (psf) ;
 397         psf->datalength = psf->filelength - psf->dataoffset ;
 398
 399         if (psf->datalength % PAF24_BLOCK_SIZE)
 400         {       if (psf->mode == SFM_READ)
 401                         psf_log_printf (psf, "*** Warning : file seems to be truncated.\n") ;
 402                 ppaf24->max_blocks = psf->datalength / ppaf24->blocksize + 1 ;
 403                 }
 404         else
 405                 ppaf24->max_blocks = psf->datalength / ppaf24->blocksize ;
 406
 407         ppaf24->read_block = 0 ;
 408         if (psf->mode == SFM_RDWR)
 409                 ppaf24->write_block = ppaf24->max_blocks ;
 410         else
 411                 ppaf24->write_block = 0 ;
 412
 413         psf->sf.frames = ppaf24->samplesperblock * ppaf24->max_blocks ;
 414         ppaf24->sample_count = psf->sf.frames ;
 415
 416         return 0 ;
 417 } /* paf24_init */
 418
 419 static sf_count_t
 420 paf24_seek (SF_PRIVATE *psf, int mode, sf_count_t offset)
 421 {       PAF24_PRIVATE   *ppaf24 ;
 422         int                             newblock, newsample ;
 423
 424         if (psf->fdata == NULL)
 425         {       psf->error = SFE_INTERNAL ;
 426                 return PSF_SEEK_ERROR ;
 427                 } ;
 428
 429         ppaf24 = (PAF24_PRIVATE*) psf->fdata ;
 430
 431         if (mode == SFM_READ && ppaf24->write_count > 0)
 432                 paf24_write_block (psf, ppaf24) ;
 433
 434         newblock        = offset / ppaf24->samplesperblock ;
 435         newsample       = offset % ppaf24->samplesperblock ;
 436
 437         switch (mode)
 438         {       case SFM_READ :
 439                                 if (psf->last_op == SFM_WRITE && ppaf24->write_count)
 440                                         paf24_write_block (psf, ppaf24) ;
 441
 442                                 psf_fseek (psf, psf->dataoffset + newblock * ppaf24->blocksize, SEEK_SET) ;
 443                                 ppaf24->read_block = newblock ;
 444                                 paf24_read_block (psf, ppaf24) ;
 445                                 ppaf24->read_count = newsample ;
 446                                 break ;
 447
 448                 case SFM_WRITE :
 449                                 if (offset > ppaf24->sample_count)
 450                                 {       psf->error = SFE_BAD_SEEK ;
 451                                         return PSF_SEEK_ERROR ;
 452                                         } ;
 453
 454                                 if (psf->last_op == SFM_WRITE && ppaf24->write_count)
 455                                         paf24_write_block (psf, ppaf24) ;
 456
 457                                 psf_fseek (psf, psf->dataoffset + newblock * ppaf24->blocksize, SEEK_SET) ;
 458                                 ppaf24->write_block = newblock ;
 459                                 paf24_read_block (psf, ppaf24) ;
 460                                 ppaf24->write_count = newsample ;
 461                                 break ;
 462
 463                 default :
 464                                 psf->error = SFE_BAD_SEEK ;
 465                                 return PSF_SEEK_ERROR ;
 466                 } ;
 467
 468         return newblock * ppaf24->samplesperblock + newsample ;
 469 } /* paf24_seek */
 470
 471 static int
 472 paf24_close (SF_PRIVATE *psf)
 473 {       PAF24_PRIVATE *ppaf24 ;
 474
 475         if (psf->fdata == NULL)
 476                 return 0 ;
 477
 478         ppaf24 = (PAF24_PRIVATE*) psf->fdata ;
 479
 480         if (psf->mode == SFM_WRITE || psf->mode == SFM_RDWR)
 481         {       if (ppaf24->write_count > 0)
 482                         paf24_write_block (psf, ppaf24) ;
 483                 } ;
 484
 485         return 0 ;
 486 } /* paf24_close */
 487
 488 /*---------------------------------------------------------------------------
 489 */
 490 static int
 491 paf24_read_block (SF_PRIVATE *psf, PAF24_PRIVATE *ppaf24)
 492 {       int                             k, channel ;
 493         unsigned char   *cptr ;
 494
 495         ppaf24->read_block ++ ;
 496         ppaf24->read_count = 0 ;
 497
 498         if (ppaf24->read_block * ppaf24->samplesperblock > ppaf24->sample_count)
 499         {       memset (ppaf24->samples, 0, ppaf24->samplesperblock * ppaf24->channels) ;
 500                 return 1 ;
 501                 } ;
 502
 503         /* Read the block. */
 504         if ((k = psf_fread (ppaf24->block, 1, ppaf24->blocksize, psf)) != ppaf24->blocksize)
 505                 psf_log_printf (psf, "*** Warning : short read (%d != %d).\n", k, ppaf24->blocksize) ;
 506
 507
 508         if (CPU_IS_LITTLE_ENDIAN)
 509         {       /* Do endian swapping if necessary. */
 510                 if (psf->endian == SF_ENDIAN_BIG)
 511                         endswap_int_array       (ppaf24->data, 8 * ppaf24->channels) ;
 512
 513                 /* Unpack block. */
 514                 for (k = 0 ; k < PAF24_SAMPLES_PER_BLOCK * ppaf24->channels ; k++)
 515                 {       channel = k % ppaf24->channels ;
 516                         cptr = ppaf24->block + PAF24_BLOCK_SIZE * channel + 3 * (k / ppaf24->channels) ;
 517                         ppaf24->samples [k] = (cptr [0] << 8) | (cptr [1] << 16) | (cptr [2] << 24) ;
 518                         } ;
 519                 }
 520         else
 521         {       /* Do endian swapping if necessary. */
 522                 if (psf->endian == SF_ENDIAN_BIG)
 523                         endswap_int_array       (ppaf24->data, 8 * ppaf24->channels) ;
 524
 525                 /* Unpack block. */
 526                 for (k = 0 ; k < PAF24_SAMPLES_PER_BLOCK * ppaf24->channels ; k++)
 527                 {       channel = k % ppaf24->channels ;
 528                         cptr = ppaf24->block + PAF24_BLOCK_SIZE * channel + 3 * (k / ppaf24->channels) ;
 529                         ppaf24->samples [k] = (cptr [0] << 8) | (cptr [1] << 16) | (cptr [2] << 24) ;
 530                         } ;
 531                 } ;
 532
 533         return 1 ;
 534 } /* paf24_read_block */
 535
 536 static int
 537 paf24_read (SF_PRIVATE *psf, PAF24_PRIVATE *ppaf24, int *ptr, int len)
 538 {       int     count, total = 0 ;
 539
 540         while (total < len)
 541         {       if (ppaf24->read_block * ppaf24->samplesperblock >= ppaf24->sample_count)
 542                 {       memset (&(ptr [total]), 0, (len - total) * sizeof (int)) ;
 543                         return total ;
 544                         } ;
 545
 546                 if (ppaf24->read_count >= ppaf24->samplesperblock)
 547                         paf24_read_block (psf, ppaf24) ;
 548
 549                 count = (ppaf24->samplesperblock - ppaf24->read_count) * ppaf24->channels ;
 550                 count = (len - total > count) ? count : len - total ;
 551
 552                 memcpy (&(ptr [total]), &(ppaf24->samples [ppaf24->read_count * ppaf24->channels]), count * sizeof (int)) ;
 553                 total += count ;
 554                 ppaf24->read_count += count / ppaf24->channels ;
 555                 } ;
 556
 557         return total ;
 558 } /* paf24_read */
 559
 560 static sf_count_t
 561 paf24_read_s (SF_PRIVATE *psf, short *ptr, sf_count_t len)
 562 {       PAF24_PRIVATE   *ppaf24 ;
 563         int                             *iptr ;
 564         int                             k, bufferlen, readcount, count ;
 565         sf_count_t              total = 0 ;
 566
 567         if (psf->fdata == NULL)
 568                 return 0 ;
 569         ppaf24 = (PAF24_PRIVATE*) psf->fdata ;
 570
 571         iptr = psf->u.ibuf ;
 572         bufferlen = ARRAY_LEN (psf->u.ibuf) ;
 573         while (len > 0)
 574         {       readcount = (len >= bufferlen) ? bufferlen : len ;
 575                 count = paf24_read (psf, ppaf24, iptr, readcount) ;
 576                 for (k = 0 ; k < readcount ; k++)
 577                         ptr [total + k] = iptr [k] >> 16 ;
 578                 total += count ;
 579                 len -= readcount ;
 580                 } ;
 581         return total ;
 582 } /* paf24_read_s */
 583
 584 static sf_count_t
 585 paf24_read_i (SF_PRIVATE *psf, int *ptr, sf_count_t len)
 586 {       PAF24_PRIVATE *ppaf24 ;
 587         int                             total ;
 588
 589         if (psf->fdata == NULL)
 590                 return 0 ;
 591         ppaf24 = (PAF24_PRIVATE*) psf->fdata ;
 592
 593         total = paf24_read (psf, ppaf24, ptr, len) ;
 594
 595         return total ;
 596 } /* paf24_read_i */
 597
 598 static sf_count_t
 599 paf24_read_f (SF_PRIVATE *psf, float *ptr, sf_count_t len)
 600 {       PAF24_PRIVATE   *ppaf24 ;
 601         int                             *iptr ;
 602         int                             k, bufferlen, readcount, count ;
 603         sf_count_t              total = 0 ;
 604         float                   normfact ;
 605
 606         if (psf->fdata == NULL)
 607                 return 0 ;
 608         ppaf24 = (PAF24_PRIVATE*) psf->fdata ;
 609
 610         normfact = (psf->norm_float == SF_TRUE) ? (1.0 / 0x80000000) : (1.0 / 0x100) ;
 611
 612         iptr = psf->u.ibuf ;
 613         bufferlen = ARRAY_LEN (psf->u.ibuf) ;
 614         while (len > 0)
 615         {       readcount = (len >= bufferlen) ? bufferlen : len ;
 616                 count = paf24_read (psf, ppaf24, iptr, readcount) ;
 617                 for (k = 0 ; k < readcount ; k++)
 618                         ptr [total + k] = normfact * iptr [k] ;
 619                 total += count ;
 620                 len -= readcount ;
 621                 } ;
 622         return total ;
 623 } /* paf24_read_f */
 624
 625 static sf_count_t
 626 paf24_read_d (SF_PRIVATE *psf, double *ptr, sf_count_t len)
 627 {       PAF24_PRIVATE   *ppaf24 ;
 628         int                             *iptr ;
 629         int                             k, bufferlen, readcount, count ;
 630         sf_count_t              total = 0 ;
 631         double                  normfact ;
 632
 633         if (psf->fdata == NULL)
 634                 return 0 ;
 635         ppaf24 = (PAF24_PRIVATE*) psf->fdata ;
 636
 637         normfact = (psf->norm_double == SF_TRUE) ? (1.0 / 0x80000000) : (1.0 / 0x100) ;
 638
 639         iptr = psf->u.ibuf ;
 640         bufferlen = ARRAY_LEN (psf->u.ibuf) ;
 641         while (len > 0)
 642         {       readcount = (len >= bufferlen) ? bufferlen : len ;
 643                 count = paf24_read (psf, ppaf24, iptr, readcount) ;
 644                 for (k = 0 ; k < readcount ; k++)
 645                         ptr [total + k] = normfact * iptr [k] ;
 646                 total += count ;
 647                 len -= readcount ;
 648                 } ;
 649         return total ;
 650 } /* paf24_read_d */
 651
 652 /*---------------------------------------------------------------------------
 653 */
 654
 655 static int
 656 paf24_write_block (SF_PRIVATE *psf, PAF24_PRIVATE *ppaf24)
 657 {       int                             k, nextsample, channel ;
 658         unsigned char   *cptr ;
 659
 660         /* First pack block. */
 661
 662         if (CPU_IS_LITTLE_ENDIAN)
 663         {       for (k = 0 ; k < PAF24_SAMPLES_PER_BLOCK * ppaf24->channels ; k++)
 664                 {       channel = k % ppaf24->channels ;
 665                         cptr = ppaf24->block + PAF24_BLOCK_SIZE * channel + 3 * (k / ppaf24->channels) ;
 666                         nextsample = ppaf24->samples [k] >> 8 ;
 667                         cptr [0] = nextsample ;
 668                         cptr [1] = nextsample >> 8 ;
 669                         cptr [2] = nextsample >> 16 ;
 670                         } ;
 671
 672                 /* Do endian swapping if necessary. */
 673                 if (psf->endian == SF_ENDIAN_BIG)
 674                         endswap_int_array (ppaf24->data, 8 * ppaf24->channels) ;
 675                 }
 676         else if (CPU_IS_BIG_ENDIAN)
 677         {       /* This is correct. */
 678                 for (k = 0 ; k < PAF24_SAMPLES_PER_BLOCK * ppaf24->channels ; k++)
 679                 {       channel = k % ppaf24->channels ;
 680                         cptr = ppaf24->block + PAF24_BLOCK_SIZE * channel + 3 * (k / ppaf24->channels) ;
 681                         nextsample = ppaf24->samples [k] >> 8 ;
 682                         cptr [0] = nextsample ;
 683                         cptr [1] = nextsample >> 8 ;
 684                         cptr [2] = nextsample >> 16 ;
 685                         } ;
 686                 if (psf->endian == SF_ENDIAN_BIG)
 687                         endswap_int_array (ppaf24->data, 8 * ppaf24->channels) ;
 688                 } ;
 689
 690         /* Write block to disk. */
 691         if ((k = psf_fwrite (ppaf24->block, 1, ppaf24->blocksize, psf)) != ppaf24->blocksize)
 692                 psf_log_printf (psf, "*** Warning : short write (%d != %d).\n", k, ppaf24->blocksize) ;
 693
 694         if (ppaf24->sample_count < ppaf24->write_block * ppaf24->samplesperblock + ppaf24->write_count)
 695                 ppaf24->sample_count = ppaf24->write_block * ppaf24->samplesperblock + ppaf24->write_count ;
 696
 697         if (ppaf24->write_count == ppaf24->samplesperblock)
 698         {       ppaf24->write_block ++ ;
 699                 ppaf24->write_count = 0 ;
 700                 } ;
 701
 702         return 1 ;
 703 } /* paf24_write_block */
 704
 705 static int
 706 paf24_write (SF_PRIVATE *psf, PAF24_PRIVATE *ppaf24, const int *ptr, int len)
 707 {       int             count, total = 0 ;
 708
 709         while (total < len)
 710         {       count = (ppaf24->samplesperblock - ppaf24->write_count) * ppaf24->channels ;
 711
 712                 if (count > len - total)
 713                         count = len - total ;
 714
 715                 memcpy (&(ppaf24->samples [ppaf24->write_count * ppaf24->channels]), &(ptr [total]), count * sizeof (int)) ;
 716                 total += count ;
 717                 ppaf24->write_count += count / ppaf24->channels ;
 718
 719                 if (ppaf24->write_count >= ppaf24->samplesperblock)
 720                         paf24_write_block (psf, ppaf24) ;
 721                 } ;
 722
 723         return total ;
 724 } /* paf24_write */
 725
 726 static sf_count_t
 727 paf24_write_s (SF_PRIVATE *psf, const short *ptr, sf_count_t len)
 728 {       PAF24_PRIVATE   *ppaf24 ;
 729         int                             *iptr ;
 730         int                             k, bufferlen, writecount = 0, count ;
 731         sf_count_t              total = 0 ;
 732
 733         if (psf->fdata == NULL)
 734                 return 0 ;
 735         ppaf24 = (PAF24_PRIVATE*) psf->fdata ;
 736
 737         iptr = psf->u.ibuf ;
 738         bufferlen = ARRAY_LEN (psf->u.ibuf) ;
 739         while (len > 0)
 740         {       writecount = (len >= bufferlen) ? bufferlen : len ;
 741                 for (k = 0 ; k < writecount ; k++)
 742                         iptr [k] = ptr [total + k] << 16 ;
 743                 count = paf24_write (psf, ppaf24, iptr, writecount) ;
 744                 total += count ;
 745                 len -= writecount ;
 746                 if (count != writecount)
 747                         break ;
 748                 } ;
 749         return total ;
 750 } /* paf24_write_s */
 751
 752 static sf_count_t
 753 paf24_write_i (SF_PRIVATE *psf, const int *ptr, sf_count_t len)
 754 {       PAF24_PRIVATE   *ppaf24 ;
 755         int                             writecount, count ;
 756         sf_count_t              total = 0 ;
 757
 758         if (psf->fdata == NULL)
 759                 return 0 ;
 760         ppaf24 = (PAF24_PRIVATE*) psf->fdata ;
 761
 762         while (len > 0)
 763         {       writecount = (len > 0x10000000) ? 0x10000000 : (int) len ;
 764
 765                 count = paf24_write (psf, ppaf24, ptr, writecount) ;
 766
 767                 total += count ;
 768                 len -= count ;
 769                 if (count != writecount)
 770                         break ;
 771                 } ;
 772
 773         return total ;
 774 } /* paf24_write_i */
 775
 776 static sf_count_t
 777 paf24_write_f (SF_PRIVATE *psf, const float *ptr, sf_count_t len)
 778 {       PAF24_PRIVATE   *ppaf24 ;
 779         int                             *iptr ;
 780         int                             k, bufferlen, writecount = 0, count ;
 781         sf_count_t              total = 0 ;
 782         float                   normfact ;
 783
 784         if (psf->fdata == NULL)
 785                 return 0 ;
 786         ppaf24 = (PAF24_PRIVATE*) psf->fdata ;
 787
 788         normfact = (psf->norm_float == SF_TRUE) ? (1.0 * 0x7FFFFFFF) : (1.0 / 0x100) ;
 789
 790         iptr = psf->u.ibuf ;
 791         bufferlen = ARRAY_LEN (psf->u.ibuf) ;
 792         while (len > 0)
 793         {       writecount = (len >= bufferlen) ? bufferlen : len ;
 794                 for (k = 0 ; k < writecount ; k++)
 795                         iptr [k] = lrintf (normfact * ptr [total + k]) ;
 796                 count = paf24_write (psf, ppaf24, iptr, writecount) ;
 797                 total += count ;
 798                 len -= writecount ;
 799                 if (count != writecount)
 800                         break ;
 801                 } ;
 802
 803         return total ;
 804 } /* paf24_write_f */
 805
 806 static sf_count_t
 807 paf24_write_d (SF_PRIVATE *psf, const double *ptr, sf_count_t len)
 808 {       PAF24_PRIVATE   *ppaf24 ;
 809         int                             *iptr ;
 810         int                             k, bufferlen, writecount = 0, count ;
 811         sf_count_t              total = 0 ;
 812         double                  normfact ;
 813
 814         if (psf->fdata == NULL)
 815                 return 0 ;
 816         ppaf24 = (PAF24_PRIVATE*) psf->fdata ;
 817
 818         normfact = (psf->norm_double == SF_TRUE) ? (1.0 * 0x7FFFFFFF) : (1.0 / 0x100) ;
 819
 820         iptr = psf->u.ibuf ;
 821         bufferlen = ARRAY_LEN (psf->u.ibuf) ;
 822         while (len > 0)
 823         {       writecount = (len >= bufferlen) ? bufferlen : len ;
 824                 for (k = 0 ; k < writecount ; k++)
 825                         iptr [k] = lrint (normfact * ptr [total+k]) ;
 826                 count = paf24_write (psf, ppaf24, iptr, writecount) ;
 827                 total += count ;
 828                 len -= writecount ;
 829                 if (count != writecount)
 830                         break ;
 831                 } ;
 832
 833         return total ;
 834 } /* paf24_write_d */
 835
 836
 837 /*
 838 ** Do not edit or modify anything in this comment block.
 839 ** The arch-tag line is a file identity tag for the GNU Arch
 840 ** revision control system.
 841 **
 842 ** arch-tag: 477a5308-451e-4bbd-bab4-fab6caa4e884
 843 */