OpenBLAS/kernel/loongarch64/dot_lsx.S

/***************************************************************************
Copyright (c) 2023, 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.
*****************************************************************************/

#define ASSEMBLER

#include "common.h"

#define N      $r4
#define X      $r5
#define INCX   $r6
#define Y      $r7
#define INCY   $r8

#define I      $r17
#define TEMP   $r18

/* Don't change following FR unless you know the effects. */
#define s1     $f8
#define s2     $f9
#define a1     $f10
#define b1     $f11

PROLOGUE

#ifdef F_INTERFACE
    LDINT     N,      0(N)
    LDINT     INCX,   0(INCX)
    LDINT     INCY,   0(INCY)
#endif

    /* init $f8 and $f9 to zero */
    SUB       s1,     s1,      s1
    SUB       s2,     s2,      s2
    slli.d    INCX,   INCX,    BASE_SHIFT
    li.d      TEMP,   SIZE
    slli.d    INCY,   INCY,    BASE_SHIFT
    bge       $r0,    N,       .L999
    bne       INCX,   TEMP,    .L20   /* inc_x=1 */
    bne       INCY,   TEMP,    .L20   /* inc_y=1 */

    /* !((inc_x == 1) && (inc_y == 1)) */

    /* init $vr8 and $vr9 to zero */
#ifdef DOUBLE
    vldrepl.d  $vr0,  X,       0
#else
    vldrepl.w  $vr0,  X,       0
#endif
#ifdef DSDOT
    vfcvtl.d.s        $vr0,    $vr0
    vfsub.d   $vr8,   $vr0,    $vr0
    vfsub.d   $vr9,   $vr0,    $vr0
#else
    VFSUB     $vr8,   $vr0,    $vr0
    VFSUB     $vr9,   $vr0,    $vr0
#endif

#ifdef DOUBLE
    srai.d    I,      N,       3
#else
    srai.d    I,      N,       4
#endif
    bge       $r0,    I,       .L12   /* FLOAT: <16 ; DOUBLE: <8 */
    .align  3
.L11:
    /* FLOAT: 16~ ; DOUBLE: 8~ */
    vld       $vr0,   X,       0
    vld       $vr1,   X,       16
    vld       $vr2,   X,       32
    vld       $vr3,   X,       48
    vld       $vr4,   Y,       0
    vld       $vr5,   Y,       16
    vld       $vr6,   Y,       32
    vld       $vr7,   Y,       48
    addi.w    I,      I,       -1
    addi.d    X,      X,       64
    addi.d    Y,      Y,       64
#ifdef DSDOT
    vfcvtl.d.s        $vr10,   $vr0
    vfcvtl.d.s        $vr11,   $vr4
    vfcvth.d.s        $vr12,   $vr0
    vfcvth.d.s        $vr13,   $vr4
    vfmadd.d $vr8,    $vr10,   $vr12,  $vr8
    vfmadd.d $vr9,    $vr11,   $vr13,  $vr9
    vfcvtl.d.s        $vr10,   $vr1
    vfcvtl.d.s        $vr11,   $vr5
    vfcvth.d.s        $vr12,   $vr1
    vfcvth.d.s        $vr13,   $vr5
    vfmadd.d $vr8,    $vr10,   $vr12,  $vr8
    vfmadd.d $vr9,    $vr11,   $vr13,  $vr9
    vfcvtl.d.s        $vr10,   $vr2
    vfcvtl.d.s        $vr11,   $vr6
    vfcvth.d.s        $vr12,   $vr2
    vfcvth.d.s        $vr13,   $vr6
    vfmadd.d $vr8,    $vr10,   $vr12,  $vr8
    vfmadd.d $vr9,    $vr11,   $vr13,  $vr9
    vfcvtl.d.s        $vr10,   $vr3
    vfcvtl.d.s        $vr11,   $vr7
    vfcvth.d.s        $vr12,   $vr3
    vfcvth.d.s        $vr13,   $vr7
    vfmadd.d $vr8,    $vr10,   $vr12,  $vr8
    vfmadd.d $vr9,    $vr11,   $vr13,  $vr9
