ssprfs_1.c
Go to the documentation of this file.
00001 /*****************************************************************************
00002   Copyright (c) 2010, Intel Corp.
00003   All rights reserved.
00004 
00005   Redistribution and use in source and binary forms, with or without
00006   modification, are permitted provided that the following conditions are met:
00007 
00008     * Redistributions of source code must retain the above copyright notice,
00009       this list of conditions and the following disclaimer.
00010     * Redistributions in binary form must reproduce the above copyright
00011       notice, this list of conditions and the following disclaimer in the
00012       documentation and/or other materials provided with the distribution.
00013     * Neither the name of Intel Corporation nor the names of its contributors
00014       may be used to endorse or promote products derived from this software
00015       without specific prior written permission.
00016 
00017   THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
00018   AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
00019   IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
00020   ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE
00021   LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
00022   CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
00023   SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
00024   INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
00025   CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
00026   ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF
00027   THE POSSIBILITY OF SUCH DAMAGE.
00028 *****************************************************************************/
00029 /*  Contents: test routine for C interface to LAPACK
00030 *   Author: Intel Corporation
00031 *   Created in March, 2010
00032 *
00033 * Purpose
00034 *
00035 * ssprfs_1 is the test program for the C interface to LAPACK
00036 * routine ssprfs
00037 * The program doesn't require an input, the input data is hardcoded in the
00038 * test program.
00039 * The program tests the C interface in the four combinations:
00040 *   1) column-major layout, middle-level interface
00041 *   2) column-major layout, high-level interface
00042 *   3) row-major layout, middle-level interface
00043 *   4) row-major layout, high-level interface
00044 * The output of the C interface function is compared to those obtained from
00045 * the corresponiding LAPACK routine with the same input data, and the
00046 * comparison diagnostics is then printed on the standard output having PASSED
00047 * keyword if the test is passed, and FAILED keyword if the test isn't passed.
00048 *****************************************************************************/
00049 #include <stdio.h>
00050 #include "lapacke.h"
00051 #include "lapacke_utils.h"
00052 #include "test_utils.h"
00053 
00054 static void init_scalars_ssprfs( char *uplo, lapack_int *n, lapack_int *nrhs,
00055                                  lapack_int *ldb, lapack_int *ldx );
00056 static void init_ap( lapack_int size, float *ap );
00057 static void init_afp( lapack_int size, float *afp );
00058 static void init_ipiv( lapack_int size, lapack_int *ipiv );
00059 static void init_b( lapack_int size, float *b );
00060 static void init_x( lapack_int size, float *x );
00061 static void init_ferr( lapack_int size, float *ferr );
00062 static void init_berr( lapack_int size, float *berr );
00063 static void init_work( lapack_int size, float *work );
00064 static void init_iwork( lapack_int size, lapack_int *iwork );
00065 static int compare_ssprfs( float *x, float *x_i, float *ferr, float *ferr_i,
00066                            float *berr, float *berr_i, lapack_int info,
00067                            lapack_int info_i, lapack_int ldx, lapack_int nrhs );
00068 
00069 int main(void)
00070 {
00071     /* Local scalars */
00072     char uplo, uplo_i;
00073     lapack_int n, n_i;
00074     lapack_int nrhs, nrhs_i;
00075     lapack_int ldb, ldb_i;
00076     lapack_int ldb_r;
00077     lapack_int ldx, ldx_i;
00078     lapack_int ldx_r;
00079     lapack_int info, info_i;
00080     lapack_int i;
00081     int failed;
00082 
00083     /* Local arrays */
00084     float *ap = NULL, *ap_i = NULL;
00085     float *afp = NULL, *afp_i = NULL;
00086     lapack_int *ipiv = NULL, *ipiv_i = NULL;
00087     float *b = NULL, *b_i = NULL;
00088     float *x = NULL, *x_i = NULL;
00089     float *ferr = NULL, *ferr_i = NULL;
00090     float *berr = NULL, *berr_i = NULL;
00091     float *work = NULL, *work_i = NULL;
00092     lapack_int *iwork = NULL, *iwork_i = NULL;
00093     float *x_save = NULL;
00094     float *ferr_save = NULL;
00095     float *berr_save = NULL;
00096     float *ap_r = NULL;
00097     float *afp_r = NULL;
00098     float *b_r = NULL;
00099     float *x_r = NULL;
00100 
00101     /* Iniitialize the scalar parameters */
00102     init_scalars_ssprfs( &uplo, &n, &nrhs, &ldb, &ldx );
00103     ldb_r = nrhs+2;
00104     ldx_r = nrhs+2;
00105     uplo_i = uplo;
00106     n_i = n;
00107     nrhs_i = nrhs;
00108     ldb_i = ldb;
00109     ldx_i = ldx;
00110 
00111     /* Allocate memory for the LAPACK routine arrays */
00112     ap = (float *)LAPACKE_malloc( ((n*(n+1)/2)) * sizeof(float) );
00113     afp = (float *)LAPACKE_malloc( ((n*(n+1)/2)) * sizeof(float) );
00114     ipiv = (lapack_int *)LAPACKE_malloc( n * sizeof(lapack_int) );
00115     b = (float *)LAPACKE_malloc( ldb*nrhs * sizeof(float) );
00116     x = (float *)LAPACKE_malloc( ldx*nrhs * sizeof(float) );
00117     ferr = (float *)LAPACKE_malloc( nrhs * sizeof(float) );
00118     berr = (float *)LAPACKE_malloc( nrhs * sizeof(float) );
00119     work = (float *)LAPACKE_malloc( 3*n * sizeof(float) );
00120     iwork = (lapack_int *)LAPACKE_malloc( n * sizeof(lapack_int) );
00121 
00122     /* Allocate memory for the C interface function arrays */
00123     ap_i = (float *)LAPACKE_malloc( ((n*(n+1)/2)) * sizeof(float) );
00124     afp_i = (float *)LAPACKE_malloc( ((n*(n+1)/2)) * sizeof(float) );
00125     ipiv_i = (lapack_int *)LAPACKE_malloc( n * sizeof(lapack_int) );
00126     b_i = (float *)LAPACKE_malloc( ldb*nrhs * sizeof(float) );
00127     x_i = (float *)LAPACKE_malloc( ldx*nrhs * sizeof(float) );
00128     ferr_i = (float *)LAPACKE_malloc( nrhs * sizeof(float) );
00129     berr_i = (float *)LAPACKE_malloc( nrhs * sizeof(float) );
00130     work_i = (float *)LAPACKE_malloc( 3*n * sizeof(float) );
00131     iwork_i = (lapack_int *)LAPACKE_malloc( n * sizeof(lapack_int) );
00132 
00133     /* Allocate memory for the backup arrays */
00134     x_save = (float *)LAPACKE_malloc( ldx*nrhs * sizeof(float) );
00135     ferr_save = (float *)LAPACKE_malloc( nrhs * sizeof(float) );
00136     berr_save = (float *)LAPACKE_malloc( nrhs * sizeof(float) );
00137 
00138     /* Allocate memory for the row-major arrays */
00139     ap_r = (float *)LAPACKE_malloc( n*(n+1)/2 * sizeof(float) );
00140     afp_r = (float *)LAPACKE_malloc( n*(n+1)/2 * sizeof(float) );
00141     b_r = (float *)LAPACKE_malloc( n*(nrhs+2) * sizeof(float) );
00142     x_r = (float *)LAPACKE_malloc( n*(nrhs+2) * sizeof(float) );
00143 
00144     /* Initialize input arrays */
00145     init_ap( (n*(n+1)/2), ap );
00146     init_afp( (n*(n+1)/2), afp );
00147     init_ipiv( n, ipiv );
00148     init_b( ldb*nrhs, b );
00149     init_x( ldx*nrhs, x );
00150     init_ferr( nrhs, ferr );
00151     init_berr( nrhs, berr );
00152     init_work( 3*n, work );
00153     init_iwork( n, iwork );
00154 
00155     /* Backup the ouptut arrays */
00156     for( i = 0; i < ldx*nrhs; i++ ) {
00157         x_save[i] = x[i];
00158     }
00159     for( i = 0; i < nrhs; i++ ) {
00160         ferr_save[i] = ferr[i];
00161     }
00162     for( i = 0; i < nrhs; i++ ) {
00163         berr_save[i] = berr[i];
00164     }
00165 
00166     /* Call the LAPACK routine */
00167     ssprfs_( &uplo, &n, &nrhs, ap, afp, ipiv, b, &ldb, x, &ldx, ferr, berr,
00168              work, iwork, &info );
00169 
00170     /* Initialize input data, call the column-major middle-level
00171      * interface to LAPACK routine and check the results */
00172     for( i = 0; i < (n*(n+1)/2); i++ ) {
00173         ap_i[i] = ap[i];
00174     }
00175     for( i = 0; i < (n*(n+1)/2); i++ ) {
00176         afp_i[i] = afp[i];
00177     }
00178     for( i = 0; i < n; i++ ) {
00179         ipiv_i[i] = ipiv[i];
00180     }
00181     for( i = 0; i < ldb*nrhs; i++ ) {
00182         b_i[i] = b[i];
00183     }
00184     for( i = 0; i < ldx*nrhs; i++ ) {
00185         x_i[i] = x_save[i];
00186     }
00187     for( i = 0; i < nrhs; i++ ) {
00188         ferr_i[i] = ferr_save[i];
00189     }
00190     for( i = 0; i < nrhs; i++ ) {
00191         berr_i[i] = berr_save[i];
00192     }
00193     for( i = 0; i < 3*n; i++ ) {
00194         work_i[i] = work[i];
00195     }
00196     for( i = 0; i < n; i++ ) {
00197         iwork_i[i] = iwork[i];
00198     }
00199     info_i = LAPACKE_ssprfs_work( LAPACK_COL_MAJOR, uplo_i, n_i, nrhs_i, ap_i,
00200                                   afp_i, ipiv_i, b_i, ldb_i, x_i, ldx_i, ferr_i,
00201                                   berr_i, work_i, iwork_i );
00202 
00203     failed = compare_ssprfs( x, x_i, ferr, ferr_i, berr, berr_i, info, info_i,
00204                              ldx, nrhs );
00205     if( failed == 0 ) {
00206         printf( "PASSED: column-major middle-level interface to ssprfs\n" );
00207     } else {
00208         printf( "FAILED: column-major middle-level interface to ssprfs\n" );
00209     }
00210 
00211     /* Initialize input data, call the column-major high-level
00212      * interface to LAPACK routine and check the results */
00213     for( i = 0; i < (n*(n+1)/2); i++ ) {
00214         ap_i[i] = ap[i];
00215     }
00216     for( i = 0; i < (n*(n+1)/2); i++ ) {
00217         afp_i[i] = afp[i];
00218     }
00219     for( i = 0; i < n; i++ ) {
00220         ipiv_i[i] = ipiv[i];
00221     }
00222     for( i = 0; i < ldb*nrhs; i++ ) {
00223         b_i[i] = b[i];
00224     }
00225     for( i = 0; i < ldx*nrhs; i++ ) {
00226         x_i[i] = x_save[i];
00227     }
00228     for( i = 0; i < nrhs; i++ ) {
00229         ferr_i[i] = ferr_save[i];
00230     }
00231     for( i = 0; i < nrhs; i++ ) {
00232         berr_i[i] = berr_save[i];
00233     }
00234     for( i = 0; i < 3*n; i++ ) {
00235         work_i[i] = work[i];
00236     }
00237     for( i = 0; i < n; i++ ) {
00238         iwork_i[i] = iwork[i];
00239     }
00240     info_i = LAPACKE_ssprfs( LAPACK_COL_MAJOR, uplo_i, n_i, nrhs_i, ap_i, afp_i,
00241                              ipiv_i, b_i, ldb_i, x_i, ldx_i, ferr_i, berr_i );
00242 
00243     failed = compare_ssprfs( x, x_i, ferr, ferr_i, berr, berr_i, info, info_i,
00244                              ldx, nrhs );
00245     if( failed == 0 ) {
00246         printf( "PASSED: column-major high-level interface to ssprfs\n" );
00247     } else {
00248         printf( "FAILED: column-major high-level interface to ssprfs\n" );
00249     }
00250 
00251     /* Initialize input data, call the row-major middle-level
00252      * interface to LAPACK routine and check the results */
00253     for( i = 0; i < (n*(n+1)/2); i++ ) {
00254         ap_i[i] = ap[i];
00255     }
00256     for( i = 0; i < (n*(n+1)/2); i++ ) {
00257         afp_i[i] = afp[i];
00258     }
00259     for( i = 0; i < n; i++ ) {
00260         ipiv_i[i] = ipiv[i];
00261     }
00262     for( i = 0; i < ldb*nrhs; i++ ) {
00263         b_i[i] = b[i];
00264     }
00265     for( i = 0; i < ldx*nrhs; i++ ) {
00266         x_i[i] = x_save[i];
00267     }
00268     for( i = 0; i < nrhs; i++ ) {
00269         ferr_i[i] = ferr_save[i];
00270     }
00271     for( i = 0; i < nrhs; i++ ) {
00272         berr_i[i] = berr_save[i];
00273     }
00274     for( i = 0; i < 3*n; i++ ) {
00275         work_i[i] = work[i];
00276     }
00277     for( i = 0; i < n; i++ ) {
00278         iwork_i[i] = iwork[i];
00279     }
00280 
00281     LAPACKE_spp_trans( LAPACK_COL_MAJOR, uplo, n, ap_i, ap_r );
00282     LAPACKE_spp_trans( LAPACK_COL_MAJOR, uplo, n, afp_i, afp_r );
00283     LAPACKE_sge_trans( LAPACK_COL_MAJOR, n, nrhs, b_i, ldb, b_r, nrhs+2 );
00284     LAPACKE_sge_trans( LAPACK_COL_MAJOR, n, nrhs, x_i, ldx, x_r, nrhs+2 );
00285     info_i = LAPACKE_ssprfs_work( LAPACK_ROW_MAJOR, uplo_i, n_i, nrhs_i, ap_r,
00286                                   afp_r, ipiv_i, b_r, ldb_r, x_r, ldx_r, ferr_i,
00287                                   berr_i, work_i, iwork_i );
00288 
00289     LAPACKE_sge_trans( LAPACK_ROW_MAJOR, n, nrhs, x_r, nrhs+2, x_i, ldx );
00290 
00291     failed = compare_ssprfs( x, x_i, ferr, ferr_i, berr, berr_i, info, info_i,
00292                              ldx, nrhs );
00293     if( failed == 0 ) {
00294         printf( "PASSED: row-major middle-level interface to ssprfs\n" );
00295     } else {
00296         printf( "FAILED: row-major middle-level interface to ssprfs\n" );
00297     }
00298 
00299     /* Initialize input data, call the row-major high-level
00300      * interface to LAPACK routine and check the results */
00301     for( i = 0; i < (n*(n+1)/2); i++ ) {
00302         ap_i[i] = ap[i];
00303     }
00304     for( i = 0; i < (n*(n+1)/2); i++ ) {
00305         afp_i[i] = afp[i];
00306     }
00307     for( i = 0; i < n; i++ ) {
00308         ipiv_i[i] = ipiv[i];
00309     }
00310     for( i = 0; i < ldb*nrhs; i++ ) {
00311         b_i[i] = b[i];
00312     }
00313     for( i = 0; i < ldx*nrhs; i++ ) {
00314         x_i[i] = x_save[i];
00315     }
00316     for( i = 0; i < nrhs; i++ ) {
00317         ferr_i[i] = ferr_save[i];
00318     }
00319     for( i = 0; i < nrhs; i++ ) {
00320         berr_i[i] = berr_save[i];
00321     }
00322     for( i = 0; i < 3*n; i++ ) {
00323         work_i[i] = work[i];
00324     }
00325     for( i = 0; i < n; i++ ) {
00326         iwork_i[i] = iwork[i];
00327     }
00328 
00329     /* Init row_major arrays */
00330     LAPACKE_spp_trans( LAPACK_COL_MAJOR, uplo, n, ap_i, ap_r );
00331     LAPACKE_spp_trans( LAPACK_COL_MAJOR, uplo, n, afp_i, afp_r );
00332     LAPACKE_sge_trans( LAPACK_COL_MAJOR, n, nrhs, b_i, ldb, b_r, nrhs+2 );
00333     LAPACKE_sge_trans( LAPACK_COL_MAJOR, n, nrhs, x_i, ldx, x_r, nrhs+2 );
00334     info_i = LAPACKE_ssprfs( LAPACK_ROW_MAJOR, uplo_i, n_i, nrhs_i, ap_r, afp_r,
00335                              ipiv_i, b_r, ldb_r, x_r, ldx_r, ferr_i, berr_i );
00336 
00337     LAPACKE_sge_trans( LAPACK_ROW_MAJOR, n, nrhs, x_r, nrhs+2, x_i, ldx );
00338 
00339     failed = compare_ssprfs( x, x_i, ferr, ferr_i, berr, berr_i, info, info_i,
00340                              ldx, nrhs );
00341     if( failed == 0 ) {
00342         printf( "PASSED: row-major high-level interface to ssprfs\n" );
00343     } else {
00344         printf( "FAILED: row-major high-level interface to ssprfs\n" );
00345     }
00346 
00347     /* Release memory */
00348     if( ap != NULL ) {
00349         LAPACKE_free( ap );
00350     }
00351     if( ap_i != NULL ) {
00352         LAPACKE_free( ap_i );
00353     }
00354     if( ap_r != NULL ) {
00355         LAPACKE_free( ap_r );
00356     }
00357     if( afp != NULL ) {
00358         LAPACKE_free( afp );
00359     }
00360     if( afp_i != NULL ) {
00361         LAPACKE_free( afp_i );
00362     }
00363     if( afp_r != NULL ) {
00364         LAPACKE_free( afp_r );
00365     }
00366     if( ipiv != NULL ) {
00367         LAPACKE_free( ipiv );
00368     }
00369     if( ipiv_i != NULL ) {
00370         LAPACKE_free( ipiv_i );
00371     }
00372     if( b != NULL ) {
00373         LAPACKE_free( b );
00374     }
00375     if( b_i != NULL ) {
00376         LAPACKE_free( b_i );
00377     }
00378     if( b_r != NULL ) {
00379         LAPACKE_free( b_r );
00380     }
00381     if( x != NULL ) {
00382         LAPACKE_free( x );
00383     }
00384     if( x_i != NULL ) {
00385         LAPACKE_free( x_i );
00386     }
00387     if( x_r != NULL ) {
00388         LAPACKE_free( x_r );
00389     }
00390     if( x_save != NULL ) {
00391         LAPACKE_free( x_save );
00392     }
00393     if( ferr != NULL ) {
00394         LAPACKE_free( ferr );
00395     }
00396     if( ferr_i != NULL ) {
00397         LAPACKE_free( ferr_i );
00398     }
00399     if( ferr_save != NULL ) {
00400         LAPACKE_free( ferr_save );
00401     }
00402     if( berr != NULL ) {
00403         LAPACKE_free( berr );
00404     }
00405     if( berr_i != NULL ) {
00406         LAPACKE_free( berr_i );
00407     }
00408     if( berr_save != NULL ) {
00409         LAPACKE_free( berr_save );
00410     }
00411     if( work != NULL ) {
00412         LAPACKE_free( work );
00413     }
00414     if( work_i != NULL ) {
00415         LAPACKE_free( work_i );
00416     }
00417     if( iwork != NULL ) {
00418         LAPACKE_free( iwork );
00419     }
00420     if( iwork_i != NULL ) {
00421         LAPACKE_free( iwork_i );
00422     }
00423 
00424     return 0;
00425 }
00426 
00427 /* Auxiliary function: ssprfs scalar parameters initialization */
00428 static void init_scalars_ssprfs( char *uplo, lapack_int *n, lapack_int *nrhs,
00429                                  lapack_int *ldb, lapack_int *ldx )
00430 {
00431     *uplo = 'L';
00432     *n = 4;
00433     *nrhs = 2;
00434     *ldb = 8;
00435     *ldx = 8;
00436 
00437     return;
00438 }
00439 
00440 /* Auxiliary functions: ssprfs array parameters initialization */
00441 static void init_ap( lapack_int size, float *ap ) {
00442     lapack_int i;
00443     for( i = 0; i < size; i++ ) {
00444         ap[i] = 0;
00445     }
00446     ap[0] = 2.069999933e+000;
00447     ap[1] = 3.869999886e+000;
00448     ap[2] = 4.199999809e+000;
00449     ap[3] = -1.149999976e+000;
00450     ap[4] = -2.099999934e-001;
00451     ap[5] = 1.870000005e+000;
00452     ap[6] = 6.299999952e-001;
00453     ap[7] = 1.149999976e+000;
00454     ap[8] = 2.059999943e+000;
00455     ap[9] = -1.809999943e+000;
00456 }
00457 static void init_afp( lapack_int size, float *afp ) {
00458     lapack_int i;
00459     for( i = 0; i < size; i++ ) {
00460         afp[i] = 0;
00461     }
00462     afp[0] = 2.069999933e+000;
00463     afp[1] = 4.199999809e+000;
00464     afp[2] = 2.230414152e-001;
00465     afp[3] = 6.536583900e-001;
00466     afp[4] = 1.149999976e+000;
00467     afp[5] = 8.115010262e-001;
00468     afp[6] = -5.959697366e-001;
00469     afp[7] = -2.590677261e+000;
00470     afp[8] = 3.030846119e-001;
00471     afp[9] = 4.073850811e-001;
00472 }
00473 static void init_ipiv( lapack_int size, lapack_int *ipiv ) {
00474     lapack_int i;
00475     for( i = 0; i < size; i++ ) {
00476         ipiv[i] = 0;
00477     }
00478     ipiv[0] = -3;
00479     ipiv[1] = -3;
00480     ipiv[2] = 3;
00481     ipiv[3] = 4;
00482 }
00483 static void init_b( lapack_int size, float *b ) {
00484     lapack_int i;
00485     for( i = 0; i < size; i++ ) {
00486         b[i] = 0;
00487     }
00488     b[0] = -9.500000000e+000;  /* b[0,0] */
00489     b[8] = 2.785000038e+001;  /* b[0,1] */
00490     b[1] = -8.380000114e+000;  /* b[1,0] */
00491     b[9] = 9.899999619e+000;  /* b[1,1] */
00492     b[2] = -6.070000172e+000;  /* b[2,0] */
00493     b[10] = 1.925000000e+001;  /* b[2,1] */
00494     b[3] = -9.599999785e-001;  /* b[3,0] */
00495     b[11] = 3.930000067e+000;  /* b[3,1] */
00496 }
00497 static void init_x( lapack_int size, float *x ) {
00498     lapack_int i;
00499     for( i = 0; i < size; i++ ) {
00500         x[i] = 0;
00501     }
00502     x[0] = -4.000000954e+000;  /* x[0,0] */
00503     x[8] = 1.000000715e+000;  /* x[0,1] */
00504     x[1] = -1.000000119e+000;  /* x[1,0] */
00505     x[9] = 4.000000954e+000;  /* x[1,1] */
00506     x[2] = 2.000000715e+000;  /* x[2,0] */
00507     x[10] = 2.999999046e+000;  /* x[2,1] */
00508     x[3] = 5.000000954e+000;  /* x[3,0] */
00509     x[11] = 1.999998569e+000;  /* x[3,1] */
00510 }
00511 static void init_ferr( lapack_int size, float *ferr ) {
00512     lapack_int i;
00513     for( i = 0; i < size; i++ ) {
00514         ferr[i] = 0;
00515     }
00516 }
00517 static void init_berr( lapack_int size, float *berr ) {
00518     lapack_int i;
00519     for( i = 0; i < size; i++ ) {
00520         berr[i] = 0;
00521     }
00522 }
00523 static void init_work( lapack_int size, float *work ) {
00524     lapack_int i;
00525     for( i = 0; i < size; i++ ) {
00526         work[i] = 0;
00527     }
00528 }
00529 static void init_iwork( lapack_int size, lapack_int *iwork ) {
00530     lapack_int i;
00531     for( i = 0; i < size; i++ ) {
00532         iwork[i] = 0;
00533     }
00534 }
00535 
00536 /* Auxiliary function: C interface to ssprfs results check */
00537 /* Return value: 0 - test is passed, non-zero - test is failed */
00538 static int compare_ssprfs( float *x, float *x_i, float *ferr, float *ferr_i,
00539                            float *berr, float *berr_i, lapack_int info,
00540                            lapack_int info_i, lapack_int ldx, lapack_int nrhs )
00541 {
00542     lapack_int i;
00543     int failed = 0;
00544     for( i = 0; i < ldx*nrhs; i++ ) {
00545         failed += compare_floats(x[i],x_i[i]);
00546     }
00547     for( i = 0; i < nrhs; i++ ) {
00548         failed += compare_floats(ferr[i],ferr_i[i]);
00549     }
00550     for( i = 0; i < nrhs; i++ ) {
00551         failed += compare_floats(berr[i],berr_i[i]);
00552     }
00553     failed += (info == info_i) ? 0 : 1;
00554     if( info != 0 || info_i != 0 ) {
00555         printf( "info=%d, info_i=%d\n",(int)info,(int)info_i );
00556     }
00557 
00558     return failed;
00559 }


swiftnav
Author(s):
autogenerated on Sat Jun 8 2019 18:56:13