/*********************************************************************/ /* Copyright 2009, 2010 The University of Texas at Austin. */ /* 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. */ /* */ /* THIS SOFTWARE IS PROVIDED BY THE UNIVERSITY OF TEXAS AT */ /* AUSTIN ``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 UNIVERSITY OF TEXAS AT */ /* AUSTIN 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. */ /* */ /* The views and conclusions contained in the software and */ /* documentation are those of the authors and should not be */ /* interpreted as representing official policies, either expressed */ /* or implied, of The University of Texas at Austin. */ /*********************************************************************/ #define ASSEMBLER #include "common.h" #define N $4 #define X $5 #define INCX $6 #define Y $7 #define INCY $8 #define XX $9 #define YY $10 #define C $f17 #define S $f18 #define I $2 #define TEMP $3 #define a1 $f4 #define a2 $f5 #define a3 $f6 #define a4 $f7 #define b1 $f8 #define b2 $f9 #define b3 $f10 #define b4 $f11 #define t1 $f0 #define t2 $f1 #define t3 $f2 #define t4 $f3 PROLOGUE dsll INCX, INCX, ZBASE_SHIFT li TEMP, 2 * SIZE blez N, .L999 dsll INCY, INCY, ZBASE_SHIFT bne INCX, TEMP, .L20 dsra I, N, 1 bne INCY, TEMP, .L20 NOP blez I, .L15 daddiu I, I, -1 LD a1, 0 * SIZE(X) LD b1, 0 * SIZE(Y) LD a2, 1 * SIZE(X) LD b2, 1 * SIZE(Y) LD a3, 2 * SIZE(X) LD b3, 2 * SIZE(Y) MUL t1, S, b1 LD a4, 3 * SIZE(X) MUL t2, C, b1 LD b4, 3 * SIZE(Y) MUL t3, S, b2 blez I, .L13 MUL t4, C, b2 .align 3 .L12: MADD t1, t1, C, a1 LD b1, 4 * SIZE(Y) NMSUB t2, t2, S, a1 LD a1, 4 * SIZE(X) MADD t3, t3, C, a2 LD b2, 5 * SIZE(Y) NMSUB t4, t4, S, a2 LD a2, 5 * SIZE(X) ST t1, 0 * SIZE(X) MUL t1, S, b3 ST t2, 0 * SIZE(Y) MUL t2, C, b3 ST t3, 1 * SIZE(X) MUL t3, S, b4 ST t4, 1 * SIZE(Y) MUL t4, C, b4 MADD t1, t1, C, a3 LD b3, 6 * SIZE(Y) NMSUB t2, t2, S, a3 LD a3, 6 * SIZE(X) MADD t3, t3, C, a4 LD b4, 7 * SIZE(Y) NMSUB t4, t4, S, a4 LD a4, 7 * SIZE(X) ST t1, 2 * SIZE(X) MUL t1, S, b1 ST t2, 2 * SIZE(Y) MUL t2, C, b1 ST t3, 3 * SIZE(X) MUL t3, S, b2 ST t4, 3 * SIZE(Y) MUL t4, C, b2 daddiu I, I, -1 daddiu X, X, 4 * SIZE bgtz I, .L12 daddiu Y, Y, 4 * SIZE .align 3 .L13: MADD t1, t1, C, a1 NMSUB t2, t2, S, a1 MADD t3, t3, C, a2 NMSUB t4, t4, S, a2 ST t1, 0 * SIZE(X) MUL t1, S, b3 ST t2, 0 * SIZE(Y) MUL t2, C, b3 ST t3, 1 * SIZE(X) MUL t3, S, b4 ST t4, 1 * SIZE(Y) MUL t4, C, b4 MADD t1, t1, C, a3 NMSUB t2, t2, S, a3 MADD t3, t3, C, a4 daddiu X, X, 4 * SIZE NMSUB t4, t4, S, a4 daddiu Y, Y, 4 * SIZE ST t1, -2 * SIZE(X) ST t2, -2 * SIZE(Y) ST t3, -1 * SIZE(X) ST t4, -1 * SIZE(Y) .align 3 .L15: andi I, N, 1 blez I, .L999 NOP LD a1, 0 * SIZE(X) LD a2, 1 * SIZE(X) LD b1, 0 * SIZE(Y) LD b2, 1 * SIZE(Y) MUL t1, S, b1 MUL t2, C, b1 MUL t3, S, b2 MUL t4, C, b2 MADD t1, t1, C, a1 NMSUB t2, t2, S, a1 MADD t3, t3, C, a2 NMSUB t4, t4, S, a2 ST t1, 0 * SIZE(X) ST t2, 0 * SIZE(Y) ST t3, 1 * SIZE(X) j .L999 ST t4, 1 * SIZE(Y) .align 3 .L20: move XX, X move YY, Y blez I, .L25 daddiu I, I, -1 LD a1, 0 * SIZE(X) LD b1, 0 * SIZE(Y) LD a2, 1 * SIZE(X) dadd X, X, INCX LD b2, 1 * SIZE(Y) dadd Y, Y, INCY LD a3, 0 * SIZE(X) LD b3, 0 * SIZE(Y) LD a4, 1 * SIZE(X) dadd X, X, INCX MUL t1, S, b1 LD b4, 1 * SIZE(Y) MUL t2, C, b1 dadd Y, Y, INCY MUL t3, S, b2 blez I, .L23 MUL t4, C, b2 .align 3 .L22: MADD t1, t1, C, a1 LD b1, 0 * SIZE(Y) NMSUB t2, t2, S, a1 LD a1, 0 * SIZE(X) MADD t3, t3, C, a2 LD b2, 1 * SIZE(Y) dadd Y, Y, INCY NMSUB t4, t4, S, a2 LD a2, 1 * SIZE(X) dadd X, X, INCX ST t1, 0 * SIZE(XX) MUL t1, S, b3 ST t2, 0 * SIZE(YY) MUL t2, C, b3 ST t3, 1 * SIZE(XX) dadd XX, XX, INCX MUL t3, S, b4 ST t4, 1 * SIZE(YY) dadd YY, YY, INCY MUL t4, C, b4 MADD t1, t1, C, a3 LD b3, 0 * SIZE(Y) NMSUB t2, t2, S, a3 LD a3, 0 * SIZE(X) MADD t3, t3, C, a4 LD b4, 1 * SIZE(Y) dadd Y, Y, INCY NMSUB t4, t4, S, a4 LD a4, 1 * SIZE(X) dadd X, X, INCX ST t1, 0 * SIZE(XX) MUL t1, S, b1 ST t2, 0 * SIZE(YY) MUL t2, C, b1 ST t3, 1 * SIZE(XX) dadd XX, XX, INCX MUL t3, S, b2 ST t4, 1 * SIZE(YY) MUL t4, C, b2 daddiu I, I, -1 bgtz I, .L22 dadd YY, YY, INCY .align 3 .L23: MADD t1, t1, C, a1 NMSUB t2, t2, S, a1 MADD t3, t3, C, a2 NMSUB t4, t4, S, a2 ST t1, 0 * SIZE(XX) MUL t1, S, b3 ST t2, 0 * SIZE(YY) MUL t2, C, b3 ST t3, 1 * SIZE(XX) dadd XX, XX, INCX MUL t3, S, b4 ST t4, 1 * SIZE(YY) dadd YY, YY, INCY MUL t4, C, b4 MADD t1, t1, C, a3 NMSUB t2, t2, S, a3 MADD t3, t3, C, a4 NMSUB t4, t4, S, a4 ST t1, 0 * SIZE(XX) ST t2, 0 * SIZE(YY) ST t3, 1 * SIZE(XX) dadd XX, XX, INCX ST t4, 1 * SIZE(YY) dadd YY, YY, INCY .align 3 .L25: andi I, N, 1 blez I, .L999 NOP .align 3 .L26: LD a1, 0 * SIZE(X) LD a2, 1 * SIZE(X) LD b1, 0 * SIZE(Y) LD b2, 1 * SIZE(Y) MUL t1, S, b1 MUL t2, C, b1 MUL t3, S, b2 MUL t4, C, b2 MADD t1, t1, C, a1 NMSUB t2, t2, S, a1 MADD t3, t3, C, a2 NMSUB t4, t4, S, a2 ST t1, 0 * SIZE(X) ST t2, 0 * SIZE(Y) ST t3, 1 * SIZE(X) ST t4, 1 * SIZE(Y) .align 3 .L999: j $31 NOP EPILOGUE