/*********************************************************************/ /* 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 M r3 #define N r4 #define A r5 #define LDA r6 #define B r7 #define AO1 r8 #define AO2 r9 #define AO3 r10 #define AO4 r11 #define J r12 #define INC r30 #define INC2 r31 #define c01 f0 #define c02 f1 #define c03 f2 #define c04 f3 #define c05 f4 #define c06 f5 #define c07 f6 #define c08 f7 #define c09 f8 #define c10 f9 #define c11 f10 #define c12 f11 #define c13 f12 #define c14 f13 #define c15 f14 #define c16 f15 #define sel_p f16 #define sel_s f17 #define c17 f18 #define c18 f19 PROLOGUE PROFCODE li r0, -16 stfpdux f14, SP, r0 stfpdux f15, SP, r0 stfpdux f16, SP, r0 stfpdux f17, SP, r0 stfpdux f18, SP, r0 stfpdux f19, SP, r0 stwu r31, -4(SP) stwu r30, -4(SP) lis r9, 0x3f80 lis r10, 0xbf80 stwu r9, -4(SP) stwu r10, -4(SP) stwu r10, -4(SP) stwu r9, -4(SP) slwi LDA, LDA, BASE_SHIFT cmpwi cr0, M, 0 ble- .L99 cmpwi cr0, N, 0 ble- .L99 andi. r0, A, 2 * SIZE - 1 bne .L100 andi. r0, LDA, 2 * SIZE - 1 bne .L100 li r0, 8 addi SP, SP, -8 lfpsux sel_p, SP, r0 lfpsux sel_s, SP, r0 li INC, 1 * SIZE li INC2, 2 * SIZE subi A, A, 2 * SIZE subi B, B, 2 * SIZE srawi. J, N, 2 ble .L20 .align 4 .L11: mr AO1, A add AO2, A, LDA add AO3, AO2, LDA add AO4, AO3, LDA add A, AO4, LDA srawi. r0, M, 3 mtspr CTR, r0 ble .L15 .align 4 .L12: LFPDUX c01, AO1, INC2 LFXDUX c05, AO2, INC2 LFPDUX c09, AO3, INC2 LFXDUX c13, AO4, INC2 LFPDUX c02, AO1, INC2 LFXDUX c06, AO2, INC2 LFPDUX c10, AO3, INC2 LFXDUX c14, AO4, INC2 LFPDUX c03, AO1, INC2 LFXDUX c07, AO2, INC2 LFPDUX c11, AO3, INC2 LFXDUX c15, AO4, INC2 LFPDUX c04, AO1, INC2 LFXDUX c08, AO2, INC2 LFPDUX c12, AO3, INC2 LFXDUX c16, AO4, INC2 fpsel c17, sel_p, c01, c05 fpsel c18, sel_p, c09, c13 fpsel c01, sel_s, c01, c05 fpsel c05, sel_s, c09, c13 fpsel c09, sel_p, c02, c06 fpsel c13, sel_p, c10, c14 STFPDUX c17, B, INC2 fpsel c02, sel_s, c02, c06 STFPDUX c18, B, INC2 fpsel c06, sel_s, c10, c14 STFXDUX c01, B, INC2 fpsel c10, sel_p, c03, c07 STFXDUX c05, B, INC2 fpsel c14, sel_p, c11, c15 STFPDUX c09, B, INC2 fpsel c03, sel_s, c03, c07 STFPDUX c13, B, INC2 fpsel c07, sel_s, c11, c15 STFXDUX c02, B, INC2 fpsel c11, sel_p, c04, c08 STFXDUX c06, B, INC2 fpsel c15, sel_p, c12, c16 STFPDUX c10, B, INC2 fpsel c04, sel_s, c04, c08 STFPDUX c14, B, INC2 fpsel c08, sel_s, c12, c16 STFXDUX c03, B, INC2 STFXDUX c07, B, INC2 STFPDUX c11, B, INC2 STFPDUX c15, B, INC2 STFXDUX c04, B, INC2 STFXDUX c08, B, INC2 bdnz .L12 .align 4 .L15: andi. r0, M, 7 ble .L19 andi. r0, M, 4 beq .L16 LFPDUX c01, AO1, INC2 LFXDUX c05, AO2, INC2 LFPDUX c09, AO3, INC2 LFXDUX c13, AO4, INC2 LFPDUX c02, AO1, INC2 LFXDUX c06, AO2, INC2 LFPDUX c10, AO3, INC2 LFXDUX c14, AO4, INC2 fpsel c17, sel_p, c01, c05 fpsel c18, sel_p, c09, c13 fpsel c01, sel_s, c01, c05 fpsel c05, sel_s, c09, c13 fpsel c09, sel_p, c02, c06 fpsel c13, sel_p, c10, c14 STFPDUX c17, B, INC2 fpsel c02, sel_s, c02, c06 STFPDUX c18, B, INC2 fpsel c06, sel_s, c10, c14 STFXDUX c01, B, INC2 STFXDUX c05, B, INC2 STFPDUX c09, B, INC2 STFPDUX c13, B, INC2 STFXDUX c02, B, INC2 STFXDUX c06, B, INC2 .align 4 .L16: andi. r0, M, 2 beq .L17 LFPDUX c01, AO1, INC2 LFXDUX c05, AO2, INC2 LFPDUX c09, AO3, INC2 LFXDUX c13, AO4, INC2 fpsel c17, sel_p, c01, c05 fpsel c18, sel_p, c09, c13 fpsel c01, sel_s, c01, c05 fpsel c05, sel_s, c09, c13 STFPDUX c17, B, INC2 STFPDUX c18, B, INC2 STFXDUX c01, B, INC2 STFXDUX c05, B, INC2 .align 4 .L17: andi. r0, M, 1 beq .L19 LFDUX c01, AO1, INC2 LFDUX c02, AO2, INC2 LFDUX c03, AO3, INC2 LFDUX c04, AO4, INC2 fsmfp c01, c02 fsmfp c03, c04 STFPDUX c01, B, INC2 STFPDUX c03, B, INC2 .align 4 .L19: addic. J, J, -1 bgt .L11 .align 4 .L20: andi. J, N, 2 ble .L30 mr AO1, A add AO2, A, LDA add A, AO2, LDA srawi. r0, M, 3 mtspr CTR, r0 ble .L25 .align 4 .L22: LFPDUX c01, AO1, INC2 LFXDUX c05, AO2, INC2 LFPDUX c02, AO1, INC2 LFXDUX c06, AO2, INC2 LFPDUX c03, AO1, INC2 LFXDUX c07, AO2, INC2 LFPDUX c04, AO1, INC2 LFXDUX c08, AO2, INC2 fpsel c17, sel_p, c01, c05 fpsel c01, sel_s, c01, c05 fpsel c09, sel_p, c02, c06 fpsel c02, sel_s, c02, c06 fpsel c10, sel_p, c03, c07 fpsel c03, sel_s, c03, c07 STFPDUX c17, B, INC2 fpsel c11, sel_p, c04, c08 STFXDUX c01, B, INC2 fpsel c04, sel_s, c04, c08 STFPDUX c09, B, INC2 STFXDUX c02, B, INC2 STFPDUX c10, B, INC2 STFXDUX c03, B, INC2 STFPDUX c11, B, INC2 STFXDUX c04, B, INC2 bdnz .L22 .align 4 .L25: andi. r0, M, 7 ble .L30 andi. r0, M, 4 beq .L26 LFPDUX c01, AO1, INC2 LFXDUX c05, AO2, INC2 LFPDUX c02, AO1, INC2 LFXDUX c06, AO2, INC2 fpsel c17, sel_p, c01, c05 fpsel c01, sel_s, c01, c05 fpsel c09, sel_p, c02, c06 fpsel c02, sel_s, c02, c06 STFPDUX c17, B, INC2 STFXDUX c01, B, INC2 STFPDUX c09, B, INC2 STFXDUX c02, B, INC2 .align 4 .L26: andi. r0, M, 2 beq .L27 LFPDUX c01, AO1, INC2 LFXDUX c05, AO2, INC2 fpsel c17, sel_p, c01, c05 fpsel c01, sel_s, c01, c05 STFPDUX c17, B, INC2 STFXDUX c01, B, INC2 .align 4 .L27: andi. r0, M, 1 beq .L30 LFDUX c01, AO1, INC2 LFDUX c02, AO2, INC2 fsmfp c01, c02 STFPDUX c01, B, INC2 .align 4 .L30: andi. J, N, 1 ble .L99 mr AO1, A srawi. r0, M, 3 mtspr CTR, r0 ble .L35 .align 4 .L32: LFPDUX c01, AO1, INC2 LFPDUX c02, AO1, INC2 LFPDUX c03, AO1, INC2 LFPDUX c04, AO1, INC2 STFPDUX c01, B, INC2 STFPDUX c02, B, INC2 STFPDUX c03, B, INC2 STFPDUX c04, B, INC2 bdnz .L32 .align 4 .L35: andi. r0, M, 7 ble .L99 andi. r0, M, 4 beq .L36 LFPDUX c01, AO1, INC2 LFPDUX c02, AO1, INC2 STFPDUX c01, B, INC2 STFPDUX c02, B, INC2 .align 4 .L36: andi. r0, M, 2 beq .L37 LFPDUX c01, AO1, INC2 STFPDUX c01, B, INC2 .align 4 .L37: andi. r0, M, 1 beq .L99 LFDX c01, AO1, INC2 STFDX c01, B, INC2 .align 4 .L99: addi SP, SP, 4 lwzu r30, 4(SP) lwzu r31, 4(SP) subi SP, SP, 12 li r0, 16 lfpdux f19, SP, r0 lfpdux f18, SP, r0 lfpdux f17, SP, r0 lfpdux f16, SP, r0 lfpdux f15, SP, r0 lfpdux f14, SP, r0 addi SP, SP, 16 blr .align 4 .L100: li INC, 1 * SIZE li INC2, 2 * SIZE subi A, A, 1 * SIZE subi B, B, 2 * SIZE srawi. J, N, 2 ble .L120 .align 4 .L111: mr AO1, A add AO2, A, LDA add AO3, AO2, LDA add AO4, AO3, LDA add A, AO4, LDA srawi. r0, M, 3 mtspr CTR, r0 ble .L115 .align 4 .L112: LFDUX c01, AO1, INC LFDUX c02, AO1, INC LFDUX c03, AO1, INC LFDUX c04, AO1, INC LFDUX c09, AO1, INC LFDUX c10, AO1, INC LFDUX c11, AO1, INC LFDUX c12, AO1, INC LFSDUX c01, AO2, INC LFSDUX c02, AO2, INC LFSDUX c03, AO2, INC LFSDUX c04, AO2, INC LFSDUX c09, AO2, INC LFSDUX c10, AO2, INC LFSDUX c11, AO2, INC LFSDUX c12, AO2, INC LFDUX c05, AO3, INC LFDUX c06, AO3, INC LFDUX c07, AO3, INC LFDUX c08, AO3, INC LFDUX c13, AO3, INC LFDUX c14, AO3, INC LFDUX c15, AO3, INC LFDUX c16, AO3, INC LFSDUX c05, AO4, INC LFSDUX c06, AO4, INC LFSDUX c07, AO4, INC LFSDUX c08, AO4, INC LFSDUX c13, AO4, INC LFSDUX c14, AO4, INC LFSDUX c15, AO4, INC LFSDUX c16, AO4, INC STFPDUX c01, B, INC2 STFPDUX c05, B, INC2 STFPDUX c02, B, INC2 STFPDUX c06, B, INC2 STFPDUX c03, B, INC2 STFPDUX c07, B, INC2 STFPDUX c04, B, INC2 STFPDUX c08, B, INC2 STFPDUX c09, B, INC2 STFPDUX c13, B, INC2 STFPDUX c10, B, INC2 STFPDUX c14, B, INC2 STFPDUX c11, B, INC2 STFPDUX c15, B, INC2 STFPDUX c12, B, INC2 STFPDUX c16, B, INC2 bdnz .L112 .align 4 .L115: andi. r0, M, 7 ble .L119 andi. r0, M, 4 beq .L116 LFDUX c01, AO1, INC LFDUX c02, AO1, INC LFDUX c03, AO1, INC LFDUX c04, AO1, INC LFSDUX c01, AO2, INC LFSDUX c02, AO2, INC LFSDUX c03, AO2, INC LFSDUX c04, AO2, INC LFDUX c05, AO3, INC LFDUX c06, AO3, INC LFDUX c07, AO3, INC LFDUX c08, AO3, INC LFSDUX c05, AO4, INC LFSDUX c06, AO4, INC LFSDUX c07, AO4, INC LFSDUX c08, AO4, INC STFPDUX c01, B, INC2 STFPDUX c05, B, INC2 STFPDUX c02, B, INC2 STFPDUX c06, B, INC2 STFPDUX c03, B, INC2 STFPDUX c07, B, INC2 STFPDUX c04, B, INC2 STFPDUX c08, B, INC2 .align 4 .L116: andi. r0, M, 2 beq .L117 LFDUX c01, AO1, INC LFDUX c02, AO1, INC LFSDUX c01, AO2, INC LFSDUX c02, AO2, INC LFDUX c05, AO3, INC LFDUX c06, AO3, INC LFSDUX c05, AO4, INC LFSDUX c06, AO4, INC STFPDUX c01, B, INC2 STFPDUX c05, B, INC2 STFPDUX c02, B, INC2 STFPDUX c06, B, INC2 .align 4 .L117: andi. r0, M, 1 beq .L119 LFDUX c01, AO1, INC LFDUX c05, AO3, INC nop nop LFSDUX c01, AO2, INC LFSDUX c05, AO4, INC STFPDUX c01, B, INC2 STFPDUX c05, B, INC2 .align 4 .L119: addic. J, J, -1 bgt .L111 .align 4 .L120: andi. J, N, 2 ble .L130 mr AO1, A add AO2, A, LDA add A, AO2, LDA srawi. r0, M, 3 mtspr CTR, r0 ble .L125 .align 4 .L122: LFDUX c01, AO1, INC LFDUX c02, AO1, INC LFDUX c03, AO1, INC LFDUX c04, AO1, INC LFDUX c09, AO1, INC LFDUX c10, AO1, INC LFDUX c11, AO1, INC LFDUX c12, AO1, INC LFSDUX c01, AO2, INC LFSDUX c02, AO2, INC LFSDUX c03, AO2, INC LFSDUX c04, AO2, INC LFSDUX c09, AO2, INC LFSDUX c10, AO2, INC LFSDUX c11, AO2, INC LFSDUX c12, AO2, INC STFPDUX c01, B, INC2 STFPDUX c02, B, INC2 STFPDUX c03, B, INC2 STFPDUX c04, B, INC2 STFPDUX c09, B, INC2 STFPDUX c10, B, INC2 STFPDUX c11, B, INC2 STFPDUX c12, B, INC2 bdnz .L122 .align 4 .L125: andi. r0, M, 7 ble .L130 andi. r0, M, 4 beq .L126 LFDUX c01, AO1, INC LFDUX c02, AO1, INC LFDUX c03, AO1, INC LFDUX c04, AO1, INC LFSDUX c01, AO2, INC LFSDUX c02, AO2, INC LFSDUX c03, AO2, INC LFSDUX c04, AO2, INC STFPDUX c01, B, INC2 STFPDUX c02, B, INC2 STFPDUX c03, B, INC2 STFPDUX c04, B, INC2 .align 4 .L126: andi. r0, M, 2 beq .L127 LFDUX c01, AO1, INC LFDUX c02, AO1, INC LFSDUX c01, AO2, INC LFSDUX c02, AO2, INC STFPDUX c01, B, INC2 STFPDUX c02, B, INC2 .align 4 .L127: andi. r0, M, 1 beq .L130 LFDUX c01, AO1, INC LFDUX c02, AO2, INC fsmfp c01, c02 STFPDUX c01, B, INC2 .align 4 .L130: andi. J, N, 1 ble .L999 mr AO1, A srawi. r0, M, 3 mtspr CTR, r0 ble .L135 .align 4 .L132: LFDUX c01, AO1, INC LFDUX c02, AO1, INC LFDUX c03, AO1, INC LFDUX c04, AO1, INC LFDUX c05, AO1, INC LFDUX c06, AO1, INC LFDUX c07, AO1, INC LFDUX c08, AO1, INC fsmfp c01, c02 fsmfp c03, c04 fsmfp c05, c06 fsmfp c07, c08 STFPDUX c01, B, INC2 STFPDUX c03, B, INC2 STFPDUX c05, B, INC2 STFPDUX c07, B, INC2 bdnz .L132 .align 4 .L135: andi. r0, M, 7 ble .L999 andi. r0, M, 4 beq .L136 LFDUX c01, AO1, INC LFDUX c02, AO1, INC LFDUX c03, AO1, INC LFDUX c04, AO1, INC fsmfp c01, c02 fsmfp c03, c04 STFPDUX c01, B, INC2 STFPDUX c03, B, INC2 .align 4 .L136: andi. r0, M, 2 beq .L137 LFDUX c01, AO1, INC LFDUX c02, AO1, INC fsmfp c01, c02 STFPDUX c01, B, INC2 .align 4 .L137: andi. r0, M, 1 beq .L999 LFDX c01, AO1, INC STFDX c01, B, INC2 .align 4 .L999: addi SP, SP, 12 lwzu r30, 4(SP) lwzu r31, 4(SP) subi SP, SP, 12 li r0, 16 lfpdux f19, SP, r0 lfpdux f18, SP, r0 lfpdux f17, SP, r0 lfpdux f16, SP, r0 lfpdux f15, SP, r0 lfpdux f14, SP, r0 addi SP, SP, 16 blr EPILOGUE