#else
    VFMADD    $vr8,   $vr0,    $vr4,   $vr8
    VFMADD    $vr9,   $vr1,    $vr5,   $vr9
    VFMADD    $vr8,   $vr2,    $vr6,   $vr8
    VFMADD    $vr9,   $vr3,    $vr7,   $vr9
#endif
    bnez      I,      .L11
    .align  3
.L12:
#ifdef DOUBLE
    andi      I,      N,       0x7
    srai.d    I,      I,       1
#else
    andi      I,      N,       0xf
    srai.d    I,      I,       2
#endif
    bge       $r0,    I,       .L14   /* DOUBLE: <2 ; FLOAT: <4 */
    .align  3
.L13:
    /* FLOAT: 4~15 ; DOUBLE: 2~7 */
    vld       $vr0,   X,       0
    vld       $vr4,   Y,       0
    addi.w    I,      I,       -1
    addi.d    X,      X,       16
    addi.d    Y,      Y,       16
#ifdef DSDOT
    vfcvtl.d.s        $vr10,   $vr0
    vfcvtl.d.s        $vr11,   $vr4
    vfcvth.d.s        $vr12,   $vr0
    vfcvth.d.s        $vr13,   $vr4
    vfmadd.d  $vr8,   $vr10,   $vr12,  $vr8
    vfmadd.d  $vr9,   $vr11,   $vr13,  $vr9
#else
    VFMADD    $vr8,   $vr0,    $vr4,   $vr8
#endif
    bnez      I,      .L13
    .align  3
.L14:
    /* store dot in s1 $f8 */
#ifdef DSDOT
    vfadd.d   $vr8,   $vr8,    $vr9
    fsub.s    s2,     s2,      s2,  /* set s2 to 0.0 */
    vpackod.d $vr0,   $vr8,    $vr8
    vfadd.d   $vr8,   $vr8,    $vr0
#else
    VFADD     $vr8,   $vr8,    $vr9
    SUB       s2,     s2,      s2   /* set s2 to 0.0 */
    vpackod.d $vr0,   $vr8,    $vr8
#ifdef DOUBLE
    VFADD     $vr8,   $vr8,    $vr0
#else
    VFADD     $vr8,   $vr8,    $vr0
    vpackod.w $vr0,   $vr8,    $vr8
    VFADD     $vr8,   $vr8,    $vr0
#endif /* defined DOUBLE */
#endif /* defined DSDOT */
    .align  3
.L15:
#ifdef DOUBLE
    andi      I,      N,       0x1
#else
    andi      I,      N,       0x3
#endif
    bge       $r0,    I,       .L999  /* =0 */
   .align  3
.L16:
    /* DOUBLE: 1 ; FLOAT: 1~3 */
    LD        a1,     X,       0
    LD        b1,     Y,       0
#ifdef DSDOT
    fcvt.d.s  a1,     a1
    fcvt.d.s  b1,     b1
    fmadd.d   s1,     b1,      a1,     s1
#else
    MADD      s1,     b1,      a1,     s1
#endif
    addi.d    I,      I,       -1
    addi.d    X,      X,       SIZE
    addi.d    Y,      Y,       SIZE
    bnez      I,      .L16
    b         .L999
    .align 3

.L20:
/* !((inc_x == 1) && (inc_y == 1)) */
    srai.d    I,      N,       3
#ifdef F_INTERFACE
    bgez    INCX, .L21
    addi.d  TEMP, N, -1
    mult    TEMP, INCX
    mflo    TEMP
    dsub    X, X, TEMP
    .align 3

.L21:
    bgez    INCY, .L22
    addi.d  TEMP, N, -1
    mult    TEMP, INCY
    mflo    TEMP
    dsub    Y, Y, TEMP
    .align 3

.L22:
#endif
    bge $r0,    I, .L25  /* <8 */
    .align 3

