]> Cypherpunks.ru repositories - gostls13.git/blobdiff - test/cmplxdivide.c
cmd/compile/internal/inline: score call sites exposed by inlines
[gostls13.git] / test / cmplxdivide.c
index b3c6055ed2301288f273767a6b41e714707c7ada..89a2868b75bdb249b905f4243cd5c1f8485569a8 100644 (file)
@@ -1,8 +1,19 @@
-// Copyright 2010 The Go Authors.  All rights reserved.
+// Copyright 2010 The Go Authors. All rights reserved.
 // Use of this source code is governed by a BSD-style
 // license that can be found in the LICENSE file.
 
-// gcc '-std=c99' cmplxdivide.c && a.out >cmplxdivide1.go
+// This C program generates the file cmplxdivide1.go. It uses the
+// output of the operations by C99 as the reference to check
+// the implementation of complex numbers in Go.
+// The generated file, cmplxdivide1.go, is compiled along
+// with the driver cmplxdivide.go (the names are confusing
+// and unimaginative) to run the actual test. This is done by
+// the usual test runner.
+//
+// The file cmplxdivide1.go is checked in to the repository, but
+// if it needs to be regenerated, compile and run this C program
+// like this:
+//     gcc '-std=c99' cmplxdivide.c && a.out >cmplxdivide1.go
 
 #include <complex.h>
 #include <math.h>
 #define nelem(x) (sizeof(x)/sizeof((x)[0]))
 
 double f[] = {
-       0,
-       1,
-       -1,
-       2,
+       0.0,
+       -0.0,
+       1.0,
+       -1.0,
+       2.0,
        NAN,
        INFINITY,
        -INFINITY,
 };
 
-char*
-fmt(double g)
-{
+char* fmt(double g) {
        static char buf[10][30];
        static int n;
        char *p;
-       
+
        p = buf[n++];
-       if(n == 10)
+       if(n == 10) {
                n = 0;
+       }
+
        sprintf(p, "%g", g);
-       if(strcmp(p, "-0") == 0)
-               strcpy(p, "negzero");
-       return p;
-}
 
-int
-iscnan(double complex d)
-{
-       return !isinf(creal(d)) && !isinf(cimag(d)) && (isnan(creal(d)) || isnan(cimag(d)));
-}
+       if(strcmp(p, "0") == 0) {
+               strcpy(p, "zero");
+               return p;
+       }
+
+       if(strcmp(p, "-0") == 0) {
+               strcpy(p, "-zero");
+               return p;
+       }
 
-double complex zero;   // attempt to hide zero division from gcc
+       return p;
+}
 
-int
-main(void)
-{
+int main(void) {
        int i, j, k, l;
        double complex n, d, q;
-       
+
+       printf("// skip\n");
        printf("// # generated by cmplxdivide.c\n");
        printf("\n");
        printf("package main\n");
-       printf("var tests = []Test{\n");
+       printf("\n");
+       printf("import \"math\"\n");
+       printf("\n");
+       printf("var (\n");
+       printf("\tnan     = math.NaN()\n");
+       printf("\tinf     = math.Inf(1)\n");
+       printf("\tzero    = 0.0\n");
+       printf(")\n");
+       printf("\n");
+       printf("var tests = []struct {\n");
+       printf("\tf, g complex128\n");
+       printf("\tout  complex128\n");
+       printf("}{\n");
+
        for(i=0; i<nelem(f); i++)
        for(j=0; j<nelem(f); j++)
        for(k=0; k<nelem(f); k++)
@@ -62,17 +87,8 @@ main(void)
                n = f[i] + f[j]*I;
                d = f[k] + f[l]*I;
                q = n/d;
-               
-               // BUG FIX.
-               // Gcc gets the wrong answer for NaN/0 unless both sides are NaN.
-               // That is, it treats (NaN+NaN*I)/0 = NaN+NaN*I (a complex NaN)
-               // but it then computes (1+NaN*I)/0 = Inf+NaN*I (a complex infinity).
-               // Since both numerators are complex NaNs, it seems that the
-               // results should agree in kind.  Override the gcc computation in this case.
-               if(iscnan(n) && d == 0)
-                       q = (NAN+NAN*I) / zero;
 
-               printf("\tTest{cmplx(%s, %s), cmplx(%s, %s), cmplx(%s, %s)},\n",
+               printf("\t{complex(%s, %s), complex(%s, %s), complex(%s, %s)},\n",
                        fmt(creal(n)), fmt(cimag(n)),
                        fmt(creal(d)), fmt(cimag(d)),
                        fmt(creal(q)), fmt(cimag(q)));