Blame thirdparty/openblas/xianyi-OpenBLAS-e6e87a2/kernel/x86/trsm_kernel_LT_2x2.S

kusano 2b45e8
/*********************************************************************/
kusano 2b45e8
/* Copyright 2009, 2010 The University of Texas at Austin.           */
kusano 2b45e8
/* All rights reserved.                                              */
kusano 2b45e8
/*                                                                   */
kusano 2b45e8
/* Redistribution and use in source and binary forms, with or        */
kusano 2b45e8
/* without modification, are permitted provided that the following   */
kusano 2b45e8
/* conditions are met:                                               */
kusano 2b45e8
/*                                                                   */
kusano 2b45e8
/*   1. Redistributions of source code must retain the above         */
kusano 2b45e8
/*      copyright notice, this list of conditions and the following  */
kusano 2b45e8
/*      disclaimer.                                                  */
kusano 2b45e8
/*                                                                   */
kusano 2b45e8
/*   2. Redistributions in binary form must reproduce the above      */
kusano 2b45e8
/*      copyright notice, this list of conditions and the following  */
kusano 2b45e8
/*      disclaimer in the documentation and/or other materials       */
kusano 2b45e8
/*      provided with the distribution.                              */
kusano 2b45e8
/*                                                                   */
kusano 2b45e8
/*    THIS  SOFTWARE IS PROVIDED  BY THE  UNIVERSITY OF  TEXAS AT    */
kusano 2b45e8
/*    AUSTIN  ``AS IS''  AND ANY  EXPRESS OR  IMPLIED WARRANTIES,    */
kusano 2b45e8
/*    INCLUDING, BUT  NOT LIMITED  TO, THE IMPLIED  WARRANTIES OF    */
kusano 2b45e8
/*    MERCHANTABILITY  AND FITNESS FOR  A PARTICULAR  PURPOSE ARE    */
kusano 2b45e8
/*    DISCLAIMED.  IN  NO EVENT SHALL THE UNIVERSITY  OF TEXAS AT    */
kusano 2b45e8
/*    AUSTIN OR CONTRIBUTORS BE  LIABLE FOR ANY DIRECT, INDIRECT,    */
kusano 2b45e8
/*    INCIDENTAL,  SPECIAL, EXEMPLARY,  OR  CONSEQUENTIAL DAMAGES    */
kusano 2b45e8
/*    (INCLUDING, BUT  NOT LIMITED TO,  PROCUREMENT OF SUBSTITUTE    */
kusano 2b45e8
/*    GOODS  OR  SERVICES; LOSS  OF  USE,  DATA,  OR PROFITS;  OR    */
kusano 2b45e8
/*    BUSINESS INTERRUPTION) HOWEVER CAUSED  AND ON ANY THEORY OF    */
kusano 2b45e8
/*    LIABILITY, WHETHER  IN CONTRACT, STRICT  LIABILITY, OR TORT    */
kusano 2b45e8
/*    (INCLUDING NEGLIGENCE OR OTHERWISE)  ARISING IN ANY WAY OUT    */
kusano 2b45e8
/*    OF  THE  USE OF  THIS  SOFTWARE,  EVEN  IF ADVISED  OF  THE    */
kusano 2b45e8
/*    POSSIBILITY OF SUCH DAMAGE.                                    */
kusano 2b45e8
/*                                                                   */
kusano 2b45e8
/* The views and conclusions contained in the software and           */
kusano 2b45e8
/* documentation are those of the authors and should not be          */
kusano 2b45e8
/* interpreted as representing official policies, either expressed   */
kusano 2b45e8
/* or implied, of The University of Texas at Austin.                 */
kusano 2b45e8
/*********************************************************************/
kusano 2b45e8
kusano 2b45e8
#define ASSEMBLER
kusano 2b45e8
#include "common.h"
kusano 2b45e8
kusano 2b45e8
#define STACK	16
kusano 2b45e8
#define ARGS	16
kusano 2b45e8
	
kusano 2b45e8
#define J	 0 + STACK(%esp)
kusano 2b45e8
#define KK	 4 + STACK(%esp)
kusano 2b45e8
#define KKK	 8 + STACK(%esp)
kusano 2b45e8
#define AORIG	12 + STACK(%esp)
kusano 2b45e8
kusano 2b45e8
#define M	 4 + STACK + ARGS(%esp)
kusano 2b45e8
#define N	 8 + STACK + ARGS(%esp)
kusano 2b45e8
#define K	12 + STACK + ARGS(%esp)
kusano 2b45e8
#define ALPHA	16 + STACK + ARGS(%esp)
kusano 2b45e8
#ifdef DOUBLE
kusano 2b45e8
#define A	24 + STACK + ARGS(%esp)
kusano 2b45e8
#define B	28 + STACK + ARGS(%esp)
kusano 2b45e8
#define C	32 + STACK + ARGS(%esp)
kusano 2b45e8
#define LDC	36 + STACK + ARGS(%esp)
kusano 2b45e8
#define OFFSET	40 + STACK + ARGS(%esp)
kusano 2b45e8
#else
kusano 2b45e8
#define A	20 + STACK + ARGS(%esp)
kusano 2b45e8
#define B	24 + STACK + ARGS(%esp)
kusano 2b45e8
#define C	28 + STACK + ARGS(%esp)
kusano 2b45e8
#define LDC	32 + STACK + ARGS(%esp)
kusano 2b45e8
#define OFFSET	36 + STACK + ARGS(%esp)
kusano 2b45e8
#endif
kusano 2b45e8
kusano 2b45e8
#define PREFETCH_OFFSET 48
kusano 2b45e8
kusano 2b45e8
#if defined(PENTIUM3) || defined(PENTIUMM)
kusano 2b45e8
#define REP rep
kusano 2b45e8
#else
kusano 2b45e8
#define REP rep
kusano 2b45e8
#endif
kusano 2b45e8
	
kusano 2b45e8
#define AA	%edx
kusano 2b45e8
#define BB	%ecx
kusano 2b45e8
kusano 2b45e8
	PROLOGUE
kusano 2b45e8
kusano 2b45e8
	subl	$ARGS, %esp	# Generate Stack Frame
kusano 2b45e8
kusano 2b45e8
	pushl	%ebp
kusano 2b45e8
	pushl	%edi
kusano 2b45e8
	pushl	%esi
kusano 2b45e8
	pushl	%ebx
kusano 2b45e8
kusano 2b45e8
	PROFCODE
kusano 2b45e8
kusano 2b45e8
	movl	LDC, %ebp		# ldc			# MEMORY
kusano 2b45e8
	movl	B,   %ebx
kusano 2b45e8
	leal	(, %ebp, SIZE), %ebp
kusano 2b45e8
kusano 2b45e8
#ifdef LN
kusano 2b45e8
       movl	M, %eax
kusano 2b45e8
       leal	(, %eax, SIZE), %eax
kusano 2b45e8
       addl	%eax, C
kusano 2b45e8
       imull	K, %eax
kusano 2b45e8
       addl	%eax, A
kusano 2b45e8
#endif
kusano 2b45e8
kusano 2b45e8
#ifdef RT
kusano 2b45e8
       movl	N, %eax
kusano 2b45e8
       leal	(, %eax, SIZE), %eax
kusano 2b45e8
       imull	K, %eax
