// Solve A * x = b, input A and b.
void cblas_backsolver(double *A, double *x, double *b, int N)
{
// int i;
// for (i=0; i<N; i++) x[i] = b[i];
cblas_dcopy(N, b, 1, x, 1);
cblas_dtrsv(CblasRowMajor, CblasUpper, CblasNoTrans, CblasNonUnit, N, A, N, x, 1);
}
inline void trsv( const Order, const UpLo, const Trans, const Diag,
const int n, const double* a, const int lda, double* x,
const int incx ) {
cblas_dtrsv( cblas_option< Order >::value, cblas_option< UpLo >::value,
cblas_option< Trans >::value, cblas_option< Diag >::value, n, a,
lda, x, incx );
}
inline void trsv(
CBLAS_ORDER order, CBLAS_UPLO uplo,
CBLAS_TRANSPOSE transA, CBLAS_DIAG unit,
int n,
double const *A, int lda, double *b, int strideX
){
cblas_dtrsv(order, uplo, transA, unit,n, A, lda, b, strideX);
}
int main(int argc, char const *argv[])
{
cblas_dtrsv(order, uploA, transA, diagA, N, A, lda, X, incx);
printResult();
ContextDestroy();
return 0;
}
JNIEXPORT void JNICALL Java_uncomplicate_neanderthal_CBLAS_dtrsv
(JNIEnv *env, jclass clazz,
jint Order, jint Uplo, jint TransA, jint Diag,
jint N,
jobject A, jint offsetA, jint lda,
jobject X, jint offsetX, jint incX) {
double *cA = (double *) (*env)->GetDirectBufferAddress(env, A);
double *cX = (double *) (*env)->GetDirectBufferAddress(env, X);
cblas_dtrsv(Order, Uplo, TransA, Diag, N, cA + offsetA, lda, cX + offsetX, incX);
};
void VectorView::multiply_by_invtriangularmatrix(const SymmMatrixView &m,
const bool transpose)
{
CBLAS_TRANSPOSE tr;
if (transpose){
tr = CblasTrans;
} else {
tr = CblasNoTrans;
}
assert(m.length() == length());
cblas_dtrsv(CblasColMajor, CblasUpper, tr, CblasNonUnit,
m.length(), m.data(), m.stride(), data_, 1);
}
void F77_dtrsv(int *order, char *uplow, char *transp, char *diagn,
int *n, double *a, int *lda, double *x, int *incx ) {
double *A;
int i,j,LDA;
enum CBLAS_TRANSPOSE trans;
enum CBLAS_UPLO uplo;
enum CBLAS_DIAG diag;
get_transpose_type(transp,&trans);
get_uplo_type(uplow,&uplo);
get_diag_type(diagn,&diag);
if (*order == TEST_ROW_MJR) {
LDA = *n+1;
A = ( double* )malloc( (*n)*LDA*sizeof( double ) );
for( i=0; i<*n; i++ )
for( j=0; j<*n; j++ )
A[ LDA*i+j ]=a[ (*lda)*j+i ];
cblas_dtrsv(CblasRowMajor, uplo, trans, diag, *n, A, LDA, x, *incx );
free(A);
}
else
cblas_dtrsv(CblasColMajor, uplo, trans, diag, *n, a, *lda, x, *incx );
}
int
gsl_blas_dtrsv (CBLAS_UPLO_t Uplo, CBLAS_TRANSPOSE_t TransA,
CBLAS_DIAG_t Diag, const gsl_matrix * A, gsl_vector * X)
{
const size_t M = A->size1;
const size_t N = A->size2;
if (M != N)
{
GSL_ERROR ("matrix must be square", GSL_ENOTSQR);
}
else if (N != X->size)
{
GSL_ERROR ("invalid length", GSL_EBADLEN);
}
cblas_dtrsv (CblasRowMajor, Uplo, TransA, Diag, INT (N), A->data,
INT (A->tda), X->data, INT (X->stride));
return GSL_SUCCESS;
}
void solve_triangle_eq_mkl()
{
char *results_file = "solve_triangle_eq_mkl.txt";
FILE *res;
if((res=fopen(results_file, "w"))==NULL)
{
printf("Can't open file %s.\n", results_file);
exit(1);
}
for(int i = 10; i < ROW; i*=10)
{
double *A = get_vector(i*i,i);
double *X = get_vector(i,i);
double start = omp_get_wtime();
cblas_dtrsv(CblasRowMajor, CblasUpper, CblasNoTrans, CblasNonUnit, i, A, i, X, 1);
double end = omp_get_wtime();
fprintf(res, "%lf\n", end-start);
free(X);free(A);
}
fclose(res);
}
void F77_d2chke(char *rout) {
char *sf = ( rout ) ;
double A[2] = {0.0,0.0},
X[2] = {0.0,0.0},
Y[2] = {0.0,0.0},
ALPHA=0.0, BETA=0.0;
extern int cblas_info, cblas_lerr, cblas_ok;
extern int RowMajorStrg;
extern char *cblas_rout;
if (link_xerbla) /* call these first to link */
{
cblas_xerbla(cblas_info,cblas_rout,"");
F77_xerbla(cblas_rout,&cblas_info);
}
cblas_ok = TRUE ;
cblas_lerr = PASSED ;
if (strncmp( sf,"cblas_dgemv",11)==0) {
cblas_rout = "cblas_dgemv";
cblas_info = 1;
cblas_dgemv(INVALID, CblasNoTrans, 0, 0,
ALPHA, A, 1, X, 1, BETA, Y, 1 );
chkxer();
cblas_info = 2; RowMajorStrg = FALSE;
cblas_dgemv(CblasColMajor, INVALID, 0, 0,
ALPHA, A, 1, X, 1, BETA, Y, 1 );
chkxer();
cblas_info = 3; RowMajorStrg = FALSE;
cblas_dgemv(CblasColMajor, CblasNoTrans, INVALID, 0,
ALPHA, A, 1, X, 1, BETA, Y, 1 );
chkxer();
cblas_info = 4; RowMajorStrg = FALSE;
cblas_dgemv(CblasColMajor, CblasNoTrans, 0, INVALID,
ALPHA, A, 1, X, 1, BETA, Y, 1 );
chkxer();
cblas_info = 7; RowMajorStrg = FALSE;
cblas_dgemv(CblasColMajor, CblasNoTrans, 2, 0,
ALPHA, A, 1, X, 1, BETA, Y, 1 );
chkxer();
cblas_info = 9; RowMajorStrg = FALSE;
cblas_dgemv(CblasColMajor, CblasNoTrans, 0, 0,
ALPHA, A, 1, X, 0, BETA, Y, 1 );
chkxer();
cblas_info = 12; RowMajorStrg = FALSE;
cblas_dgemv(CblasColMajor, CblasNoTrans, 0, 0,
ALPHA, A, 1, X, 1, BETA, Y, 0 );
chkxer();
cblas_info = 2; RowMajorStrg = TRUE; RowMajorStrg = TRUE;
cblas_dgemv(CblasRowMajor, INVALID, 0, 0,
ALPHA, A, 1, X, 1, BETA, Y, 1 );
chkxer();
cblas_info = 3; RowMajorStrg = TRUE;
cblas_dgemv(CblasRowMajor, CblasNoTrans, INVALID, 0,
ALPHA, A, 1, X, 1, BETA, Y, 1 );
chkxer();
cblas_info = 4; RowMajorStrg = TRUE;
cblas_dgemv(CblasRowMajor, CblasNoTrans, 0, INVALID,
ALPHA, A, 1, X, 1, BETA, Y, 1 );
chkxer();
cblas_info = 7; RowMajorStrg = TRUE;
cblas_dgemv(CblasRowMajor, CblasNoTrans, 0, 2,
ALPHA, A, 1, X, 1, BETA, Y, 1 );
chkxer();
cblas_info = 9; RowMajorStrg = TRUE;
cblas_dgemv(CblasRowMajor, CblasNoTrans, 0, 0,
ALPHA, A, 1, X, 0, BETA, Y, 1 );
chkxer();
cblas_info = 12; RowMajorStrg = TRUE;
cblas_dgemv(CblasRowMajor, CblasNoTrans, 0, 0,
ALPHA, A, 1, X, 1, BETA, Y, 0 );
chkxer();
} else if (strncmp( sf,"cblas_dgbmv",11)==0) {
cblas_rout = "cblas_dgbmv";
cblas_info = 1; RowMajorStrg = FALSE;
cblas_dgbmv(INVALID, CblasNoTrans, 0, 0, 0, 0,
ALPHA, A, 1, X, 1, BETA, Y, 1 );
chkxer();
cblas_info = 2; RowMajorStrg = FALSE;
cblas_dgbmv(CblasColMajor, INVALID, 0, 0, 0, 0,
ALPHA, A, 1, X, 1, BETA, Y, 1 );
chkxer();
cblas_info = 3; RowMajorStrg = FALSE;
cblas_dgbmv(CblasColMajor, CblasNoTrans, INVALID, 0, 0, 0,
ALPHA, A, 1, X, 1, BETA, Y, 1 );
chkxer();
cblas_info = 4; RowMajorStrg = FALSE;
cblas_dgbmv(CblasColMajor, CblasNoTrans, 0, INVALID, 0, 0,
ALPHA, A, 1, X, 1, BETA, Y, 1 );
chkxer();
cblas_info = 5; RowMajorStrg = FALSE;
cblas_dgbmv(CblasColMajor, CblasNoTrans, 0, 0, INVALID, 0,
ALPHA, A, 1, X, 1, BETA, Y, 1 );
chkxer();
cblas_info = 6; RowMajorStrg = FALSE;
cblas_dgbmv(CblasColMajor, CblasNoTrans, 2, 0, 0, INVALID,
ALPHA, A, 1, X, 1, BETA, Y, 1 );
chkxer();
cblas_info = 9; RowMajorStrg = FALSE;
cblas_dgbmv(CblasColMajor, CblasNoTrans, 0, 0, 1, 0,
ALPHA, A, 1, X, 1, BETA, Y, 1 );
chkxer();
cblas_info = 11; RowMajorStrg = FALSE;
cblas_dgbmv(CblasColMajor, CblasNoTrans, 0, 0, 0, 0,
ALPHA, A, 1, X, 0, BETA, Y, 1 );
chkxer();
cblas_info = 14; RowMajorStrg = FALSE;
cblas_dgbmv(CblasColMajor, CblasNoTrans, 0, 0, 0, 0,
ALPHA, A, 1, X, 1, BETA, Y, 0 );
chkxer();
cblas_info = 2; RowMajorStrg = TRUE;
cblas_dgbmv(CblasRowMajor, INVALID, 0, 0, 0, 0,
ALPHA, A, 1, X, 1, BETA, Y, 1 );
chkxer();
cblas_info = 3; RowMajorStrg = TRUE;
cblas_dgbmv(CblasRowMajor, CblasNoTrans, INVALID, 0, 0, 0,
ALPHA, A, 1, X, 1, BETA, Y, 1 );
chkxer();
cblas_info = 4; RowMajorStrg = TRUE;
cblas_dgbmv(CblasRowMajor, CblasNoTrans, 0, INVALID, 0, 0,
ALPHA, A, 1, X, 1, BETA, Y, 1 );
chkxer();
cblas_info = 5; RowMajorStrg = TRUE;
cblas_dgbmv(CblasRowMajor, CblasNoTrans, 0, 0, INVALID, 0,
ALPHA, A, 1, X, 1, BETA, Y, 1 );
chkxer();
cblas_info = 6; RowMajorStrg = TRUE;
cblas_dgbmv(CblasRowMajor, CblasNoTrans, 2, 0, 0, INVALID,
ALPHA, A, 1, X, 1, BETA, Y, 1 );
chkxer();
cblas_info = 9; RowMajorStrg = TRUE;
cblas_dgbmv(CblasRowMajor, CblasNoTrans, 0, 0, 1, 0,
ALPHA, A, 1, X, 1, BETA, Y, 1 );
chkxer();
cblas_info = 11; RowMajorStrg = TRUE;
cblas_dgbmv(CblasRowMajor, CblasNoTrans, 0, 0, 0, 0,
ALPHA, A, 1, X, 0, BETA, Y, 1 );
chkxer();
cblas_info = 14; RowMajorStrg = TRUE;
cblas_dgbmv(CblasRowMajor, CblasNoTrans, 0, 0, 0, 0,
ALPHA, A, 1, X, 1, BETA, Y, 0 );
chkxer();
} else if (strncmp( sf,"cblas_dsymv",11)==0) {
cblas_rout = "cblas_dsymv";
cblas_info = 1; RowMajorStrg = FALSE;
cblas_dsymv(INVALID, CblasUpper, 0,
ALPHA, A, 1, X, 1, BETA, Y, 1 );
chkxer();
cblas_info = 2; RowMajorStrg = FALSE;
cblas_dsymv(CblasColMajor, INVALID, 0,
ALPHA, A, 1, X, 1, BETA, Y, 1 );
chkxer();
cblas_info = 3; RowMajorStrg = FALSE;
cblas_dsymv(CblasColMajor, CblasUpper, INVALID,
ALPHA, A, 1, X, 1, BETA, Y, 1 );
chkxer();
cblas_info = 6; RowMajorStrg = FALSE;
cblas_dsymv(CblasColMajor, CblasUpper, 2,
ALPHA, A, 1, X, 1, BETA, Y, 1 );
chkxer();
cblas_info = 8; RowMajorStrg = FALSE;
cblas_dsymv(CblasColMajor, CblasUpper, 0,
ALPHA, A, 1, X, 0, BETA, Y, 1 );
chkxer();
cblas_info = 11; RowMajorStrg = FALSE;
cblas_dsymv(CblasColMajor, CblasUpper, 0,
ALPHA, A, 1, X, 1, BETA, Y, 0 );
chkxer();
cblas_info = 2; RowMajorStrg = TRUE;
cblas_dsymv(CblasRowMajor, INVALID, 0,
ALPHA, A, 1, X, 1, BETA, Y, 1 );
chkxer();
cblas_info = 3; RowMajorStrg = TRUE;
cblas_dsymv(CblasRowMajor, CblasUpper, INVALID,
ALPHA, A, 1, X, 1, BETA, Y, 1 );
chkxer();
cblas_info = 6; RowMajorStrg = TRUE;
cblas_dsymv(CblasRowMajor, CblasUpper, 2,
ALPHA, A, 1, X, 1, BETA, Y, 1 );
chkxer();
cblas_info = 8; RowMajorStrg = TRUE;
cblas_dsymv(CblasRowMajor, CblasUpper, 0,
ALPHA, A, 1, X, 0, BETA, Y, 1 );
chkxer();
cblas_info = 11; RowMajorStrg = TRUE;
cblas_dsymv(CblasRowMajor, CblasUpper, 0,
ALPHA, A, 1, X, 1, BETA, Y, 0 );
chkxer();
} else if (strncmp( sf,"cblas_dsbmv",11)==0) {
cblas_rout = "cblas_dsbmv";
cblas_info = 1; RowMajorStrg = FALSE;
cblas_dsbmv(INVALID, CblasUpper, 0, 0,
ALPHA, A, 1, X, 1, BETA, Y, 1 );
chkxer();
cblas_info = 2; RowMajorStrg = FALSE;
cblas_dsbmv(CblasColMajor, INVALID, 0, 0,
ALPHA, A, 1, X, 1, BETA, Y, 1 );
chkxer();
cblas_info = 3; RowMajorStrg = FALSE;
cblas_dsbmv(CblasColMajor, CblasUpper, INVALID, 0,
ALPHA, A, 1, X, 1, BETA, Y, 1 );
chkxer();
cblas_info = 4; RowMajorStrg = FALSE;
cblas_dsbmv(CblasColMajor, CblasUpper, 0, INVALID,
ALPHA, A, 1, X, 1, BETA, Y, 1 );
chkxer();
cblas_info = 7; RowMajorStrg = FALSE;
cblas_dsbmv(CblasColMajor, CblasUpper, 0, 1,
ALPHA, A, 1, X, 1, BETA, Y, 1 );
chkxer();
cblas_info = 9; RowMajorStrg = FALSE;
cblas_dsbmv(CblasColMajor, CblasUpper, 0, 0,
ALPHA, A, 1, X, 0, BETA, Y, 1 );
chkxer();
cblas_info = 12; RowMajorStrg = FALSE;
cblas_dsbmv(CblasColMajor, CblasUpper, 0, 0,
ALPHA, A, 1, X, 1, BETA, Y, 0 );
chkxer();
cblas_info = 2; RowMajorStrg = TRUE;
cblas_dsbmv(CblasRowMajor, INVALID, 0, 0,
ALPHA, A, 1, X, 1, BETA, Y, 1 );
chkxer();
cblas_info = 3; RowMajorStrg = TRUE;
cblas_dsbmv(CblasRowMajor, CblasUpper, INVALID, 0,
ALPHA, A, 1, X, 1, BETA, Y, 1 );
chkxer();
cblas_info = 4; RowMajorStrg = TRUE;
cblas_dsbmv(CblasRowMajor, CblasUpper, 0, INVALID,
ALPHA, A, 1, X, 1, BETA, Y, 1 );
chkxer();
cblas_info = 7; RowMajorStrg = TRUE;
cblas_dsbmv(CblasRowMajor, CblasUpper, 0, 1,
ALPHA, A, 1, X, 1, BETA, Y, 1 );
chkxer();
cblas_info = 9; RowMajorStrg = TRUE;
cblas_dsbmv(CblasRowMajor, CblasUpper, 0, 0,
ALPHA, A, 1, X, 0, BETA, Y, 1 );
chkxer();
cblas_info = 12; RowMajorStrg = TRUE;
cblas_dsbmv(CblasRowMajor, CblasUpper, 0, 0,
ALPHA, A, 1, X, 1, BETA, Y, 0 );
chkxer();
} else if (strncmp( sf,"cblas_dspmv",11)==0) {
cblas_rout = "cblas_dspmv";
cblas_info = 1; RowMajorStrg = FALSE;
cblas_dspmv(INVALID, CblasUpper, 0,
ALPHA, A, X, 1, BETA, Y, 1 );
chkxer();
cblas_info = 2; RowMajorStrg = FALSE;
cblas_dspmv(CblasColMajor, INVALID, 0,
ALPHA, A, X, 1, BETA, Y, 1 );
chkxer();
cblas_info = 3; RowMajorStrg = FALSE;
cblas_dspmv(CblasColMajor, CblasUpper, INVALID,
ALPHA, A, X, 1, BETA, Y, 1 );
chkxer();
cblas_info = 7; RowMajorStrg = FALSE;
cblas_dspmv(CblasColMajor, CblasUpper, 0,
ALPHA, A, X, 0, BETA, Y, 1 );
chkxer();
cblas_info = 10; RowMajorStrg = FALSE;
cblas_dspmv(CblasColMajor, CblasUpper, 0,
ALPHA, A, X, 1, BETA, Y, 0 );
chkxer();
cblas_info = 2; RowMajorStrg = TRUE;
cblas_dspmv(CblasRowMajor, INVALID, 0,
ALPHA, A, X, 1, BETA, Y, 1 );
chkxer();
cblas_info = 3; RowMajorStrg = TRUE;
cblas_dspmv(CblasRowMajor, CblasUpper, INVALID,
ALPHA, A, X, 1, BETA, Y, 1 );
chkxer();
cblas_info = 7; RowMajorStrg = TRUE;
cblas_dspmv(CblasRowMajor, CblasUpper, 0,
ALPHA, A, X, 0, BETA, Y, 1 );
chkxer();
cblas_info = 10; RowMajorStrg = TRUE;
cblas_dspmv(CblasRowMajor, CblasUpper, 0,
ALPHA, A, X, 1, BETA, Y, 0 );
chkxer();
} else if (strncmp( sf,"cblas_dtrmv",11)==0) {
cblas_rout = "cblas_dtrmv";
cblas_info = 1; RowMajorStrg = FALSE;
cblas_dtrmv(INVALID, CblasUpper, CblasNoTrans,
CblasNonUnit, 0, A, 1, X, 1 );
chkxer();
cblas_info = 2; RowMajorStrg = FALSE;
cblas_dtrmv(CblasColMajor, INVALID, CblasNoTrans,
CblasNonUnit, 0, A, 1, X, 1 );
chkxer();
cblas_info = 3; RowMajorStrg = FALSE;
cblas_dtrmv(CblasColMajor, CblasUpper, INVALID,
CblasNonUnit, 0, A, 1, X, 1 );
chkxer();
cblas_info = 4; RowMajorStrg = FALSE;
cblas_dtrmv(CblasColMajor, CblasUpper, CblasNoTrans,
INVALID, 0, A, 1, X, 1 );
chkxer();
cblas_info = 5; RowMajorStrg = FALSE;
cblas_dtrmv(CblasColMajor, CblasUpper, CblasNoTrans,
CblasNonUnit, INVALID, A, 1, X, 1 );
chkxer();
cblas_info = 7; RowMajorStrg = FALSE;
cblas_dtrmv(CblasColMajor, CblasUpper, CblasNoTrans,
CblasNonUnit, 2, A, 1, X, 1 );
chkxer();
cblas_info = 9; RowMajorStrg = FALSE;
cblas_dtrmv(CblasColMajor, CblasUpper, CblasNoTrans,
CblasNonUnit, 0, A, 1, X, 0 );
chkxer();
cblas_info = 2; RowMajorStrg = TRUE;
cblas_dtrmv(CblasRowMajor, INVALID, CblasNoTrans,
CblasNonUnit, 0, A, 1, X, 1 );
chkxer();
cblas_info = 3; RowMajorStrg = TRUE;
cblas_dtrmv(CblasRowMajor, CblasUpper, INVALID,
CblasNonUnit, 0, A, 1, X, 1 );
chkxer();
cblas_info = 4; RowMajorStrg = TRUE;
cblas_dtrmv(CblasRowMajor, CblasUpper, CblasNoTrans,
INVALID, 0, A, 1, X, 1 );
chkxer();
cblas_info = 5; RowMajorStrg = TRUE;
cblas_dtrmv(CblasRowMajor, CblasUpper, CblasNoTrans,
CblasNonUnit, INVALID, A, 1, X, 1 );
chkxer();
cblas_info = 7; RowMajorStrg = TRUE;
cblas_dtrmv(CblasRowMajor, CblasUpper, CblasNoTrans,
CblasNonUnit, 2, A, 1, X, 1 );
chkxer();
cblas_info = 9; RowMajorStrg = TRUE;
cblas_dtrmv(CblasRowMajor, CblasUpper, CblasNoTrans,
CblasNonUnit, 0, A, 1, X, 0 );
chkxer();
} else if (strncmp( sf,"cblas_dtbmv",11)==0) {
cblas_rout = "cblas_dtbmv";
cblas_info = 1; RowMajorStrg = FALSE;
cblas_dtbmv(INVALID, CblasUpper, CblasNoTrans,
CblasNonUnit, 0, 0, A, 1, X, 1 );
chkxer();
cblas_info = 2; RowMajorStrg = FALSE;
cblas_dtbmv(CblasColMajor, INVALID, CblasNoTrans,
CblasNonUnit, 0, 0, A, 1, X, 1 );
chkxer();
cblas_info = 3; RowMajorStrg = FALSE;
cblas_dtbmv(CblasColMajor, CblasUpper, INVALID,
CblasNonUnit, 0, 0, A, 1, X, 1 );
chkxer();
cblas_info = 4; RowMajorStrg = FALSE;
cblas_dtbmv(CblasColMajor, CblasUpper, CblasNoTrans,
INVALID, 0, 0, A, 1, X, 1 );
chkxer();
cblas_info = 5; RowMajorStrg = FALSE;
cblas_dtbmv(CblasColMajor, CblasUpper, CblasNoTrans,
CblasNonUnit, INVALID, 0, A, 1, X, 1 );
chkxer();
cblas_info = 6; RowMajorStrg = FALSE;
cblas_dtbmv(CblasColMajor, CblasUpper, CblasNoTrans,
CblasNonUnit, 0, INVALID, A, 1, X, 1 );
chkxer();
cblas_info = 8; RowMajorStrg = FALSE;
cblas_dtbmv(CblasColMajor, CblasUpper, CblasNoTrans,
CblasNonUnit, 0, 1, A, 1, X, 1 );
chkxer();
cblas_info = 10; RowMajorStrg = FALSE;
cblas_dtbmv(CblasColMajor, CblasUpper, CblasNoTrans,
CblasNonUnit, 0, 0, A, 1, X, 0 );
chkxer();
cblas_info = 2; RowMajorStrg = TRUE;
cblas_dtbmv(CblasRowMajor, INVALID, CblasNoTrans,
CblasNonUnit, 0, 0, A, 1, X, 1 );
chkxer();
cblas_info = 3; RowMajorStrg = TRUE;
cblas_dtbmv(CblasRowMajor, CblasUpper, INVALID,
CblasNonUnit, 0, 0, A, 1, X, 1 );
chkxer();
cblas_info = 4; RowMajorStrg = TRUE;
cblas_dtbmv(CblasRowMajor, CblasUpper, CblasNoTrans,
INVALID, 0, 0, A, 1, X, 1 );
chkxer();
cblas_info = 5; RowMajorStrg = TRUE;
cblas_dtbmv(CblasRowMajor, CblasUpper, CblasNoTrans,
CblasNonUnit, INVALID, 0, A, 1, X, 1 );
chkxer();
cblas_info = 6; RowMajorStrg = TRUE;
cblas_dtbmv(CblasRowMajor, CblasUpper, CblasNoTrans,
CblasNonUnit, 0, INVALID, A, 1, X, 1 );
chkxer();
cblas_info = 8; RowMajorStrg = TRUE;
cblas_dtbmv(CblasRowMajor, CblasUpper, CblasNoTrans,
CblasNonUnit, 0, 1, A, 1, X, 1 );
chkxer();
cblas_info = 10; RowMajorStrg = TRUE;
cblas_dtbmv(CblasRowMajor, CblasUpper, CblasNoTrans,
CblasNonUnit, 0, 0, A, 1, X, 0 );
chkxer();
} else if (strncmp( sf,"cblas_dtpmv",11)==0) {
cblas_rout = "cblas_dtpmv";
cblas_info = 1; RowMajorStrg = FALSE;
cblas_dtpmv(INVALID, CblasUpper, CblasNoTrans,
CblasNonUnit, 0, A, X, 1 );
chkxer();
cblas_info = 2; RowMajorStrg = FALSE;
cblas_dtpmv(CblasColMajor, INVALID, CblasNoTrans,
CblasNonUnit, 0, A, X, 1 );
chkxer();
cblas_info = 3; RowMajorStrg = FALSE;
cblas_dtpmv(CblasColMajor, CblasUpper, INVALID,
CblasNonUnit, 0, A, X, 1 );
chkxer();
cblas_info = 4; RowMajorStrg = FALSE;
cblas_dtpmv(CblasColMajor, CblasUpper, CblasNoTrans,
INVALID, 0, A, X, 1 );
chkxer();
cblas_info = 5; RowMajorStrg = FALSE;
cblas_dtpmv(CblasColMajor, CblasUpper, CblasNoTrans,
CblasNonUnit, INVALID, A, X, 1 );
chkxer();
cblas_info = 8; RowMajorStrg = FALSE;
cblas_dtpmv(CblasColMajor, CblasUpper, CblasNoTrans,
CblasNonUnit, 0, A, X, 0 );
chkxer();
cblas_info = 2; RowMajorStrg = TRUE;
cblas_dtpmv(CblasRowMajor, INVALID, CblasNoTrans,
CblasNonUnit, 0, A, X, 1 );
chkxer();
cblas_info = 3; RowMajorStrg = TRUE;
cblas_dtpmv(CblasRowMajor, CblasUpper, INVALID,
CblasNonUnit, 0, A, X, 1 );
chkxer();
cblas_info = 4; RowMajorStrg = TRUE;
cblas_dtpmv(CblasRowMajor, CblasUpper, CblasNoTrans,
INVALID, 0, A, X, 1 );
chkxer();
cblas_info = 5; RowMajorStrg = TRUE;
cblas_dtpmv(CblasRowMajor, CblasUpper, CblasNoTrans,
CblasNonUnit, INVALID, A, X, 1 );
chkxer();
cblas_info = 8; RowMajorStrg = TRUE;
cblas_dtpmv(CblasRowMajor, CblasUpper, CblasNoTrans,
CblasNonUnit, 0, A, X, 0 );
chkxer();
} else if (strncmp( sf,"cblas_dtrsv",11)==0) {
cblas_rout = "cblas_dtrsv";
cblas_info = 1; RowMajorStrg = FALSE;
cblas_dtrsv(INVALID, CblasUpper, CblasNoTrans,
CblasNonUnit, 0, A, 1, X, 1 );
chkxer();
cblas_info = 2; RowMajorStrg = FALSE;
cblas_dtrsv(CblasColMajor, INVALID, CblasNoTrans,
CblasNonUnit, 0, A, 1, X, 1 );
chkxer();
cblas_info = 3; RowMajorStrg = FALSE;
cblas_dtrsv(CblasColMajor, CblasUpper, INVALID,
CblasNonUnit, 0, A, 1, X, 1 );
chkxer();
cblas_info = 4; RowMajorStrg = FALSE;
cblas_dtrsv(CblasColMajor, CblasUpper, CblasNoTrans,
INVALID, 0, A, 1, X, 1 );
chkxer();
cblas_info = 5; RowMajorStrg = FALSE;
cblas_dtrsv(CblasColMajor, CblasUpper, CblasNoTrans,
CblasNonUnit, INVALID, A, 1, X, 1 );
chkxer();
cblas_info = 7; RowMajorStrg = FALSE;
cblas_dtrsv(CblasColMajor, CblasUpper, CblasNoTrans,
CblasNonUnit, 2, A, 1, X, 1 );
chkxer();
cblas_info = 9; RowMajorStrg = FALSE;
cblas_dtrsv(CblasColMajor, CblasUpper, CblasNoTrans,
CblasNonUnit, 0, A, 1, X, 0 );
chkxer();
cblas_info = 2; RowMajorStrg = TRUE;
cblas_dtrsv(CblasRowMajor, INVALID, CblasNoTrans,
CblasNonUnit, 0, A, 1, X, 1 );
chkxer();
cblas_info = 3; RowMajorStrg = TRUE;
cblas_dtrsv(CblasRowMajor, CblasUpper, INVALID,
CblasNonUnit, 0, A, 1, X, 1 );
chkxer();
cblas_info = 4; RowMajorStrg = TRUE;
cblas_dtrsv(CblasRowMajor, CblasUpper, CblasNoTrans,
INVALID, 0, A, 1, X, 1 );
chkxer();
cblas_info = 5; RowMajorStrg = TRUE;
cblas_dtrsv(CblasRowMajor, CblasUpper, CblasNoTrans,
CblasNonUnit, INVALID, A, 1, X, 1 );
chkxer();
cblas_info = 7; RowMajorStrg = TRUE;
cblas_dtrsv(CblasRowMajor, CblasUpper, CblasNoTrans,
CblasNonUnit, 2, A, 1, X, 1 );
chkxer();
cblas_info = 9; RowMajorStrg = TRUE;
cblas_dtrsv(CblasRowMajor, CblasUpper, CblasNoTrans,
CblasNonUnit, 0, A, 1, X, 0 );
chkxer();
} else if (strncmp( sf,"cblas_dtbsv",11)==0) {
cblas_rout = "cblas_dtbsv";
cblas_info = 1; RowMajorStrg = FALSE;
cblas_dtbsv(INVALID, CblasUpper, CblasNoTrans,
CblasNonUnit, 0, 0, A, 1, X, 1 );
chkxer();
cblas_info = 2; RowMajorStrg = FALSE;
cblas_dtbsv(CblasColMajor, INVALID, CblasNoTrans,
CblasNonUnit, 0, 0, A, 1, X, 1 );
chkxer();
cblas_info = 3; RowMajorStrg = FALSE;
cblas_dtbsv(CblasColMajor, CblasUpper, INVALID,
CblasNonUnit, 0, 0, A, 1, X, 1 );
chkxer();
cblas_info = 4; RowMajorStrg = FALSE;
cblas_dtbsv(CblasColMajor, CblasUpper, CblasNoTrans,
INVALID, 0, 0, A, 1, X, 1 );
chkxer();
cblas_info = 5; RowMajorStrg = FALSE;
cblas_dtbsv(CblasColMajor, CblasUpper, CblasNoTrans,
CblasNonUnit, INVALID, 0, A, 1, X, 1 );
chkxer();
cblas_info = 6; RowMajorStrg = FALSE;
cblas_dtbsv(CblasColMajor, CblasUpper, CblasNoTrans,
CblasNonUnit, 0, INVALID, A, 1, X, 1 );
chkxer();
cblas_info = 8; RowMajorStrg = FALSE;
cblas_dtbsv(CblasColMajor, CblasUpper, CblasNoTrans,
CblasNonUnit, 0, 1, A, 1, X, 1 );
chkxer();
cblas_info = 10; RowMajorStrg = FALSE;
cblas_dtbsv(CblasColMajor, CblasUpper, CblasNoTrans,
CblasNonUnit, 0, 0, A, 1, X, 0 );
chkxer();
cblas_info = 2; RowMajorStrg = TRUE;
cblas_dtbsv(CblasRowMajor, INVALID, CblasNoTrans,
CblasNonUnit, 0, 0, A, 1, X, 1 );
chkxer();
cblas_info = 3; RowMajorStrg = TRUE;
cblas_dtbsv(CblasRowMajor, CblasUpper, INVALID,
CblasNonUnit, 0, 0, A, 1, X, 1 );
chkxer();
cblas_info = 4; RowMajorStrg = TRUE;
cblas_dtbsv(CblasRowMajor, CblasUpper, CblasNoTrans,
INVALID, 0, 0, A, 1, X, 1 );
chkxer();
cblas_info = 5; RowMajorStrg = TRUE;
cblas_dtbsv(CblasRowMajor, CblasUpper, CblasNoTrans,
CblasNonUnit, INVALID, 0, A, 1, X, 1 );
chkxer();
cblas_info = 6; RowMajorStrg = TRUE;
cblas_dtbsv(CblasRowMajor, CblasUpper, CblasNoTrans,
CblasNonUnit, 0, INVALID, A, 1, X, 1 );
chkxer();
cblas_info = 8; RowMajorStrg = TRUE;
cblas_dtbsv(CblasRowMajor, CblasUpper, CblasNoTrans,
CblasNonUnit, 0, 1, A, 1, X, 1 );
chkxer();
cblas_info = 10; RowMajorStrg = TRUE;
cblas_dtbsv(CblasRowMajor, CblasUpper, CblasNoTrans,
CblasNonUnit, 0, 0, A, 1, X, 0 );
chkxer();
} else if (strncmp( sf,"cblas_dtpsv",11)==0) {
cblas_rout = "cblas_dtpsv";
cblas_info = 1; RowMajorStrg = FALSE;
cblas_dtpsv(INVALID, CblasUpper, CblasNoTrans,
CblasNonUnit, 0, A, X, 1 );
chkxer();
cblas_info = 2; RowMajorStrg = FALSE;
cblas_dtpsv(CblasColMajor, INVALID, CblasNoTrans,
CblasNonUnit, 0, A, X, 1 );
chkxer();
cblas_info = 3; RowMajorStrg = FALSE;
cblas_dtpsv(CblasColMajor, CblasUpper, INVALID,
CblasNonUnit, 0, A, X, 1 );
chkxer();
cblas_info = 4; RowMajorStrg = FALSE;
cblas_dtpsv(CblasColMajor, CblasUpper, CblasNoTrans,
INVALID, 0, A, X, 1 );
chkxer();
cblas_info = 5; RowMajorStrg = FALSE;
cblas_dtpsv(CblasColMajor, CblasUpper, CblasNoTrans,
CblasNonUnit, INVALID, A, X, 1 );
chkxer();
cblas_info = 8; RowMajorStrg = FALSE;
cblas_dtpsv(CblasColMajor, CblasUpper, CblasNoTrans,
CblasNonUnit, 0, A, X, 0 );
chkxer();
cblas_info = 2; RowMajorStrg = TRUE;
cblas_dtpsv(CblasRowMajor, INVALID, CblasNoTrans,
CblasNonUnit, 0, A, X, 1 );
chkxer();
cblas_info = 3; RowMajorStrg = TRUE;
cblas_dtpsv(CblasRowMajor, CblasUpper, INVALID,
CblasNonUnit, 0, A, X, 1 );
chkxer();
cblas_info = 4; RowMajorStrg = TRUE;
cblas_dtpsv(CblasRowMajor, CblasUpper, CblasNoTrans,
INVALID, 0, A, X, 1 );
chkxer();
cblas_info = 5; RowMajorStrg = TRUE;
cblas_dtpsv(CblasRowMajor, CblasUpper, CblasNoTrans,
CblasNonUnit, INVALID, A, X, 1 );
chkxer();
cblas_info = 8; RowMajorStrg = TRUE;
cblas_dtpsv(CblasRowMajor, CblasUpper, CblasNoTrans,
CblasNonUnit, 0, A, X, 0 );
chkxer();
} else if (strncmp( sf,"cblas_dger",10)==0) {
cblas_rout = "cblas_dger";
cblas_info = 1; RowMajorStrg = FALSE;
cblas_dger(INVALID, 0, 0, ALPHA, X, 1, Y, 1, A, 1 );
chkxer();
cblas_info = 2; RowMajorStrg = FALSE;
cblas_dger(CblasColMajor, INVALID, 0, ALPHA, X, 1, Y, 1, A, 1 );
chkxer();
cblas_info = 3; RowMajorStrg = FALSE;
cblas_dger(CblasColMajor, 0, INVALID, ALPHA, X, 1, Y, 1, A, 1 );
chkxer();
cblas_info = 6; RowMajorStrg = FALSE;
cblas_dger(CblasColMajor, 0, 0, ALPHA, X, 0, Y, 1, A, 1 );
chkxer();
cblas_info = 8; RowMajorStrg = FALSE;
cblas_dger(CblasColMajor, 0, 0, ALPHA, X, 1, Y, 0, A, 1 );
chkxer();
cblas_info = 10; RowMajorStrg = FALSE;
cblas_dger(CblasColMajor, 2, 0, ALPHA, X, 1, Y, 1, A, 1 );
chkxer();
cblas_info = 2; RowMajorStrg = TRUE;
cblas_dger(CblasRowMajor, INVALID, 0, ALPHA, X, 1, Y, 1, A, 1 );
chkxer();
cblas_info = 3; RowMajorStrg = TRUE;
cblas_dger(CblasRowMajor, 0, INVALID, ALPHA, X, 1, Y, 1, A, 1 );
chkxer();
cblas_info = 6; RowMajorStrg = TRUE;
cblas_dger(CblasRowMajor, 0, 0, ALPHA, X, 0, Y, 1, A, 1 );
chkxer();
cblas_info = 8; RowMajorStrg = TRUE;
cblas_dger(CblasRowMajor, 0, 0, ALPHA, X, 1, Y, 0, A, 1 );
chkxer();
cblas_info = 10; RowMajorStrg = TRUE;
cblas_dger(CblasRowMajor, 0, 2, ALPHA, X, 1, Y, 1, A, 1 );
chkxer();
} else if (strncmp( sf,"cblas_dsyr2",11)==0) {
cblas_rout = "cblas_dsyr2";
cblas_info = 1; RowMajorStrg = FALSE;
cblas_dsyr2(INVALID, CblasUpper, 0, ALPHA, X, 1, Y, 1, A, 1 );
chkxer();
cblas_info = 2; RowMajorStrg = FALSE;
cblas_dsyr2(CblasColMajor, INVALID, 0, ALPHA, X, 1, Y, 1, A, 1 );
chkxer();
cblas_info = 3; RowMajorStrg = FALSE;
cblas_dsyr2(CblasColMajor, CblasUpper, INVALID, ALPHA, X, 1, Y, 1, A, 1 );
chkxer();
cblas_info = 6; RowMajorStrg = FALSE;
cblas_dsyr2(CblasColMajor, CblasUpper, 0, ALPHA, X, 0, Y, 1, A, 1 );
chkxer();
cblas_info = 8; RowMajorStrg = FALSE;
cblas_dsyr2(CblasColMajor, CblasUpper, 0, ALPHA, X, 1, Y, 0, A, 1 );
chkxer();
cblas_info = 10; RowMajorStrg = FALSE;
cblas_dsyr2(CblasColMajor, CblasUpper, 2, ALPHA, X, 1, Y, 1, A, 1 );
chkxer();
cblas_info = 2; RowMajorStrg = TRUE;
cblas_dsyr2(CblasRowMajor, INVALID, 0, ALPHA, X, 1, Y, 1, A, 1 );
chkxer();
cblas_info = 3; RowMajorStrg = TRUE;
cblas_dsyr2(CblasRowMajor, CblasUpper, INVALID, ALPHA, X, 1, Y, 1, A, 1 );
chkxer();
cblas_info = 6; RowMajorStrg = TRUE;
cblas_dsyr2(CblasRowMajor, CblasUpper, 0, ALPHA, X, 0, Y, 1, A, 1 );
chkxer();
cblas_info = 8; RowMajorStrg = TRUE;
cblas_dsyr2(CblasRowMajor, CblasUpper, 0, ALPHA, X, 1, Y, 0, A, 1 );
chkxer();
cblas_info = 10; RowMajorStrg = TRUE;
cblas_dsyr2(CblasRowMajor, CblasUpper, 2, ALPHA, X, 1, Y, 1, A, 1 );
chkxer();
} else if (strncmp( sf,"cblas_dspr2",11)==0) {
cblas_rout = "cblas_dspr2";
cblas_info = 1; RowMajorStrg = FALSE;
cblas_dspr2(INVALID, CblasUpper, 0, ALPHA, X, 1, Y, 1, A );
chkxer();
cblas_info = 2; RowMajorStrg = FALSE;
cblas_dspr2(CblasColMajor, INVALID, 0, ALPHA, X, 1, Y, 1, A );
chkxer();
cblas_info = 3; RowMajorStrg = FALSE;
cblas_dspr2(CblasColMajor, CblasUpper, INVALID, ALPHA, X, 1, Y, 1, A );
chkxer();
cblas_info = 6; RowMajorStrg = FALSE;
cblas_dspr2(CblasColMajor, CblasUpper, 0, ALPHA, X, 0, Y, 1, A );
chkxer();
cblas_info = 8; RowMajorStrg = FALSE;
cblas_dspr2(CblasColMajor, CblasUpper, 0, ALPHA, X, 1, Y, 0, A );
chkxer();
cblas_info = 2; RowMajorStrg = TRUE;
cblas_dspr2(CblasRowMajor, INVALID, 0, ALPHA, X, 1, Y, 1, A );
chkxer();
cblas_info = 3; RowMajorStrg = TRUE;
cblas_dspr2(CblasRowMajor, CblasUpper, INVALID, ALPHA, X, 1, Y, 1, A );
chkxer();
cblas_info = 6; RowMajorStrg = TRUE;
cblas_dspr2(CblasRowMajor, CblasUpper, 0, ALPHA, X, 0, Y, 1, A );
chkxer();
cblas_info = 8; RowMajorStrg = TRUE;
cblas_dspr2(CblasRowMajor, CblasUpper, 0, ALPHA, X, 1, Y, 0, A );
chkxer();
} else if (strncmp( sf,"cblas_dsyr",10)==0) {
cblas_rout = "cblas_dsyr";
cblas_info = 1; RowMajorStrg = FALSE;
cblas_dsyr(INVALID, CblasUpper, 0, ALPHA, X, 1, A, 1 );
chkxer();
cblas_info = 2; RowMajorStrg = FALSE;
cblas_dsyr(CblasColMajor, INVALID, 0, ALPHA, X, 1, A, 1 );
chkxer();
cblas_info = 3; RowMajorStrg = FALSE;
cblas_dsyr(CblasColMajor, CblasUpper, INVALID, ALPHA, X, 1, A, 1 );
chkxer();
cblas_info = 6; RowMajorStrg = FALSE;
cblas_dsyr(CblasColMajor, CblasUpper, 0, ALPHA, X, 0, A, 1 );
chkxer();
cblas_info = 8; RowMajorStrg = FALSE;
cblas_dsyr(CblasColMajor, CblasUpper, 2, ALPHA, X, 1, A, 1 );
chkxer();
cblas_info = 2; RowMajorStrg = TRUE;
cblas_dsyr(CblasRowMajor, INVALID, 0, ALPHA, X, 1, A, 1 );
chkxer();
cblas_info = 3; RowMajorStrg = TRUE;
cblas_dsyr(CblasRowMajor, CblasUpper, INVALID, ALPHA, X, 1, A, 1 );
chkxer();
cblas_info = 6; RowMajorStrg = TRUE;
cblas_dsyr(CblasRowMajor, CblasUpper, 0, ALPHA, X, 0, A, 1 );
chkxer();
cblas_info = 8; RowMajorStrg = TRUE;
cblas_dsyr(CblasRowMajor, CblasUpper, 2, ALPHA, X, 1, A, 1 );
chkxer();
} else if (strncmp( sf,"cblas_dspr",10)==0) {
cblas_rout = "cblas_dspr";
cblas_info = 1; RowMajorStrg = FALSE;
cblas_dspr(INVALID, CblasUpper, 0, ALPHA, X, 1, A );
chkxer();
cblas_info = 2; RowMajorStrg = FALSE;
cblas_dspr(CblasColMajor, INVALID, 0, ALPHA, X, 1, A );
chkxer();
cblas_info = 3; RowMajorStrg = FALSE;
cblas_dspr(CblasColMajor, CblasUpper, INVALID, ALPHA, X, 1, A );
chkxer();
cblas_info = 6; RowMajorStrg = FALSE;
cblas_dspr(CblasColMajor, CblasUpper, 0, ALPHA, X, 0, A );
chkxer();
cblas_info = 2; RowMajorStrg = FALSE;
cblas_dspr(CblasColMajor, INVALID, 0, ALPHA, X, 1, A );
chkxer();
cblas_info = 3; RowMajorStrg = FALSE;
cblas_dspr(CblasColMajor, CblasUpper, INVALID, ALPHA, X, 1, A );
chkxer();
cblas_info = 6; RowMajorStrg = FALSE;
cblas_dspr(CblasColMajor, CblasUpper, 0, ALPHA, X, 0, A );
chkxer();
}
if (cblas_ok == TRUE)
printf(" %-12s PASSED THE TESTS OF ERROR-EXITS\n", cblas_rout);
else
printf("******* %s FAILED THE TESTS OF ERROR-EXITS *******\n",cblas_rout);
}
void My_dtrsv(const enum CBLAS_UPLO Uplo, const enum CBLAS_TRANSPOSE TransA, const enum CBLAS_DIAG Diag,
const gsl_matrix * A, gsl_vector* b)
{
int N = b->size;
cblas_dtrsv(CblasRowMajor, Uplo , TransA ,Diag, N, A->data, A->tda, b->data, b->stride);
}
void
test_trsv (void) {
const double flteps = 1e-4, dbleps = 1e-6;
{
int order = 101;
int trans = 111;
int uplo = 121;
int diag = 131;
int N = 1;
int lda = 1;
float A[] = { 0.995f };
float X[] = { 0.348f };
int incX = -1;
float x_expected[] = { 0.349749f };
cblas_strsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], flteps, "strsv(case 1150)");
}
};
};
{
int order = 101;
int trans = 111;
int uplo = 121;
int diag = 132;
int N = 1;
int lda = 1;
float A[] = { 0.995f };
float X[] = { 0.348f };
int incX = -1;
float x_expected[] = { 0.348f };
cblas_strsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], flteps, "strsv(case 1151)");
}
};
};
{
int order = 101;
int trans = 111;
int uplo = 122;
int diag = 131;
int N = 1;
int lda = 1;
float A[] = { 0.995f };
float X[] = { 0.348f };
int incX = -1;
float x_expected[] = { 0.349749f };
cblas_strsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], flteps, "strsv(case 1152)");
}
};
};
{
int order = 101;
int trans = 111;
int uplo = 122;
int diag = 132;
int N = 1;
int lda = 1;
float A[] = { 0.995f };
float X[] = { 0.348f };
int incX = -1;
float x_expected[] = { 0.348f };
cblas_strsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], flteps, "strsv(case 1153)");
}
};
};
{
int order = 102;
int trans = 111;
int uplo = 121;
int diag = 131;
int N = 1;
int lda = 1;
float A[] = { 0.995f };
float X[] = { 0.348f };
int incX = -1;
float x_expected[] = { 0.349749f };
cblas_strsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], flteps, "strsv(case 1154)");
}
};
};
{
int order = 102;
int trans = 111;
int uplo = 121;
int diag = 132;
int N = 1;
int lda = 1;
float A[] = { 0.995f };
float X[] = { 0.348f };
int incX = -1;
float x_expected[] = { 0.348f };
cblas_strsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], flteps, "strsv(case 1155)");
}
};
};
{
int order = 102;
int trans = 111;
int uplo = 122;
int diag = 131;
int N = 1;
int lda = 1;
float A[] = { 0.995f };
float X[] = { 0.348f };
int incX = -1;
float x_expected[] = { 0.349749f };
cblas_strsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], flteps, "strsv(case 1156)");
}
};
};
{
int order = 102;
int trans = 111;
int uplo = 122;
int diag = 132;
int N = 1;
int lda = 1;
float A[] = { 0.995f };
float X[] = { 0.348f };
int incX = -1;
float x_expected[] = { 0.348f };
cblas_strsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], flteps, "strsv(case 1157)");
}
};
};
{
int order = 101;
int trans = 112;
int uplo = 121;
int diag = 131;
int N = 1;
int lda = 1;
float A[] = { 0.793f };
float X[] = { 0.338f };
int incX = -1;
float x_expected[] = { 0.42623f };
cblas_strsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], flteps, "strsv(case 1158)");
}
};
};
{
int order = 101;
int trans = 112;
int uplo = 121;
int diag = 132;
int N = 1;
int lda = 1;
float A[] = { 0.793f };
float X[] = { 0.338f };
int incX = -1;
float x_expected[] = { 0.338f };
cblas_strsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], flteps, "strsv(case 1159)");
}
};
};
{
int order = 101;
int trans = 112;
int uplo = 122;
int diag = 131;
int N = 1;
int lda = 1;
float A[] = { 0.793f };
float X[] = { 0.338f };
int incX = -1;
float x_expected[] = { 0.42623f };
cblas_strsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], flteps, "strsv(case 1160)");
}
};
};
{
int order = 101;
int trans = 112;
int uplo = 122;
int diag = 132;
int N = 1;
int lda = 1;
float A[] = { 0.793f };
float X[] = { 0.338f };
int incX = -1;
float x_expected[] = { 0.338f };
cblas_strsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], flteps, "strsv(case 1161)");
}
};
};
{
int order = 102;
int trans = 112;
int uplo = 121;
int diag = 131;
int N = 1;
int lda = 1;
float A[] = { 0.793f };
float X[] = { 0.338f };
int incX = -1;
float x_expected[] = { 0.42623f };
cblas_strsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], flteps, "strsv(case 1162)");
}
};
};
{
int order = 102;
int trans = 112;
int uplo = 121;
int diag = 132;
int N = 1;
int lda = 1;
float A[] = { 0.793f };
float X[] = { 0.338f };
int incX = -1;
float x_expected[] = { 0.338f };
cblas_strsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], flteps, "strsv(case 1163)");
}
};
};
{
int order = 102;
int trans = 112;
int uplo = 122;
int diag = 131;
int N = 1;
int lda = 1;
float A[] = { 0.793f };
float X[] = { 0.338f };
int incX = -1;
float x_expected[] = { 0.42623f };
cblas_strsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], flteps, "strsv(case 1164)");
}
};
};
{
int order = 102;
int trans = 112;
int uplo = 122;
int diag = 132;
int N = 1;
int lda = 1;
float A[] = { 0.793f };
float X[] = { 0.338f };
int incX = -1;
float x_expected[] = { 0.338f };
cblas_strsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], flteps, "strsv(case 1165)");
}
};
};
{
int order = 101;
int trans = 111;
int uplo = 121;
int diag = 131;
int N = 1;
int lda = 1;
double A[] = { -0.21 };
double X[] = { 0.473 };
int incX = -1;
double x_expected[] = { -2.25238095238 };
cblas_dtrsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], dbleps, "dtrsv(case 1166)");
}
};
};
{
int order = 101;
int trans = 111;
int uplo = 121;
int diag = 132;
int N = 1;
int lda = 1;
double A[] = { -0.21 };
double X[] = { 0.473 };
int incX = -1;
double x_expected[] = { 0.473 };
cblas_dtrsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], dbleps, "dtrsv(case 1167)");
}
};
};
{
int order = 101;
int trans = 111;
int uplo = 122;
int diag = 131;
int N = 1;
int lda = 1;
double A[] = { -0.21 };
double X[] = { 0.473 };
int incX = -1;
double x_expected[] = { -2.25238095238 };
cblas_dtrsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], dbleps, "dtrsv(case 1168)");
}
};
};
{
int order = 101;
int trans = 111;
int uplo = 122;
int diag = 132;
int N = 1;
int lda = 1;
double A[] = { -0.21 };
double X[] = { 0.473 };
int incX = -1;
double x_expected[] = { 0.473 };
cblas_dtrsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], dbleps, "dtrsv(case 1169)");
}
};
};
{
int order = 102;
int trans = 111;
int uplo = 121;
int diag = 131;
int N = 1;
int lda = 1;
double A[] = { -0.21 };
double X[] = { 0.473 };
int incX = -1;
double x_expected[] = { -2.25238095238 };
cblas_dtrsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], dbleps, "dtrsv(case 1170)");
}
};
};
{
int order = 102;
int trans = 111;
int uplo = 121;
int diag = 132;
int N = 1;
int lda = 1;
double A[] = { -0.21 };
double X[] = { 0.473 };
int incX = -1;
double x_expected[] = { 0.473 };
cblas_dtrsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], dbleps, "dtrsv(case 1171)");
}
};
};
{
int order = 102;
int trans = 111;
int uplo = 122;
int diag = 131;
int N = 1;
int lda = 1;
double A[] = { -0.21 };
double X[] = { 0.473 };
int incX = -1;
double x_expected[] = { -2.25238095238 };
cblas_dtrsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], dbleps, "dtrsv(case 1172)");
}
};
};
{
int order = 102;
int trans = 111;
int uplo = 122;
int diag = 132;
int N = 1;
int lda = 1;
double A[] = { -0.21 };
double X[] = { 0.473 };
int incX = -1;
double x_expected[] = { 0.473 };
cblas_dtrsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], dbleps, "dtrsv(case 1173)");
}
};
};
{
int order = 101;
int trans = 112;
int uplo = 121;
int diag = 131;
int N = 1;
int lda = 1;
double A[] = { 0.748 };
double X[] = { 0.979 };
int incX = -1;
double x_expected[] = { 1.30882352941 };
cblas_dtrsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], dbleps, "dtrsv(case 1174)");
}
};
};
{
int order = 101;
int trans = 112;
int uplo = 121;
int diag = 132;
int N = 1;
int lda = 1;
double A[] = { 0.748 };
double X[] = { 0.979 };
int incX = -1;
double x_expected[] = { 0.979 };
cblas_dtrsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], dbleps, "dtrsv(case 1175)");
}
};
};
{
int order = 101;
int trans = 112;
int uplo = 122;
int diag = 131;
int N = 1;
int lda = 1;
double A[] = { 0.748 };
double X[] = { 0.979 };
int incX = -1;
double x_expected[] = { 1.30882352941 };
cblas_dtrsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], dbleps, "dtrsv(case 1176)");
}
};
};
{
int order = 101;
int trans = 112;
int uplo = 122;
int diag = 132;
int N = 1;
int lda = 1;
double A[] = { 0.748 };
double X[] = { 0.979 };
int incX = -1;
double x_expected[] = { 0.979 };
cblas_dtrsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], dbleps, "dtrsv(case 1177)");
}
};
};
{
int order = 102;
int trans = 112;
int uplo = 121;
int diag = 131;
int N = 1;
int lda = 1;
double A[] = { 0.748 };
double X[] = { 0.979 };
int incX = -1;
double x_expected[] = { 1.30882352941 };
cblas_dtrsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], dbleps, "dtrsv(case 1178)");
}
};
};
{
int order = 102;
int trans = 112;
int uplo = 121;
int diag = 132;
int N = 1;
int lda = 1;
double A[] = { 0.748 };
double X[] = { 0.979 };
int incX = -1;
double x_expected[] = { 0.979 };
cblas_dtrsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], dbleps, "dtrsv(case 1179)");
}
};
};
{
int order = 102;
int trans = 112;
int uplo = 122;
int diag = 131;
int N = 1;
int lda = 1;
double A[] = { 0.748 };
double X[] = { 0.979 };
int incX = -1;
double x_expected[] = { 1.30882352941 };
cblas_dtrsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], dbleps, "dtrsv(case 1180)");
}
};
};
{
int order = 102;
int trans = 112;
int uplo = 122;
int diag = 132;
int N = 1;
int lda = 1;
double A[] = { 0.748 };
double X[] = { 0.979 };
int incX = -1;
double x_expected[] = { 0.979 };
cblas_dtrsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], dbleps, "dtrsv(case 1181)");
}
};
};
{
int order = 101;
int trans = 111;
int uplo = 121;
int diag = 131;
int N = 1;
int lda = 1;
float A[] = { 0.529f, -0.348f };
float X[] = { -0.95f, 0.343f };
int incX = -1;
float x_expected[] = { -1.55112f, -0.372004f };
cblas_ctrsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctrsv(case 1182) real");
gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctrsv(case 1182) imag");
};
};
};
{
int order = 101;
int trans = 111;
int uplo = 121;
int diag = 132;
int N = 1;
int lda = 1;
float A[] = { 0.529f, -0.348f };
float X[] = { -0.95f, 0.343f };
int incX = -1;
float x_expected[] = { -0.95f, 0.343f };
cblas_ctrsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctrsv(case 1183) real");
gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctrsv(case 1183) imag");
};
};
};
{
int order = 101;
int trans = 111;
int uplo = 122;
int diag = 131;
int N = 1;
int lda = 1;
float A[] = { 0.529f, -0.348f };
float X[] = { -0.95f, 0.343f };
int incX = -1;
float x_expected[] = { -1.55112f, -0.372004f };
cblas_ctrsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctrsv(case 1184) real");
gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctrsv(case 1184) imag");
};
};
};
{
int order = 101;
int trans = 111;
int uplo = 122;
int diag = 132;
int N = 1;
int lda = 1;
float A[] = { 0.529f, -0.348f };
float X[] = { -0.95f, 0.343f };
int incX = -1;
float x_expected[] = { -0.95f, 0.343f };
cblas_ctrsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctrsv(case 1185) real");
gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctrsv(case 1185) imag");
};
};
};
{
int order = 102;
int trans = 111;
int uplo = 121;
int diag = 131;
int N = 1;
int lda = 1;
float A[] = { 0.529f, -0.348f };
float X[] = { -0.95f, 0.343f };
int incX = -1;
float x_expected[] = { -1.55112f, -0.372004f };
cblas_ctrsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctrsv(case 1186) real");
gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctrsv(case 1186) imag");
};
};
};
{
int order = 102;
int trans = 111;
int uplo = 121;
int diag = 132;
int N = 1;
int lda = 1;
float A[] = { 0.529f, -0.348f };
float X[] = { -0.95f, 0.343f };
int incX = -1;
float x_expected[] = { -0.95f, 0.343f };
cblas_ctrsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctrsv(case 1187) real");
gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctrsv(case 1187) imag");
};
};
};
{
int order = 102;
int trans = 111;
int uplo = 122;
int diag = 131;
int N = 1;
int lda = 1;
float A[] = { 0.529f, -0.348f };
float X[] = { -0.95f, 0.343f };
int incX = -1;
float x_expected[] = { -1.55112f, -0.372004f };
cblas_ctrsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctrsv(case 1188) real");
gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctrsv(case 1188) imag");
};
};
};
{
int order = 102;
int trans = 111;
int uplo = 122;
int diag = 132;
int N = 1;
int lda = 1;
float A[] = { 0.529f, -0.348f };
float X[] = { -0.95f, 0.343f };
int incX = -1;
float x_expected[] = { -0.95f, 0.343f };
cblas_ctrsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctrsv(case 1189) real");
gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctrsv(case 1189) imag");
};
};
};
{
int order = 101;
int trans = 112;
int uplo = 121;
int diag = 131;
int N = 1;
int lda = 1;
float A[] = { 0.6f, 0.041f };
float X[] = { 0.896f, -0.447f };
int incX = -1;
float x_expected[] = { 1.43572f, -0.843108f };
cblas_ctrsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctrsv(case 1190) real");
gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctrsv(case 1190) imag");
};
};
};
{
int order = 101;
int trans = 112;
int uplo = 121;
int diag = 132;
int N = 1;
int lda = 1;
float A[] = { 0.6f, 0.041f };
float X[] = { 0.896f, -0.447f };
int incX = -1;
float x_expected[] = { 0.896f, -0.447f };
cblas_ctrsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctrsv(case 1191) real");
gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctrsv(case 1191) imag");
};
};
};
{
int order = 101;
int trans = 112;
int uplo = 122;
int diag = 131;
int N = 1;
int lda = 1;
float A[] = { 0.6f, 0.041f };
float X[] = { 0.896f, -0.447f };
int incX = -1;
float x_expected[] = { 1.43572f, -0.843108f };
cblas_ctrsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctrsv(case 1192) real");
gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctrsv(case 1192) imag");
};
};
};
{
int order = 101;
int trans = 112;
int uplo = 122;
int diag = 132;
int N = 1;
int lda = 1;
float A[] = { 0.6f, 0.041f };
float X[] = { 0.896f, -0.447f };
int incX = -1;
float x_expected[] = { 0.896f, -0.447f };
cblas_ctrsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctrsv(case 1193) real");
gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctrsv(case 1193) imag");
};
};
};
{
int order = 102;
int trans = 112;
int uplo = 121;
int diag = 131;
int N = 1;
int lda = 1;
float A[] = { 0.6f, 0.041f };
float X[] = { 0.896f, -0.447f };
int incX = -1;
float x_expected[] = { 1.43572f, -0.843108f };
cblas_ctrsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctrsv(case 1194) real");
gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctrsv(case 1194) imag");
};
};
};
{
int order = 102;
int trans = 112;
int uplo = 121;
int diag = 132;
int N = 1;
int lda = 1;
float A[] = { 0.6f, 0.041f };
float X[] = { 0.896f, -0.447f };
int incX = -1;
float x_expected[] = { 0.896f, -0.447f };
cblas_ctrsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctrsv(case 1195) real");
gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctrsv(case 1195) imag");
};
};
};
{
int order = 102;
int trans = 112;
int uplo = 122;
int diag = 131;
int N = 1;
int lda = 1;
float A[] = { 0.6f, 0.041f };
float X[] = { 0.896f, -0.447f };
int incX = -1;
float x_expected[] = { 1.43572f, -0.843108f };
cblas_ctrsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctrsv(case 1196) real");
gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctrsv(case 1196) imag");
};
};
};
{
int order = 102;
int trans = 112;
int uplo = 122;
int diag = 132;
int N = 1;
int lda = 1;
float A[] = { 0.6f, 0.041f };
float X[] = { 0.896f, -0.447f };
int incX = -1;
float x_expected[] = { 0.896f, -0.447f };
cblas_ctrsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctrsv(case 1197) real");
gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctrsv(case 1197) imag");
};
};
};
{
int order = 101;
int trans = 113;
int uplo = 121;
int diag = 131;
int N = 1;
int lda = 1;
float A[] = { 0.397f, 0.683f };
float X[] = { 0.765f, 0.18f };
int incX = -1;
float x_expected[] = { 0.289642f, 0.951701f };
cblas_ctrsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctrsv(case 1198) real");
gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctrsv(case 1198) imag");
};
};
};
{
int order = 101;
int trans = 113;
int uplo = 121;
int diag = 132;
int N = 1;
int lda = 1;
float A[] = { 0.397f, 0.683f };
float X[] = { 0.765f, 0.18f };
int incX = -1;
float x_expected[] = { 0.765f, 0.18f };
cblas_ctrsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctrsv(case 1199) real");
gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctrsv(case 1199) imag");
};
};
};
{
int order = 101;
int trans = 113;
int uplo = 122;
int diag = 131;
int N = 1;
int lda = 1;
float A[] = { 0.397f, 0.683f };
float X[] = { 0.765f, 0.18f };
int incX = -1;
float x_expected[] = { 0.289642f, 0.951701f };
cblas_ctrsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctrsv(case 1200) real");
gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctrsv(case 1200) imag");
};
};
};
{
int order = 101;
int trans = 113;
int uplo = 122;
int diag = 132;
int N = 1;
int lda = 1;
float A[] = { 0.397f, 0.683f };
float X[] = { 0.765f, 0.18f };
int incX = -1;
float x_expected[] = { 0.765f, 0.18f };
cblas_ctrsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctrsv(case 1201) real");
gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctrsv(case 1201) imag");
};
};
};
{
int order = 102;
int trans = 113;
int uplo = 121;
int diag = 131;
int N = 1;
int lda = 1;
float A[] = { 0.397f, 0.683f };
float X[] = { 0.765f, 0.18f };
int incX = -1;
float x_expected[] = { 0.289642f, 0.951701f };
cblas_ctrsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctrsv(case 1202) real");
gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctrsv(case 1202) imag");
};
};
};
{
int order = 102;
int trans = 113;
int uplo = 121;
int diag = 132;
int N = 1;
int lda = 1;
float A[] = { 0.397f, 0.683f };
float X[] = { 0.765f, 0.18f };
int incX = -1;
float x_expected[] = { 0.765f, 0.18f };
cblas_ctrsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctrsv(case 1203) real");
gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctrsv(case 1203) imag");
};
};
};
{
int order = 102;
int trans = 113;
int uplo = 122;
int diag = 131;
int N = 1;
int lda = 1;
float A[] = { 0.397f, 0.683f };
float X[] = { 0.765f, 0.18f };
int incX = -1;
float x_expected[] = { 0.289642f, 0.951701f };
cblas_ctrsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctrsv(case 1204) real");
gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctrsv(case 1204) imag");
};
};
};
{
int order = 102;
int trans = 113;
int uplo = 122;
int diag = 132;
int N = 1;
int lda = 1;
float A[] = { 0.397f, 0.683f };
float X[] = { 0.765f, 0.18f };
int incX = -1;
float x_expected[] = { 0.765f, 0.18f };
cblas_ctrsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctrsv(case 1205) real");
gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctrsv(case 1205) imag");
};
};
};
{
int order = 101;
int trans = 111;
int uplo = 121;
int diag = 131;
int N = 1;
int lda = 1;
double A[] = { 0.977, -0.955 };
double X[] = { -0.627, 0.281 };
int incX = -1;
double x_expected[] = { -0.471957414573, -0.173714770642 };
cblas_ztrsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztrsv(case 1206) real");
gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztrsv(case 1206) imag");
};
};
};
{
int order = 101;
int trans = 111;
int uplo = 121;
int diag = 132;
int N = 1;
int lda = 1;
double A[] = { 0.977, -0.955 };
double X[] = { -0.627, 0.281 };
int incX = -1;
double x_expected[] = { -0.627, 0.281 };
cblas_ztrsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztrsv(case 1207) real");
gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztrsv(case 1207) imag");
};
};
};
{
int order = 101;
int trans = 111;
int uplo = 122;
int diag = 131;
int N = 1;
int lda = 1;
double A[] = { 0.977, -0.955 };
double X[] = { -0.627, 0.281 };
int incX = -1;
double x_expected[] = { -0.471957414573, -0.173714770642 };
cblas_ztrsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztrsv(case 1208) real");
gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztrsv(case 1208) imag");
};
};
};
{
int order = 101;
int trans = 111;
int uplo = 122;
int diag = 132;
int N = 1;
int lda = 1;
double A[] = { 0.977, -0.955 };
double X[] = { -0.627, 0.281 };
int incX = -1;
double x_expected[] = { -0.627, 0.281 };
cblas_ztrsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztrsv(case 1209) real");
gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztrsv(case 1209) imag");
};
};
};
{
int order = 102;
int trans = 111;
int uplo = 121;
int diag = 131;
int N = 1;
int lda = 1;
double A[] = { 0.977, -0.955 };
double X[] = { -0.627, 0.281 };
int incX = -1;
double x_expected[] = { -0.471957414573, -0.173714770642 };
cblas_ztrsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztrsv(case 1210) real");
gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztrsv(case 1210) imag");
};
};
};
{
int order = 102;
int trans = 111;
int uplo = 121;
int diag = 132;
int N = 1;
int lda = 1;
double A[] = { 0.977, -0.955 };
double X[] = { -0.627, 0.281 };
int incX = -1;
double x_expected[] = { -0.627, 0.281 };
cblas_ztrsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztrsv(case 1211) real");
gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztrsv(case 1211) imag");
};
};
};
{
int order = 102;
int trans = 111;
int uplo = 122;
int diag = 131;
int N = 1;
int lda = 1;
double A[] = { 0.977, -0.955 };
double X[] = { -0.627, 0.281 };
int incX = -1;
double x_expected[] = { -0.471957414573, -0.173714770642 };
cblas_ztrsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztrsv(case 1212) real");
gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztrsv(case 1212) imag");
};
};
};
{
int order = 102;
int trans = 111;
int uplo = 122;
int diag = 132;
int N = 1;
int lda = 1;
double A[] = { 0.977, -0.955 };
double X[] = { -0.627, 0.281 };
int incX = -1;
double x_expected[] = { -0.627, 0.281 };
cblas_ztrsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztrsv(case 1213) real");
gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztrsv(case 1213) imag");
};
};
};
{
int order = 101;
int trans = 112;
int uplo = 121;
int diag = 131;
int N = 1;
int lda = 1;
double A[] = { 0.076, -0.16 };
double X[] = { 0.3, -0.874 };
int incX = -1;
double x_expected[] = { 5.18357980622, -0.587200407955 };
cblas_ztrsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztrsv(case 1214) real");
gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztrsv(case 1214) imag");
};
};
};
{
int order = 101;
int trans = 112;
int uplo = 121;
int diag = 132;
int N = 1;
int lda = 1;
double A[] = { 0.076, -0.16 };
double X[] = { 0.3, -0.874 };
int incX = -1;
double x_expected[] = { 0.3, -0.874 };
cblas_ztrsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztrsv(case 1215) real");
gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztrsv(case 1215) imag");
};
};
};
{
int order = 101;
int trans = 112;
int uplo = 122;
int diag = 131;
int N = 1;
int lda = 1;
double A[] = { 0.076, -0.16 };
double X[] = { 0.3, -0.874 };
int incX = -1;
double x_expected[] = { 5.18357980622, -0.587200407955 };
cblas_ztrsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztrsv(case 1216) real");
gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztrsv(case 1216) imag");
};
};
};
{
int order = 101;
int trans = 112;
int uplo = 122;
int diag = 132;
int N = 1;
int lda = 1;
double A[] = { 0.076, -0.16 };
double X[] = { 0.3, -0.874 };
int incX = -1;
double x_expected[] = { 0.3, -0.874 };
cblas_ztrsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztrsv(case 1217) real");
gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztrsv(case 1217) imag");
};
};
};
{
int order = 102;
int trans = 112;
int uplo = 121;
int diag = 131;
int N = 1;
int lda = 1;
double A[] = { 0.076, -0.16 };
double X[] = { 0.3, -0.874 };
int incX = -1;
double x_expected[] = { 5.18357980622, -0.587200407955 };
cblas_ztrsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztrsv(case 1218) real");
gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztrsv(case 1218) imag");
};
};
};
{
int order = 102;
int trans = 112;
int uplo = 121;
int diag = 132;
int N = 1;
int lda = 1;
double A[] = { 0.076, -0.16 };
double X[] = { 0.3, -0.874 };
int incX = -1;
double x_expected[] = { 0.3, -0.874 };
cblas_ztrsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztrsv(case 1219) real");
gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztrsv(case 1219) imag");
};
};
};
{
int order = 102;
int trans = 112;
int uplo = 122;
int diag = 131;
int N = 1;
int lda = 1;
double A[] = { 0.076, -0.16 };
double X[] = { 0.3, -0.874 };
int incX = -1;
double x_expected[] = { 5.18357980622, -0.587200407955 };
cblas_ztrsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztrsv(case 1220) real");
gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztrsv(case 1220) imag");
};
};
};
{
int order = 102;
int trans = 112;
int uplo = 122;
int diag = 132;
int N = 1;
int lda = 1;
double A[] = { 0.076, -0.16 };
double X[] = { 0.3, -0.874 };
int incX = -1;
double x_expected[] = { 0.3, -0.874 };
cblas_ztrsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztrsv(case 1221) real");
gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztrsv(case 1221) imag");
};
};
};
{
int order = 101;
int trans = 113;
int uplo = 121;
int diag = 131;
int N = 1;
int lda = 1;
double A[] = { 0.372, -0.745 };
double X[] = { -0.085, -0.303 };
int incX = -1;
double x_expected[] = { -0.371144591432, -0.0712292456544 };
cblas_ztrsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztrsv(case 1222) real");
gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztrsv(case 1222) imag");
};
};
};
{
int order = 101;
int trans = 113;
int uplo = 121;
int diag = 132;
int N = 1;
int lda = 1;
double A[] = { 0.372, -0.745 };
double X[] = { -0.085, -0.303 };
int incX = -1;
double x_expected[] = { -0.085, -0.303 };
cblas_ztrsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztrsv(case 1223) real");
gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztrsv(case 1223) imag");
};
};
};
{
int order = 101;
int trans = 113;
int uplo = 122;
int diag = 131;
int N = 1;
int lda = 1;
double A[] = { 0.372, -0.745 };
double X[] = { -0.085, -0.303 };
int incX = -1;
double x_expected[] = { -0.371144591432, -0.0712292456544 };
cblas_ztrsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztrsv(case 1224) real");
gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztrsv(case 1224) imag");
};
};
};
{
int order = 101;
int trans = 113;
int uplo = 122;
int diag = 132;
int N = 1;
int lda = 1;
double A[] = { 0.372, -0.745 };
double X[] = { -0.085, -0.303 };
int incX = -1;
double x_expected[] = { -0.085, -0.303 };
cblas_ztrsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztrsv(case 1225) real");
gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztrsv(case 1225) imag");
};
};
};
{
int order = 102;
int trans = 113;
int uplo = 121;
int diag = 131;
int N = 1;
int lda = 1;
double A[] = { 0.372, -0.745 };
double X[] = { -0.085, -0.303 };
int incX = -1;
double x_expected[] = { -0.371144591432, -0.0712292456544 };
cblas_ztrsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztrsv(case 1226) real");
gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztrsv(case 1226) imag");
};
};
};
{
int order = 102;
int trans = 113;
int uplo = 121;
int diag = 132;
int N = 1;
int lda = 1;
double A[] = { 0.372, -0.745 };
double X[] = { -0.085, -0.303 };
int incX = -1;
double x_expected[] = { -0.085, -0.303 };
cblas_ztrsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztrsv(case 1227) real");
gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztrsv(case 1227) imag");
};
};
};
{
int order = 102;
int trans = 113;
int uplo = 122;
int diag = 131;
int N = 1;
int lda = 1;
double A[] = { 0.372, -0.745 };
double X[] = { -0.085, -0.303 };
int incX = -1;
double x_expected[] = { -0.371144591432, -0.0712292456544 };
cblas_ztrsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztrsv(case 1228) real");
gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztrsv(case 1228) imag");
};
};
};
{
int order = 102;
int trans = 113;
int uplo = 122;
int diag = 132;
int N = 1;
int lda = 1;
double A[] = { 0.372, -0.745 };
double X[] = { -0.085, -0.303 };
int incX = -1;
double x_expected[] = { -0.085, -0.303 };
cblas_ztrsv(order, uplo, trans, diag, N, A, lda, X, incX);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztrsv(case 1229) real");
gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztrsv(case 1229) imag");
};
};
};
}
__cminpack_attr__
void __cminpack_func__(lmpar)(int n, real *r, int ldr,
const int *ipvt, const real *diag, const real *qtb, real delta,
real *par, real *x, real *sdiag, real *wa1,
real *wa2)
{
/* Initialized data */
#define p1 .1
#define p001 .001
/* System generated locals */
real d1, d2;
/* Local variables */
int j, l;
real fp;
real parc, parl;
int iter;
real temp, paru, dwarf;
int nsing;
real gnorm;
real dxnorm;
/* ********** */
/* subroutine lmpar */
/* given an m by n matrix a, an n by n nonsingular diagonal */
/* matrix d, an m-vector b, and a positive number delta, */
/* the problem is to determine a value for the parameter */
/* par such that if x solves the system */
/* a*x = b , sqrt(par)*d*x = 0 , */
/* in the least squares sense, and dxnorm is the euclidean */
/* norm of d*x, then either par is zero and */
/* (dxnorm-delta) .le. 0.1*delta , */
/* or par is positive and */
/* abs(dxnorm-delta) .le. 0.1*delta . */
/* this subroutine completes the solution of the problem */
/* if it is provided with the necessary information from the */
/* qr factorization, with column pivoting, of a. that is, if */
/* a*p = q*r, where p is a permutation matrix, q has orthogonal */
/* columns, and r is an upper triangular matrix with diagonal */
/* elements of nonincreasing magnitude, then lmpar expects */
/* the full upper triangle of r, the permutation matrix p, */
/* and the first n components of (q transpose)*b. on output */
/* lmpar also provides an upper triangular matrix s such that */
/* t t t */
/* p *(a *a + par*d*d)*p = s *s . */
/* s is employed within lmpar and may be of separate interest. */
/* only a few iterations are generally needed for convergence */
/* of the algorithm. if, however, the limit of 10 iterations */
/* is reached, then the output par will contain the best */
/* value obtained so far. */
/* the subroutine statement is */
/* subroutine lmpar(n,r,ldr,ipvt,diag,qtb,delta,par,x,sdiag, */
/* wa1,wa2) */
/* where */
/* n is a positive integer input variable set to the order of r. */
/* r is an n by n array. on input the full upper triangle */
/* must contain the full upper triangle of the matrix r. */
/* on output the full upper triangle is unaltered, and the */
/* strict lower triangle contains the strict upper triangle */
/* (transposed) of the upper triangular matrix s. */
/* ldr is a positive integer input variable not less than n */
/* which specifies the leading dimension of the array r. */
/* ipvt is an integer input array of length n which defines the */
/* permutation matrix p such that a*p = q*r. column j of p */
/* is column ipvt(j) of the identity matrix. */
/* diag is an input array of length n which must contain the */
/* diagonal elements of the matrix d. */
/* qtb is an input array of length n which must contain the first */
/* n elements of the vector (q transpose)*b. */
/* delta is a positive input variable which specifies an upper */
/* bound on the euclidean norm of d*x. */
/* par is a nonnegative variable. on input par contains an */
/* initial estimate of the levenberg-marquardt parameter. */
/* on output par contains the final estimate. */
/* x is an output array of length n which contains the least */
/* squares solution of the system a*x = b, sqrt(par)*d*x = 0, */
/* for the output par. */
/* sdiag is an output array of length n which contains the */
/* diagonal elements of the upper triangular matrix s. */
/* wa1 and wa2 are work arrays of length n. */
/* subprograms called */
/* minpack-supplied ... dpmpar,enorm,qrsolv */
/* fortran-supplied ... dabs,dmax1,dmin1,dsqrt */
/* argonne national laboratory. minpack project. march 1980. */
/* burton s. garbow, kenneth e. hillstrom, jorge j. more */
/* ********** */
/* dwarf is the smallest positive magnitude. */
dwarf = __cminpack_func__(dpmpar)(2);
/* compute and store in x the gauss-newton direction. if the */
/* jacobian is rank-deficient, obtain a least squares solution. */
nsing = n;
for (j = 0; j < n; ++j) {
wa1[j] = qtb[j];
if (r[j + j * ldr] == 0. && nsing == n) {
nsing = j;
}
if (nsing < n) {
wa1[j] = 0.;
}
}
# ifdef USE_CBLAS
cblas_dtrsv(CblasColMajor, CblasUpper, CblasNoTrans, CblasNonUnit, nsing, r, ldr, wa1, 1);
# else
if (nsing >= 1) {
int k;
for (k = 1; k <= nsing; ++k) {
j = nsing - k;
wa1[j] /= r[j + j * ldr];
temp = wa1[j];
if (j >= 1) {
int i;
for (i = 0; i < j; ++i) {
wa1[i] -= r[i + j * ldr] * temp;
}
}
}
}
# endif
for (j = 0; j < n; ++j) {
l = ipvt[j]-1;
x[l] = wa1[j];
}
/* initialize the iteration counter. */
/* evaluate the function at the origin, and test */
/* for acceptance of the gauss-newton direction. */
iter = 0;
for (j = 0; j < n; ++j) {
wa2[j] = diag[j] * x[j];
}
dxnorm = __cminpack_enorm__(n, wa2);
fp = dxnorm - delta;
if (fp <= p1 * delta) {
goto TERMINATE;
}
/* if the jacobian is not rank deficient, the newton */
/* step provides a lower bound, parl, for the zero of */
/* the function. otherwise set this bound to zero. */
parl = 0.;
if (nsing >= n) {
for (j = 0; j < n; ++j) {
l = ipvt[j]-1;
wa1[j] = diag[l] * (wa2[l] / dxnorm);
}
# ifdef USE_CBLAS
cblas_dtrsv(CblasColMajor, CblasUpper, CblasTrans, CblasNonUnit, n, r, ldr, wa1, 1);
# else
for (j = 0; j < n; ++j) {
real sum = 0.;
if (j >= 1) {
int i;
for (i = 0; i < j; ++i) {
sum += r[i + j * ldr] * wa1[i];
}
}
wa1[j] = (wa1[j] - sum) / r[j + j * ldr];
}
# endif
temp = __cminpack_enorm__(n, wa1);
parl = fp / delta / temp / temp;
}
/* calculate an upper bound, paru, for the zero of the function. */
for (j = 0; j < n; ++j) {
real sum;
# ifdef USE_CBLAS
sum = cblas_ddot(j+1, &r[j*ldr], 1, qtb, 1);
# else
sum = 0.;
int i;
for (i = 0; i <= j; ++i) {
sum += r[i + j * ldr] * qtb[i];
}
# endif
l = ipvt[j]-1;
wa1[j] = sum / diag[l];
}
gnorm = __cminpack_enorm__(n, wa1);
paru = gnorm / delta;
if (paru == 0.) {
paru = dwarf / min(delta,(real)p1) /* / p001 ??? */;
}
/* if the input par lies outside of the interval (parl,paru), */
/* set par to the closer endpoint. */
*par = max(*par,parl);
*par = min(*par,paru);
if (*par == 0.) {
*par = gnorm / dxnorm;
}
/* beginning of an iteration. */
for (;;) {
++iter;
/* evaluate the function at the current value of par. */
if (*par == 0.) {
/* Computing MAX */
d1 = dwarf, d2 = p001 * paru;
*par = max(d1,d2);
}
temp = sqrt(*par);
for (j = 0; j < n; ++j) {
wa1[j] = temp * diag[j];
}
__cminpack_func__(qrsolv)(n, r, ldr, ipvt, wa1, qtb, x, sdiag, wa2);
for (j = 0; j < n; ++j) {
wa2[j] = diag[j] * x[j];
}
dxnorm = __cminpack_enorm__(n, wa2);
temp = fp;
fp = dxnorm - delta;
/* if the function is small enough, accept the current value */
/* of par. also test for the exceptional cases where parl */
/* is zero or the number of iterations has reached 10. */
if (fabs(fp) <= p1 * delta || (parl == 0. && fp <= temp && temp < 0.) || iter == 10) {
goto TERMINATE;
}
/* compute the newton correction. */
# ifdef USE_CBLAS
for (j = 0; j < nsing; ++j) {
l = ipvt[j]-1;
wa1[j] = diag[l] * (wa2[l] / dxnorm);
}
for (j = nsing; j < n; ++j) {
wa1[j] = 0.;
}
/* exchange the diagonal of r with sdiag */
cblas_dswap(n, r, ldr+1, sdiag, 1);
/* solve lower(r).x = wa1, result id put in wa1 */
cblas_dtrsv(CblasColMajor, CblasLower, CblasNoTrans, CblasNonUnit, nsing, r, ldr, wa1, 1);
/* exchange the diagonal of r with sdiag */
cblas_dswap( n, r, ldr+1, sdiag, 1);
# else /* !USE_CBLAS */
for (j = 0; j < n; ++j) {
l = ipvt[j]-1;
wa1[j] = diag[l] * (wa2[l] / dxnorm);
}
for (j = 0; j < n; ++j) {
wa1[j] /= sdiag[j];
temp = wa1[j];
if (n > j+1) {
int i;
for (i = j+1; i < n; ++i) {
wa1[i] -= r[i + j * ldr] * temp;
}
}
}
# endif /* !USE_CBLAS */
temp = __cminpack_enorm__(n, wa1);
parc = fp / delta / temp / temp;
/* depending on the sign of the function, update parl or paru. */
if (fp > 0.) {
parl = max(parl,*par);
}
if (fp < 0.) {
paru = min(paru,*par);
}
/* compute an improved estimate for par. */
/* Computing MAX */
d1 = parl, d2 = *par + parc;
*par = max(d1,d2);
/* end of an iteration. */
}
TERMINATE:
/* termination. */
if (iter == 0) {
*par = 0.;
}
/* last card of subroutine lmpar. */
} /* lmpar_ */
void HPL_dtrsv
(
const enum HPL_ORDER ORDER,
const enum HPL_UPLO UPLO,
const enum HPL_TRANS TRANS,
const enum HPL_DIAG DIAG,
const int N,
const double * A,
const int LDA,
double * X,
const int INCX
)
{
/*
* Purpose
* =======
*
* HPL_dtrsv solves one of the systems of equations
*
* A * x = b, or A^T * x = b,
*
* where b and x are n-element vectors and A is an n by n non-unit, or
* unit, upper or lower triangular matrix.
*
* No test for singularity or near-singularity is included in this
* routine. Such tests must be performed before calling this routine.
*
* Arguments
* =========
*
* ORDER (local input) const enum HPL_ORDER
* On entry, ORDER specifies the storage format of the operands
* as follows:
* ORDER = HplRowMajor,
* ORDER = HplColumnMajor.
*
* UPLO (local input) const enum HPL_UPLO
* On entry, UPLO specifies whether the upper or lower
* triangular part of the array A is to be referenced. When
* UPLO==HplUpper, only the upper triangular part of A is to be
* referenced, otherwise only the lower triangular part of A is
* to be referenced.
*
* TRANS (local input) const enum HPL_TRANS
* On entry, TRANS specifies the equations to be solved as
* follows:
* TRANS==HplNoTrans A * x = b,
* TRANS==HplTrans A^T * x = b.
*
* DIAG (local input) const enum HPL_DIAG
* On entry, DIAG specifies whether A is unit triangular or
* not. When DIAG==HplUnit, A is assumed to be unit triangular,
* and otherwise, A is not assumed to be unit triangular.
*
* N (local input) const int
* On entry, N specifies the order of the matrix A. N must be at
* least zero.
*
* A (local input) const double *
* On entry, A points to an array of size equal to or greater
* than LDA * n. Before entry with UPLO==HplUpper, the leading
* n by n upper triangular part of the array A must contain the
* upper triangular matrix and the strictly lower triangular
* part of A is not referenced. When UPLO==HplLower on entry,
* the leading n by n lower triangular part of the array A must
* contain the lower triangular matrix and the strictly upper
* triangular part of A is not referenced.
*
* Note that when DIAG==HplUnit, the diagonal elements of A
* not referenced either, but are assumed to be unity.
*
* LDA (local input) const int
* On entry, LDA specifies the leading dimension of A as
* declared in the calling (sub) program. LDA must be at
* least MAX(1,n).
*
* X (local input/output) double *
* On entry, X is an incremented array of dimension at least
* ( 1 + ( n - 1 ) * abs( INCX ) ) that contains the vector x.
* Before entry, the incremented array X must contain the n
* element right-hand side vector b. On exit, X is overwritten
* with the solution vector x.
*
* INCX (local input) const int
* On entry, INCX specifies the increment for the elements of X.
* INCX must not be zero.
*
* ---------------------------------------------------------------------
*/
START_TRACE( DTRSV )
cblas_dtrsv( ORDER, UPLO, TRANS, DIAG, N, A, LDA, X, INCX );
END_TRACE
/*
* End of HPL_dtrsv
*/
}
void gmres(double *A, double *D, double *x, double *b, int N, int max_restart, int max_iter, double tol)
{
int i, j, k, l, m, N2;
double resid, *normb, *beta, *temp_nrm, temp, *r, *q, *Aq, *qA, *Dq, *w, *cs, *sn, *s, *y, *Q, *H, *res;
normb = (double *) malloc(1*sizeof(double));
beta = (double *) malloc(1*sizeof(double));
temp_nrm = (double *) malloc(1*sizeof(double));
Q = (double *) malloc(N*N*(max_iter+1)*sizeof(double));
H = (double *) malloc((N+1)*max_iter*sizeof(double));
r = (double *) malloc(N*N*sizeof(double));
q = (double *) malloc(N*N*sizeof(double));
Aq = (double *) malloc(N*N*sizeof(double));
qA = (double *) malloc(N*N*sizeof(double));
Dq = (double *) malloc(N*N*sizeof(double));
w = (double *) malloc(N*N*sizeof(double));
cs = (double *) malloc((max_iter+1)*sizeof(double));
sn = (double *) malloc((max_iter+1)*sizeof(double));
s = (double *) malloc((max_iter+1)*sizeof(double));
y = (double *) malloc((max_iter+1)*sizeof(double));
res = (double *) malloc(max_iter*sizeof(double));
N2 = N*N;
norm(b, normb, N2);
for (k=0; k<N2; k++) r[k] = b[k];
norm(r, beta, N2);
if ((resid = *beta / *normb) <= tol)
{
tol = resid;
max_iter = 0;
}
for (m=0; m<max_restart; m++)
{
for (i=0; i<N2; i++) Q[i] = r[i] / *beta;
for (i=0; i<max_iter; i++) s[i+1] = 0.0;
s[0] = *beta;
for (i = 0; i<max_iter; i++)
{
q_subQ(q, Q, N2, i);
matrix_matrix(A, q, Aq, N);
matrix_matrix(q, A, qA, N);
matrix_matrix(D, q, Dq, N);
for (k=0; k<N2; k++) w[k] = Aq[k] + qA[k] + Dq[k];
for (k=0; k<=i; k++)
{
q_subQ(q, Q, N2, k);
H[max_iter*k+i] = inner_product(q, w, N2);
w_shift(w, q, H[max_iter*k+i], N2);
}
/* for (k=0; k<=i; k++)
{
H[max_iter*k+i] = 0.0;
for (j=0; j<N2; j++) H[max_iter*k+i] += Q[N2*k+j]*w[j];
}
for (k=0; k<=i; k++)
{
for (j=0; j<N2; j++) w[j] = w[j] - H[max_iter*k+i]*Q[N2*k+j];
}
*/
norm(w, temp_nrm, N2);
H[max_iter*(i+1)+i] = *temp_nrm;
subQ_v(Q, w, N2, i+1, H[max_iter*(i+1)+i]);
for (k = 0; k < i; k++)
{
//ApplyPlaneRotation(H(k,i), H(k+1,i), cs(k), sn(k))
temp = cs[k]*H[max_iter*k+i] + sn[k]*H[max_iter*(k+1)+i];
H[max_iter*(k+1)+i] = -1.0*sn[k]*H[max_iter*k+i] + cs[k]*H[max_iter*(k+1)+i];
H[max_iter*k+i] = temp;
}
GeneratePlaneRotation(H[max_iter*i+i], H[max_iter*(i+1)+i], cs, sn, i);
//ApplyPlaneRotation(H(i,i), H(i+1,i), cs(i), sn(i))
H[max_iter*i+i] = cs[i]*H[max_iter*i+i] + sn[i]*H[max_iter*(i+1)+i];
H[max_iter*(i+1)+i] = 0.0;
//ApplyPlaneRotation(s(i), s(i+1), cs(i), sn(i));
temp = cs[i]*s[i];
s[i+1] = -1.0*sn[i]*s[i];
s[i] = temp;
resid = fabs(s[i+1] / *beta);
res[i] = resid;
if (resid < tol)
{
// backsolve(H, s, y, N, max_iter, i);
for (k=0; k<max_iter+1; k++) y[k] = s[k];
cblas_dtrsv(CblasRowMajor, CblasUpper, CblasNoTrans, CblasNonUnit, i, H, max_iter, y, 1);
for(j=0; j<N; j++)
{
for (l=0; l<N; l++)
{
for(k=0; k<=i; k++)
{
x[N*j+l] += Q[N2*k+N*j+l]*y[k];
}
}
}
break;
}
}//end inside for
if (resid < tol)
{
printf(" resid = %e \n", resid);
printf(" Converges at %d cycle %d step. \n", m, i+1);
break;
}
// Caution : i = i + 1.
i = i - 1;
backsolve(H, s, y, N, max_iter, i);
for(j=0; j<N; j++)
{
for (l=0; l<N; l++)
{
for(k=0; k<=i; k++)
{
x[N*j+l] += Q[N2*k+N*j+l]*y[k];
}
}
}
matrix_matrix(A, x, Aq, N);
matrix_matrix(x, A, qA, N);
matrix_matrix(D, x, Dq, N);
for (j=0; j<N2; j++) r[j] = b[j] - (Aq[j] + qA[j] + Dq[j]);
norm(r, beta, N2);
s[i+1] = *beta;
resid = s[i+1] / *normb;
if ( resid < tol)
{
printf(" resid = %e \n", resid);
printf(" Converges at %d cycle %d step. \n", m, i);
break;
}
}//end outside for
}
