gcc/config/i386/crtfastmath.c

   1 /*
   2  * Copyright (C) 2005 Free Software Foundation, Inc.
   3  *
   4  * This file is free software; you can redistribute it and/or modify it
   5  * under the terms of the GNU General Public License as published by the
   6  * Free Software Foundation; either version 2, or (at your option) any
   7  * later version.
   8  *
   9  * In addition to the permissions in the GNU General Public License, the
  10  * Free Software Foundation gives you unlimited permission to link the
  11  * compiled version of this file with other programs, and to distribute
  12  * those programs without any restriction coming from the use of this
  13  * file.  (The General Public License restrictions do apply in other
  14  * respects; for example, they cover modification of the file, and
  15  * distribution when not linked into another program.)
  16  *
  17  * This file is distributed in the hope that it will be useful, but
  18  * WITHOUT ANY WARRANTY; without even the implied warranty of
  19  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  20  * General Public License for more details.
  21  *
  22  * You should have received a copy of the GNU General Public License
  23  * along with this program; see the file COPYING.  If not, write to
  24  * the Free Software Foundation, 59 Temple Place - Suite 330,
  25  * Boston, MA 02111-1307, USA.
  26  *
  27  *    As a special exception, if you link this library with files
  28  *    compiled with GCC to produce an executable, this does not cause
  29  *    the resulting executable to be covered by the GNU General Public License.
  30  *    This exception does not however invalidate any other reasons why
  31  *    the executable file might be covered by the GNU General Public License.
  32  */
  33
  34 #define MXCSR_DAZ (1 << 6)      /* Enable denormals are zero mode */
  35 #define MXCSR_FTZ (1 << 15)     /* Enable flush to zero mode */
  36
  37 static void __attribute__((constructor))
  38 set_fast_math (void)
  39 {
  40 #ifndef __x86_64__
  41   /* SSE is the part of 64bit. Only need to check it for 32bit.  */
  42   unsigned int eax, ebx, ecx, edx;
  43
  44   /* See if we can use cpuid.  */
  45   asm volatile ("pushfl; pushfl; popl %0; movl %0,%1; xorl %2,%0;"
  46                 "pushl %0; popfl; pushfl; popl %0; popfl"
  47                 : "=&r" (eax), "=&r" (ebx)
  48                 : "i" (0x00200000));
  49
  50   if (((eax ^ ebx) & 0x00200000) == 0)
  51     return;
  52
  53   /* Check the highest input value for eax.  */
  54   asm volatile ("xchgl %%ebx, %1; cpuid; xchgl %%ebx, %1"
  55                 : "=a" (eax), "=r" (ebx), "=c" (ecx), "=d" (edx)
  56                 : "0" (0));
  57
  58   if (eax == 0)
  59     return;
  60
  61   asm volatile ("xchgl %%ebx, %1; cpuid; xchgl %%ebx, %1"
  62                 : "=a" (eax), "=r" (ebx), "=c" (ecx), "=d" (edx)
  63                 : "0" (1));
  64
  65   if (edx & (1 << 25))
  66 #endif
  67     {
  68       unsigned int mxcsr = __builtin_ia32_stmxcsr ();
  69       mxcsr |= MXCSR_DAZ | MXCSR_FTZ;
  70       __builtin_ia32_ldmxcsr (mxcsr);
  71     }
  72 }