int clapack_sgesv(const enum CBLAS_ORDER Order, const int N, const int NRHS,
float *A, const int lda, int *ipiv,
float *B, const int ldb)
/*
* clapack_sgesv computes the solution to a system of linear equations
* A * X = B,
* where A is an N-by-N matrix and X and B are N-by-NRHS matrices.
* The LU factorization used to factor A is dependent on the Order parameter,
* as detailed in the leading comments of clapack_sgetrf.
* The factored form of A is then used solve the system of equations A * X = B.
* A is overwritten with the appropriate LU factorization, and B, which
* contains B on input, is overwritten with the solution X on output.
*/
{
int ierr = 0;
if (Order != CblasRowMajor && Order != CblasColMajor)
{
ierr = -1;
cblas_xerbla(1, "clapack_sgesv",
"Order must be %d or %d, but is set to %d.\n",
CblasRowMajor, CblasColMajor, Order);
}
if (N < 0)
{
ierr = -2;
cblas_xerbla(2, "clapack_sgesv",
"N cannot be less than zero 0,; is set to %d.\n", N);
}
if (NRHS < 0)
{
ierr = -3;
cblas_xerbla(3, "clapack_sgesv",
"NRHS cannot be less than zero 0,; is set to %d.\n", NRHS);
}
if (lda < N || lda < 1)
{
ierr = -5;
cblas_xerbla(5, "clapack_sgesv",
"lda must be >= MAX(N,1): lda=%d N=%d\n", lda, N);
}
if (ldb < N || ldb < 1)
{
ierr = -8;
cblas_xerbla(8, "clapack_sgesv",
"ldb must be >= MAX(N,1): ldb=%d N=%d\n", ldb, N);
}
if (!ierr) ierr = ATL_sgetrf(Order, N, N, A, lda, ipiv);
if (!ierr) ATL_sgetrs(Order, CblasNoTrans, N, NRHS, A, lda, ipiv, B, ldb);
return(ierr);
}
int clapack_sgetrs
(const enum CBLAS_ORDER Order, const enum CBLAS_TRANSPOSE Trans,
const int N, const int NRHS, const float *A, const int lda,
const int *ipiv, float *B, const int ldb)
{
int ierr=0;
if (Order != CblasRowMajor && Order != CblasColMajor)
{
ierr = -1;
cblas_xerbla(1, "clapack_sgetrs",
"Order must be %d or %d, but is set to %d\n",
CblasRowMajor, CblasColMajor, Order);
}
if (Trans != CblasNoTrans && Trans != CblasTrans && Trans != CblasConjTrans)
{
ierr = -2;
cblas_xerbla(2, "clapack_sgetrs",
"Trans must be %d, %d, or %d, but is set to %d\n",
CblasNoTrans, CblasTrans, CblasConjTrans);
}
if (N < 0)
{
ierr = -3;
cblas_xerbla(3, "clapack_sgetrs",
"N cannot be less than zero 0,; is set to %d.\n", N);
}
if (NRHS < 0)
{
ierr = -4;
cblas_xerbla(4, "clapack_sgetrs",
"NRHS cannot be less than zero 0,; is set to %d.\n", NRHS);
}
if (lda < N || lda < 1)
{
ierr = -6;
cblas_xerbla(6, "clapack_sgetrs",
"lda must be >= MAX(N,1): lda=%d N=%d\n", lda, N);
}
if (ldb < N || ldb < 1)
{
ierr = -9;
cblas_xerbla(9, "clapack_sgetrs",
"ldb must be >= MAX(N,1): lda=%d N=%d\n", lda, N);
}
if (!ierr) ATL_sgetrs(Order, Trans, N, NRHS, A, lda, ipiv, B, ldb);
return(ierr);
}