kusano 2b45e8
       addl	%eax, %ebx
kusano 2b45e8
kusano 2b45e8
       movl	N,    %eax
kusano 2b45e8
       imull	%ebp, %eax
kusano 2b45e8
       addl	%eax, C
kusano 2b45e8
#endif
kusano 2b45e8
kusano 2b45e8
#ifdef RN
kusano 2b45e8
       movl	OFFSET, %eax
kusano 2b45e8
       negl	%eax
kusano 2b45e8
       movl	%eax, KK
kusano 2b45e8
#endif	
kusano 2b45e8
kusano 2b45e8
#ifdef RT
kusano 2b45e8
       movl	N, %eax
kusano 2b45e8
       subl	OFFSET, %eax
kusano 2b45e8
       movl	%eax, KK
kusano 2b45e8
#endif
kusano 2b45e8
kusano 2b45e8
	movl	N,   %eax		# j = (n >> 1)		# MEMORY
kusano 2b45e8
	sarl	$1,  %eax		 
kusano 2b45e8
	movl	%eax, J			# j = (n >> 1)		# MEMORY
kusano 2b45e8
	je	.L8
kusano 2b45e8
	ALIGN_4
kusano 2b45e8
kusano 2b45e8
.L34:
kusano 2b45e8
#if defined(LT) || defined(RN)
kusano 2b45e8
	movl	A, AA
kusano 2b45e8
#else
kusano 2b45e8
	movl	A, %eax
kusano 2b45e8
	movl	%eax, AORIG
kusano 2b45e8
#endif
kusano 2b45e8
kusano 2b45e8
#ifdef RT
kusano 2b45e8
	movl	K, %eax
kusano 2b45e8
	sall	$1 + BASE_SHIFT, %eax
kusano 2b45e8
	subl	%eax, %ebx
kusano 2b45e8
#endif
kusano 2b45e8
	lea	(, %ebp, 2), %eax
kusano 2b45e8
kusano 2b45e8
#ifdef RT
kusano 2b45e8
	subl	%eax, C
kusano 2b45e8
#endif
kusano 2b45e8
	movl	C, %edi
kusano 2b45e8
#ifndef RT
kusano 2b45e8
	addl	%eax, C
kusano 2b45e8
#endif
kusano 2b45e8
kusano 2b45e8
#ifdef LN
kusano 2b45e8
	movl	OFFSET, %eax
kusano 2b45e8
	addl	M, %eax
kusano 2b45e8
	movl	%eax, KK
kusano 2b45e8
#endif	
kusano 2b45e8
kusano 2b45e8
#ifdef LT
kusano 2b45e8
	movl	OFFSET, %eax
kusano 2b45e8
	movl	%eax, KK
kusano 2b45e8
#endif
kusano 2b45e8
kusano 2b45e8
	movl	M,  %esi
kusano 2b45e8
	sarl	$1, %esi
kusano 2b45e8
	je	.L12
kusano 2b45e8
	ALIGN_4
kusano 2b45e8
kusano 2b45e8
.MainHead:
kusano 2b45e8
#ifdef LN
kusano 2b45e8
       movl	K, %eax
kusano 2b45e8
       sall	$1 + BASE_SHIFT, %eax
kusano 2b45e8
       subl	%eax, AORIG
kusano 2b45e8
#endif
kusano 2b45e8
kusano 2b45e8
#if defined(LN) || defined(RT)
kusano 2b45e8
	movl	KK, %eax
kusano 2b45e8
	leal	(, %eax, SIZE), %eax
kusano 2b45e8
	movl	AORIG, AA
kusano 2b45e8
	leal	(AA,   %eax, 2), AA
kusano 2b45e8
	leal	(%ebx, %eax, 2), BB
kusano 2b45e8
#else
kusano 2b45e8
	movl	%ebx, BB
kusano 2b45e8
#endif	
kusano 2b45e8
kusano 2b45e8
	fldz
kusano 2b45e8
	fldz
kusano 2b45e8
	fldz
kusano 2b45e8
	fldz
kusano 2b45e8
kusano 2b45e8
	FLD	 4 * SIZE(BB)		# b5
kusano 2b45e8
	FLD	 4 * SIZE(AA)		# a5
kusano 2b45e8
	FLD	 0 * SIZE(BB)		# b1
kusano 2b45e8
	FLD	 0 * SIZE(AA)		# a1
kusano 2b45e8
kusano 2b45e8
#if   defined(HAVE_3DNOW)
kusano 2b45e8
	prefetchw	2 * SIZE(%edi)
kusano 2b45e8
 	prefetchw	2 * SIZE(%edi, %ebp, 1)
kusano 2b45e8
#elif defined(HAVE_SSE)
kusano 2b45e8
	prefetchnta	2 * SIZE(%edi)
kusano 2b45e8
 	prefetchnta	2 * SIZE(%edi, %ebp, 1)
kusano 2b45e8
#endif
kusano 2b45e8
kusano 2b45e8
#if defined(LT) || defined(RN)
kusano 2b45e8
	movl	KK, %eax
kusano 2b45e8
#else
kusano 2b45e8
	movl	K,  %eax
kusano 2b45e8
	subl	KK, %eax
kusano 2b45e8
#endif
kusano 2b45e8
	sarl	$2, %eax
kusano 2b45e8
 	je	.L16
kusano 2b45e8
	ALIGN_4
kusano 2b45e8
kusano 2b45e8
.MainLoop:
kusano 2b45e8
#if defined(HAVE_3DNOW)
kusano 2b45e8
	prefetch	(PREFETCH_OFFSET) * SIZE(BB)
kusano 2b45e8
	nop
kusano 2b45e8
#elif defined(HAVE_SSE)
kusano 2b45e8
	prefetchnta	(PREFETCH_OFFSET) * SIZE(BB)
kusano 2b45e8
#if (L2_SIZE == 524288)
kusano 2b45e8
	prefetcht0	(PREFETCH_OFFSET) * SIZE(AA)
kusano 2b45e8
#endif
kusano 2b45e8
#endif
kusano 2b45e8
kusano 2b45e8
	fmul	%st, %st(1)
kusano 2b45e8
	FMUL	 1 * SIZE(BB)
kusano 2b45e8
	fxch	%st(1)
kusano 2b45e8
	faddp	%st, %st(4)
kusano 2b45e8
	FLD	 0 * SIZE(BB)
kusano 2b45e8
	fxch	%st(1)
kusano 2b45e8
	faddp	%st, %st(5)
kusano 2b45e8
	FLD	 1 * SIZE(AA)
kusano 2b45e8
	fmul	%st, %st(1)
kusano 2b45e8
	FMUL	 1 * SIZE(BB)
kusano 2b45e8
	fxch	%st(1)
kusano 2b45e8
	faddp	%st, %st(6)
kusano 2b45e8
	FLD	 2 * SIZE(BB)
kusano 2b45e8
	fxch	%st(1)
kusano 2b45e8
	faddp	%st, %st(7)
kusano 2b45e8
	FLD	 2 * SIZE(AA)
kusano 2b45e8
kusano 2b45e8
	fmul	%st, %st(1)
kusano 2b45e8
	FMUL	 3 * SIZE(BB)
kusano 2b45e8
	fxch	%st(1)
kusano 2b45e8
	faddp	%st, %st(4)
kusano 2b45e8
	FLD	 2 * SIZE(BB)