.L23:
    LD a1,  X,   0 * SIZE
    add.d   X, X, INCX
    LD b1,  Y,   0 * SIZE
    add.d   Y, Y, INCY
#ifdef DSDOT
    fcvt.d.s  a1, a1
    fcvt.d.s  b1, b1
    fmadd.d  s1, b1, a1, s1
#else
    MADD  s1, b1, a1, s1
#endif

    LD a1,  X,   0 * SIZE
    add.d   X, X, INCX
    LD b1,  Y,   0 * SIZE
    add.d   Y, Y, INCY
#ifdef DSDOT
    fcvt.d.s  a1, a1
    fcvt.d.s  b1, b1
    fmadd.d  s2, b1, a1, s2
#else
    MADD  s2, b1, a1, s2
#endif

    LD a1,  X,   0 * SIZE
    add.d   X, X, INCX
    LD b1,  Y,   0 * SIZE
    add.d   Y, Y, INCY
#ifdef DSDOT
    fcvt.d.s  a1, a1
    fcvt.d.s  b1, b1
    fmadd.d  s1, b1, a1, s1
#else
    MADD  s1, b1, a1, s1
#endif

    LD a1,  X,   0 * SIZE
    add.d   X, X, INCX
    LD b1,  Y,   0 * SIZE
    add.d   Y, Y, INCY
#ifdef DSDOT
    fcvt.d.s  a1, a1
    fcvt.d.s  b1, b1
    fmadd.d  s2, b1, a1, s2
#else
    MADD  s2, b1, a1, s2
#endif

    LD a1,  X,   0 * SIZE
    add.d   X, X, INCX
    LD b1,  Y,   0 * SIZE
    add.d   Y, Y, INCY
#ifdef DSDOT
    fcvt.d.s  a1, a1
    fcvt.d.s  b1, b1
    fmadd.d  s1, b1, a1, s1
#else
    MADD  s1, b1, a1, s1
#endif

    LD a1,  X,   0 * SIZE
    add.d   X, X, INCX
    LD b1,  Y,   0 * SIZE
    add.d   Y, Y, INCY
#ifdef DSDOT
    fcvt.d.s  a1, a1
    fcvt.d.s  b1, b1
    fmadd.d  s2, b1, a1, s2
#else
    MADD  s2, b1, a1, s2
#endif

    LD a1,  X,   0 * SIZE
    add.d   X, X, INCX
    LD b1,  Y,   0 * SIZE
    add.d   Y, Y, INCY
#ifdef DSDOT
    fcvt.d.s  a1, a1
    fcvt.d.s  b1, b1
    fmadd.d  s1, b1, a1, s1
#else
    MADD  s1, b1, a1, s1
#endif

    LD a1,  X,   0 * SIZE
    add.d   X, X, INCX
    LD b1,  Y,   0 * SIZE
    add.d   Y, Y, INCY
    addi.d  I, I, -1
#ifdef DSDOT
    fcvt.d.s  a1, a1
    fcvt.d.s  b1, b1
    fmadd.d  s2, b1, a1, s2
#else
    MADD  s2, b1, a1, s2
#endif
    blt $r0,    I, .L23
    .align 3

.L25:
    andi    I,  N, 7
    bge $r0,    I, .L999
    .align  3

.L26:
    LD a1,  X,   0 * SIZE
    add.d   X, X, INCX
    LD b1,  Y,   0 * SIZE
    add.d   Y, Y, INCY
    addi.d  I, I, -1
#ifdef DSDOT
    fcvt.d.s  a1, a1
    fcvt.d.s  b1, b1
    fmadd.d  s1, b1, a1, s1
#else
    MADD  s1, b1, a1, s1
#endif
    blt $r0,    I, .L26
    .align 3

.L999:
#ifdef DSDOT
    fadd.d    $f0,    s1,      s2
#else
    ADD       $f0,    s1,      s2
#endif
    move      $r4,    $r17
    jirl      $r0,    $r1,     0x0

EPILOGUE