/* * Copyright (c) 2017, 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. Oracle designates this * particular file as subject to the "Classpath" exception as provided * by Oracle in the LICENSE file that accompanied this code. * * 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 * questions. */ package jdk.incubator.vector; import java.nio.ByteBuffer; import java.nio.ByteOrder; import java.nio.DoubleBuffer; import java.nio.ReadOnlyBufferException; import java.util.Arrays; import java.util.Objects; import java.util.function.IntUnaryOperator; import jdk.internal.misc.Unsafe; import jdk.internal.vm.annotation.ForceInline; import static jdk.incubator.vector.VectorIntrinsics.*; @SuppressWarnings("cast") final class Double128Vector extends DoubleVector { private static final VectorSpecies SPECIES = DoubleVector.SPECIES_128; static final Double128Vector ZERO = new Double128Vector(); static final int LENGTH = SPECIES.length(); // Index vector species private static final IntVector.IntSpecies INDEX_SPECIES; static { int bitSize = Vector.bitSizeForVectorLength(int.class, LENGTH); INDEX_SPECIES = (IntVector.IntSpecies) IntVector.species(VectorShape.forBitSize(bitSize)); } private final double[] vec; // Don't access directly, use getElements() instead. private double[] getElements() { return VectorIntrinsics.maybeRebox(this).vec; } Double128Vector() { vec = new double[SPECIES.length()]; } Double128Vector(double[] v) { vec = v; } @Override public int length() { return LENGTH; } // Unary operator @Override Double128Vector uOp(FUnOp f) { double[] vec = getElements(); double[] res = new double[length()]; for (int i = 0; i < length(); i++) { res[i] = f.apply(i, vec[i]); } return new Double128Vector(res); } @Override Double128Vector uOp(VectorMask o, FUnOp f) { double[] vec = getElements(); double[] res = new double[length()]; boolean[] mbits = ((Double128Mask)o).getBits(); for (int i = 0; i < length(); i++) { res[i] = mbits[i] ? f.apply(i, vec[i]) : vec[i]; } return new Double128Vector(res); } // Binary operator @Override Double128Vector bOp(Vector o, FBinOp f) { double[] res = new double[length()]; double[] vec1 = this.getElements(); double[] vec2 = ((Double128Vector)o).getElements(); for (int i = 0; i < length(); i++) { res[i] = f.apply(i, vec1[i], vec2[i]); } return new Double128Vector(res); } @Override Double128Vector bOp(Vector o1, VectorMask o2, FBinOp f) { double[] res = new double[length()]; double[] vec1 = this.getElements(); double[] vec2 = ((Double128Vector)o1).getElements(); boolean[] mbits = ((Double128Mask)o2).getBits(); for (int i = 0; i < length(); i++) { res[i] = mbits[i] ? f.apply(i, vec1[i], vec2[i]) : vec1[i]; } return new Double128Vector(res); } // Trinary operator @Override Double128Vector tOp(Vector o1, Vector o2, FTriOp f) { double[] res = new double[length()]; double[] vec1 = this.getElements(); double[] vec2 = ((Double128Vector)o1).getElements(); double[] vec3 = ((Double128Vector)o2).getElements(); for (int i = 0; i < length(); i++) { res[i] = f.apply(i, vec1[i], vec2[i], vec3[i]); } return new Double128Vector(res); } @Override Double128Vector tOp(Vector o1, Vector o2, VectorMask o3, FTriOp f) { double[] res = new double[length()]; double[] vec1 = getElements(); double[] vec2 = ((Double128Vector)o1).getElements(); double[] vec3 = ((Double128Vector)o2).getElements(); boolean[] mbits = ((Double128Mask)o3).getBits(); for (int i = 0; i < length(); i++) { res[i] = mbits[i] ? f.apply(i, vec1[i], vec2[i], vec3[i]) : vec1[i]; } return new Double128Vector(res); } @Override double rOp(double v, FBinOp f) { double[] vec = getElements(); for (int i = 0; i < length(); i++) { v = f.apply(i, v, vec[i]); } return v; } @Override @ForceInline public Vector cast(VectorSpecies s) { Objects.requireNonNull(s); if (s.length() != LENGTH) throw new IllegalArgumentException("Vector length this species length differ"); return VectorIntrinsics.cast( Double128Vector.class, double.class, LENGTH, s.vectorType(), s.elementType(), LENGTH, this, s, (species, vector) -> vector.castDefault(species) ); } @SuppressWarnings("unchecked") @ForceInline private Vector castDefault(VectorSpecies s) { int limit = s.length(); Class stype = s.elementType(); if (stype == byte.class) { byte[] a = new byte[limit]; for (int i = 0; i < limit; i++) { a[i] = (byte) this.lane(i); } return (Vector) ByteVector.fromArray((VectorSpecies) s, a, 0); } else if (stype == short.class) { short[] a = new short[limit]; for (int i = 0; i < limit; i++) { a[i] = (short) this.lane(i); } return (Vector) ShortVector.fromArray((VectorSpecies) s, a, 0); } else if (stype == int.class) { int[] a = new int[limit]; for (int i = 0; i < limit; i++) { a[i] = (int) this.lane(i); } return (Vector) IntVector.fromArray((VectorSpecies) s, a, 0); } else if (stype == long.class) { long[] a = new long[limit]; for (int i = 0; i < limit; i++) { a[i] = (long) this.lane(i); } return (Vector) LongVector.fromArray((VectorSpecies) s, a, 0); } else if (stype == float.class) { float[] a = new float[limit]; for (int i = 0; i < limit; i++) { a[i] = (float) this.lane(i); } return (Vector) FloatVector.fromArray((VectorSpecies) s, a, 0); } else if (stype == double.class) { double[] a = new double[limit]; for (int i = 0; i < limit; i++) { a[i] = (double) this.lane(i); } return (Vector) DoubleVector.fromArray((VectorSpecies) s, a, 0); } else { throw new UnsupportedOperationException("Bad lane type for casting."); } } @Override @ForceInline @SuppressWarnings("unchecked") public Vector reinterpret(VectorSpecies s) { Objects.requireNonNull(s); if(s.elementType().equals(double.class)) { return (Vector) reshape((VectorSpecies)s); } if(s.bitSize() == bitSize()) { return reinterpretType(s); } return defaultReinterpret(s); } @ForceInline private Vector reinterpretType(VectorSpecies s) { Objects.requireNonNull(s); Class stype = s.elementType(); if (stype == byte.class) { return VectorIntrinsics.reinterpret( Double128Vector.class, double.class, LENGTH, Byte128Vector.class, byte.class, Byte128Vector.LENGTH, this, s, (species, vector) -> vector.defaultReinterpret(species) ); } else if (stype == short.class) { return VectorIntrinsics.reinterpret( Double128Vector.class, double.class, LENGTH, Short128Vector.class, short.class, Short128Vector.LENGTH, this, s, (species, vector) -> vector.defaultReinterpret(species) ); } else if (stype == int.class) { return VectorIntrinsics.reinterpret( Double128Vector.class, double.class, LENGTH, Int128Vector.class, int.class, Int128Vector.LENGTH, this, s, (species, vector) -> vector.defaultReinterpret(species) ); } else if (stype == long.class) { return VectorIntrinsics.reinterpret( Double128Vector.class, double.class, LENGTH, Long128Vector.class, long.class, Long128Vector.LENGTH, this, s, (species, vector) -> vector.defaultReinterpret(species) ); } else if (stype == float.class) { return VectorIntrinsics.reinterpret( Double128Vector.class, double.class, LENGTH, Float128Vector.class, float.class, Float128Vector.LENGTH, this, s, (species, vector) -> vector.defaultReinterpret(species) ); } else if (stype == double.class) { return VectorIntrinsics.reinterpret( Double128Vector.class, double.class, LENGTH, Double128Vector.class, double.class, Double128Vector.LENGTH, this, s, (species, vector) -> vector.defaultReinterpret(species) ); } else { throw new UnsupportedOperationException("Bad lane type for casting."); } } @Override @ForceInline public DoubleVector reshape(VectorSpecies s) { Objects.requireNonNull(s); if (s.bitSize() == 64 && (s.vectorType() == Double64Vector.class)) { return VectorIntrinsics.reinterpret( Double128Vector.class, double.class, LENGTH, Double64Vector.class, double.class, Double64Vector.LENGTH, this, s, (species, vector) -> (DoubleVector) vector.defaultReinterpret(species) ); } else if (s.bitSize() == 128 && (s.vectorType() == Double128Vector.class)) { return VectorIntrinsics.reinterpret( Double128Vector.class, double.class, LENGTH, Double128Vector.class, double.class, Double128Vector.LENGTH, this, s, (species, vector) -> (DoubleVector) vector.defaultReinterpret(species) ); } else if (s.bitSize() == 256 && (s.vectorType() == Double256Vector.class)) { return VectorIntrinsics.reinterpret( Double128Vector.class, double.class, LENGTH, Double256Vector.class, double.class, Double256Vector.LENGTH, this, s, (species, vector) -> (DoubleVector) vector.defaultReinterpret(species) ); } else if (s.bitSize() == 512 && (s.vectorType() == Double512Vector.class)) { return VectorIntrinsics.reinterpret( Double128Vector.class, double.class, LENGTH, Double512Vector.class, double.class, Double512Vector.LENGTH, this, s, (species, vector) -> (DoubleVector) vector.defaultReinterpret(species) ); } else if ((s.bitSize() > 0) && (s.bitSize() <= 2048) && (s.bitSize() % 128 == 0) && (s.vectorType() == DoubleMaxVector.class)) { return VectorIntrinsics.reinterpret( Double128Vector.class, double.class, LENGTH, DoubleMaxVector.class, double.class, DoubleMaxVector.LENGTH, this, s, (species, vector) -> (DoubleVector) vector.defaultReinterpret(species) ); } else { throw new InternalError("Unimplemented size"); } } // Binary operations with scalars @Override @ForceInline public DoubleVector add(double o) { return add((Double128Vector)DoubleVector.broadcast(SPECIES, o)); } @Override @ForceInline public DoubleVector add(double o, VectorMask m) { return add((Double128Vector)DoubleVector.broadcast(SPECIES, o), m); } @Override @ForceInline public DoubleVector sub(double o) { return sub((Double128Vector)DoubleVector.broadcast(SPECIES, o)); } @Override @ForceInline public DoubleVector sub(double o, VectorMask m) { return sub((Double128Vector)DoubleVector.broadcast(SPECIES, o), m); } @Override @ForceInline public DoubleVector mul(double o) { return mul((Double128Vector)DoubleVector.broadcast(SPECIES, o)); } @Override @ForceInline public DoubleVector mul(double o, VectorMask m) { return mul((Double128Vector)DoubleVector.broadcast(SPECIES, o), m); } @Override @ForceInline public DoubleVector min(double o) { return min((Double128Vector)DoubleVector.broadcast(SPECIES, o)); } @Override @ForceInline public DoubleVector max(double o) { return max((Double128Vector)DoubleVector.broadcast(SPECIES, o)); } @Override @ForceInline public VectorMask equal(double o) { return equal((Double128Vector)DoubleVector.broadcast(SPECIES, o)); } @Override @ForceInline public VectorMask notEqual(double o) { return notEqual((Double128Vector)DoubleVector.broadcast(SPECIES, o)); } @Override @ForceInline public VectorMask lessThan(double o) { return lessThan((Double128Vector)DoubleVector.broadcast(SPECIES, o)); } @Override @ForceInline public VectorMask lessThanEq(double o) { return lessThanEq((Double128Vector)DoubleVector.broadcast(SPECIES, o)); } @Override @ForceInline public VectorMask greaterThan(double o) { return greaterThan((Double128Vector)DoubleVector.broadcast(SPECIES, o)); } @Override @ForceInline public VectorMask greaterThanEq(double o) { return greaterThanEq((Double128Vector)DoubleVector.broadcast(SPECIES, o)); } @Override @ForceInline public DoubleVector blend(double o, VectorMask m) { return blend((Double128Vector)DoubleVector.broadcast(SPECIES, o), m); } @Override @ForceInline public DoubleVector div(double o) { return div((Double128Vector)DoubleVector.broadcast(SPECIES, o)); } @Override @ForceInline public DoubleVector div(double o, VectorMask m) { return div((Double128Vector)DoubleVector.broadcast(SPECIES, o), m); } @Override @ForceInline public Double128Vector div(Vector v, VectorMask m) { return blend(div(v), m); } @Override @ForceInline public DoubleVector atan2(double o) { return atan2((Double128Vector)DoubleVector.broadcast(SPECIES, o)); } @Override @ForceInline public DoubleVector atan2(double o, VectorMask m) { return atan2((Double128Vector)DoubleVector.broadcast(SPECIES, o), m); } @Override @ForceInline public DoubleVector pow(double o) { return pow((Double128Vector)DoubleVector.broadcast(SPECIES, o)); } @Override @ForceInline public DoubleVector pow(double o, VectorMask m) { return pow((Double128Vector)DoubleVector.broadcast(SPECIES, o), m); } @Override @ForceInline public DoubleVector fma(double o1, double o2) { return fma((Double128Vector)DoubleVector.broadcast(SPECIES, o1), (Double128Vector)DoubleVector.broadcast(SPECIES, o2)); } @Override @ForceInline public DoubleVector fma(double o1, double o2, VectorMask m) { return fma((Double128Vector)DoubleVector.broadcast(SPECIES, o1), (Double128Vector)DoubleVector.broadcast(SPECIES, o2), m); } @Override @ForceInline public DoubleVector hypot(double o) { return hypot((Double128Vector)DoubleVector.broadcast(SPECIES, o)); } @Override @ForceInline public DoubleVector hypot(double o, VectorMask m) { return hypot((Double128Vector)DoubleVector.broadcast(SPECIES, o), m); } // Unary operations @ForceInline @Override public Double128Vector neg(VectorMask m) { return blend(neg(), m); } @Override @ForceInline public Double128Vector abs() { return VectorIntrinsics.unaryOp( VECTOR_OP_ABS, Double128Vector.class, double.class, LENGTH, this, v1 -> v1.uOp((i, a) -> (double) Math.abs(a))); } @ForceInline @Override public Double128Vector abs(VectorMask m) { return blend(abs(), m); } @Override @ForceInline public Double128Vector neg() { return VectorIntrinsics.unaryOp( VECTOR_OP_NEG, Double128Vector.class, double.class, LENGTH, this, v1 -> v1.uOp((i, a) -> (double) -a)); } @Override @ForceInline public Double128Vector div(Vector o) { Objects.requireNonNull(o); Double128Vector v = (Double128Vector)o; return VectorIntrinsics.binaryOp( VECTOR_OP_DIV, Double128Vector.class, double.class, LENGTH, this, v, (v1, v2) -> v1.bOp(v2, (i, a, b) -> (double)(a / b))); } @Override @ForceInline public Double128Vector sqrt() { return VectorIntrinsics.unaryOp( VECTOR_OP_SQRT, Double128Vector.class, double.class, LENGTH, this, v1 -> v1.uOp((i, a) -> (double) Math.sqrt((double) a))); } @Override @ForceInline public Double128Vector exp() { return (Double128Vector) VectorIntrinsics.unaryOp( VECTOR_OP_EXP, Double128Vector.class, double.class, LENGTH, this, v1 -> ((Double128Vector)v1).uOp((i, a) -> (double) Math.exp((double) a))); } @Override @ForceInline public Double128Vector log1p() { return (Double128Vector) VectorIntrinsics.unaryOp( VECTOR_OP_LOG1P, Double128Vector.class, double.class, LENGTH, this, v1 -> ((Double128Vector)v1).uOp((i, a) -> (double) Math.log1p((double) a))); } @Override @ForceInline public Double128Vector log() { return (Double128Vector) VectorIntrinsics.unaryOp( VECTOR_OP_LOG, Double128Vector.class, double.class, LENGTH, this, v1 -> ((Double128Vector)v1).uOp((i, a) -> (double) Math.log((double) a))); } @Override @ForceInline public Double128Vector log10() { return (Double128Vector) VectorIntrinsics.unaryOp( VECTOR_OP_LOG10, Double128Vector.class, double.class, LENGTH, this, v1 -> ((Double128Vector)v1).uOp((i, a) -> (double) Math.log10((double) a))); } @Override @ForceInline public Double128Vector expm1() { return (Double128Vector) VectorIntrinsics.unaryOp( VECTOR_OP_EXPM1, Double128Vector.class, double.class, LENGTH, this, v1 -> ((Double128Vector)v1).uOp((i, a) -> (double) Math.expm1((double) a))); } @Override @ForceInline public Double128Vector cbrt() { return (Double128Vector) VectorIntrinsics.unaryOp( VECTOR_OP_CBRT, Double128Vector.class, double.class, LENGTH, this, v1 -> ((Double128Vector)v1).uOp((i, a) -> (double) Math.cbrt((double) a))); } @Override @ForceInline public Double128Vector sin() { return (Double128Vector) VectorIntrinsics.unaryOp( VECTOR_OP_SIN, Double128Vector.class, double.class, LENGTH, this, v1 -> ((Double128Vector)v1).uOp((i, a) -> (double) Math.sin((double) a))); } @Override @ForceInline public Double128Vector cos() { return (Double128Vector) VectorIntrinsics.unaryOp( VECTOR_OP_COS, Double128Vector.class, double.class, LENGTH, this, v1 -> ((Double128Vector)v1).uOp((i, a) -> (double) Math.cos((double) a))); } @Override @ForceInline public Double128Vector tan() { return (Double128Vector) VectorIntrinsics.unaryOp( VECTOR_OP_TAN, Double128Vector.class, double.class, LENGTH, this, v1 -> ((Double128Vector)v1).uOp((i, a) -> (double) Math.tan((double) a))); } @Override @ForceInline public Double128Vector asin() { return (Double128Vector) VectorIntrinsics.unaryOp( VECTOR_OP_ASIN, Double128Vector.class, double.class, LENGTH, this, v1 -> ((Double128Vector)v1).uOp((i, a) -> (double) Math.asin((double) a))); } @Override @ForceInline public Double128Vector acos() { return (Double128Vector) VectorIntrinsics.unaryOp( VECTOR_OP_ACOS, Double128Vector.class, double.class, LENGTH, this, v1 -> ((Double128Vector)v1).uOp((i, a) -> (double) Math.acos((double) a))); } @Override @ForceInline public Double128Vector atan() { return (Double128Vector) VectorIntrinsics.unaryOp( VECTOR_OP_ATAN, Double128Vector.class, double.class, LENGTH, this, v1 -> ((Double128Vector)v1).uOp((i, a) -> (double) Math.atan((double) a))); } @Override @ForceInline public Double128Vector sinh() { return (Double128Vector) VectorIntrinsics.unaryOp( VECTOR_OP_SINH, Double128Vector.class, double.class, LENGTH, this, v1 -> ((Double128Vector)v1).uOp((i, a) -> (double) Math.sinh((double) a))); } @Override @ForceInline public Double128Vector cosh() { return (Double128Vector) VectorIntrinsics.unaryOp( VECTOR_OP_COSH, Double128Vector.class, double.class, LENGTH, this, v1 -> ((Double128Vector)v1).uOp((i, a) -> (double) Math.cosh((double) a))); } @Override @ForceInline public Double128Vector tanh() { return (Double128Vector) VectorIntrinsics.unaryOp( VECTOR_OP_TANH, Double128Vector.class, double.class, LENGTH, this, v1 -> ((Double128Vector)v1).uOp((i, a) -> (double) Math.tanh((double) a))); } @Override @ForceInline public Double128Vector pow(Vector o) { Objects.requireNonNull(o); Double128Vector v = (Double128Vector)o; return (Double128Vector) VectorIntrinsics.binaryOp( VECTOR_OP_POW, Double128Vector.class, double.class, LENGTH, this, v, (v1, v2) -> ((Double128Vector)v1).bOp(v2, (i, a, b) -> (double)(Math.pow(a,b)))); } @Override @ForceInline public Double128Vector hypot(Vector o) { Objects.requireNonNull(o); Double128Vector v = (Double128Vector)o; return (Double128Vector) VectorIntrinsics.binaryOp( VECTOR_OP_HYPOT, Double128Vector.class, double.class, LENGTH, this, v, (v1, v2) -> ((Double128Vector)v1).bOp(v2, (i, a, b) -> (double)(Math.hypot(a,b)))); } @Override @ForceInline public Double128Vector atan2(Vector o) { Objects.requireNonNull(o); Double128Vector v = (Double128Vector)o; return (Double128Vector) VectorIntrinsics.binaryOp( VECTOR_OP_ATAN2, Double128Vector.class, double.class, LENGTH, this, v, (v1, v2) -> ((Double128Vector)v1).bOp(v2, (i, a, b) -> (double)(Math.atan2(a,b)))); } // Binary operations @Override @ForceInline public Double128Vector add(Vector o) { Objects.requireNonNull(o); Double128Vector v = (Double128Vector)o; return VectorIntrinsics.binaryOp( VECTOR_OP_ADD, Double128Vector.class, double.class, LENGTH, this, v, (v1, v2) -> v1.bOp(v2, (i, a, b) -> (double)(a + b))); } @Override @ForceInline public Double128Vector add(Vector v, VectorMask m) { return blend(add(v), m); } @Override @ForceInline public Double128Vector sub(Vector o) { Objects.requireNonNull(o); Double128Vector v = (Double128Vector)o; return VectorIntrinsics.binaryOp( VECTOR_OP_SUB, Double128Vector.class, double.class, LENGTH, this, v, (v1, v2) -> v1.bOp(v2, (i, a, b) -> (double)(a - b))); } @Override @ForceInline public Double128Vector sub(Vector v, VectorMask m) { return blend(sub(v), m); } @Override @ForceInline public Double128Vector mul(Vector o) { Objects.requireNonNull(o); Double128Vector v = (Double128Vector)o; return VectorIntrinsics.binaryOp( VECTOR_OP_MUL, Double128Vector.class, double.class, LENGTH, this, v, (v1, v2) -> v1.bOp(v2, (i, a, b) -> (double)(a * b))); } @Override @ForceInline public Double128Vector mul(Vector v, VectorMask m) { return blend(mul(v), m); } @Override @ForceInline public Double128Vector min(Vector o) { Objects.requireNonNull(o); Double128Vector v = (Double128Vector)o; return (Double128Vector) VectorIntrinsics.binaryOp( VECTOR_OP_MIN, Double128Vector.class, double.class, LENGTH, this, v, (v1, v2) -> v1.bOp(v2, (i, a, b) -> (double) Math.min(a, b))); } @Override @ForceInline public Double128Vector min(Vector v, VectorMask m) { return blend(min(v), m); } @Override @ForceInline public Double128Vector max(Vector o) { Objects.requireNonNull(o); Double128Vector v = (Double128Vector)o; return VectorIntrinsics.binaryOp( VECTOR_OP_MAX, Double128Vector.class, double.class, LENGTH, this, v, (v1, v2) -> v1.bOp(v2, (i, a, b) -> (double) Math.max(a, b))); } @Override @ForceInline public Double128Vector max(Vector v, VectorMask m) { return blend(max(v), m); } // Ternary operations @Override @ForceInline public Double128Vector fma(Vector o1, Vector o2) { Objects.requireNonNull(o1); Objects.requireNonNull(o2); Double128Vector v1 = (Double128Vector)o1; Double128Vector v2 = (Double128Vector)o2; return VectorIntrinsics.ternaryOp( VECTOR_OP_FMA, Double128Vector.class, double.class, LENGTH, this, v1, v2, (w1, w2, w3) -> w1.tOp(w2, w3, (i, a, b, c) -> Math.fma(a, b, c))); } // Type specific horizontal reductions @Override @ForceInline public double addLanes() { long bits = (long) VectorIntrinsics.reductionCoerced( VECTOR_OP_ADD, Double128Vector.class, double.class, LENGTH, this, v -> { double r = v.rOp((double) 0, (i, a, b) -> (double) (a + b)); return (long)Double.doubleToLongBits(r); }); return Double.longBitsToDouble(bits); } @Override @ForceInline public double mulLanes() { long bits = (long) VectorIntrinsics.reductionCoerced( VECTOR_OP_MUL, Double128Vector.class, double.class, LENGTH, this, v -> { double r = v.rOp((double) 1, (i, a, b) -> (double) (a * b)); return (long)Double.doubleToLongBits(r); }); return Double.longBitsToDouble(bits); } @Override @ForceInline public double minLanes() { long bits = (long) VectorIntrinsics.reductionCoerced( VECTOR_OP_MIN, Double128Vector.class, double.class, LENGTH, this, v -> { double r = v.rOp(Double.POSITIVE_INFINITY , (i, a, b) -> (double) Math.min(a, b)); return (long)Double.doubleToLongBits(r); }); return Double.longBitsToDouble(bits); } @Override @ForceInline public double maxLanes() { long bits = (long) VectorIntrinsics.reductionCoerced( VECTOR_OP_MAX, Double128Vector.class, double.class, LENGTH, this, v -> { double r = v.rOp(Double.NEGATIVE_INFINITY, (i, a, b) -> (double) Math.max(a, b)); return (long)Double.doubleToLongBits(r); }); return Double.longBitsToDouble(bits); } @Override @ForceInline public double addLanes(VectorMask m) { return DoubleVector.broadcast(SPECIES, (double) 0).blend(this, m).addLanes(); } @Override @ForceInline public double mulLanes(VectorMask m) { return DoubleVector.broadcast(SPECIES, (double) 1).blend(this, m).mulLanes(); } @Override @ForceInline public double minLanes(VectorMask m) { return DoubleVector.broadcast(SPECIES, Double.MAX_VALUE).blend(this, m).minLanes(); } @Override @ForceInline public double maxLanes(VectorMask m) { return DoubleVector.broadcast(SPECIES, Double.MIN_VALUE).blend(this, m).maxLanes(); } @Override @ForceInline public VectorShuffle toShuffle() { double[] a = toArray(); int[] sa = new int[a.length]; for (int i = 0; i < a.length; i++) { sa[i] = (int) a[i]; } return VectorShuffle.fromArray(SPECIES, sa, 0); } // Memory operations private static final int ARRAY_SHIFT = 31 - Integer.numberOfLeadingZeros(Unsafe.ARRAY_DOUBLE_INDEX_SCALE); private static final int BOOLEAN_ARRAY_SHIFT = 31 - Integer.numberOfLeadingZeros(Unsafe.ARRAY_BOOLEAN_INDEX_SCALE); @Override @ForceInline public void intoArray(double[] a, int ix) { Objects.requireNonNull(a); ix = VectorIntrinsics.checkIndex(ix, a.length, LENGTH); VectorIntrinsics.store(Double128Vector.class, double.class, LENGTH, a, (((long) ix) << ARRAY_SHIFT) + Unsafe.ARRAY_DOUBLE_BASE_OFFSET, this, a, ix, (arr, idx, v) -> v.forEach((i, e) -> arr[idx + i] = e)); } @Override @ForceInline public final void intoArray(double[] a, int ax, VectorMask m) { DoubleVector oldVal = DoubleVector.fromArray(SPECIES, a, ax); DoubleVector newVal = oldVal.blend(this, m); newVal.intoArray(a, ax); } @Override @ForceInline public void intoArray(double[] a, int ix, int[] b, int iy) { Objects.requireNonNull(a); Objects.requireNonNull(b); // Index vector: vix[0:n] = i -> ix + indexMap[iy + i] IntVector vix = IntVector.fromArray(INDEX_SPECIES, b, iy).add(ix); vix = VectorIntrinsics.checkIndex(vix, a.length); VectorIntrinsics.storeWithMap(Double128Vector.class, double.class, LENGTH, Int64Vector.class, a, Unsafe.ARRAY_DOUBLE_BASE_OFFSET, vix, this, a, ix, b, iy, (arr, idx, v, indexMap, idy) -> v.forEach((i, e) -> arr[idx+indexMap[idy+i]] = e)); } @Override @ForceInline public final void intoArray(double[] a, int ax, VectorMask m, int[] b, int iy) { // @@@ This can result in out of bounds errors for unset mask lanes DoubleVector oldVal = DoubleVector.fromArray(SPECIES, a, ax, b, iy); DoubleVector newVal = oldVal.blend(this, m); newVal.intoArray(a, ax, b, iy); } @Override @ForceInline public void intoByteArray(byte[] a, int ix) { Objects.requireNonNull(a); ix = VectorIntrinsics.checkIndex(ix, a.length, bitSize() / Byte.SIZE); VectorIntrinsics.store(Double128Vector.class, double.class, LENGTH, a, ((long) ix) + Unsafe.ARRAY_BYTE_BASE_OFFSET, this, a, ix, (c, idx, v) -> { ByteBuffer bbc = ByteBuffer.wrap(c, idx, c.length - idx).order(ByteOrder.nativeOrder()); DoubleBuffer tb = bbc.asDoubleBuffer(); v.forEach((i, e) -> tb.put(e)); }); } @Override @ForceInline public final void intoByteArray(byte[] a, int ix, VectorMask m) { Double128Vector oldVal = (Double128Vector) DoubleVector.fromByteArray(SPECIES, a, ix); Double128Vector newVal = oldVal.blend(this, m); newVal.intoByteArray(a, ix); } @Override @ForceInline public void intoByteBuffer(ByteBuffer bb, int ix) { if (bb.order() != ByteOrder.nativeOrder()) { throw new IllegalArgumentException(); } if (bb.isReadOnly()) { throw new ReadOnlyBufferException(); } ix = VectorIntrinsics.checkIndex(ix, bb.limit(), bitSize() / Byte.SIZE); VectorIntrinsics.store(Double128Vector.class, double.class, LENGTH, U.getReference(bb, BYTE_BUFFER_HB), ix + U.getLong(bb, BUFFER_ADDRESS), this, bb, ix, (c, idx, v) -> { ByteBuffer bbc = c.duplicate().position(idx).order(ByteOrder.nativeOrder()); DoubleBuffer tb = bbc.asDoubleBuffer(); v.forEach((i, e) -> tb.put(e)); }); } @Override @ForceInline public void intoByteBuffer(ByteBuffer bb, int ix, VectorMask m) { Double128Vector oldVal = (Double128Vector) DoubleVector.fromByteBuffer(SPECIES, bb, ix); Double128Vector newVal = oldVal.blend(this, m); newVal.intoByteBuffer(bb, ix); } // @Override public String toString() { return Arrays.toString(getElements()); } @Override public boolean equals(Object o) { if (this == o) return true; if (o == null || this.getClass() != o.getClass()) return false; Double128Vector that = (Double128Vector) o; return this.equal(that).allTrue(); } @Override public int hashCode() { return Arrays.hashCode(vec); } // Binary test @Override Double128Mask bTest(Vector o, FBinTest f) { double[] vec1 = getElements(); double[] vec2 = ((Double128Vector)o).getElements(); boolean[] bits = new boolean[length()]; for (int i = 0; i < length(); i++){ bits[i] = f.apply(i, vec1[i], vec2[i]); } return new Double128Mask(bits); } // Comparisons @Override @ForceInline public Double128Mask equal(Vector o) { Objects.requireNonNull(o); Double128Vector v = (Double128Vector)o; return VectorIntrinsics.compare( BT_eq, Double128Vector.class, Double128Mask.class, double.class, LENGTH, this, v, (v1, v2) -> v1.bTest(v2, (i, a, b) -> a == b)); } @Override @ForceInline public Double128Mask notEqual(Vector o) { Objects.requireNonNull(o); Double128Vector v = (Double128Vector)o; return VectorIntrinsics.compare( BT_ne, Double128Vector.class, Double128Mask.class, double.class, LENGTH, this, v, (v1, v2) -> v1.bTest(v2, (i, a, b) -> a != b)); } @Override @ForceInline public Double128Mask lessThan(Vector o) { Objects.requireNonNull(o); Double128Vector v = (Double128Vector)o; return VectorIntrinsics.compare( BT_lt, Double128Vector.class, Double128Mask.class, double.class, LENGTH, this, v, (v1, v2) -> v1.bTest(v2, (i, a, b) -> a < b)); } @Override @ForceInline public Double128Mask lessThanEq(Vector o) { Objects.requireNonNull(o); Double128Vector v = (Double128Vector)o; return VectorIntrinsics.compare( BT_le, Double128Vector.class, Double128Mask.class, double.class, LENGTH, this, v, (v1, v2) -> v1.bTest(v2, (i, a, b) -> a <= b)); } @Override @ForceInline public Double128Mask greaterThan(Vector o) { Objects.requireNonNull(o); Double128Vector v = (Double128Vector)o; return (Double128Mask) VectorIntrinsics.compare( BT_gt, Double128Vector.class, Double128Mask.class, double.class, LENGTH, this, v, (v1, v2) -> v1.bTest(v2, (i, a, b) -> a > b)); } @Override @ForceInline public Double128Mask greaterThanEq(Vector o) { Objects.requireNonNull(o); Double128Vector v = (Double128Vector)o; return VectorIntrinsics.compare( BT_ge, Double128Vector.class, Double128Mask.class, double.class, LENGTH, this, v, (v1, v2) -> v1.bTest(v2, (i, a, b) -> a >= b)); } // Foreach @Override void forEach(FUnCon f) { double[] vec = getElements(); for (int i = 0; i < length(); i++) { f.apply(i, vec[i]); } } @Override void forEach(VectorMask o, FUnCon f) { boolean[] mbits = ((Double128Mask)o).getBits(); forEach((i, a) -> { if (mbits[i]) { f.apply(i, a); } }); } Long128Vector toBits() { double[] vec = getElements(); long[] res = new long[this.species().length()]; for(int i = 0; i < this.species().length(); i++){ res[i] = Double.doubleToLongBits(vec[i]); } return new Long128Vector(res); } @Override @ForceInline public Double128Vector rotateLanesLeft(int j) { int L = length(); if (j < 0) { throw new IllegalArgumentException("Index " + j + " must be zero or positive"); } else { j = j & (L-1); VectorShuffle PermMask = VectorShuffle.shuffleIota(SPECIES, L - j); return this.rearrange(PermMask); } } @Override @ForceInline public Double128Vector rotateLanesRight(int j) { int L = length(); if (j < 0) { throw new IllegalArgumentException("Index " + j + " must be zero or positive"); } else { j = j & (L-1); VectorShuffle PermMask = VectorShuffle.shuffleIota(SPECIES, j); return this.rearrange(PermMask); } } @Override @ForceInline @SuppressWarnings("unchecked") public Double128Vector shiftLanesLeft(int j) { int L = length(); if (j < 0) { throw new IllegalArgumentException("Index " + j + " must be zero or positive"); } else if ( j >= L ) { return ZERO; } else { Double128Shuffle Iota = (Double128Shuffle)(VectorShuffle.shuffleIota(SPECIES, L-j)); VectorMask BlendMask = Iota.toVector().lessThan(Double128Vector.broadcast(SPECIES, (double)(L-j))); Iota = (Double128Shuffle)(VectorShuffle.shuffleIota(SPECIES, L -j)); return ZERO.blend(this.rearrange(Iota),BlendMask); } } @Override @ForceInline @SuppressWarnings("unchecked") public Double128Vector shiftLanesRight(int j) { int L = length(); if (j < 0) { throw new IllegalArgumentException("Index " + j + " must be zero or positive"); } else if ( j >= L ) { return ZERO; } else { Double128Shuffle Iota = (Double128Shuffle)(VectorShuffle.shuffleIota(SPECIES, j)); VectorMask BlendMask = Iota.toVector().greaterThanEq(Double128Vector.broadcast(SPECIES, (double)(j))); Iota = (Double128Shuffle)(VectorShuffle.shuffleIota(SPECIES, j)); return ZERO.blend(this.rearrange(Iota),BlendMask); } } @Override @ForceInline public Double128Vector rearrange(Vector v, VectorShuffle s, VectorMask m) { return this.rearrange(s).blend(v.rearrange(s), m); } @Override @ForceInline public Double128Vector rearrange(VectorShuffle o1) { Objects.requireNonNull(o1); Double128Shuffle s = (Double128Shuffle)o1; return VectorIntrinsics.rearrangeOp( Double128Vector.class, Double128Shuffle.class, double.class, LENGTH, this, s, (v1, s_) -> v1.uOp((i, a) -> { int ei = s_.lane(i); return v1.lane(ei); })); } @Override @ForceInline public Double128Vector blend(Vector o1, VectorMask o2) { Objects.requireNonNull(o1); Objects.requireNonNull(o2); Double128Vector v = (Double128Vector)o1; Double128Mask m = (Double128Mask)o2; return VectorIntrinsics.blend( Double128Vector.class, Double128Mask.class, double.class, LENGTH, this, v, m, (v1, v2, m_) -> v1.bOp(v2, (i, a, b) -> m_.lane(i) ? b : a)); } // Accessors @Override public double lane(int i) { if (i < 0 || i >= LENGTH) { throw new IllegalArgumentException("Index " + i + " must be zero or positive, and less than " + LENGTH); } long bits = (long) VectorIntrinsics.extract( Double128Vector.class, double.class, LENGTH, this, i, (vec, ix) -> { double[] vecarr = vec.getElements(); return (long)Double.doubleToLongBits(vecarr[ix]); }); return Double.longBitsToDouble(bits); } @Override public Double128Vector with(int i, double e) { if (i < 0 || i >= LENGTH) { throw new IllegalArgumentException("Index " + i + " must be zero or positive, and less than " + LENGTH); } return VectorIntrinsics.insert( Double128Vector.class, double.class, LENGTH, this, i, (long)Double.doubleToLongBits(e), (v, ix, bits) -> { double[] res = v.getElements().clone(); res[ix] = Double.longBitsToDouble((long)bits); return new Double128Vector(res); }); } // Mask static final class Double128Mask extends AbstractMask { static final Double128Mask TRUE_MASK = new Double128Mask(true); static final Double128Mask FALSE_MASK = new Double128Mask(false); private final boolean[] bits; // Don't access directly, use getBits() instead. public Double128Mask(boolean[] bits) { this(bits, 0); } public Double128Mask(boolean[] bits, int offset) { boolean[] a = new boolean[species().length()]; for (int i = 0; i < a.length; i++) { a[i] = bits[offset + i]; } this.bits = a; } public Double128Mask(boolean val) { boolean[] bits = new boolean[species().length()]; Arrays.fill(bits, val); this.bits = bits; } boolean[] getBits() { return VectorIntrinsics.maybeRebox(this).bits; } @Override Double128Mask uOp(MUnOp f) { boolean[] res = new boolean[species().length()]; boolean[] bits = getBits(); for (int i = 0; i < species().length(); i++) { res[i] = f.apply(i, bits[i]); } return new Double128Mask(res); } @Override Double128Mask bOp(VectorMask o, MBinOp f) { boolean[] res = new boolean[species().length()]; boolean[] bits = getBits(); boolean[] mbits = ((Double128Mask)o).getBits(); for (int i = 0; i < species().length(); i++) { res[i] = f.apply(i, bits[i], mbits[i]); } return new Double128Mask(res); } @Override public VectorSpecies species() { return SPECIES; } @Override public Double128Vector toVector() { double[] res = new double[species().length()]; boolean[] bits = getBits(); for (int i = 0; i < species().length(); i++) { // -1 will result in the most significant bit being set in // addition to some or all other bits res[i] = (double) (bits[i] ? -1 : 0); } return new Double128Vector(res); } @Override @ForceInline @SuppressWarnings("unchecked") public VectorMask cast(VectorSpecies species) { if (length() != species.length()) throw new IllegalArgumentException("VectorMask length and species length differ"); Class stype = species.elementType(); boolean [] maskArray = toArray(); if (stype == byte.class) { return (VectorMask ) new Byte128Vector.Byte128Mask(maskArray); } else if (stype == short.class) { return (VectorMask ) new Short128Vector.Short128Mask(maskArray); } else if (stype == int.class) { return (VectorMask ) new Int128Vector.Int128Mask(maskArray); } else if (stype == long.class) { return (VectorMask ) new Long128Vector.Long128Mask(maskArray); } else if (stype == float.class) { return (VectorMask ) new Float128Vector.Float128Mask(maskArray); } else if (stype == double.class) { return (VectorMask ) new Double128Vector.Double128Mask(maskArray); } else { throw new UnsupportedOperationException("Bad lane type for casting."); } } // Unary operations @Override @ForceInline public Double128Mask not() { return (Double128Mask) VectorIntrinsics.unaryOp( VECTOR_OP_NOT, Double128Mask.class, long.class, LENGTH, this, (m1) -> m1.uOp((i, a) -> !a)); } // Binary operations @Override @ForceInline public Double128Mask and(VectorMask o) { Objects.requireNonNull(o); Double128Mask m = (Double128Mask)o; return VectorIntrinsics.binaryOp(VECTOR_OP_AND, Double128Mask.class, long.class, LENGTH, this, m, (m1, m2) -> m1.bOp(m2, (i, a, b) -> a & b)); } @Override @ForceInline public Double128Mask or(VectorMask o) { Objects.requireNonNull(o); Double128Mask m = (Double128Mask)o; return VectorIntrinsics.binaryOp(VECTOR_OP_OR, Double128Mask.class, long.class, LENGTH, this, m, (m1, m2) -> m1.bOp(m2, (i, a, b) -> a | b)); } // Reductions @Override @ForceInline public boolean anyTrue() { return VectorIntrinsics.test(BT_ne, Double128Mask.class, long.class, LENGTH, this, this, (m, __) -> anyTrueHelper(((Double128Mask)m).getBits())); } @Override @ForceInline public boolean allTrue() { return VectorIntrinsics.test(BT_overflow, Double128Mask.class, long.class, LENGTH, this, VectorMask.maskAllTrue(species()), (m, __) -> allTrueHelper(((Double128Mask)m).getBits())); } } // Shuffle static final class Double128Shuffle extends AbstractShuffle { Double128Shuffle(byte[] reorder) { super(reorder); } public Double128Shuffle(int[] reorder) { super(reorder); } public Double128Shuffle(int[] reorder, int i) { super(reorder, i); } public Double128Shuffle(IntUnaryOperator f) { super(f); } @Override public VectorSpecies species() { return SPECIES; } private DoubleVector toVector_helper() { double[] va = new double[SPECIES.length()]; for (int i = 0; i < va.length; i++) { va[i] = (double) lane(i); } return DoubleVector.fromArray(SPECIES, va, 0); } @Override @ForceInline public DoubleVector toVector() { return VectorIntrinsics.shuffleToVector(Double128Vector.class, double.class, Double128Shuffle.class, this, SPECIES.length(), (s) -> (((Double128Shuffle)(s)).toVector_helper())); } @Override @ForceInline @SuppressWarnings("unchecked") public VectorShuffle cast(VectorSpecies species) { if (length() != species.length()) throw new IllegalArgumentException("Shuffle length and species length differ"); Class stype = species.elementType(); int [] shuffleArray = toArray(); if (stype == byte.class) { return (VectorShuffle) new Byte128Vector.Byte128Shuffle(shuffleArray); } else if (stype == short.class) { return (VectorShuffle) new Short128Vector.Short128Shuffle(shuffleArray); } else if (stype == int.class) { return (VectorShuffle) new Int128Vector.Int128Shuffle(shuffleArray); } else if (stype == long.class) { return (VectorShuffle) new Long128Vector.Long128Shuffle(shuffleArray); } else if (stype == float.class) { return (VectorShuffle) new Float128Vector.Float128Shuffle(shuffleArray); } else if (stype == double.class) { return (VectorShuffle) new Double128Vector.Double128Shuffle(shuffleArray); } else { throw new UnsupportedOperationException("Bad lane type for casting."); } } @Override public Double128Shuffle rearrange(VectorShuffle o) { Double128Shuffle s = (Double128Shuffle) o; byte[] r = new byte[reorder.length]; for (int i = 0; i < reorder.length; i++) { r[i] = reorder[s.reorder[i]]; } return new Double128Shuffle(r); } } // VectorSpecies @Override public VectorSpecies species() { return SPECIES; } }