kusano 2b45e8
	fxch	%st(1)
kusano 2b45e8
	faddp	%st, %st(5)
kusano 2b45e8
	FLD	 3 * SIZE(AA)
kusano 2b45e8
	fmul	%st, %st(1)
kusano 2b45e8
	FMUL	 3 * SIZE(BB)
kusano 2b45e8
	fxch	%st(1)
kusano 2b45e8
	faddp	%st, %st(6)
kusano 2b45e8
	FLD	 8 * SIZE(BB)
kusano 2b45e8
	fxch	%st(1)
kusano 2b45e8
	faddp	%st, %st(7)
kusano 2b45e8
	FLD	 8 * SIZE(AA)
kusano 2b45e8
	fxch	%st(2)
kusano 2b45e8
kusano 2b45e8
#if !defined(HAVE_3DNOW) && defined(HAVE_SSE)  && defined(DOUBLE)
kusano 2b45e8
	prefetchnta	(PREFETCH_OFFSET + 4) * SIZE(BB)
kusano 2b45e8
#if (L2_SIZE == 524288)
kusano 2b45e8
	prefetcht0	(PREFETCH_OFFSET + 4) * SIZE(AA)
kusano 2b45e8
#endif
kusano 2b45e8
#endif
kusano 2b45e8
kusano 2b45e8
	fmul	%st, %st(3)
kusano 2b45e8
	FMUL	 5 * SIZE(BB)
kusano 2b45e8
	fxch	%st(3)
kusano 2b45e8
	faddp	%st, %st(4)
kusano 2b45e8
	FLD	 4 * SIZE(BB)
kusano 2b45e8
	fxch	%st(3)
kusano 2b45e8
	faddp	%st, %st(5)
kusano 2b45e8
	FLD	 5 * SIZE(AA)
kusano 2b45e8
	fmul	%st, %st(3)
kusano 2b45e8
	FMUL	 5 * SIZE(BB)
kusano 2b45e8
	fxch	%st(3)
kusano 2b45e8
	faddp	%st, %st(6)
kusano 2b45e8
	FLD	 6 * SIZE(BB)
kusano 2b45e8
	fxch	%st(3)
kusano 2b45e8
	faddp	%st, %st(7)
kusano 2b45e8
	FLD	 6 * SIZE(AA)
kusano 2b45e8
kusano 2b45e8
	fmul	%st, %st(3)
kusano 2b45e8
	FMUL	 7 * SIZE(BB)
kusano 2b45e8
	fxch	%st(3)
kusano 2b45e8
	faddp	%st, %st(4)
kusano 2b45e8
	FLD	 6 * SIZE(BB)
kusano 2b45e8
	fxch	%st(3)
kusano 2b45e8
	faddp	%st, %st(5)
kusano 2b45e8
	FLD	 7 * SIZE(AA)
kusano 2b45e8
	fmul	%st, %st(3)
kusano 2b45e8
	FMUL	 7 * SIZE(BB)
kusano 2b45e8
	fxch	%st(3)
kusano 2b45e8
	faddp	%st, %st(6)
kusano 2b45e8
	FLD	12 * SIZE(BB)
kusano 2b45e8
	fxch	%st(3)
kusano 2b45e8
	faddp	%st, %st(7)
kusano 2b45e8
	FLD	12 * SIZE(AA)
kusano 2b45e8
	fxch	%st(2)
kusano 2b45e8
kusano 2b45e8
	subl	$-8 * SIZE, BB
kusano 2b45e8
	subl	$-8 * SIZE, AA
kusano 2b45e8
	decl	%eax			# l --
kusano 2b45e8
	jne	.MainLoop
kusano 2b45e8
	ALIGN_4
kusano 2b45e8
kusano 2b45e8
.L16:
kusano 2b45e8
#if defined(LT) || defined(RN)
kusano 2b45e8
	movl	KK, %eax
kusano 2b45e8
#else
kusano 2b45e8
	movl	K,  %eax
kusano 2b45e8
	subl	KK, %eax
kusano 2b45e8
#endif
kusano 2b45e8
	and	$3,  %eax
kusano 2b45e8
	je	.L21
kusano 2b45e8
	ALIGN_4
kusano 2b45e8
kusano 2b45e8
.SubLoop:
kusano 2b45e8
	fmul	%st, %st(1)
kusano 2b45e8
	FMUL	 1 * SIZE(BB)
kusano 2b45e8
	fxch	%st(1)
kusano 2b45e8
	faddp	%st, %st(4)
kusano 2b45e8
	FLD	 0 * SIZE(BB)
kusano 2b45e8
	fxch	%st(1)
kusano 2b45e8
	faddp	%st, %st(5)
kusano 2b45e8
	FLD	 1 * SIZE(AA)
kusano 2b45e8
	fmul	%st, %st(1)
kusano 2b45e8
	FMUL	 1 * SIZE(BB)
kusano 2b45e8
	fxch	%st(1)
kusano 2b45e8
	faddp	%st, %st(6)
kusano 2b45e8
	FLD	 2 * SIZE(BB)
kusano 2b45e8
	fxch	%st(1)
kusano 2b45e8
	faddp	%st, %st(7)
kusano 2b45e8
	FLD	 2 * SIZE(AA)
kusano 2b45e8
kusano 2b45e8
	addl	$2 * SIZE,BB
kusano 2b45e8
	addl	$2 * SIZE,AA
kusano 2b45e8
	decl	%eax
kusano 2b45e8
	jne	 .SubLoop
kusano 2b45e8
	ALIGN_4
kusano 2b45e8
kusano 2b45e8
.L21:
kusano 2b45e8
	ffreep	%st(0)
kusano 2b45e8
	ffreep	%st(0)
kusano 2b45e8
	ffreep	%st(0)
kusano 2b45e8
	ffreep	%st(0)
kusano 2b45e8
kusano 2b45e8
#if defined(LN) || defined(RT)
kusano 2b45e8
	movl	KK, %eax
kusano 2b45e8
#ifdef LN
kusano 2b45e8
	subl	$2, %eax
kusano 2b45e8
#else
kusano 2b45e8
	subl	$2, %eax
kusano 2b45e8
#endif
kusano 2b45e8
kusano 2b45e8
	leal	(, %eax, SIZE), %eax
kusano 2b45e8
kusano 2b45e8
	movl	AORIG, AA
kusano 2b45e8
	leal	(AA,   %eax, 2), AA
kusano 2b45e8
	leal	(%ebx, %eax, 2), BB
kusano 2b45e8
#endif
kusano 2b45e8
kusano 2b45e8
#if defined(LN) || defined(LT)
kusano 2b45e8
	FLD	0 * SIZE(BB)
kusano 2b45e8
	fsubp	%st, %st(1)
kusano 2b45e8
	FLD	1 * SIZE(BB)
kusano 2b45e8
	fsubp	%st, %st(2)
kusano 2b45e8
	FLD	2 * SIZE(BB)
kusano 2b45e8
	fsubp	%st, %st(3)
kusano 2b45e8
	FLD	3 * SIZE(BB)
kusano 2b45e8
	fsubp	%st, %st(4)
kusano 2b45e8
#else
kusano 2b45e8
	FLD	0 * SIZE(AA)
kusano 2b45e8
	fsubp	%st, %st(1)
