Print this page
patch mul-ste


  22 /*
  23  * Copyright 2008 Sun Microsystems, Inc.  All rights reserved.
  24  * Use is subject to license terms.
  25  */
  26 
  27 #pragma ident   "%Z%%M% %I%     %E% SMI"
  28 
  29 #include "lint.h"
  30 #include "base_conversion.h"
  31 #include <sys/isa_defs.h>
  32 
  33 #define CSR_DEFAULT 0x1f80
  34 
  35 /* The following should be coded as inline expansion templates.  */
  36 
  37 /*
  38  * Multiplies two normal or subnormal doubles, returns result and exceptions.
  39  */
  40 double
  41 __mul_set(double x, double y, int *pe) {
  42         extern void _putmxcsr(), _getmxcsr();
  43         int csr;
  44         double z;
  45 
  46         _putmxcsr(CSR_DEFAULT);
  47         z = x * y;
  48         _getmxcsr(&csr);




  49         if ((csr & 0x3f) == 0) {
  50                 *pe = 0;
  51         } else {
  52                 /* Result may not be exact. */
  53                 *pe = 1;
  54         }
  55         return (z);
  56 }
  57 
  58 /*
  59  * Divides two normal or subnormal doubles x/y, returns result and exceptions.
  60  */
  61 double
  62 __div_set(double x, double y, int *pe) {
  63         extern void _putmxcsr(), _getmxcsr();
  64         int csr;
  65         double z;
  66 
  67         _putmxcsr(CSR_DEFAULT);
  68         z = x / y;




  22 /*
  23  * Copyright 2008 Sun Microsystems, Inc.  All rights reserved.
  24  * Use is subject to license terms.
  25  */
  26 
  27 #pragma ident   "%Z%%M% %I%     %E% SMI"
  28 
  29 #include "lint.h"
  30 #include "base_conversion.h"
  31 #include <sys/isa_defs.h>
  32 
  33 #define CSR_DEFAULT 0x1f80
  34 
  35 /* The following should be coded as inline expansion templates.  */
  36 
  37 /*
  38  * Multiplies two normal or subnormal doubles, returns result and exceptions.
  39  */
  40 double
  41 __mul_set(double x, double y, int *pe) {
  42         int csr = CSR_DEFAULT;

  43         double z;
  44 
  45         __asm__ __volatile__(
  46             "ldmxcsr    %4\n"
  47             "mulsd      %2,%3\n"
  48             "stmxcsr    %0\n"
  49             : "=m" (csr), "=C" (z)
  50             : "C" (x), "C" (y), "m" (csr));
  51 
  52         if ((csr & 0x3f) == 0) {
  53                 *pe = 0;
  54         } else {
  55                 /* Result may not be exact. */
  56                 *pe = 1;
  57         }
  58         return (z);
  59 }
  60 
  61 /*
  62  * Divides two normal or subnormal doubles x/y, returns result and exceptions.
  63  */
  64 double
  65 __div_set(double x, double y, int *pe) {
  66         extern void _putmxcsr(), _getmxcsr();
  67         int csr;
  68         double z;
  69 
  70         _putmxcsr(CSR_DEFAULT);
  71         z = x / y;