void ATL_zreftrmmRLTN
(
const int M,
const int N,
const double * ALPHA,
const double * A,
const int LDA,
double * B,
const int LDB
)
{
/*
* Purpose
* =======
*
* ATL_zreftrmmRLTN( ... )
*
* <=>
*
* ATL_zreftrmm
* ( AtlasRight, AtlasLower, AtlasTrans, AtlasNonUnit, ... )
*
* See ATL_zreftrmm for details.
*
* ---------------------------------------------------------------------
*/
/*
* .. Local Variables ..
*/
double t0_i, t0_r;
int i, iajk, ibij, ibik, j, jak, jbj, jbk, k,
lda2 = ( LDA << 1 ), ldb2 = ( LDB << 1 );
/* ..
* .. Executable Statements ..
*
*/
for( k = N-1, jak = (N-1)*lda2, jbk = (N-1)*ldb2;
k >= 0; k--, jak -= lda2, jbk -= ldb2 )
{
for( j = k+1, iajk = ( (k+1) << 1 )+jak, jbj = (k+1)*ldb2;
j < N; j++, iajk += 2, jbj += ldb2 )
{
Mmul( ALPHA[0], ALPHA[1], A[iajk], A[iajk+1], t0_r, t0_i );
for( i = 0, ibij = jbj, ibik = jbk;
i < M; i++, ibij += 2, ibik += 2 )
{ Mmla( t0_r, t0_i, B[ibik], B[ibik+1], B[ibij], B[ibij+1] ); }
}
Mset( ALPHA[0], ALPHA[1], t0_r, t0_i );
i = ( k << 1 ) + jak; Mdscl( A[i], A[i+1], t0_r, t0_i );
for( i = 0, ibik = jbk; i < M; i++, ibik += 2 )
{ Mdscl( t0_r, t0_i, B[ibik], B[ibik+1] ); }
}
/*
* End of ATL_zreftrmmRLTN
*/
}
void ATL_zreftrmmRUNN
(
const int M,
const int N,
const double * ALPHA,
const double * A,
const int LDA,
double * B,
const int LDB
)
{
/*
* Purpose
* =======
*
* ATL_zreftrmmRUNN( ... )
*
* <=>
*
* ATL_zreftrmm
* ( AtlasRight, AtlasUpper, AtlasNoTrans, AtlasNonUnit, ... )
*
* See ATL_zreftrmm for details.
*
* ---------------------------------------------------------------------
*/
/*
* .. Local Variables ..
*/
double t0_i, t0_r;
int i, iakj, ibij, ibik, j, jaj, jbj, jbk, k,
lda2 = ( LDA << 1 ), ldb2 = ( LDB << 1 );
/* ..
* .. Executable Statements ..
*
*/
for( j = N-1, jaj = (N-1)*lda2, jbj = (N-1)*ldb2;
j >= 0; j--, jaj -= lda2, jbj -= ldb2 )
{
Mset( ALPHA[0], ALPHA[1], t0_r, t0_i );
i = ( j << 1 ) + jaj; Mdscl( A[i], A[i+1], t0_r, t0_i );
for( i = 0, ibij = jbj; i < M; i++, ibij += 2 )
{ Mdscl( t0_r, t0_i, B[ibij], B[ibij+1] ); }
for( k = 0, iakj = jaj, jbk = 0; k < j; k++, iakj += 2, jbk += ldb2 )
{
Mmul( ALPHA[0], ALPHA[1], A[iakj], A[iakj+1], t0_r, t0_i );
for( i = 0, ibij = jbj, ibik = jbk; i < M; i++, ibij += 2, ibik += 2 )
{ Mmla( t0_r, t0_i, B[ibik], B[ibik+1], B[ibij], B[ibij+1] ); }
}
}
/*
* End of ATL_zreftrmmRUNN
*/
}
void ATL_zreftrsmRLCN
(
const int M,
const int N,
const double * ALPHA,
const double * A,
const int LDA,
double * B,
const int LDB
)
{
/*
* Purpose
* =======
*
* ATL_zreftrsmRLCN( ... )
*
* <=>
*
* ATL_zreftrsm
* ( AtlasRight, AtlasLower, AtlasConjTrans, AtlasNonUnit, ... )
*
* See ATL_zreftrsm for details.
*
* ---------------------------------------------------------------------
*/
/*
* .. Local Variables ..
*/
register double t0_i, t0_r;
int i, iajk, ibij, ibik, j, jak, jbj, jbk, k,
lda2 = ( LDA << 1 ), ldb2 = ( LDB << 1 );
/* ..
* .. Executable Statements ..
*
*/
for( k = 0, jak = 0, jbk = 0; k < N; k++, jak += lda2, jbk += ldb2 )
{
for( i = 0, ibik = jbk; i < M; i++, ibik += 2 )
{ j = ( k << 1 ) + jak; Mddiv( A[j], -A[j+1], B[ibik], B[ibik+1] ); }
for( j = k+1, iajk = ((k+1) << 1)+jak, jbj = (k+1)*ldb2;
j < N; j++, iajk += 2, jbj += ldb2 )
{
Mset( A[iajk], -A[iajk+1], t0_r, t0_i );
for( i = 0, ibij = jbj, ibik = jbk; i < M; i++, ibij += 2, ibik += 2 )
{ Mmls( t0_r, t0_i, B[ibik], B[ibik+1], B[ibij], B[ibij+1] ); }
}
for( i = 0, ibik = jbk; i < M; i++, ibik += 2 )
{ Mdscl( ALPHA[0], ALPHA[1], B[ibik], B[ibik+1] ); }
}
/*
* End of ATL_zreftrsmRLCN
*/
}
void ATL_zreftrsmLLNU
(
const int M,
const int N,
const double * ALPHA,
const double * A,
const int LDA,
double * B,
const int LDB
)
{
/*
* Purpose
* =======
*
* ATL_zreftrsmLLNU( ... )
*
* <=>
*
* ATL_zreftrsm
* ( AtlasLeft, AtlasLower, AtlasNoTrans, AtlasUnit, ... )
*
* See ATL_zreftrsm for details.
*
* ---------------------------------------------------------------------
*/
/*
* .. Local Variables ..
*/
int i, iaik, ibij, ibkj, j, jak, jbj, k,
lda2 = ( LDA << 1 ), ldb2 = ( LDB << 1 );
/* ..
* .. Executable Statements ..
*
*/
for( j = 0, jbj = 0; j < N; j++, jbj += ldb2 )
{
for( i = 0, ibij = jbj; i < M; i++, ibij += 2 )
{
Mdscl( ALPHA[0], ALPHA[1], B[ibij], B[ibij+1] );
}
for( k = 0, jak = 0, ibkj = jbj; k < M; k++, jak += lda2, ibkj += 2 )
{
for( i = k+1, iaik = ((k+1) << 1)+jak, ibij = ((k+1) << 1)+jbj;
i < M; i++, iaik += 2, ibij += 2 )
{
Mmls( A[iaik], A[iaik+1], B[ibkj], B[ibkj+1], B[ibij], B[ibij+1] );
}
}
}
/*
* End of ATL_zreftrsmLLNU
*/
}
void ATL_zreftrmmRUTU
(
const int M,
const int N,
const double * ALPHA,
const double * A,
const int LDA,
double * B,
const int LDB
)
{
/*
* Purpose
* =======
*
* ATL_zreftrmmRUTU( ... )
*
* <=>
*
* ATL_zreftrmm
* ( AtlasRight, AtlasUpper, AtlasTrans, AtlasUnit, ... )
*
* See ATL_zreftrmm for details.
*
* ---------------------------------------------------------------------
*/
/*
* .. Local Variables ..
*/
double t0_i, t0_r;
int i, iajk, ibij, ibik, j, jak, jbj, jbk, k,
lda2 = ( LDA << 1 ), ldb2 = ( LDB << 1 );
/* ..
* .. Executable Statements ..
*
*/
for( k = 0, jak = 0, jbk = 0; k < N; k++, jak += lda2, jbk += ldb2 )
{
for( j = 0, iajk = jak, jbj = 0; j < k; j++, iajk += 2, jbj += ldb2 )
{
Mmul( ALPHA[0], ALPHA[1], A[iajk], A[iajk+1], t0_r, t0_i );
for( i = 0, ibij = jbj, ibik = jbk;
i < M; i++, ibij += 2, ibik += 2 )
{ Mmla( t0_r, t0_i, B[ibik], B[ibik+1], B[ibij], B[ibij+1] ); }
}
Mset( ALPHA[0], ALPHA[1], t0_r, t0_i );
for( i = 0, ibik = jbk; i < M; i++, ibik += 2 )
{ Mdscl( t0_r, t0_i, B[ibik], B[ibik+1] ); }
}
/*
* End of ATL_zreftrmmRUTU
*/
}
void ATL_zreftrsmRLNN
(
const int M,
const int N,
const double * ALPHA,
const double * A,
const int LDA,
double * B,
const int LDB
)
{
/*
* Purpose
* =======
*
* ATL_zreftrsmRLNN( ... )
*
* <=>
*
* ATL_zreftrsm
* ( AtlasRight, AtlasLower, AtlasNoTrans, AtlasNonUnit, ... )
*
* See ATL_zreftrsm for details.
*
* ---------------------------------------------------------------------
*/
/*
* .. Local Variables ..
*/
int i, iakj, ibij, ibik, j, jaj, jbj, jbk, k,
lda2 = ( LDA << 1 ), ldb2 = ( LDB << 1 );
/* ..
* .. Executable Statements ..
*
*/
for( j = N-1, jaj = (N-1)*lda2, jbj = (N-1)*ldb2;
j >= 0; j--, jaj -= lda2, jbj -= ldb2 )
{
for( i = 0, ibij = jbj; i < M; i++, ibij += 2 )
{ Mdscl( ALPHA[0], ALPHA[1], B[ibij], B[ibij+1] ); }
for( k = j+1, iakj = ((j+1) << 1)+jaj, jbk = (j+1)*ldb2;
k < N; k++, iakj += 2, jbk += ldb2 )
{
for( i = 0, ibij = jbj, ibik = jbk; i < M; i++, ibij += 2, ibik += 2 )
{ Mmls( A[iakj], A[iakj+1], B[ibik], B[ibik+1], B[ibij], B[ibij+1] ); }
}
for( i = 0, ibij = jbj; i < M; i++, ibij += 2 )
{ k = ( j << 1 ) + jaj; Mddiv( A[k], A[k+1], B[ibij], B[ibij+1] ); }
}
/*
* End of ATL_zreftrsmRLNN
*/
}
void ATL_zreftpmvLCN
(
const int N,
const double * A,
const int LDA,
double * X,
const int INCX
)
{
/*
* Purpose
* =======
*
* ATL_zreftpmvLCN( ... )
*
* <=>
*
* ATL_zreftpmv( AtlasLower, AtlasConj, AtlasNonUnit, ... )
*
* See ATL_zreftpmv for details.
*
* ---------------------------------------------------------------------
*/
/*
* .. Local Variables ..
*/
register double t0_i, t0_r;
int i, iaij, incx2 = 2 * INCX, ix, j, jaj, jx,
lda2 = ( ( LDA - N + 1 ) << 1 );
/* ..
* .. Executable Statements ..
*
*/
for( j = N-1, jaj = (((LDA+1)*(N-1)) << 1) - (N-1)*N, jx = (N-1)*incx2;
j >= 0; j--, jx -= incx2 )
{
Mset( X[jx], X[jx+1], t0_r, t0_i );
Mdscl( A[jaj], -A[jaj+1], X[jx], X[jx+1] );
for( i = j+1, iaij = jaj+2, ix = jx+incx2;
i < N; i++, iaij += 2, ix += incx2 )
{ Mmla( A[iaij], -A[iaij+1], t0_r, t0_i, X[ix] , X[ix+1] ); }
lda2 += 2; jaj -= lda2;
}
/*
* End of ATL_zreftpmvLCN
*/
}
