atl_dmm2x1x24_5pabc.c

来自「基于Blas CLapck的.用过的人知道是干啥的」· C语言 代码 · 共 3,370 行 · 第 1/4 页

	fmulp %st,%st(7)
	fldl	-56(pB0)
	fxch %st(5)
	faddp %st,%st(3)
	fldl	-56(pA0)
	fmul %st(5),%st
	fxch %st(2)
	faddp %st,%st(6)
	fldl	-56(pA1)
	fmulp %st,%st(5)
	fldl	-48(pB0)
	fxch %st(3)
	faddp %st,%st(1)
	fldl	-48(pA0)
	fmul %st(3),%st
	fxch %st(6)
	faddp %st,%st(4)
	fldl	-48(pA1)
	fmulp %st,%st(3)
	fldl	-40(pB0)
	fxch %st(1)
	faddp %st,%st(7)
	fldl	-40(pA0)
	fmul %st(1),%st
	fxch %st(4)
	faddp %st,%st(2)
	fldl	-40(pA1)
	fmulp %st,%st(1)
	fldl	-32(pB0)
	fxch %st(7)
	faddp %st,%st(5)
	fldl	-32(pA0)
	fmul %st(7),%st
	fxch %st(2)
	faddp %st,%st(6)
	fldl	-32(pA1)
	fmulp %st,%st(7)
	fldl	-24(pB0)
	fxch %st(5)
	faddp %st,%st(3)
	fldl	-24(pA0)
	fmul %st(5),%st
	fxch %st(6)
	faddp %st,%st(4)
	fldl	-24(pA1)
	fmulp %st,%st(5)
	fldl	-16(pB0)
	fxch %st(3)
	faddp %st,%st(1)
	fldl	-16(pA0)
	fmul %st(3),%st
	fxch %st(4)
	faddp %st,%st(2)
	fldl	-16(pA1)
	fmulp %st,%st(3)
	fldl	-8(pB0)
	fxch %st(1)
	faddp %st,%st(7)
	fldl	-8(pA0)
	fmul %st(1),%st
	fxch %st(2)
	faddp %st,%st(6)
	fldl	-8(pA1)
	fmulp %st,%st(1)
	fldl	(pB0)
	fxch %st(7)
	faddp %st,%st(5)
	fldl	(pA0)
	fmul %st(7),%st
	fxch %st(6)
	faddp %st,%st(4)
	fldl	(pA1)
	fmulp %st,%st(7)
	fldl	8(pB0)
	fxch %st(5)
	faddp %st,%st(3)
	fldl	8(pA0)
	fmul %st(5),%st
	fxch %st(4)
	faddp %st,%st(2)
	fldl	8(pA1)
	fmulp %st,%st(5)
	fldl	16(pB0)
	fxch %st(3)
	faddp %st,%st(1)
	fldl	16(pA0)
	fmul %st(3),%st
	fxch %st(2)
	faddp %st,%st(6)
	fldl	16(pA1)
	fmulp %st,%st(3)
	fldl	24(pB0)
	fxch %st(1)
	faddp %st,%st(7)
	fldl	24(pA0)
	fmul %st(1),%st
	fxch %st(6)
	faddp %st,%st(4)
	fldl	24(pA1)
	fmulp %st,%st(1)
	fldl	32(pB0)
	fxch %st(7)
	faddp %st,%st(5)
	fldl	32(pA0)
	fmul %st(7),%st
	fxch %st(4)
	faddp %st,%st(2)
	fldl	32(pA1)
	fmulp %st,%st(7)
	fldl	40(pB0)
	fxch %st(5)
	faddp %st,%st(3)
	fldl	40(pA0)
	fmul %st(5),%st
	fxch %st(2)
	faddp %st,%st(6)
	fldl	40(pA1)
	fmulp %st,%st(5)
	fldl	48(pB0)
	fxch %st(3)
	faddp %st,%st(1)
	fldl	48(pA0)
	fmul %st(3),%st
	fxch %st(6)
	faddp %st,%st(4)
	fldl	48(pA1)
	fmulp %st,%st(3)
	fldl	56(pB0)
	fxch %st(1)
	faddp %st,%st(7)
	fldl	56(pA0)
	fmul %st(1),%st
	fxch %st(4)
	faddp %st,%st(2)
	fldl	56(pA1)
	fmulp %st,%st(1)
	fldl	64(pB0)
	fxch %st(7)
	faddp %st,%st(5)
	fldl	64(pA0)
	fmul %st(7),%st
	fxch %st(2)
	faddp %st,%st(6)
	fldl	64(pA1)
	fmulp %st,%st(7)
	fldl	72(pB0)
	fxch %st(5)
	faddp %st,%st(3)
	fldl	72(pA0)
	fmul %st(5),%st
	fxch %st(6)
	faddp %st,%st(4)
	fldl	72(pA1)
	fmulp %st,%st(5)
	fldl	80(pB0)
	fxch %st(3)
	faddp %st,%st(1)
	fldl	80(pA0)
	fmul %st(3),%st
	fxch %st(4)
	faddp %st,%st(2)
	fldl	80(pA1)
	fmulp %st,%st(3)
	fldl	88(pB0)
	fxch %st(1)
	faddp %st,%st(7)
	fldl	88(pA0)
	fmul %st(1),%st
	fxch %st(2)
	faddp %st,%st(6)
	fldl	88(pA1)
	fmulp %st,%st(1)
	fxch %st(6)
	addl	$NB2so, pA0
	faddp %st,%st(4)
	fxch %st(4)
	addl	$NB2so, pA1
	faddp %st,%st(2)
	faddp %st,%st(2)
	faddp %st,%st(2)
	faddp %st,%st(2)