kusano 2b45e8
	FLD	1 * SIZE(AA)
kusano 2b45e8
	fsubp	%st, %st(3)
kusano 2b45e8
	FLD	2 * SIZE(AA)
kusano 2b45e8
	fsubp	%st, %st(2)
kusano 2b45e8
	FLD	3 * SIZE(AA)
kusano 2b45e8
	fsubp	%st, %st(4)
kusano 2b45e8
#endif
kusano 2b45e8
kusano 2b45e8
#ifdef LN
kusano 2b45e8
       FLD	3 * SIZE(AA)
kusano 2b45e8
       fmul	%st, %st(3)
kusano 2b45e8
       fmulp	%st, %st(4)
kusano 2b45e8
kusano 2b45e8
       FLD	2 * SIZE(AA)
kusano 2b45e8
       fmul	%st(3), %st
kusano 2b45e8
       FLD	2 * SIZE(AA)
kusano 2b45e8
       fmul	%st(5), %st
kusano 2b45e8
kusano 2b45e8
       fsubrp	%st, %st(3)
kusano 2b45e8
       fsubrp	%st, %st(1)
kusano 2b45e8
kusano 2b45e8
       FLD	0 * SIZE(AA)
kusano 2b45e8
       fmul	%st, %st(1)
kusano 2b45e8
       fmulp	%st, %st(2)
kusano 2b45e8
#endif
kusano 2b45e8
kusano 2b45e8
#ifdef LT
kusano 2b45e8
       FLD	0 * SIZE(AA)
kusano 2b45e8
       fmul	%st, %st(1)
kusano 2b45e8
       fmulp	%st, %st(2)
kusano 2b45e8
kusano 2b45e8
       FLD	1 * SIZE(AA)
kusano 2b45e8
       fmul	%st(1), %st
kusano 2b45e8
       FLD	1 * SIZE(AA)
kusano 2b45e8
       fmul	%st(3), %st
kusano 2b45e8
kusano 2b45e8
       fsubrp	%st, %st(5)
kusano 2b45e8
       fsubrp	%st, %st(3)
kusano 2b45e8
kusano 2b45e8
       FLD	3 * SIZE(AA)
kusano 2b45e8
       fmul	%st, %st(3)
kusano 2b45e8
       fmulp	%st, %st(4)
kusano 2b45e8
#endif
kusano 2b45e8
kusano 2b45e8
#ifdef RN
kusano 2b45e8
       FLD	0 * SIZE(BB)
kusano 2b45e8
       fmul	%st, %st(1)
kusano 2b45e8
       fmulp	%st, %st(3)
kusano 2b45e8
kusano 2b45e8
       FLD	1 * SIZE(BB)
kusano 2b45e8
       fmul	%st(1), %st
kusano 2b45e8
       FLD	1 * SIZE(BB)
kusano 2b45e8
       fmul	%st(4), %st
kusano 2b45e8
kusano 2b45e8
       fsubrp	%st, %st(5)
kusano 2b45e8
       fsubrp	%st, %st(2)
kusano 2b45e8
kusano 2b45e8
       FLD	3 * SIZE(BB)
kusano 2b45e8
       fmul	%st, %st(2)
kusano 2b45e8
       fmulp	%st, %st(4)
kusano 2b45e8
#endif
kusano 2b45e8
kusano 2b45e8
#ifdef RT
kusano 2b45e8
       FLD	3 * SIZE(BB)
kusano 2b45e8
       fmul	%st, %st(2)
kusano 2b45e8
       fmulp	%st, %st(4)
kusano 2b45e8
kusano 2b45e8
       FLD	2 * SIZE(BB)
kusano 2b45e8
       fmul	%st(2), %st
kusano 2b45e8
       FLD	2 * SIZE(BB)
kusano 2b45e8
       fmul	%st(5), %st
kusano 2b45e8
kusano 2b45e8
       fsubrp	%st, %st(4)
kusano 2b45e8
       fsubrp	%st, %st(1)
kusano 2b45e8
kusano 2b45e8
       FLD	0 * SIZE(BB)
kusano 2b45e8
       fmul	%st, %st(1)
kusano 2b45e8
       fmulp	%st, %st(3)
kusano 2b45e8
#endif
kusano 2b45e8
kusano 2b45e8
#ifdef LN
kusano 2b45e8
	subl	$2 * SIZE, %edi
kusano 2b45e8
#endif
kusano 2b45e8
kusano 2b45e8
#if defined(LN) || defined(LT)
kusano 2b45e8
	FSTU	0 * SIZE(BB)
kusano 2b45e8
	fxch	%st(1)
kusano 2b45e8
	FSTU	1 * SIZE(BB)
kusano 2b45e8
	fxch	%st(2)
kusano 2b45e8
	FSTU	2 * SIZE(BB)
kusano 2b45e8
	fxch	%st(3)
kusano 2b45e8
	FSTU	3 * SIZE(BB)
kusano 2b45e8
kusano 2b45e8
	FST	1 * SIZE(%edi,%ebp)
kusano 2b45e8
	FST	0 * SIZE(%edi)
kusano 2b45e8
	FST	0 * SIZE(%edi,%ebp)
kusano 2b45e8
	FST	1 * SIZE(%edi)
kusano 2b45e8
#else
kusano 2b45e8
	FSTU	0 * SIZE(AA)
kusano 2b45e8
	fxch	%st(2)
kusano 2b45e8
	FSTU	1 * SIZE(AA)
kusano 2b45e8
	fxch	%st(1)
kusano 2b45e8
	FSTU	2 * SIZE(AA)
kusano 2b45e8
	fxch	%st(3)
kusano 2b45e8
	FSTU	3 * SIZE(AA)
kusano 2b45e8
kusano 2b45e8
	FST	1 * SIZE(%edi,%ebp)
kusano 2b45e8
	FST	1 * SIZE(%edi)
kusano 2b45e8
	FST	0 * SIZE(%edi)
kusano 2b45e8
	FST	0 * SIZE(%edi,%ebp)
kusano 2b45e8
#endif
kusano 2b45e8
kusano 2b45e8
#ifndef LN
kusano 2b45e8
	addl	$2 * SIZE, %edi
kusano 2b45e8
#endif
kusano 2b45e8
kusano 2b45e8
#if defined(LT) || defined(RN)
kusano 2b45e8
	movl	K,  %eax
kusano 2b45e8
	subl	KK, %eax
kusano 2b45e8
	leal	(,%eax, SIZE), %eax
kusano 2b45e8
	leal	(AA, %eax, 2), AA
kusano 2b45e8
	leal	(BB, %eax, 2), BB
kusano 2b45e8
#endif
kusano 2b45e8
kusano 2b45e8
#ifdef LN
kusano 2b45e8
	subl	$2, KK
kusano 2b45e8
#endif
kusano 2b45e8
kusano 2b45e8
#ifdef LT
kusano 2b45e8
	addl	$2, KK
kusano 2b45e8
#endif
kusano 2b45e8
kusano 2b45e8
#ifdef RT
kusano 2b45e8
       movl	K, %eax
kusano 2b45e8
       sall	$1 + BASE_SHIFT, %eax
kusano 2b45e8
       addl	%eax, AORIG
kusano 2b45e8
#endif
kusano 2b45e8
kusano 2b45e8
	decl	%esi			# i --
kusano 2b45e8
	jne	.MainHead
kusano 2b45e8
	ALIGN_4
kusano 2b45e8
kusano 2b45e8
.L12:
kusano 2b45e8
	movl	 M, %eax		# m			# MEMORY
kusano 2b45e8
	andl	$1, %eax
kusano 2b45e8
	je	.L27
kusano 2b45e8
kusano 2b45e8
#ifdef LN
kusano 2b45e8
       movl	K, %eax
kusano 2b45e8
       sall	$0 + BASE_SHIFT, %eax
kusano 2b45e8
       subl	%eax, AORIG
kusano 2b45e8
#endif
kusano 2b45e8
kusano 2b45e8
#if defined(LN) || defined(RT)
kusano 2b45e8
	movl	KK, %eax
kusano 2b45e8
	leal	(, %eax, SIZE), %eax
kusano 2b45e8
	movl	AORIG, AA
kusano 2b45e8
	leal	(AA,   %eax, 1), AA
kusano 2b45e8
	leal	(%ebx, %eax, 2), BB
kusano 2b45e8
#else
kusano 2b45e8
	movl	%ebx, BB
kusano 2b45e8
#endif	
kusano 2b45e8
kusano 2b45e8
	fldz
kusano 2b45e8
	fldz
kusano 2b45e8
kusano 2b45e8
	FLD	0 * SIZE(AA)		# temp1 = *(aoffset + 0)
kusano 2b45e8
kusano 2b45e8
#if defined(LT) || defined(RN)
kusano 2b45e8
	movl	KK, %eax
kusano 2b45e8
#else
kusano 2b45e8
	movl	K,  %eax
kusano 2b45e8
	subl	KK, %eax
kusano 2b45e8
#endif
kusano 2b45e8
	sarl	$1,%eax			# k >> 1		# MEMORY
kusano 2b45e8
	je	 .L54
kusano 2b45e8
	ALIGN_4
kusano 2b45e8
kusano 2b45e8
.L55:
kusano 2b45e8
	FLD	0 * SIZE(BB)		# temp2 = *(boffset + 0)
kusano 2b45e8
	rep
kusano 2b45e8
	fmul	%st(1), %st
kusano 2b45e8
	faddp	%st, %st(2)
kusano 2b45e8
kusano 2b45e8
	FMUL	1 * SIZE(BB)		# temp2 = *(boffset + 0)
kusano 2b45e8
	faddp	%st, %st(2)
kusano 2b45e8
	FLD	1 * SIZE(AA)		# temp1 = *(aoffset + 0)
kusano 2b45e8
kusano 2b45e8
	FLD	2 * SIZE(BB)		# temp2 = *(boffset + 0)
kusano 2b45e8
	rep
kusano 2b45e8
	fmul	%st(1), %st
kusano 2b45e8
	faddp	%st, %st(2)
kusano 2b45e8
kusano 2b45e8
	FMUL	3 * SIZE(BB)		# temp2 = *(boffset + 0)
kusano 2b45e8
	faddp	%st, %st(2)
kusano 2b45e8
	FLD	2 * SIZE(AA)		# temp1 = *(aoffset + 0)
kusano 2b45e8
kusano 2b45e8
	addl	$2 * SIZE, AA
kusano 2b45e8
	addl	$4 * SIZE, BB
kusano 2b45e8
	decl	%eax
kusano 2b45e8
	jne	.L55
kusano 2b45e8
	ALIGN_4
kusano 2b45e8
kusano 2b45e8
.L54:
kusano 2b45e8
#if defined(LT) || defined(RN)
kusano 2b45e8
	movl	KK, %eax
kusano 2b45e8
#else
kusano 2b45e8
	movl	K,  %eax
kusano 2b45e8
	subl	KK, %eax
kusano 2b45e8
#endif
kusano 2b45e8
	andl	$1,%eax			# k & 1
kusano 2b45e8
	je	.L33
kusano 2b45e8
	ALIGN_4
kusano 2b45e8
kusano 2b45e8
	FLD	0 * SIZE(BB)		# temp2 = *(boffset + 0)
kusano 2b45e8
	rep
kusano 2b45e8
	fmul	%st(1), %st
kusano 2b45e8
	faddp	%st, %st(2)
kusano 2b45e8
kusano 2b45e8
	FMUL	1 * SIZE(BB)		# temp2 = *(boffset + 0)
kusano 2b45e8
	faddp	%st, %st(2)
kusano 2b45e8
	FLD	1 * SIZE(AA)		# temp1 = *(aoffset + 0)
kusano 2b45e8
kusano 2b45e8
	addl	$1 * SIZE, AA
kusano 2b45e8
	addl	$2 * SIZE, BB
kusano 2b45e8
	ALIGN_4
kusano 2b45e8
kusano 2b45e8
.L33:
kusano 2b45e8
	ffreep	%st(0)
kusano 2b45e8
kusano 2b45e8
#if defined(LN) || defined(RT)
kusano 2b45e8
	movl	KK, %eax
kusano 2b45e8
#ifdef LN
kusano 2b45e8
	subl	$1, %eax
kusano 2b45e8
#else
kusano 2b45e8
	subl	$2, %eax
kusano 2b45e8
#endif
kusano 2b45e8
kusano 2b45e8
	leal	(, %eax, SIZE), %eax
kusano 2b45e8
kusano 2b45e8
	movl	AORIG, AA
kusano 2b45e8
	leal	(AA,   %eax, 1), AA
kusano 2b45e8
	leal	(%ebx, %eax, 2), BB
kusano 2b45e8
#endif
kusano 2b45e8
kusano 2b45e8
#if defined(LN) || defined(LT)
kusano 2b45e8
	FLD	0 * SIZE(BB)
kusano 2b45e8
	fsubp	%st, %st(1)
kusano 2b45e8
	FLD	1 * SIZE(BB)
kusano 2b45e8
	fsubp	%st, %st(2)
kusano 2b45e8
#else
kusano 2b45e8
	FLD	0 * SIZE(AA)
kusano 2b45e8
	fsubp	%st, %st(1)
kusano 2b45e8
	FLD	1 * SIZE(AA)
kusano 2b45e8
	fsubp	%st, %st(2)
kusano 2b45e8
#endif
kusano 2b45e8
kusano 2b45e8
#if defined(LN) || defined(LT)
kusano 2b45e8
       FLD	0 * SIZE(AA)
kusano 2b45e8
       fmul	%st, %st(1)
kusano 2b45e8
       fmulp	%st, %st(2)
kusano 2b45e8
#endif
kusano 2b45e8
kusano 2b45e8
#ifdef RN
kusano 2b45e8
       FLD	0 * SIZE(BB)
kusano 2b45e8
       fmulp	%st, %st(1)
kusano 2b45e8
kusano 2b45e8
       FLD	1 * SIZE(BB)
kusano 2b45e8
       fmul	%st(1), %st
kusano 2b45e8
kusano 2b45e8
       fsubrp	%st, %st(2)
kusano 2b45e8
kusano 2b45e8
       FLD	3 * SIZE(BB)
kusano 2b45e8
       fmulp	%st, %st(2)
kusano 2b45e8
#endif
kusano 2b45e8
kusano 2b45e8
#ifdef RT
kusano 2b45e8
       FLD	3 * SIZE(BB)
kusano 2b45e8
       fmulp	%st, %st(2)
