LAPACK 3.11.0
LAPACK: Linear Algebra PACKage
Loading...
Searching...
No Matches

◆ LAPACKE_sporfsx_work()

lapack_int LAPACKE_sporfsx_work ( int  matrix_layout,
char  uplo,
char  equed,
lapack_int  n,
lapack_int  nrhs,
const float *  a,
lapack_int  lda,
const float *  af,
lapack_int  ldaf,
const float *  s,
const float *  b,
lapack_int  ldb,
float *  x,
lapack_int  ldx,
float *  rcond,
float *  berr,
lapack_int  n_err_bnds,
float *  err_bnds_norm,
float *  err_bnds_comp,
lapack_int  nparams,
float *  params,
float *  work,
lapack_int iwork 
)

Definition at line 35 of file lapacke_sporfsx_work.c.

44{
45 lapack_int info = 0;
46 if( matrix_layout == LAPACK_COL_MAJOR ) {
47 /* Call LAPACK function and adjust info */
48 LAPACK_sporfsx( &uplo, &equed, &n, &nrhs, a, &lda, af, &ldaf, s, b,
49 &ldb, x, &ldx, rcond, berr, &n_err_bnds, err_bnds_norm,
50 err_bnds_comp, &nparams, params, work, iwork, &info );
51 if( info < 0 ) {
52 info = info - 1;
53 }
54 } else if( matrix_layout == LAPACK_ROW_MAJOR ) {
55 lapack_int lda_t = MAX(1,n);
56 lapack_int ldaf_t = MAX(1,n);
57 lapack_int ldb_t = MAX(1,n);
58 lapack_int ldx_t = MAX(1,n);
59 float* a_t = NULL;
60 float* af_t = NULL;
61 float* b_t = NULL;
62 float* x_t = NULL;
63 float* err_bnds_norm_t = NULL;
64 float* err_bnds_comp_t = NULL;
65 /* Check leading dimension(s) */
66 if( lda < n ) {
67 info = -7;
68 LAPACKE_xerbla( "LAPACKE_sporfsx_work", info );
69 return info;
70 }
71 if( ldaf < n ) {
72 info = -9;
73 LAPACKE_xerbla( "LAPACKE_sporfsx_work", info );
74 return info;
75 }
76 if( ldb < nrhs ) {
77 info = -12;
78 LAPACKE_xerbla( "LAPACKE_sporfsx_work", info );
79 return info;
80 }
81 if( ldx < nrhs ) {
82 info = -14;
83 LAPACKE_xerbla( "LAPACKE_sporfsx_work", info );
84 return info;
85 }
86 /* Allocate memory for temporary array(s) */
87 a_t = (float*)LAPACKE_malloc( sizeof(float) * lda_t * MAX(1,n) );
88 if( a_t == NULL ) {
90 goto exit_level_0;
91 }
92 af_t = (float*)LAPACKE_malloc( sizeof(float) * ldaf_t * MAX(1,n) );
93 if( af_t == NULL ) {
95 goto exit_level_1;
96 }
97 b_t = (float*)LAPACKE_malloc( sizeof(float) * ldb_t * MAX(1,nrhs) );
98 if( b_t == NULL ) {
100 goto exit_level_2;
101 }
102 x_t = (float*)LAPACKE_malloc( sizeof(float) * ldx_t * MAX(1,nrhs) );
103 if( x_t == NULL ) {
105 goto exit_level_3;
106 }
107 err_bnds_norm_t = (float*)
108 LAPACKE_malloc( sizeof(float) * nrhs * MAX(1,n_err_bnds) );
109 if( err_bnds_norm_t == NULL ) {
111 goto exit_level_4;
112 }
113 err_bnds_comp_t = (float*)
114 LAPACKE_malloc( sizeof(float) * nrhs * MAX(1,n_err_bnds) );
115 if( err_bnds_comp_t == NULL ) {
117 goto exit_level_5;
118 }
119 /* Transpose input matrices */
120 LAPACKE_ssy_trans( matrix_layout, uplo, n, a, lda, a_t, lda_t );
121 LAPACKE_ssy_trans( matrix_layout, uplo, n, af, ldaf, af_t, ldaf_t );
122 LAPACKE_sge_trans( matrix_layout, n, nrhs, b, ldb, b_t, ldb_t );
123 LAPACKE_sge_trans( matrix_layout, n, nrhs, x, ldx, x_t, ldx_t );
124 /* Call LAPACK function and adjust info */
125 LAPACK_sporfsx( &uplo, &equed, &n, &nrhs, a_t, &lda_t, af_t, &ldaf_t, s,
126 b_t, &ldb_t, x_t, &ldx_t, rcond, berr, &n_err_bnds,
127 err_bnds_norm_t, err_bnds_comp_t, &nparams, params,
128 work, iwork, &info );
129 if( info < 0 ) {
130 info = info - 1;
131 }
132 /* Transpose output matrices */
133 LAPACKE_sge_trans( LAPACK_COL_MAJOR, n, nrhs, x_t, ldx_t, x, ldx );
134 LAPACKE_sge_trans( LAPACK_COL_MAJOR, nrhs, n_err_bnds, err_bnds_norm_t,
135 nrhs, err_bnds_norm, nrhs );
136 LAPACKE_sge_trans( LAPACK_COL_MAJOR, nrhs, n_err_bnds, err_bnds_comp_t,
137 nrhs, err_bnds_comp, nrhs );
138 /* Release memory and exit */
139 LAPACKE_free( err_bnds_comp_t );
140exit_level_5:
141 LAPACKE_free( err_bnds_norm_t );
142exit_level_4:
143 LAPACKE_free( x_t );
144exit_level_3:
145 LAPACKE_free( b_t );
146exit_level_2:
147 LAPACKE_free( af_t );
148exit_level_1:
149 LAPACKE_free( a_t );
150exit_level_0:
151 if( info == LAPACK_TRANSPOSE_MEMORY_ERROR ) {
152 LAPACKE_xerbla( "LAPACKE_sporfsx_work", info );
153 }
154 } else {
155 info = -1;
156 LAPACKE_xerbla( "LAPACKE_sporfsx_work", info );
157 }
158 return info;
159}
#define LAPACK_sporfsx(...)
Definition: lapack.h:12952
#define lapack_int
Definition: lapack.h:87
#define LAPACK_COL_MAJOR
Definition: lapacke.h:53
#define LAPACKE_free(p)
Definition: lapacke.h:46
#define LAPACK_ROW_MAJOR
Definition: lapacke.h:52
#define LAPACKE_malloc(size)
Definition: lapacke.h:43
#define LAPACK_TRANSPOSE_MEMORY_ERROR
Definition: lapacke.h:56
void LAPACKE_xerbla(const char *name, lapack_int info)
void LAPACKE_ssy_trans(int matrix_layout, char uplo, lapack_int n, const float *in, lapack_int ldin, float *out, lapack_int ldout)
void LAPACKE_sge_trans(int matrix_layout, lapack_int m, lapack_int n, const float *in, lapack_int ldin, float *out, lapack_int ldout)
#define MAX(x, y)
Definition: lapacke_utils.h:46
Here is the call graph for this function:
Here is the caller graph for this function: