diff --git a/kernel/power/KERNEL.POWER8 b/kernel/power/KERNEL.POWER8 index 99bb38096..8fbceb052 100644 --- a/kernel/power/KERNEL.POWER8 +++ b/kernel/power/KERNEL.POWER8 @@ -135,15 +135,15 @@ DROTKERNEL = drot.c #CROTKERNEL = ../arm/zrot.c #ZROTKERNEL = ../arm/zrot.c # -#SSCALKERNEL = ../arm/scal.c +SSCALKERNEL = sscal.c DSCALKERNEL = dscal.c #CSCALKERNEL = ../arm/zscal.c -#ZSCALKERNEL = zscal.c +ZSCALKERNEL = zscal.c # SSWAPKERNEL = sswap.c DSWAPKERNEL = dswap.c -#CSWAPKERNEL = ../arm/zswap.c -#ZSWAPKERNEL = ../arm/zswap.c +CSWAPKERNEL = cswap.c +ZSWAPKERNEL = zswap.c # #SGEMVNKERNEL = ../arm/gemv_n.c diff --git a/kernel/power/cswap.c b/kernel/power/cswap.c new file mode 100644 index 000000000..da97c896e --- /dev/null +++ b/kernel/power/cswap.c @@ -0,0 +1,175 @@ +/*************************************************************************** +Copyright (c) 2013-2016, The OpenBLAS Project +All rights reserved. +Redistribution and use in source and binary forms, with or without +modification, are permitted provided that the following conditions are +met: +1. Redistributions of source code must retain the above copyright +notice, this list of conditions and the following disclaimer. +2. Redistributions in binary form must reproduce the above copyright +notice, this list of conditions and the following disclaimer in +the documentation and/or other materials provided with the +distribution. +3. Neither the name of the OpenBLAS project nor the names of +its contributors may be used to endorse or promote products +derived from this software without specific prior written permission. +THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" +AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE +IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE +ARE DISCLAIMED. IN NO EVENT SHALL THE OPENBLAS PROJECT OR CONTRIBUTORS BE +LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL +DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR +SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER +CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, +OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE +USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE. +*****************************************************************************/ + +/************************************************************************************** +* 2016/03/27 Werner Saar (wernsaar@googlemail.com) +* BLASTEST : OK +* CTEST : OK +* TEST : OK +* LAPACK-TEST : OK +**************************************************************************************/ + +#include "common.h" + + +#if defined(POWER8) +#include "cswap_microk_power8.c" +#endif + + +#ifndef HAVE_KERNEL_32 + +static void cswap_kernel_32(BLASLONG n, FLOAT *x, FLOAT *y) +{ + + BLASLONG i=0; + FLOAT f0, f1, f2, f3, f4, f5, f6, f7; + FLOAT g0, g1, g2, g3, g4, g5, g6, g7; + FLOAT *x1=x; + FLOAT *y1=y; + + while ( i 0 ) + { + cswap_kernel_32(n1, x, y); + i=n1; + ix = 2* n1; + iy = 2* n1; + } + + while(i < n) + { + + temp[0] = x[ix] ; + temp[1] = x[ix+1] ; + x[ix] = y[iy] ; + x[ix+1] = y[iy+1] ; + y[iy] = temp[0] ; + y[iy+1] = temp[1] ; + + ix += 2 ; + iy += 2 ; + i++ ; + + + } + + + } + else + { + + inc_x2 = 2 * inc_x; + inc_y2 = 2 * inc_y; + + while(i < n) + { + + temp[0] = x[ix] ; + temp[1] = x[ix+1] ; + x[ix] = y[iy] ; + x[ix+1] = y[iy+1] ; + y[iy] = temp[0] ; + y[iy+1] = temp[1] ; + + ix += inc_x2 ; + iy += inc_y2 ; + i++ ; + + } + + } + return(0); + + +} + + diff --git a/kernel/power/cswap_microk_power8.c b/kernel/power/cswap_microk_power8.c new file mode 100644 index 000000000..90ab59c54 --- /dev/null +++ b/kernel/power/cswap_microk_power8.c @@ -0,0 +1,180 @@ +/*************************************************************************** +Copyright (c) 2013-2016, The OpenBLAS Project +All rights reserved. +Redistribution and use in source and binary forms, with or without +modification, are permitted provided that the following conditions are +met: +1. Redistributions of source code must retain the above copyright +notice, this list of conditions and the following disclaimer. +2. Redistributions in binary form must reproduce the above copyright +notice, this list of conditions and the following disclaimer in +the documentation and/or other materials provided with the +distribution. +3. Neither the name of the OpenBLAS project nor the names of +its contributors may be used to endorse or promote products +derived from this software without specific prior written permission. +THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" +AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE +IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE +ARE DISCLAIMED. IN NO EVENT SHALL THE OPENBLAS PROJECT OR CONTRIBUTORS BE +LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL +DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR +SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER +CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, +OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE +USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE. +*****************************************************************************/ + +/************************************************************************************** +* 2016/03/27 Werner Saar (wernsaar@googlemail.com) +* BLASTEST : OK +* CTEST : OK +* TEST : OK +* LAPACK-TEST : OK +**************************************************************************************/ + +#define HAVE_KERNEL_32 1 + +static void cswap_kernel_32( BLASLONG n, FLOAT *x, FLOAT *y) __attribute__ ((noinline)); + +static void cswap_kernel_32( BLASLONG n, FLOAT *x, FLOAT *y) +{ + + + BLASLONG i = n; + BLASLONG o16 = 16; + BLASLONG o32 = 32; + BLASLONG o48 = 48; + BLASLONG o64 = 64; + BLASLONG o80 = 80; + BLASLONG o96 = 96; + BLASLONG o112 = 112; + FLOAT *x1=x; + FLOAT *y1=y; + FLOAT *x2=x+1; + FLOAT *y2=y+1; + BLASLONG pre = 384; + BLASLONG alpha=0; + + __asm__ __volatile__ + ( + + "addi %3, %3, -4 \n\t" + "addi %4, %4, -4 \n\t" + + ".align 5 \n\t" + "1: \n\t" + + "lxvw4x 32, 0, %2 \n\t" + "lxvw4x 33, %5, %2 \n\t" + "lxvw4x 34, %6, %2 \n\t" + "lxvw4x 35, %7, %2 \n\t" + "lxvw4x 36, %8, %2 \n\t" + "lxvw4x 37, %9, %2 \n\t" + "lxvw4x 38, %10, %2 \n\t" + "lxvw4x 39, %11, %2 \n\t" + + "addi %2, %2, 128 \n\t" + + "lxvw4x 40, 0, %2 \n\t" + "lxvw4x 41, %5, %2 \n\t" + "lxvw4x 42, %6, %2 \n\t" + "lxvw4x 43, %7, %2 \n\t" + "lxvw4x 44, %8, %2 \n\t" + "lxvw4x 45, %9, %2 \n\t" + "lxvw4x 46, %10, %2 \n\t" + "lxvw4x 47, %11, %2 \n\t" + + "addi %2, %2, 128 \n\t" + + "lxvw4x 48, 0, %1 \n\t" + "lxvw4x 49, %5, %1 \n\t" + "lxvw4x 50, %6, %1 \n\t" + "lxvw4x 51, %7, %1 \n\t" + "lxvw4x 52, %8, %1 \n\t" + "lxvw4x 53, %9, %1 \n\t" + "lxvw4x 54, %10, %1 \n\t" + "lxvw4x 55, %11, %1 \n\t" + + "addi %1, %1, 128 \n\t" + + "lxvw4x 56, 0, %1 \n\t" + "lxvw4x 57, %5, %1 \n\t" + "lxvw4x 58, %6, %1 \n\t" + "lxvw4x 59, %7, %1 \n\t" + "lxvw4x 60, %8, %1 \n\t" + "lxvw4x 61, %9, %1 \n\t" + "lxvw4x 62, %10, %1 \n\t" + "lxvw4x 63, %11, %1 \n\t" + + "addi %1, %1, 128 \n\t" + + "stxvw4x 32, 0, %3 \n\t" + "stxvw4x 33, %5, %3 \n\t" + "stxvw4x 34, %6, %3 \n\t" + "stxvw4x 35, %7, %3 \n\t" + "stxvw4x 36, %8, %3 \n\t" + "stxvw4x 37, %9, %3 \n\t" + "stxvw4x 38, %10, %3 \n\t" + "stxvw4x 39, %11, %3 \n\t" + + "addi %3, %3, 128 \n\t" + + "stxvw4x 40, 0, %3 \n\t" + "stxvw4x 41, %5, %3 \n\t" + "stxvw4x 42, %6, %3 \n\t" + "stxvw4x 43, %7, %3 \n\t" + "stxvw4x 44, %8, %3 \n\t" + "stxvw4x 45, %9, %3 \n\t" + "stxvw4x 46, %10, %3 \n\t" + "stxvw4x 47, %11, %3 \n\t" + + "addi %3, %3, 128 \n\t" + + "stxvw4x 48, 0, %4 \n\t" + "stxvw4x 49, %5, %4 \n\t" + "stxvw4x 50, %6, %4 \n\t" + "stxvw4x 51, %7, %4 \n\t" + "stxvw4x 52, %8, %4 \n\t" + "stxvw4x 53, %9, %4 \n\t" + "stxvw4x 54, %10, %4 \n\t" + "stxvw4x 55, %11, %4 \n\t" + + "addi %4, %4, 128 \n\t" + + "stxvw4x 56, 0, %4 \n\t" + "stxvw4x 57, %5, %4 \n\t" + "stxvw4x 58, %6, %4 \n\t" + "stxvw4x 59, %7, %4 \n\t" + "stxvw4x 60, %8, %4 \n\t" + "stxvw4x 61, %9, %4 \n\t" + "stxvw4x 62, %10, %4 \n\t" + "stxvw4x 63, %11, %4 \n\t" + + "addi %4, %4, 128 \n\t" + + "addic. %0 , %0 , -32 \n\t" + "bgt 1b \n\t" + + "2: \n\t" + + : + : + "r" (i), // 0 + "r" (y1), // 1 + "r" (x1), // 2 + "r" (y2), // 3 + "r" (x2), // 4 + "r" (o16), // 5 + "r" (o32), // 6 + "r" (o48), // 7 + "r" (o64), // 8 + "r" (o80), // 9 + "r" (o96), // 10 + "r" (o112) // 11 + : "cr0", "%0", "%2" , "%1", "%3", "%4", "memory" + ); + +} + + diff --git a/kernel/power/sscal.c b/kernel/power/sscal.c new file mode 100644 index 000000000..c6ef5e969 --- /dev/null +++ b/kernel/power/sscal.c @@ -0,0 +1,179 @@ +/*************************************************************************** +Copyright (c) 2013-2016, The OpenBLAS Project +All rights reserved. +Redistribution and use in source and binary forms, with or without +modification, are permitted provided that the following conditions are +met: +1. Redistributions of source code must retain the above copyright +notice, this list of conditions and the following disclaimer. +2. Redistributions in binary form must reproduce the above copyright +notice, this list of conditions and the following disclaimer in +the documentation and/or other materials provided with the +distribution. +3. Neither the name of the OpenBLAS project nor the names of +its contributors may be used to endorse or promote products +derived from this software without specific prior written permission. +THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" +AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE +IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE +ARE DISCLAIMED. IN NO EVENT SHALL THE OPENBLAS PROJECT OR CONTRIBUTORS BE +LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL +DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR +SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER +CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, +OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE +USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE. +*****************************************************************************/ + +/************************************************************************************** +* 2016/03/27 Werner Saar (wernsaar@googlemail.com) +* BLASTEST : OK +* CTEST : OK +* TEST : OK +* LAPACK-TEST : OK +**************************************************************************************/ + +#include "common.h" + +#if defined(POWER8) +#include "sscal_microk_power8.c" +#endif + + +#if !defined(HAVE_KERNEL_16) + +static void sscal_kernel_16( BLASLONG n, FLOAT *da , FLOAT *x ) +{ + + BLASLONG i; + FLOAT alpha = *da; + + for( i=0; i 0 ) + { + alpha[0]=da; + alpha[1]=da; + alpha[2]=da; + alpha[3]=da; + sscal_kernel_16_zero(n1 , alpha , x); + j=n1; + } + + while(j < n) + { + + x[j]=0.0; + j++; + } + + } + else + { + + BLASLONG n1 = n & -32; + if ( n1 > 0 ) + { + alpha[0]=da; + alpha[1]=da; + alpha[2]=da; + alpha[3]=da; + sscal_kernel_16(n1 , alpha , x); + j=n1; + } + while(j < n) + { + + x[j] = da * x[j] ; + j++; + } + } + + + } + else + { + + if ( da == 0.0 ) + { + + while(j < n) + { + + x[i]=0.0; + i += inc_x ; + j++; + } + + } + else + { + + while(j < n) + { + + x[i] = da * x[i] ; + i += inc_x ; + j++; + } + } + + } + return 0; + +} + + diff --git a/kernel/power/sscal_microk_power8.c b/kernel/power/sscal_microk_power8.c new file mode 100644 index 000000000..963cec777 --- /dev/null +++ b/kernel/power/sscal_microk_power8.c @@ -0,0 +1,218 @@ +/*************************************************************************** +Copyright (c) 2013-2016, The OpenBLAS Project +All rights reserved. +Redistribution and use in source and binary forms, with or without +modification, are permitted provided that the following conditions are +met: +1. Redistributions of source code must retain the above copyright +notice, this list of conditions and the following disclaimer. +2. Redistributions in binary form must reproduce the above copyright +notice, this list of conditions and the following disclaimer in +the documentation and/or other materials provided with the +distribution. +3. Neither the name of the OpenBLAS project nor the names of +its contributors may be used to endorse or promote products +derived from this software without specific prior written permission. +THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" +AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE +IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE +ARE DISCLAIMED. IN NO EVENT SHALL THE OPENBLAS PROJECT OR CONTRIBUTORS BE +LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL +DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR +SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER +CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, +OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE +USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE. +*****************************************************************************/ + +/************************************************************************************** +* 2016/03/27 Werner Saar (wernsaar@googlemail.com) +* BLASTEST : OK +* CTEST : OK +* TEST : OK +* LAPACK-TEST : OK +**************************************************************************************/ + +#define HAVE_KERNEL_16 1 + +static void sscal_kernel_16( BLASLONG n, FLOAT *alpha, FLOAT *x) __attribute__ ((noinline)); + +static void sscal_kernel_16( BLASLONG n, FLOAT *alpha, FLOAT *x) +{ + + + BLASLONG i = n; + BLASLONG o16 = 16; + BLASLONG o32 = 32; + BLASLONG o48 = 48; + BLASLONG o64 = 64; + BLASLONG o80 = 80; + BLASLONG o96 = 96; + BLASLONG o112 = 112; + FLOAT *x1=x; + FLOAT *x2=x+1; + BLASLONG pre = 384; + + __asm__ __volatile__ + ( + + "lxvw4x 32, 0, %3 \n\t" + "addi %1, %1, -4 \n\t" + + "dcbt %2, %4 \n\t" + + "lxvw4x 40, 0, %2 \n\t" + "lxvw4x 41, %5, %2 \n\t" + "lxvw4x 42, %6, %2 \n\t" + "lxvw4x 43, %7, %2 \n\t" + "lxvw4x 44, %8, %2 \n\t" + "lxvw4x 45, %9, %2 \n\t" + "lxvw4x 46, %10, %2 \n\t" + "lxvw4x 47, %11, %2 \n\t" + + "addi %2, %2, 128 \n\t" + + "addic. %0 , %0 , -32 \n\t" + "ble 2f \n\t" + + ".align 5 \n\t" + "1: \n\t" + + "dcbt %2, %4 \n\t" + + "xvmulsp 48, 40, 32 \n\t" + "xvmulsp 49, 41, 32 \n\t" + "lxvw4x 40, 0, %2 \n\t" + "lxvw4x 41, %5, %2 \n\t" + "xvmulsp 50, 42, 32 \n\t" + "xvmulsp 51, 43, 32 \n\t" + "lxvw4x 42, %6, %2 \n\t" + "lxvw4x 43, %7, %2 \n\t" + "xvmulsp 52, 44, 32 \n\t" + "xvmulsp 53, 45, 32 \n\t" + "lxvw4x 44, %8, %2 \n\t" + "lxvw4x 45, %9, %2 \n\t" + "xvmulsp 54, 46, 32 \n\t" + "xvmulsp 55, 47, 32 \n\t" + "lxvw4x 46, %10, %2 \n\t" + "lxvw4x 47, %11, %2 \n\t" + + "stxvw4x 48, 0, %1 \n\t" + "stxvw4x 49, %5, %1 \n\t" + "stxvw4x 50, %6, %1 \n\t" + "stxvw4x 51, %7, %1 \n\t" + "stxvw4x 52, %8, %1 \n\t" + "stxvw4x 53, %9, %1 \n\t" + "stxvw4x 54, %10, %1 \n\t" + "stxvw4x 55, %11, %1 \n\t" + + "addi %1, %1, 128 \n\t" + "addi %2, %2, 128 \n\t" + + "addic. %0 , %0 , -32 \n\t" + "bgt 1b \n\t" + + "2: \n\t" + + "xvmulsp 48, 40, 32 \n\t" + "xvmulsp 49, 41, 32 \n\t" + "xvmulsp 50, 42, 32 \n\t" + "xvmulsp 51, 43, 32 \n\t" + "xvmulsp 52, 44, 32 \n\t" + "xvmulsp 53, 45, 32 \n\t" + "xvmulsp 54, 46, 32 \n\t" + "xvmulsp 55, 47, 32 \n\t" + + "stxvw4x 48, 0, %1 \n\t" + "stxvw4x 49, %5, %1 \n\t" + "stxvw4x 50, %6, %1 \n\t" + "stxvw4x 51, %7, %1 \n\t" + "stxvw4x 52, %8, %1 \n\t" + "stxvw4x 53, %9, %1 \n\t" + "stxvw4x 54, %10, %1 \n\t" + "stxvw4x 55, %11, %1 \n\t" + + : + : + "r" (i), // 0 + "r" (x2), // 1 + "r" (x1), // 2 + "r" (alpha), // 3 + "r" (pre), // 4 + "r" (o16), // 5 + "r" (o32), // 6 + "r" (o48), // 7 + "r" (o64), // 8 + "r" (o80), // 9 + "r" (o96), // 10 + "r" (o112) // 11 + : "cr0", "%0", "%2" , "%1", "memory" + ); + +} + + +static void sscal_kernel_16_zero( BLASLONG n, FLOAT *alpha, FLOAT *x) __attribute__ ((noinline)); + +static void sscal_kernel_16_zero( BLASLONG n, FLOAT *alpha, FLOAT *x) +{ + + + BLASLONG i = n; + BLASLONG o16 = 16; + BLASLONG o32 = 32; + BLASLONG o48 = 48; + BLASLONG o64 = 64; + BLASLONG o80 = 80; + BLASLONG o96 = 96; + BLASLONG o112 = 112; + FLOAT *x1=x; + FLOAT *x2=x+1; + BLASLONG pre = 384; + + __asm__ __volatile__ + ( + + "xxlxor 32 , 32 , 32 \n\t" + "addi %1, %1, -4 \n\t" + + + ".align 5 \n\t" + "1: \n\t" + + "stxvw4x 32, 0, %1 \n\t" + "stxvw4x 32, %5, %1 \n\t" + "stxvw4x 32, %6, %1 \n\t" + "stxvw4x 32, %7, %1 \n\t" + "stxvw4x 32, %8, %1 \n\t" + "stxvw4x 32, %9, %1 \n\t" + "stxvw4x 32, %10, %1 \n\t" + "stxvw4x 32, %11, %1 \n\t" + + "addi %1, %1, 128 \n\t" + + "addic. %0 , %0 , -32 \n\t" + "bgt 1b \n\t" + + "2: \n\t" + + : + : + "r" (i), // 0 + "r" (x2), // 1 + "r" (x1), // 2 + "r" (alpha), // 3 + "r" (pre), // 4 + "r" (o16), // 5 + "r" (o32), // 6 + "r" (o48), // 7 + "r" (o64), // 8 + "r" (o80), // 9 + "r" (o96), // 10 + "r" (o112) // 11 + : "cr0", "%0", "%2" , "%1", "memory" + ); + +} + + diff --git a/kernel/power/zscal.c b/kernel/power/zscal.c new file mode 100644 index 000000000..213839a8f --- /dev/null +++ b/kernel/power/zscal.c @@ -0,0 +1,176 @@ +/*************************************************************************** +Copyright (c) 2013-2016, The OpenBLAS Project +All rights reserved. +Redistribution and use in source and binary forms, with or without +modification, are permitted provided that the following conditions are +met: +1. Redistributions of source code must retain the above copyright +notice, this list of conditions and the following disclaimer. +2. Redistributions in binary form must reproduce the above copyright +notice, this list of conditions and the following disclaimer in +the documentation and/or other materials provided with the +distribution. +3. Neither the name of the OpenBLAS project nor the names of +its contributors may be used to endorse or promote products +derived from this software without specific prior written permission. +THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" +AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE +IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE +ARE DISCLAIMED. IN NO EVENT SHALL THE OPENBLAS PROJECT OR CONTRIBUTORS BE +LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL +DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR +SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER +CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, +OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE +USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE. +*****************************************************************************/ + +/************************************************************************************** +* 2016/03/27 Werner Saar (wernsaar@googlemail.com) +* BLASTEST : OK +* CTEST : OK +* TEST : OK +* LAPACK-TEST : OK +**************************************************************************************/ + + +#include "common.h" + +#pragma GCC optimize "O1" + +#if defined(POWER8) +#include "zscal_microk_power8.c" +#endif + + +#ifndef HAVE_KERNEL_8 + +static void zscal_kernel_8(BLASLONG n, FLOAT *x, FLOAT *alpha) +{ + + BLASLONG i=0; + FLOAT *x1=x; + FLOAT alpha_r1=alpha[0]; + FLOAT alpha_r2=alpha[1]; + FLOAT alpha_i1=alpha[2]; + FLOAT alpha_i2=alpha[3]; + FLOAT temp00, temp01, temp10, temp11, temp20, temp21, temp30, temp31; + FLOAT x0_r, x0_i, x1_r, x1_i, x2_r, x2_i, x3_r, x3_i; + + while ( i 0 ) + { + alpha[0] = da_r; + alpha[1] = da_r; + alpha[2] = -da_i; + alpha[3] = da_i; + zscal_kernel_8(n1, x, alpha); + i=n1; + ip = n1 * 2; + + } + + while ( i < n ) + { + + temp = da_r * x[ip] - da_i * x[ip+1] ; + x[ip+1] = da_r * x[ip+1] + da_i * x[ip] ; + x[ip] = temp; + ip += 2; + i++; + } + + } + else + { + + inc_x2 = 2 * inc_x; + + while ( i < n ) + { + + temp = da_r * x[ip] - da_i * x[ip+1] ; + x[ip+1] = da_r * x[ip+1] + da_i * x[ip] ; + x[ip] = temp; + ip += inc_x2; + i++; + } + + + } + + return(0); + +} + + diff --git a/kernel/power/zscal_microk_power8.c b/kernel/power/zscal_microk_power8.c new file mode 100644 index 000000000..5e09d8d79 --- /dev/null +++ b/kernel/power/zscal_microk_power8.c @@ -0,0 +1,224 @@ +/*************************************************************************** +Copyright (c) 2013-2016, The OpenBLAS Project +All rights reserved. +Redistribution and use in source and binary forms, with or without +modification, are permitted provided that the following conditions are +met: +1. Redistributions of source code must retain the above copyright +notice, this list of conditions and the following disclaimer. +2. Redistributions in binary form must reproduce the above copyright +notice, this list of conditions and the following disclaimer in +the documentation and/or other materials provided with the +distribution. +3. Neither the name of the OpenBLAS project nor the names of +its contributors may be used to endorse or promote products +derived from this software without specific prior written permission. +THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" +AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE +IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE +ARE DISCLAIMED. IN NO EVENT SHALL THE OPENBLAS PROJECT OR CONTRIBUTORS BE +LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL +DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR +SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER +CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, +OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE +USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE. +*****************************************************************************/ + +/************************************************************************************** +* 2016/03/25 Werner Saar (wernsaar@googlemail.com) +* +* I don't use fused multipy-add ( lapack precision problems ) +* +* BLASTEST : OK +* CTEST : OK +* TEST : OK +* LAPACK-TEST : OK +**************************************************************************************/ + +#define HAVE_KERNEL_8 1 + +static void zscal_kernel_8( BLASLONG n, FLOAT *x, FLOAT *alpha) __attribute__ ((noinline)); + +static void zscal_kernel_8( BLASLONG n, FLOAT *x, FLOAT *alpha) +{ + + + BLASLONG i = n; + BLASLONG o16 = 16; + BLASLONG o32 = 32; + BLASLONG o48 = 48; + BLASLONG o64 = 64; + BLASLONG o80 = 80; + BLASLONG o96 = 96; + BLASLONG o112 = 112; + FLOAT *x1=x; + FLOAT *x2=x+1; + BLASLONG pre = 384; + + __asm__ __volatile__ + ( + + "lxvd2x 32, 0, %3 \n\t" // alpha_r , alpha_r + "lxvd2x 33, %5, %3 \n\t" // -alpha_i , alpha_i + "addi %1, %1, -8 \n\t" + + "dcbt %2, %4 \n\t" + + "lxvd2x 40, 0, %2 \n\t" // x0_r, x0_i + "lxvd2x 41, %5, %2 \n\t" + "lxvd2x 42, %6, %2 \n\t" + "lxvd2x 43, %7, %2 \n\t" + "lxvd2x 44, %8, %2 \n\t" + "lxvd2x 45, %9, %2 \n\t" + "lxvd2x 46, %10, %2 \n\t" + "lxvd2x 47, %11, %2 \n\t" + + "addi %2, %2, 128 \n\t" + + "addic. %0 , %0 , -8 \n\t" + "ble 2f \n\t" + + ".align 5 \n\t" + "1: \n\t" + + "dcbt %2, %4 \n\t" + + "xvmuldp 48, 40, 32 \n\t" // x0_r * alpha_r, x0_i * alpha_r + "xvmuldp 49, 41, 32 \n\t" + "xvmuldp 50, 42, 32 \n\t" + "xvmuldp 51, 43, 32 \n\t" + "xvmuldp 52, 44, 32 \n\t" + "xvmuldp 53, 45, 32 \n\t" + "xvmuldp 54, 46, 32 \n\t" + "xvmuldp 55, 47, 32 \n\t" + + "xxswapd 56, 40 \n\t" + "xxswapd 57, 41 \n\t" + "xxswapd 58, 42 \n\t" + "xxswapd 59, 43 \n\t" + "xxswapd 60, 44 \n\t" + "xxswapd 61, 45 \n\t" + "xxswapd 62, 46 \n\t" + "xxswapd 63, 47 \n\t" + + "xvmuldp 56, 56, 33 \n\t" // x0_i * -alpha_i, x0_r * alpha_i + "xvmuldp 57, 57, 33 \n\t" + + "lxvd2x 40, 0, %2 \n\t" // x0_r, x0_i + "lxvd2x 41, %5, %2 \n\t" + + "xvmuldp 58, 58, 33 \n\t" + "xvmuldp 59, 59, 33 \n\t" + + "lxvd2x 42, %6, %2 \n\t" + "lxvd2x 43, %7, %2 \n\t" + + "xvmuldp 60, 60, 33 \n\t" + "xvmuldp 61, 61, 33 \n\t" + + "lxvd2x 44, %8, %2 \n\t" + "lxvd2x 45, %9, %2 \n\t" + + "xvmuldp 62, 62, 33 \n\t" + "xvmuldp 63, 63, 33 \n\t" + + "lxvd2x 46, %10, %2 \n\t" + "lxvd2x 47, %11, %2 \n\t" + + "xvadddp 48, 48 , 56 \n\t" + "xvadddp 49, 49 , 57 \n\t" + "xvadddp 50, 50 , 58 \n\t" + "xvadddp 51, 51 , 59 \n\t" + + "stxvd2x 48, 0, %1 \n\t" + "stxvd2x 49, %5, %1 \n\t" + + "xvadddp 52, 52 , 60 \n\t" + "xvadddp 53, 53 , 61 \n\t" + + "stxvd2x 50, %6, %1 \n\t" + "stxvd2x 51, %7, %1 \n\t" + + "xvadddp 54, 54 , 62 \n\t" + "xvadddp 55, 55 , 63 \n\t" + + "stxvd2x 52, %8, %1 \n\t" + "stxvd2x 53, %9, %1 \n\t" + "stxvd2x 54, %10, %1 \n\t" + "stxvd2x 55, %11, %1 \n\t" + + "addi %1, %1, 128 \n\t" + "addi %2, %2, 128 \n\t" + + "addic. %0 , %0 , -8 \n\t" + "bgt 1b \n\t" + + "2: \n\t" + + "xvmuldp 48, 40, 32 \n\t" // x0_r * alpha_r, x0_i * alpha_r + "xvmuldp 49, 41, 32 \n\t" + "xvmuldp 50, 42, 32 \n\t" + "xvmuldp 51, 43, 32 \n\t" + "xvmuldp 52, 44, 32 \n\t" + "xvmuldp 53, 45, 32 \n\t" + "xvmuldp 54, 46, 32 \n\t" + "xvmuldp 55, 47, 32 \n\t" + + "xxswapd 56, 40 \n\t" + "xxswapd 57, 41 \n\t" + "xxswapd 58, 42 \n\t" + "xxswapd 59, 43 \n\t" + "xxswapd 60, 44 \n\t" + "xxswapd 61, 45 \n\t" + "xxswapd 62, 46 \n\t" + "xxswapd 63, 47 \n\t" + + "xvmuldp 56, 56, 33 \n\t" // x0_i * -alpha_i, x0_r * alpha_i + "xvmuldp 57, 57, 33 \n\t" + "xvmuldp 58, 58, 33 \n\t" + "xvmuldp 59, 59, 33 \n\t" + "xvmuldp 60, 60, 33 \n\t" + "xvmuldp 61, 61, 33 \n\t" + "xvmuldp 62, 62, 33 \n\t" + "xvmuldp 63, 63, 33 \n\t" + + "xvadddp 48, 48 , 56 \n\t" + "xvadddp 49, 49 , 57 \n\t" + "xvadddp 50, 50 , 58 \n\t" + "xvadddp 51, 51 , 59 \n\t" + "xvadddp 52, 52 , 60 \n\t" + "xvadddp 53, 53 , 61 \n\t" + "xvadddp 54, 54 , 62 \n\t" + "xvadddp 55, 55 , 63 \n\t" + + "stxvd2x 48, 0, %1 \n\t" + "stxvd2x 49, %5, %1 \n\t" + "stxvd2x 50, %6, %1 \n\t" + "stxvd2x 51, %7, %1 \n\t" + "stxvd2x 52, %8, %1 \n\t" + "stxvd2x 53, %9, %1 \n\t" + "stxvd2x 54, %10, %1 \n\t" + "stxvd2x 55, %11, %1 \n\t" + + + : + : + "r" (i), // 0 + "r" (x2), // 1 + "r" (x1), // 2 + "r" (alpha), // 3 + "r" (pre), // 4 + "r" (o16), // 5 + "r" (o32), // 6 + "r" (o48), // 7 + "r" (o64), // 8 + "r" (o80), // 9 + "r" (o96), // 10 + "r" (o112) // 11 + : "cr0", "%0", "%2" , "%1", "memory" + ); + +} + + diff --git a/kernel/power/zswap.c b/kernel/power/zswap.c new file mode 100644 index 000000000..5ec1eee2e --- /dev/null +++ b/kernel/power/zswap.c @@ -0,0 +1,175 @@ +/*************************************************************************** +Copyright (c) 2013-2016, The OpenBLAS Project +All rights reserved. +Redistribution and use in source and binary forms, with or without +modification, are permitted provided that the following conditions are +met: +1. Redistributions of source code must retain the above copyright +notice, this list of conditions and the following disclaimer. +2. Redistributions in binary form must reproduce the above copyright +notice, this list of conditions and the following disclaimer in +the documentation and/or other materials provided with the +distribution. +3. Neither the name of the OpenBLAS project nor the names of +its contributors may be used to endorse or promote products +derived from this software without specific prior written permission. +THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" +AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE +IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE +ARE DISCLAIMED. IN NO EVENT SHALL THE OPENBLAS PROJECT OR CONTRIBUTORS BE +LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL +DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR +SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER +CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, +OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE +USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE. +*****************************************************************************/ + +/************************************************************************************** +* 2016/03/27 Werner Saar (wernsaar@googlemail.com) +* BLASTEST : OK +* CTEST : OK +* TEST : OK +* LAPACK-TEST : OK +**************************************************************************************/ + +#include "common.h" + + +#if defined(POWER8) +#include "zswap_microk_power8.c" +#endif + + +#ifndef HAVE_KERNEL_16 + +static void zswap_kernel_16(BLASLONG n, FLOAT *x, FLOAT *y) +{ + + BLASLONG i=0; + FLOAT f0, f1, f2, f3, f4, f5, f6, f7; + FLOAT g0, g1, g2, g3, g4, g5, g6, g7; + FLOAT *x1=x; + FLOAT *y1=y; + + while ( i 0 ) + { + zswap_kernel_16(n1, x, y); + i=n1; + ix = 2* n1; + iy = 2* n1; + } + + while(i < n) + { + + temp[0] = x[ix] ; + temp[1] = x[ix+1] ; + x[ix] = y[iy] ; + x[ix+1] = y[iy+1] ; + y[iy] = temp[0] ; + y[iy+1] = temp[1] ; + + ix += 2 ; + iy += 2 ; + i++ ; + + + } + + + } + else + { + + inc_x2 = 2 * inc_x; + inc_y2 = 2 * inc_y; + + while(i < n) + { + + temp[0] = x[ix] ; + temp[1] = x[ix+1] ; + x[ix] = y[iy] ; + x[ix+1] = y[iy+1] ; + y[iy] = temp[0] ; + y[iy+1] = temp[1] ; + + ix += inc_x2 ; + iy += inc_y2 ; + i++ ; + + } + + } + return(0); + + +} + + diff --git a/kernel/power/zswap_microk_power8.c b/kernel/power/zswap_microk_power8.c new file mode 100644 index 000000000..9e5623752 --- /dev/null +++ b/kernel/power/zswap_microk_power8.c @@ -0,0 +1,180 @@ +/*************************************************************************** +Copyright (c) 2013-2016, The OpenBLAS Project +All rights reserved. +Redistribution and use in source and binary forms, with or without +modification, are permitted provided that the following conditions are +met: +1. Redistributions of source code must retain the above copyright +notice, this list of conditions and the following disclaimer. +2. Redistributions in binary form must reproduce the above copyright +notice, this list of conditions and the following disclaimer in +the documentation and/or other materials provided with the +distribution. +3. Neither the name of the OpenBLAS project nor the names of +its contributors may be used to endorse or promote products +derived from this software without specific prior written permission. +THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" +AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE +IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE +ARE DISCLAIMED. IN NO EVENT SHALL THE OPENBLAS PROJECT OR CONTRIBUTORS BE +LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL +DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR +SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER +CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, +OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE +USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE. +*****************************************************************************/ + +/************************************************************************************** +* 2016/03/27 Werner Saar (wernsaar@googlemail.com) +* BLASTEST : OK +* CTEST : OK +* TEST : OK +* LAPACK-TEST : OK +**************************************************************************************/ + +#define HAVE_KERNEL_16 1 + +static void zswap_kernel_16( BLASLONG n, FLOAT *x, FLOAT *y) __attribute__ ((noinline)); + +static void zswap_kernel_16( BLASLONG n, FLOAT *x, FLOAT *y) +{ + + + BLASLONG i = n; + BLASLONG o16 = 16; + BLASLONG o32 = 32; + BLASLONG o48 = 48; + BLASLONG o64 = 64; + BLASLONG o80 = 80; + BLASLONG o96 = 96; + BLASLONG o112 = 112; + FLOAT *x1=x; + FLOAT *y1=y; + FLOAT *x2=x+1; + FLOAT *y2=y+1; + BLASLONG pre = 384; + BLASLONG alpha=0; + + __asm__ __volatile__ + ( + + "addi %3, %3, -8 \n\t" + "addi %4, %4, -8 \n\t" + + ".align 5 \n\t" + "1: \n\t" + + "lxvd2x 32, 0, %2 \n\t" + "lxvd2x 33, %5, %2 \n\t" + "lxvd2x 34, %6, %2 \n\t" + "lxvd2x 35, %7, %2 \n\t" + "lxvd2x 36, %8, %2 \n\t" + "lxvd2x 37, %9, %2 \n\t" + "lxvd2x 38, %10, %2 \n\t" + "lxvd2x 39, %11, %2 \n\t" + + "addi %2, %2, 128 \n\t" + + "lxvd2x 40, 0, %2 \n\t" + "lxvd2x 41, %5, %2 \n\t" + "lxvd2x 42, %6, %2 \n\t" + "lxvd2x 43, %7, %2 \n\t" + "lxvd2x 44, %8, %2 \n\t" + "lxvd2x 45, %9, %2 \n\t" + "lxvd2x 46, %10, %2 \n\t" + "lxvd2x 47, %11, %2 \n\t" + + "addi %2, %2, 128 \n\t" + + "lxvd2x 48, 0, %1 \n\t" + "lxvd2x 49, %5, %1 \n\t" + "lxvd2x 50, %6, %1 \n\t" + "lxvd2x 51, %7, %1 \n\t" + "lxvd2x 52, %8, %1 \n\t" + "lxvd2x 53, %9, %1 \n\t" + "lxvd2x 54, %10, %1 \n\t" + "lxvd2x 55, %11, %1 \n\t" + + "addi %1, %1, 128 \n\t" + + "lxvd2x 56, 0, %1 \n\t" + "lxvd2x 57, %5, %1 \n\t" + "lxvd2x 58, %6, %1 \n\t" + "lxvd2x 59, %7, %1 \n\t" + "lxvd2x 60, %8, %1 \n\t" + "lxvd2x 61, %9, %1 \n\t" + "lxvd2x 62, %10, %1 \n\t" + "lxvd2x 63, %11, %1 \n\t" + + "addi %1, %1, 128 \n\t" + + "stxvd2x 32, 0, %3 \n\t" + "stxvd2x 33, %5, %3 \n\t" + "stxvd2x 34, %6, %3 \n\t" + "stxvd2x 35, %7, %3 \n\t" + "stxvd2x 36, %8, %3 \n\t" + "stxvd2x 37, %9, %3 \n\t" + "stxvd2x 38, %10, %3 \n\t" + "stxvd2x 39, %11, %3 \n\t" + + "addi %3, %3, 128 \n\t" + + "stxvd2x 40, 0, %3 \n\t" + "stxvd2x 41, %5, %3 \n\t" + "stxvd2x 42, %6, %3 \n\t" + "stxvd2x 43, %7, %3 \n\t" + "stxvd2x 44, %8, %3 \n\t" + "stxvd2x 45, %9, %3 \n\t" + "stxvd2x 46, %10, %3 \n\t" + "stxvd2x 47, %11, %3 \n\t" + + "addi %3, %3, 128 \n\t" + + "stxvd2x 48, 0, %4 \n\t" + "stxvd2x 49, %5, %4 \n\t" + "stxvd2x 50, %6, %4 \n\t" + "stxvd2x 51, %7, %4 \n\t" + "stxvd2x 52, %8, %4 \n\t" + "stxvd2x 53, %9, %4 \n\t" + "stxvd2x 54, %10, %4 \n\t" + "stxvd2x 55, %11, %4 \n\t" + + "addi %4, %4, 128 \n\t" + + "stxvd2x 56, 0, %4 \n\t" + "stxvd2x 57, %5, %4 \n\t" + "stxvd2x 58, %6, %4 \n\t" + "stxvd2x 59, %7, %4 \n\t" + "stxvd2x 60, %8, %4 \n\t" + "stxvd2x 61, %9, %4 \n\t" + "stxvd2x 62, %10, %4 \n\t" + "stxvd2x 63, %11, %4 \n\t" + + "addi %4, %4, 128 \n\t" + + "addic. %0 , %0 , -16 \n\t" + "bgt 1b \n\t" + + "2: \n\t" + + : + : + "r" (i), // 0 + "r" (y1), // 1 + "r" (x1), // 2 + "r" (y2), // 3 + "r" (x2), // 4 + "r" (o16), // 5 + "r" (o32), // 6 + "r" (o48), // 7 + "r" (o64), // 8 + "r" (o80), // 9 + "r" (o96), // 10 + "r" (o112) // 11 + : "cr0", "%0", "%2" , "%1", "%3", "%4", "memory" + ); + +} + +