/*

 * Copyright (c) 2014, Oracle and/or its affiliates. All rights reserved.

 * DO NOT ALTER OR REMOVE COPYRIGHT NOTICES OR THIS FILE HEADER.

 *

 * This code is free software; you can redistribute it and/or modify it

 * under the terms of the GNU General Public License version 2 only, as

 * published by the Free Software Foundation.

 *

 * This code is distributed in the hope that it will be useful, but WITHOUT

 * ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or

 * FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License

 * version 2 for more details (a copy is included in the LICENSE file that

 * accompanied this code).

 *

 * You should have received a copy of the GNU General Public License version

 * 2 along with this work; if not, write to the Free Software Foundation,

 * Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA.

 *

 * Please contact Oracle, 500 Oracle Parkway, Redwood Shores, CA 94065 USA

 * or visit www.oracle.com if you need additional information or have any

 * questions.

 */

/**

 * @test

 * @bug 8015396

 * @summary double a%b returns NaN for some (a,b) (|a| < inf, |b|>0) (on Core i7 980X)

 *

 * @run main ModNaN

 */

public class ModNaN {

    /* This bug was seen in the field for a particular version of the VM,

     * but never reproduced internally, and the reason was never known,

     * nor were the exact circumstances of the failure.

     */

    /*

     * Failed on Windows 7/Core i7 980X/1.6.0_38 (64-bit):

     *

     * NaNs (i.e., when it fails, this is what we expect to see):

     *

     * 8.98846567431158E307 % 1.295163E-318 = NaN

     * (0x7FE0000000000000L % 0x0000000000040000L)

     *

     * 1.7976931348623157E308 % 2.59032E-318 = NaN

     * (0x7FEFFFFFFFFFFFFFL % 0x000000000007FFFFL)

     *

     * 1.7976931348623157E308 % 1.060997895E-314 = NaN

     * (0x7FEFFFFFFFFFFFFFL % 0x000000007FFFFFFFL)

     *

     * 1.7976931348623157E308 % 6.767486E-317 = NaN

     * (0x7FEFFFFFFFFFFFFFL % 0x0000000000d10208L)

     *

     * 1.7976931348623157E308 % 7.528725E-318 = NaN

     * (0x7FEFFFFFFFFFFFFFL % 0x0000000000174077L)

     *

     * These cases did not fail, even when the previous five did:

     * 8.98846567431158E307 % 1.29516E-318 = 2.53E-321

     * (0x7fe0000000000000L % 0x000000000003ffffL)

     *

     * 1.7976931348623157E308 % 2.590327E-318 = 0.0

     * (0x7fefffffffffffffL % 0x0000000000080000L)

     *

     * 1.7976931348623157E308 % 1.060965516E-314 = 9.35818525E-315

     * (0x7fefffffffffffffL % 0x000000007ffeffffL)

     *

     */

    static double[][] bad = new double[][] {

            /*

             * These hex numbers correspond to the base-10 doubles in the

             * comment above; this can be checked by observing the output

             * of testWithPrint.

             */

            new double[] { Double.longBitsToDouble(0x7FE0000000000000L),

                    Double.longBitsToDouble(0x0000000000040000L) },

            new double[] { Double.longBitsToDouble(0x7FEFFFFFFFFFFFFFL),

                    Double.longBitsToDouble(0x000000000007FFFFL) },

            new double[] { Double.longBitsToDouble(0x7FEFFFFFFFFFFFFFL),

                    Double.longBitsToDouble(0x000000007FFFFFFFL) },

            new double[] { Double.longBitsToDouble(0x7FEFFFFFFFFFFFFFL),

                    6.767486E-317 },

            new double[] { Double.longBitsToDouble(0x7FEFFFFFFFFFFFFFL),

                    7.528725E-318 }, };

    static double[][] good = new double[][] {

            new double[] { Double.longBitsToDouble(0x7FE0000000000000L),

                    Double.longBitsToDouble(0x000000000003FFFFL) },

            new double[] { Double.longBitsToDouble(0x7FEFFFFFFFFFFFFFL),

                    Double.longBitsToDouble(0x0000000000080000L) },

            new double[] { Double.longBitsToDouble(0x7FEFFFFFFFFFFFFFL),

                    Double.longBitsToDouble(0x000000007FFEFFFFL) }, };

    public static void main(String[] args) throws InterruptedException {

        int N = 10000;

        testWithPrint();

        for (int i = 0; i < N; i++)

            testStrict();

        for (int i = 0; i < N; i++)

            test();

        Thread.sleep(1000); // pause to let the compiler work

        for (int i = 0; i < 10; i++)

            testStrict();

        for (int i = 0; i < 10; i++)

            test();

    }

    public strictfp static void testWithPrint() {

        for (double[] ab : bad) {

            double a = ab[0];

            double b = ab[1];

            double mod = a % b;

            System.out.println("" + a + "("+toHexRep(a)+") mod " +

                                    b + "("+toHexRep(b)+") yields " +

                                    mod + "("+toHexRep(mod)+")");

        }

        for (double[] ab : good) {

            double a = ab[0];

            double b = ab[1];

            double mod = a % b;

            System.out.println("" + a + "("+toHexRep(a)+") mod " +

                                    b + "("+toHexRep(b)+") yields " +

                                    mod + "("+toHexRep(mod)+")");

        }

    }

    public strictfp static void testStrict() {

        for (double[] ab : bad) {

            double a = ab[0];

            double b = ab[1];

            double mod = a % b;

            check(mod);

        }

        for (double[] ab : good) {

            double a = ab[0];

            double b = ab[1];

            double mod = a % b;

            check(mod);

        }

    }

    public static void test() {

        for (double[] ab : bad) {

            double a = ab[0];

            double b = ab[1];

            double mod = a % b;

            check(mod);

        }

        for (double[] ab : good) {

            double a = ab[0];

            double b = ab[1];

            double mod = a % b;

            check(mod);

        }

    }

    static String toHexRep(double d) {

        return "0x" + Long.toHexString(Double.doubleToRawLongBits(d)) + "L";

    }

    static void check(double mod) {

        if (Double.isNaN(mod)) {

            throw new Error("Saw a NaN, fail");

        }

    }

}