/*
 *      Write results back to C
 */
   #ifdef BETAX
        fldl    (pC)
        fldl    OFF(pC)
        fldl    BETAOFF(%esp)
        fmul    %st, %st(1)
        fmulp   %st, %st(2)
        faddp   %st, %st(3)
        faddp   %st, %st(1)
   #endif
	fstpl (pC)
	fstpl OFF(pC)
/*
 *      pC += 2;  pA += 2*NB; pB -= NB;
 */
	addl	$incCm, pC
	fldl	-96(pB0)
	fldl	-96(pA0)
	fmul %st(1),%st
	fldl	-96(pA1)
	fmulp %st,%st(2)
	fldl	-88(pB0)
	fldl	-88(pA0)
	fmul %st(1),%st
	fldl	-88(pA1)
	fmulp %st,%st(2)
	fldl	-80(pB0)
	fldl	-80(pA0)
	fmul %st(1),%st
   #if defined(BETA0) || defined(BETAX)
        fldz
   #else
	fldl (pC)
   #endif
	faddp %st,%st(5)
	fldl	-80(pA1)
	fmulp %st,%st(2)
	fldl	-72(pB0)
   #if defined(BETA0) || defined(BETAX)
        fldz
   #else
	fldl OFF(pC)
   #endif
	faddp %st,%st(7)
	fldl	-72(pA0)
	fmul %st(1),%st
	fxch %st(6)
	faddp %st,%st(4)
	fldl	-72(pA1)
	fmulp %st,%st(1)
	fldl	-64(pB0)
	fxch %st(7)
	faddp %st,%st(5)
	fldl	-64(pA0)
	fmul %st(7),%st
	fxch %st(4)
	faddp %st,%st(2)
	fldl	-64(pA1)
	fmulp %st,%st(7)
	fldl	-56(pB0)
	fxch %st(5)
	faddp %st,%st(3)
	fldl	-56(pA0)
	fmul %st(5),%st
	fxch %st(2)
	faddp %st,%st(6)
	fldl	-56(pA1)
	fmulp %st,%st(5)
	fldl	-48(pB0)
	fxch %st(3)
	faddp %st,%st(1)
	fldl	-48(pA0)
	fmul %st(3),%st
	fxch %st(6)
	faddp %st,%st(4)
	fldl	-48(pA1)
	fmulp %st,%st(3)
	fldl	-40(pB0)
	fxch %st(1)
	faddp %st,%st(7)
	fldl	-40(pA0)
	fmul %st(1),%st
	fxch %st(4)
	faddp %st,%st(2)
	fldl	-40(pA1)
	fmulp %st,%st(1)
	fldl	-32(pB0)
	fxch %st(7)
	faddp %st,%st(5)
	fldl	-32(pA0)
	fmul %st(7),%st
	fxch %st(2)
	faddp %st,%st(6)
	fldl	-32(pA1)
	fmulp %st,%st(7)
	fldl	-24(pB0)
	fxch %st(5)
	faddp %st,%st(3)
	fldl	-24(pA0)
	fmul %st(5),%st
	fxch %st(6)
	faddp %st,%st(4)
	fldl	-24(pA1)
	fmulp %st,%st(5)
	fldl	-16(pB0)
	fxch %st(3)
	faddp %st,%st(1)
	fldl	-16(pA0)
	fmul %st(3),%st
	fxch %st(4)
	faddp %st,%st(2)
	fldl	-16(pA1)
	fmulp %st,%st(3)
	fldl	-8(pB0)
	fxch %st(1)
	faddp %st,%st(7)
	fldl	-8(pA0)
	fmul %st(1),%st
	fxch %st(2)
	faddp %st,%st(6)
	fldl	-8(pA1)
	fmulp %st,%st(1)
	fldl	(pB0)
	fxch %st(7)
	faddp %st,%st(5)
	fldl	(pA0)
	fmul %st(7),%st
	fxch %st(6)
	faddp %st,%st(4)
	fldl	(pA1)
	fmulp %st,%st(7)
	fldl	8(pB0)
	fxch %st(5)
	faddp %st,%st(3)
	fldl	8(pA0)
	fmul %st(5),%st
	fxch %st(4)
	faddp %st,%st(2)
	fldl	8(pA1)
	fmulp %st,%st(5)
	fldl	16(pB0)
	fxch %st(3)
	faddp %st,%st(1)
	fldl	16(pA0)
	fmul %st(3),%st
	fxch %st(2)
	faddp %st,%st(6)
	fldl	16(pA1)
	fmulp %st,%st(3)
	fldl	24(pB0)
	fxch %st(1)
	faddp %st,%st(7)
	fldl	24(pA0)
	fmul %st(1),%st
	fxch %st(6)
	faddp %st,%st(4)
	fldl	24(pA1)
	fmulp %st,%st(1)
	fldl	32(pB0)
	fxch %st(7)
	faddp %st,%st(5)
	fldl	32(pA0)
	fmul %st(7),%st
	fxch %st(4)
	faddp %st,%st(2)
	fldl	32(pA1)
	fmulp %st,%st(7)
	fldl	40(pB0)
	fxch %st(5)
	faddp %st,%st(3)
	fldl	40(pA0)
	fmul %st(5),%st
	fxch %st(2)
	faddp %st,%st(6)
	fldl	40(pA1)
	fmulp %st,%st(5)
	fldl	48(pB0)
	fxch %st(3)
	faddp %st,%st(1)
	fldl	48(pA0)
	fmul %st(3),%st
	fxch %st(6)
	faddp %st,%st(4)
	fldl	48(pA1)
	fmulp %st,%st(3)
	fldl	56(pB0)
	fxch %st(1)
	faddp %st,%st(7)
	fldl	56(pA0)
	fmul %st(1),%st
	fxch %st(4)
	faddp %st,%st(2)
	fldl	56(pA1)
	fmulp %st,%st(1)
	fldl	64(pB0)
	fxch %st(7)
	faddp %st,%st(5)
	fldl	64(pA0)
	fmul %st(7),%st
	fxch %st(2)
	faddp %st,%st(6)
	fldl	64(pA1)
	fmulp %st,%st(7)
	fldl	72(pB0)
	fxch %st(5)
	faddp %st,%st(3)
	fldl	72(pA0)
	fmul %st(5),%st
	fxch %st(6)
	faddp %st,%st(4)
	fldl	72(pA1)
	fmulp %st,%st(5)
	fldl	80(pB0)
	fxch %st(3)
	faddp %st,%st(1)
	fldl	80(pA0)
	fmul %st(3),%st
	fxch %st(4)
	faddp %st,%st(2)
	fldl	80(pA1)
	fmulp %st,%st(3)
	fldl	88(pB0)
	fxch %st(1)
	faddp %st,%st(7)
	fldl	88(pA0)
	fmul %st(1),%st
	fxch %st(2)
	faddp %st,%st(6)
	fldl	88(pA1)
	fmulp %st,%st(1)
	fxch %st(6)
	addl	$NB2so, pA0
	faddp %st,%st(4)
	fxch %st(4)
	addl	$NB2so, pA1
	faddp %st,%st(2)
	faddp %st,%st(2)
	faddp %st,%st(2)
	faddp %st,%st(2)
