/*********************************************************************/
/* 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