kusano 2b45e8
kusano 2b45e8
       FLD	2 * SIZE(BB)
kusano 2b45e8
       fmul	%st(2), %st
kusano 2b45e8
kusano 2b45e8
       fsubrp	%st, %st(1)
kusano 2b45e8
kusano 2b45e8
       FLD	0 * SIZE(BB)
kusano 2b45e8
       fmulp	%st, %st(1)
kusano 2b45e8
#endif
kusano 2b45e8
kusano 2b45e8
#ifdef LN
kusano 2b45e8
	subl	$1 * SIZE, %edi
kusano 2b45e8
#endif
kusano 2b45e8
kusano 2b45e8
#if defined(LN) || defined(LT)
kusano 2b45e8
	FSTU	0 * SIZE(BB)
kusano 2b45e8
	fxch	%st(1)
kusano 2b45e8
	FSTU	1 * SIZE(BB)
kusano 2b45e8
#else
kusano 2b45e8
	FSTU	0 * SIZE(AA)
kusano 2b45e8
	fxch	%st(1)
kusano 2b45e8
	FSTU	1 * SIZE(AA)
kusano 2b45e8
#endif
kusano 2b45e8
kusano 2b45e8
	FST	0 * SIZE(%edi,%ebp)
kusano 2b45e8
	FST	0 * SIZE(%edi)
kusano 2b45e8
kusano 2b45e8
#ifndef LN
kusano 2b45e8
	addl	$1 * SIZE, %edi
kusano 2b45e8
#endif
kusano 2b45e8
kusano 2b45e8
#if defined(LT) || defined(RN)
kusano 2b45e8
	movl	K,  %eax
kusano 2b45e8
	subl	KK, %eax
kusano 2b45e8
	leal	(,%eax, SIZE), %eax
kusano 2b45e8
	leal	(AA, %eax, 1), AA
kusano 2b45e8
	leal	(BB, %eax, 2), BB
kusano 2b45e8
#endif
kusano 2b45e8
kusano 2b45e8
#ifdef LN
kusano 2b45e8
	subl	$1, KK
kusano 2b45e8
#endif
kusano 2b45e8
kusano 2b45e8
#ifdef LT
kusano 2b45e8
	addl	$1, KK
kusano 2b45e8
#endif
kusano 2b45e8
kusano 2b45e8
#ifdef RT
kusano 2b45e8
       movl	K, %eax
kusano 2b45e8
       sall	$0 + BASE_SHIFT, %eax
kusano 2b45e8
       addl	%eax, AORIG
kusano 2b45e8
#endif
kusano 2b45e8
	ALIGN_4
kusano 2b45e8
kusano 2b45e8
.L27:
kusano 2b45e8
#ifdef LN
kusano 2b45e8
       movl	K, %eax
kusano 2b45e8
       leal	(    , %eax, SIZE), %eax
kusano 2b45e8
       leal	(%ebx, %eax, 2), %ebx
kusano 2b45e8
#endif
kusano 2b45e8
#if defined(LT) || defined(RN)
kusano 2b45e8
	movl	BB, %ebx
kusano 2b45e8
#endif
kusano 2b45e8
kusano 2b45e8
#ifdef RN
kusano 2b45e8
	addl	$2, KK
kusano 2b45e8
#endif
kusano 2b45e8
kusano 2b45e8
#ifdef RT
kusano 2b45e8
	subl	$2, KK
kusano 2b45e8
#endif
kusano 2b45e8
kusano 2b45e8
	decl	J			# j--			# MEMORY
kusano 2b45e8
	jne	.L34
kusano 2b45e8
	ALIGN_4
kusano 2b45e8
kusano 2b45e8
.L8:
kusano 2b45e8
	movl	N,  %eax		# n			# MEMORY
kusano 2b45e8
	andl	$1, %eax
kusano 2b45e8
	je	.End
kusano 2b45e8
	
kusano 2b45e8
#if defined(LT) || defined(RN)
kusano 2b45e8
	movl	A, AA
kusano 2b45e8
#else
kusano 2b45e8
	movl	A, %eax
kusano 2b45e8
	movl	%eax, AORIG
kusano 2b45e8
#endif
kusano 2b45e8
kusano 2b45e8
#ifdef RT
kusano 2b45e8
	movl	K, %eax
kusano 2b45e8
	sall	$0 + BASE_SHIFT, %eax
kusano 2b45e8
	subl	%eax, %ebx
kusano 2b45e8
#endif
kusano 2b45e8
kusano 2b45e8
#ifdef RT
kusano 2b45e8
	subl	%ebp, C
kusano 2b45e8
#endif
kusano 2b45e8
	movl	C, %edi			# c			# MEMORY
kusano 2b45e8
#ifndef RT
kusano 2b45e8
	addl	%ebp, C
kusano 2b45e8
#endif
kusano 2b45e8
kusano 2b45e8
#ifdef LN
kusano 2b45e8
	movl	OFFSET, %eax
kusano 2b45e8
	addl	M, %eax
kusano 2b45e8
	movl	%eax, KK
kusano 2b45e8
#endif	
kusano 2b45e8
kusano 2b45e8
#ifdef LT
kusano 2b45e8
	movl	OFFSET, %eax
kusano 2b45e8
	movl	%eax, KK
kusano 2b45e8
#endif
kusano 2b45e8
kusano 2b45e8
	movl	M,  %esi		# m			# MEMORY
kusano 2b45e8
	sarl	$1, %esi		# m >> 1
kusano 2b45e8
	je	.L36
kusano 2b45e8
	ALIGN_4
kusano 2b45e8
	
kusano 2b45e8
.L46:
kusano 2b45e8
#ifdef LN
kusano 2b45e8
       movl	K, %eax
kusano 2b45e8
       sall	$1 + BASE_SHIFT, %eax
kusano 2b45e8
       subl	%eax, AORIG
kusano 2b45e8
#endif
kusano 2b45e8
kusano 2b45e8
#if defined(LN) || defined(RT)
kusano 2b45e8
	movl	KK, %eax
kusano 2b45e8
	leal	(, %eax, SIZE), %eax
kusano 2b45e8
	movl	AORIG, AA
kusano 2b45e8
	leal	(AA,   %eax, 2), AA
kusano 2b45e8
	leal	(%ebx, %eax, 1), BB
kusano 2b45e8
#else
kusano 2b45e8
	movl	%ebx, BB
kusano 2b45e8
#endif	
kusano 2b45e8
kusano 2b45e8
	fldz
kusano 2b45e8
	fldz
kusano 2b45e8
	FLD	0 * SIZE(BB)		# temp1 = *(boffset + 0)
kusano 2b45e8
kusano 2b45e8
#if defined(LT) || defined(RN)
kusano 2b45e8
	movl	KK, %eax
kusano 2b45e8
#else
kusano 2b45e8
	movl	K,  %eax
kusano 2b45e8
	subl	KK, %eax
kusano 2b45e8
#endif
kusano 2b45e8
	sarl	$1, %eax
kusano 2b45e8
	je	.L56
kusano 2b45e8
	ALIGN_4
kusano 2b45e8
kusano 2b45e8
.L57:
kusano 2b45e8
	FLD	0 * SIZE(AA)		# temp2 = *(aoffset + 0)
kusano 2b45e8
	fmul	%st(1), %st