/*
 *      Write results back to C
 */
   #ifdef BETAX
        fldl    (pC)
        fldl    OFF(pC)
        fldl    BETAOFF(%esp)
        fmul    %st, %st(1)
        fmulp   %st, %st(2)
        faddp   %st, %st(3)
        faddp   %st, %st(1)
   #endif
	fstpl (pC)
	fstpl OFF(pC)
/*
 *      pC += 2;  pA += 2*NB; pB -= NB;
 */
	addl	$incCm, pC
	fldl	-96(pB0)
	fldl	-96(pA0)
	fmul %st(1),%st
	fldl	-96(pA1)
	fmulp %st,%st(2)
	fldl	-88(pB0)
	fldl	-88(pA0)
	fmul %st(1),%st
	fldl	-88(pA1)
	fmulp %st,%st(2)
	fldl	-80(pB0)
	fldl	-80(pA0)
	fmul %st(1),%st
   #if defined(BETA0) || defined(BETAX)
        fldz
   #else
	fldl (pC)
   #endif
	faddp %st,%st(5)
	fldl	-80(pA1)
	fmulp %st,%st(2)
	fldl	-72(pB0)
   #if defined(BETA0) || defined(BETAX)
        fldz
   #else
	fldl OFF(pC)
   #endif
	faddp %st,%st(7)
	fldl	-72(pA0)
	fmul %st(1),%st
	fxch %st(6)
	faddp %st,%st(4)
	fldl	-72(pA1)
	fmulp %st,%st(1)
	fldl	-64(pB0)
	fxch %st(7)
	faddp %st,%st(5)
	fldl	-64(pA0)
	fmul %st(7),%st
	fxch %st(4)
	faddp %st,%st(2)
	fldl	-64(pA1)
	fmulp %st,%st(7)
	fldl	-56(pB0)
	fxch %st(5)
	faddp %st,%st(3)
	fldl	-56(pA0)
	fmul %st(5),%st
	fxch %st(2)
	faddp %st,%st(6)
	fldl	-56(pA1)
	fmulp %st,%st(5)
	fldl	-48(pB0)
	fxch %st(3)
	faddp %st,%st(1)
	fldl	-48(pA0)
	fmul %st(3),%st
	fxch %st(6)
	faddp %st,%st(4)
	fldl	-48(pA1)
	fmulp %st,%st(3)
	fldl	-40(pB0)
	fxch %st(1)
	faddp %st,%st(7)
	fldl	-40(pA0)
	fmul %st(1),%st
	fxch %st(4)
	faddp %st,%st(2)
	fldl	-40(pA1)
	fmulp %st,%st(1)
	fldl	-32(pB0)
	fxch %st(7)
	faddp %st,%st(5)
	fldl	-32(pA0)
	fmul %st(7),%st
	fxch %st(2)
	faddp %st,%st(6)
	fldl	-32(pA1)
	fmulp %st,%st(7)
	fldl	-24(pB0)
	fxch %st(5)
	faddp %st,%st(3)
	fldl	-24(pA0)
	fmul %st(5),%st
	fxch %st(6)
	faddp %st,%st(4)
	fldl	-24(pA1)
	fmulp %st,%st(5)
	fldl	-16(pB0)
	fxch %st(3)
	faddp %st,%st(1)
	fldl	-16(pA0)
	fmul %st(3),%st
	fxch %st(4)
	faddp %st,%st(2)
	fldl	-16(pA1)
	fmulp %st,%st(3)
	fldl	-8(pB0)
	fxch %st(1)
	faddp %st,%st(7)
	fldl	-8(pA0)
	fmul %st(1),%st
	fxch %st(2)
	faddp %st,%st(6)
	fldl	-8(pA1)
	fmulp %st,%st(1)
	fldl	(pB0)
	fxch %st(7)
	faddp %st,%st(5)
	fldl	(pA0)
	fmul %st(7),%st
	fxch %st(6)
	faddp %st,%st(4)
	fldl	(pA1)
	fmulp %st,%st(7)
	fldl	8(pB0)
	fxch %st(5)
	faddp %st,%st(3)
	fldl	8(pA0)
	fmul %st(5),%st
	fxch %st(4)
	faddp %st,%st(2)
	fldl	8(pA1)
	fmulp %st,%st(5)
	fldl	16(pB0)
	fxch %st(3)
	faddp %st,%st(1)
	fldl	16(pA0)
	fmul %st(3),%st
	fxch %st(2)
	faddp %st,%st(6)
	fldl	16(pA1)
	fmulp %st,%st(3)
	fldl	24(pB0)
	fxch %st(1)
	faddp %st,%st(7)
	fldl	24(pA0)
	fmul %st(1),%st
	fxch %st(6)
	faddp %st,%st(4)
	fldl	24(pA1)
	fmulp %st,%st(1)
	fldl	32(pB0)
	fxch %st(7)
	faddp %st,%st(5)
	fldl	32(pA0)
	fmul %st(7),%st
	fxch %st(4)
	faddp %st,%st(2)
	fldl	32(pA1)
	fmulp %st,%st(7)
	fldl	40(pB0)
	fxch %st(5)
	faddp %st,%st(3)
	fldl	40(pA0)
	fmul %st(5),%st
	fxch %st(2)
	faddp %st,%st(6)
	fldl	40(pA1)
	fmulp %st,%st(5)
	fldl	48(pB0)
	fxch %st(3)
	faddp %st,%st(1)
	fldl	48(pA0)
	fmul %st(3),%st
	fxch %st(6)
	faddp %st,%st(4)
	fldl	48(pA1)
	fmulp %st,%st(3)
	fldl	56(pB0)
	fxch %st(1)
	faddp %st,%st(7)
	fldl	56(pA0)
	fmul %st(1),%st
	fxch %st(4)
	faddp %st,%st(2)
	fldl	56(pA1)
	fmulp %st,%st(1)
	fldl	64(pB0)
	fxch %st(7)
	faddp %st,%st(5)
	fldl	64(pA0)
	fmul %st(7),%st
	fxch %st(2)
	faddp %st,%st(6)
	fldl	64(pA1)
	fmulp %st,%st(7)
	fldl	72(pB0)
	fxch %st(5)
	faddp %st,%st(3)
	fldl	72(pA0)
	fmul %st(5),%st
	fxch %st(6)
	faddp %st,%st(4)
	fldl	72(pA1)
	fmulp %st,%st(5)
	fldl	80(pB0)
	fxch %st(3)
	faddp %st,%st(1)
	fldl	80(pA0)
	fmul %st(3),%st
	fxch %st(4)
	faddp %st,%st(2)
	fldl	80(pA1)
	fmulp %st,%st(3)
	fldl	88(pB0)
	fxch %st(1)
	faddp %st,%st(7)
	fldl	88(pA0)
	fmul %st(1),%st
	fxch %st(2)
	faddp %st,%st(6)
	fldl	88(pA1)
	fmulp %st,%st(1)
	fxch %st(6)
	addl	$NB2so, pA0
	faddp %st,%st(4)
	fxch %st(4)
	addl	$NB2so, pA1
	faddp %st,%st(2)
	faddp %st,%st(2)
	faddp %st,%st(2)
	faddp %st,%st(2)
