cannam@95: /*
cannam@95:  * Copyright (c) 2003, 2007-11 Matteo Frigo
cannam@95:  * Copyright (c) 2003, 2007-11 Massachusetts Institute of Technology
cannam@95:  *
cannam@95:  * This program is free software; you can redistribute it and/or modify
cannam@95:  * it under the terms of the GNU General Public License as published by
cannam@95:  * the Free Software Foundation; either version 2 of the License, or
cannam@95:  * (at your option) any later version.
cannam@95:  *
cannam@95:  * This program is distributed in the hope that it will be useful,
cannam@95:  * but WITHOUT ANY WARRANTY; without even the implied warranty of
cannam@95:  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
cannam@95:  * GNU General Public License for more details.
cannam@95:  *
cannam@95:  * You should have received a copy of the GNU General Public License
cannam@95:  * along with this program; if not, write to the Free Software
cannam@95:  * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301  USA
cannam@95:  *
cannam@95:  */
cannam@95: 
cannam@95: 
cannam@95: #include "ifftw.h"
cannam@95: 
cannam@95: /* ``kernel'' malloc(), with proper memory alignment */
cannam@95: 
cannam@95: #if defined(HAVE_DECL_MEMALIGN) && !HAVE_DECL_MEMALIGN
cannam@95: #  if defined(HAVE_MALLOC_H)
cannam@95: #    include <malloc.h>
cannam@95: #  else
cannam@95: extern void *memalign(size_t, size_t);
cannam@95: #  endif
cannam@95: #endif
cannam@95: 
cannam@95: #if defined(HAVE_DECL_POSIX_MEMALIGN) && !HAVE_DECL_POSIX_MEMALIGN
cannam@95: extern int posix_memalign(void **, size_t, size_t);
cannam@95: #endif
cannam@95: 
cannam@95: #if defined(macintosh) /* MacOS 9 */
cannam@95: #  include <Multiprocessing.h>
cannam@95: #endif
cannam@95: 
cannam@95: #define real_free free /* memalign and malloc use ordinary free */
cannam@95: 
cannam@95: #define IS_POWER_OF_TWO(n) (((n) > 0) && (((n) & ((n) - 1)) == 0))
cannam@95: #if defined(WITH_OUR_MALLOC) && (MIN_ALIGNMENT >= 8) && IS_POWER_OF_TWO(MIN_ALIGNMENT)
cannam@95: /* Our own MIN_ALIGNMENT-aligned malloc/free.  Assumes sizeof(void*) is a
cannam@95:    power of two <= 8 and that malloc is at least sizeof(void*)-aligned.
cannam@95: 
cannam@95:    The main reason for this routine is that, as of this writing,
cannam@95:    Windows does not include any aligned allocation routines in its
cannam@95:    system libraries, and instead provides an implementation with a
cannam@95:    Visual C++ "Processor Pack" that you have to statically link into
cannam@95:    your program.  We do not want to require users to have VC++
cannam@95:    (e.g. gcc/MinGW should be fine).  Our code should be at least as good
cannam@95:    as the MS _aligned_malloc, in any case, according to second-hand
cannam@95:    reports of the algorithm it employs (also based on plain malloc). */
cannam@95: static void *our_malloc(size_t n)
cannam@95: {
cannam@95:      void *p0, *p;
cannam@95:      if (!(p0 = malloc(n + MIN_ALIGNMENT))) return (void *) 0;
cannam@95:      p = (void *) (((uintptr_t) p0 + MIN_ALIGNMENT) & (~((uintptr_t) (MIN_ALIGNMENT - 1))));
cannam@95:      *((void **) p - 1) = p0;
cannam@95:      return p;
cannam@95: }
cannam@95: static void our_free(void *p)
cannam@95: {
cannam@95:      if (p) free(*((void **) p - 1));
cannam@95: }
cannam@95: #endif
cannam@95: 
cannam@95: void *X(kernel_malloc)(size_t n)
cannam@95: {
cannam@95:      void *p;
cannam@95: 
cannam@95: #if defined(MIN_ALIGNMENT)
cannam@95: 
cannam@95: #  if defined(WITH_OUR_MALLOC)
cannam@95:      p = our_malloc(n);
cannam@95: #    undef real_free
cannam@95: #    define real_free our_free
cannam@95: 
cannam@95: #  elif defined(__FreeBSD__) && (MIN_ALIGNMENT <= 16)
cannam@95:      /* FreeBSD does not have memalign, but its malloc is 16-byte aligned. */
cannam@95:      p = malloc(n);
cannam@95: 
cannam@95: #  elif (defined(__MACOSX__) || defined(__APPLE__)) && (MIN_ALIGNMENT <= 16)
cannam@95:      /* MacOS X malloc is already 16-byte aligned */
cannam@95:      p = malloc(n);
cannam@95: 
cannam@95: #  elif defined(HAVE_MEMALIGN)
cannam@95:      p = memalign(MIN_ALIGNMENT, n);
cannam@95: 
cannam@95: #  elif defined(HAVE_POSIX_MEMALIGN)
cannam@95:      /* note: posix_memalign is broken in glibc 2.2.5: it constrains
cannam@95: 	the size, not the alignment, to be (power of two) * sizeof(void*).
cannam@95:         The bug seems to have been fixed as of glibc 2.3.1. */
cannam@95:      if (posix_memalign(&p, MIN_ALIGNMENT, n))
cannam@95: 	  p = (void*) 0;
cannam@95: 
cannam@95: #  elif defined(__ICC) || defined(__INTEL_COMPILER) || defined(HAVE__MM_MALLOC)
cannam@95:      /* Intel's C compiler defines _mm_malloc and _mm_free intrinsics */
cannam@95:      p = (void *) _mm_malloc(n, MIN_ALIGNMENT);
cannam@95: #    undef real_free
cannam@95: #    define real_free _mm_free
cannam@95: 
cannam@95: #  elif defined(_MSC_VER)
cannam@95:      /* MS Visual C++ 6.0 with a "Processor Pack" supports SIMD
cannam@95: 	and _aligned_malloc/free (uses malloc.h) */
cannam@95:      p = (void *) _aligned_malloc(n, MIN_ALIGNMENT);
cannam@95: #    undef real_free
cannam@95: #    define real_free _aligned_free
cannam@95: 
cannam@95: #  elif defined(macintosh) /* MacOS 9 */
cannam@95:      p = (void *) MPAllocateAligned(n,
cannam@95: #    if MIN_ALIGNMENT == 8
cannam@95: 				    kMPAllocate8ByteAligned,
cannam@95: #    elif MIN_ALIGNMENT == 16
cannam@95: 				    kMPAllocate16ByteAligned,
cannam@95: #    elif MIN_ALIGNMENT == 32
cannam@95: 				    kMPAllocate32ByteAligned,
cannam@95: #    else
cannam@95: #      error "Unknown alignment for MPAllocateAligned"
cannam@95: #    endif
cannam@95: 				    0);
cannam@95: #    undef real_free
cannam@95: #    define real_free MPFree
cannam@95: 
cannam@95: #  else
cannam@95:      /* Add your machine here and send a patch to fftw@fftw.org 
cannam@95:         or (e.g. for Windows) configure --with-our-malloc */
cannam@95: #    error "Don't know how to malloc() aligned memory ... try configuring --with-our-malloc"
cannam@95: #  endif
cannam@95: 
cannam@95: #else /* !defined(MIN_ALIGNMENT) */
cannam@95:      p = malloc(n);
cannam@95: #endif
cannam@95: 
cannam@95:      return p;
cannam@95: }
cannam@95: 
cannam@95: void X(kernel_free)(void *p)
cannam@95: {
cannam@95:      real_free(p);
cannam@95: }