kusano 2b45e8
	faddp	%st, %st(2)
kusano 2b45e8
kusano 2b45e8
	FMUL	1 * SIZE(AA)		# temp2 = *(aoffset + 0)
kusano 2b45e8
	faddp	%st, %st(2)
kusano 2b45e8
	FLD	1 * SIZE(BB)		# temp1 = *(boffset + 0)
kusano 2b45e8
kusano 2b45e8
	FLD	2 * SIZE(AA)		# temp2 = *(aoffset + 0)
kusano 2b45e8
	fmul	%st(1), %st
kusano 2b45e8
	faddp	%st, %st(2)
kusano 2b45e8
kusano 2b45e8
	FMUL	3 * SIZE(AA)		# temp2 = *(aoffset + 0)
kusano 2b45e8
	faddp	%st, %st(2)
kusano 2b45e8
	FLD	2 * SIZE(BB)		# temp1 = *(boffset + 0)
kusano 2b45e8
kusano 2b45e8
	addl	$4 * SIZE,AA
kusano 2b45e8
	addl	$2 * SIZE,BB
kusano 2b45e8
	dec	%eax
kusano 2b45e8
	jne	.L57
kusano 2b45e8
	ALIGN_4
kusano 2b45e8
kusano 2b45e8
.L56:
kusano 2b45e8
#if defined(LT) || defined(RN)
kusano 2b45e8
	movl	KK, %eax
kusano 2b45e8
#else
kusano 2b45e8
	movl	K,  %eax
kusano 2b45e8
	subl	KK, %eax
kusano 2b45e8
#endif
kusano 2b45e8
	andl	$1, %eax
kusano 2b45e8
	je	.L45
kusano 2b45e8
	ALIGN_4
kusano 2b45e8
kusano 2b45e8
	FLD	0 * SIZE(AA)		# temp2 = *(aoffset + 0)
kusano 2b45e8
	fmul	%st(1), %st
kusano 2b45e8
	faddp	%st, %st(2)
kusano 2b45e8
kusano 2b45e8
	FMUL	1 * SIZE(AA)		# temp2 = *(aoffset + 0)
kusano 2b45e8
	faddp	%st, %st(2)
kusano 2b45e8
	FLD	3 * SIZE(BB)		# temp1 = *(boffset + 0)
kusano 2b45e8
kusano 2b45e8
	addl	$2 * SIZE,AA
kusano 2b45e8
	addl	$1 * SIZE,BB
kusano 2b45e8
	ALIGN_4
kusano 2b45e8
kusano 2b45e8
.L45:
kusano 2b45e8
	ffreep	%st(0)
kusano 2b45e8
kusano 2b45e8
#if defined(LN) || defined(RT)
kusano 2b45e8
	movl	KK, %eax
kusano 2b45e8
#ifdef LN
kusano 2b45e8
	subl	$2, %eax
kusano 2b45e8
#else
kusano 2b45e8
	subl	$1, %eax
kusano 2b45e8
#endif
kusano 2b45e8
kusano 2b45e8
	leal	(, %eax, SIZE), %eax
kusano 2b45e8
kusano 2b45e8
	movl	AORIG, AA
kusano 2b45e8
	leal	(AA,   %eax, 2), AA
kusano 2b45e8
	leal	(%ebx, %eax, 1), BB
kusano 2b45e8
#endif
kusano 2b45e8
kusano 2b45e8
#if defined(LN) || defined(LT)
kusano 2b45e8
	FLD	0 * SIZE(BB)
kusano 2b45e8
	fsubp	%st, %st(1)
kusano 2b45e8
	FLD	1 * SIZE(BB)
kusano 2b45e8
	fsubp	%st, %st(2)
kusano 2b45e8
#else
kusano 2b45e8
	FLD	0 * SIZE(AA)
kusano 2b45e8
	fsubp	%st, %st(1)
kusano 2b45e8
	FLD	1 * SIZE(AA)
kusano 2b45e8
	fsubp	%st, %st(2)
kusano 2b45e8
#endif
kusano 2b45e8
kusano 2b45e8
#ifdef LN
kusano 2b45e8
       FLD	3 * SIZE(AA)
kusano 2b45e8
       fmulp	%st, %st(2)
kusano 2b45e8
kusano 2b45e8
       FLD	2 * SIZE(AA)
kusano 2b45e8
       fmul	%st(2), %st
kusano 2b45e8
kusano 2b45e8
       fsubrp	%st, %st(1)
kusano 2b45e8
       FLD	0 * SIZE(AA)
kusano 2b45e8
       fmulp	%st, %st(1)
kusano 2b45e8
#endif
kusano 2b45e8
kusano 2b45e8
#ifdef LT
kusano 2b45e8
       FLD	0 * SIZE(AA)
kusano 2b45e8
       fmulp	%st, %st(1)
kusano 2b45e8
kusano 2b45e8
       FLD	1 * SIZE(AA)
kusano 2b45e8
       fmul	%st(1), %st
kusano 2b45e8
kusano 2b45e8
       fsubrp	%st, %st(2)
kusano 2b45e8
kusano 2b45e8
       FLD	3 * SIZE(AA)
kusano 2b45e8
       fmulp	%st, %st(2)
kusano 2b45e8
#endif
kusano 2b45e8
kusano 2b45e8
#ifdef RN
kusano 2b45e8
       FLD	0 * SIZE(BB)
kusano 2b45e8
       fmul	%st, %st(1)
kusano 2b45e8
       fmulp	%st, %st(2)
kusano 2b45e8
#endif
kusano 2b45e8
kusano 2b45e8
#ifdef RT
kusano 2b45e8
       FLD	0 * SIZE(BB)
kusano 2b45e8
       fmul	%st, %st(1)
kusano 2b45e8
       fmulp	%st, %st(2)
kusano 2b45e8
#endif
kusano 2b45e8
kusano 2b45e8
#ifdef LN
kusano 2b45e8
	subl	$2 * SIZE, %edi
kusano 2b45e8
#endif
kusano 2b45e8
kusano 2b45e8
#if defined(LN) || defined(LT)
kusano 2b45e8
	FSTU	0 * SIZE(BB)
kusano 2b45e8
	fxch	%st(1)
kusano 2b45e8
	FSTU	1 * SIZE(BB)
kusano 2b45e8
#else
kusano 2b45e8
	FSTU	0 * SIZE(AA)
kusano 2b45e8
	fxch	%st(1)
kusano 2b45e8
	FSTU	1 * SIZE(AA)
kusano 2b45e8
#endif
kusano 2b45e8
kusano 2b45e8
	FST	1 * SIZE(%edi)
kusano 2b45e8
	FST	0 * SIZE(%edi)
kusano 2b45e8
kusano 2b45e8
#ifndef LN
kusano 2b45e8
	addl	$2 * SIZE, %edi
kusano 2b45e8
#endif
kusano 2b45e8
kusano 2b45e8
#if defined(LT) || defined(RN)
kusano 2b45e8
	movl	K,  %eax
kusano 2b45e8
	subl	KK, %eax
kusano 2b45e8
	leal	(,%eax, SIZE), %eax
kusano 2b45e8
	leal	(AA, %eax, 2), AA
kusano 2b45e8
	leal	(BB, %eax, 1), BB
kusano 2b45e8
#endif
kusano 2b45e8
kusano 2b45e8
#ifdef LN
kusano 2b45e8
	subl	$2, KK