/*
 *      Write results back to C
 */
   #ifdef BETAX
        fldl    (pC)
        fldl    OFF(pC)
        fldl    BETAOFF(%esp)
        fmul    %st, %st(1)
        fmulp   %st, %st(2)
        faddp   %st, %st(3)
        faddp   %st, %st(1)
   #endif
	fstpl (pC)
	fstpl OFF(pC)
/*
 *      pC += 2;  pA += 2*NB; pB -= NB;
 */
	addl	$incCm, pC
	fldl	-96(pB0)
	fldl	-96(pA0)
	fmul %st(1),%st
	fldl	-96(pA1)
	fmulp %st,%st(2)
	fldl	-88(pB0)
	fldl	-88(pA0)
	fmul %st(1),%st
	fldl	-88(pA1)
	fmulp %st,%st(2)
	fldl	-80(pB0)
	fldl	-80(pA0)
	fmul %st(1),%st
   #if defined(BETA0) || defined(BETAX)
        fldz
   #else
	fldl (pC)
   #endif
	faddp %st,%st(5)
	fldl	-80(pA1)
	fmulp %st,%st(2)
	fldl	-72(pB0)
   #if defined(BETA0) || defined(BETAX)
        fldz
   #else
	fldl OFF(pC)
   #endif
	faddp %st,%st(7)
	fldl	-72(pA0)
	fmul %st(1),%st
	fxch %st(6)
	faddp %st,%st(4)
	fldl	-72(pA1)
	fmulp %st,%st(1)
	fldl	-64(pB0)
	fxch %st(7)
	faddp %st,%st(5)
	fldl	-64(pA0)
	fmul %st(7),%st
	fxch %st(4)
	faddp %st,%st(2)
	fldl	-64(pA1)
	fmulp %st,%st(7)
	fldl	-56(pB0)
	fxch %st(5)
	faddp %st,%st(3)
	fldl	-56(pA0)
	fmul %st(5),%st
	fxch %st(2)
	faddp %st,%st(6)
	fldl	-56(pA1)
	fmulp %st,%st(5)
	fldl	-48(pB0)
	fxch %st(3)
	faddp %st,%st(1)
	fldl	-48(pA0)
	fmul %st(3),%st
	fxch %st(6)
	faddp %st,%st(4)
	fldl	-48(pA1)
	fmulp %st,%st(3)
	fldl	-40(pB0)
	fxch %st(1)
	faddp %st,%st(7)
	fldl	-40(pA0)
	fmul %st(1),%st
	fxch %st(4)
	faddp %st,%st(2)
	fldl	-40(pA1)
	fmulp %st,%st(1)
	fldl	-32(pB0)
	fxch %st(7)
	faddp %st,%st(5)
	fldl	-32(pA0)
	fmul %st(7),%st
	fxch %st(2)
	faddp %st,%st(6)
	fldl	-32(pA1)
	fmulp %st,%st(7)
	fldl	-24(pB0)
	fxch %st(5)
	faddp %st,%st(3)
	fldl	-24(pA0)
	fmul %st(5),%st
	fxch %st(6)
	faddp %st,%st(4)
	fldl	-24(pA1)
	fmulp %st,%st(5)
	fldl	-16(pB0)
	fxch %st(3)
	faddp %st,%st(1)
	fldl	-16(pA0)
	fmul %st(3),%st
	fxch %st(4)
	faddp %st,%st(2)
	fldl	-16(pA1)
	fmulp %st,%st(3)
	fldl	-8(pB0)
	fxch %st(1)
	faddp %st,%st(7)
	fldl	-8(pA0)
	fmul %st(1),%st
	fxch %st(2)
	faddp %st,%st(6)
	fldl	-8(pA1)
	fmulp %st,%st(1)
	fldl	(pB0)
	fxch %st(7)
	faddp %st,%st(5)
	fldl	(pA0)
	fmul %st(7),%st
	fxch %st(6)
	faddp %st,%st(4)
	fldl	(pA1)
	fmulp %st,%st(7)
	fldl	8(pB0)
	fxch %st(5)
	faddp %st,%st(3)
	fldl	8(pA0)
	fmul %st(5),%st
	fxch %st(4)
	faddp %st,%st(2)
	fldl	8(pA1)
	fmulp %st,%st(5)
	fldl	16(pB0)
	fxch %st(3)
	faddp %st,%st(1)
	fldl	16(pA0)
	fmul %st(3),%st
	fxch %st(2)
	faddp %st,%st(6)
	fldl	16(pA1)
	fmulp %st,%st(3)
	fldl	24(pB0)
	fxch %st(1)
	faddp %st,%st(7)
	fldl	24(pA0)
	fmul %st(1),%st
	fxch %st(6)
	faddp %st,%st(4)
	fldl	24(pA1)
	fmulp %st,%st(1)
	fldl	32(pB0)
	fxch %st(7)
	faddp %st,%st(5)
	fldl	32(pA0)
	fmul %st(7),%st
	fxch %st(4)
	faddp %st,%st(2)
	fldl	32(pA1)
	fmulp %st,%st(7)
	fldl	40(pB0)
	fxch %st(5)
	faddp %st,%st(3)
	fldl	40(pA0)
	fmul %st(5),%st
	fxch %st(2)
	faddp %st,%st(6)
	fldl	40(pA1)
	fmulp %st,%st(5)
	fldl	48(pB0)
	fxch %st(3)