kusano 2b45e8
#endif
kusano 2b45e8
kusano 2b45e8
#ifdef LT
kusano 2b45e8
	addl	$2, KK
kusano 2b45e8
#endif
kusano 2b45e8
kusano 2b45e8
#ifdef RT
kusano 2b45e8
       movl	K, %eax
kusano 2b45e8
       sall	$1 + BASE_SHIFT, %eax
kusano 2b45e8
       addl	%eax, AORIG
kusano 2b45e8
#endif
kusano 2b45e8
kusano 2b45e8
	decl	%esi			# i --
kusano 2b45e8
	jne	.L46
kusano 2b45e8
	ALIGN_4
kusano 2b45e8
kusano 2b45e8
.L36:
kusano 2b45e8
	movl	M,  %eax		# m			# MEMORY
kusano 2b45e8
	andl	$1, %eax		# m & 1
kusano 2b45e8
	je	.L99
kusano 2b45e8
kusano 2b45e8
#ifdef LN
kusano 2b45e8
       movl	K, %eax
kusano 2b45e8
       sall	$0 + BASE_SHIFT, %eax
kusano 2b45e8
       subl	%eax, AORIG
kusano 2b45e8
#endif
kusano 2b45e8
kusano 2b45e8
#if defined(LN) || defined(RT)
kusano 2b45e8
	movl	KK, %eax
kusano 2b45e8
	leal	(, %eax, SIZE), %eax
kusano 2b45e8
	movl	AORIG, AA
kusano 2b45e8
	leal	(AA,   %eax, 1), AA
kusano 2b45e8
	leal	(%ebx, %eax, 1), BB
kusano 2b45e8
#else
kusano 2b45e8
	movl	%ebx, BB
kusano 2b45e8
#endif	
kusano 2b45e8
kusano 2b45e8
	fldz
kusano 2b45e8
kusano 2b45e8
#if defined(LT) || defined(RN)
kusano 2b45e8
	movl	KK, %eax
kusano 2b45e8
#else
kusano 2b45e8
	movl	K,  %eax
kusano 2b45e8
	subl	KK, %eax
kusano 2b45e8
#endif
kusano 2b45e8
	test	%eax, %eax
kusano 2b45e8
	jle	.L52
kusano 2b45e8
	ALIGN_3
kusano 2b45e8
kusano 2b45e8
.L51:
kusano 2b45e8
	FLD	(AA)
kusano 2b45e8
	FMUL	(BB)
kusano 2b45e8
	addl	$1 * SIZE,AA
kusano 2b45e8
	addl	$1 * SIZE,BB
kusano 2b45e8
	faddp	%st,%st(1)
kusano 2b45e8
	decl	%eax
kusano 2b45e8
	jne	.L51
kusano 2b45e8
	ALIGN_4
kusano 2b45e8
kusano 2b45e8
.L52:
kusano 2b45e8
kusano 2b45e8
#if defined(LN) || defined(RT)
kusano 2b45e8
	movl	KK, %eax
kusano 2b45e8
#ifdef LN
kusano 2b45e8
	subl	$1, %eax
kusano 2b45e8
#else
kusano 2b45e8
	subl	$1, %eax
kusano 2b45e8
#endif
kusano 2b45e8
kusano 2b45e8
	leal	(, %eax, SIZE), %eax
kusano 2b45e8
kusano 2b45e8
	movl	AORIG, AA
kusano 2b45e8
	leal	(AA,   %eax, 1), AA
kusano 2b45e8
	leal	(%ebx, %eax, 1), BB
kusano 2b45e8
#endif
kusano 2b45e8
kusano 2b45e8
#if defined(LN) || defined(LT)
kusano 2b45e8
	FLD	0 * SIZE(BB)
kusano 2b45e8
	fsubp	%st, %st(1)
kusano 2b45e8
#else
kusano 2b45e8
	FLD	0 * SIZE(AA)
kusano 2b45e8
	fsubp	%st, %st(1)
kusano 2b45e8
#endif
kusano 2b45e8
kusano 2b45e8
#if defined(LN) || defined(LT)
kusano 2b45e8
       FMUL	0 * SIZE(AA)
kusano 2b45e8
#else
kusano 2b45e8
       FMUL	0 * SIZE(BB)
kusano 2b45e8
#endif
kusano 2b45e8
kusano 2b45e8
#ifdef LN
kusano 2b45e8
	subl	$1 * SIZE, %edi
kusano 2b45e8
#endif
kusano 2b45e8
kusano 2b45e8
#if defined(LN) || defined(LT)
kusano 2b45e8
	FSTU	0 * SIZE(BB)
kusano 2b45e8
#else
kusano 2b45e8
	FSTU	0 * SIZE(AA)
kusano 2b45e8
#endif
kusano 2b45e8
kusano 2b45e8
	FST	0 * SIZE(%edi)
kusano 2b45e8
kusano 2b45e8
#ifndef LN
kusano 2b45e8
	addl	$1 * SIZE, %edi
kusano 2b45e8
#endif
kusano 2b45e8
kusano 2b45e8
#if defined(LT) || defined(RN)
kusano 2b45e8
	movl	K,  %eax
kusano 2b45e8
	subl	KK, %eax
kusano 2b45e8
	leal	(,%eax, SIZE), %eax
kusano 2b45e8
	leal	(AA, %eax, 1), AA
kusano 2b45e8
	leal	(BB, %eax, 1), BB
kusano 2b45e8
#endif
kusano 2b45e8
kusano 2b45e8
#ifdef LN
kusano 2b45e8
	subl	$1, KK
kusano 2b45e8
#endif
kusano 2b45e8
kusano 2b45e8
#ifdef LT
kusano 2b45e8
	addl	$1, KK
kusano 2b45e8
#endif
kusano 2b45e8
kusano 2b45e8
#ifdef RT
kusano 2b45e8
       movl	K, %eax
kusano 2b45e8
       sall	$0 + BASE_SHIFT, %eax
kusano 2b45e8
       addl	%eax, AORIG
kusano 2b45e8
#endif
kusano 2b45e8
	ALIGN_4
kusano 2b45e8
kusano 2b45e8
.L99:
kusano 2b45e8
#ifdef LN
kusano 2b45e8
       movl	K, %eax
kusano 2b45e8
       leal	(%ebx, %eax, SIZE), %ebx
kusano 2b45e8
#endif
kusano 2b45e8
#if defined(LT) || defined(RN)
kusano 2b45e8
	movl	BB, %ebx
kusano 2b45e8
#endif
kusano 2b45e8
kusano 2b45e8
#ifdef RN
kusano 2b45e8
	addl	$1, KK
kusano 2b45e8
#endif
kusano 2b45e8
kusano 2b45e8
#ifdef RT
kusano 2b45e8
	subl	$1, KK
kusano 2b45e8
#endif
kusano 2b45e8
	ALIGN_4
kusano 2b45e8
kusano 2b45e8
.End:
kusano 2b45e8
	popl	%ebx
kusano 2b45e8
	popl	%esi
kusano 2b45e8
	popl	%edi
kusano 2b45e8
	popl	%ebp
kusano 2b45e8
	addl	$ARGS, %esp
kusano 2b45e8
	ret
kusano 2b45e8
kusano 2b45e8
	EPILOGUE