lapack_int LAPACKE_zlarnv( lapack_int idist, lapack_int* iseed, lapack_int n,
lapack_complex_double* x );
+lapack_int LAPACKE_slascl( int matrix_layout, char type, lapack_int kl,
+ lapack_int ku, float cfrom, float cto,
+ lapack_int m, lapack_int n, float* a,
+ lapack_int lda );
+lapack_int LAPACKE_dlascl( int matrix_layout, char type, lapack_int kl,
+ lapack_int ku, double cfrom, double cto,
+ lapack_int m, lapack_int n, double* a,
+ lapack_int lda );
+lapack_int LAPACKE_clascl( int matrix_layout, char type, lapack_int kl,
+ lapack_int ku, float cfrom, float cto,
+ lapack_int m, lapack_int n, lapack_complex_float* a,
+ lapack_int lda );
+lapack_int LAPACKE_zlascl( int matrix_layout, char type, lapack_int kl,
+ lapack_int ku, double cfrom, double cto,
+ lapack_int m, lapack_int n, lapack_complex_double* a,
+ lapack_int lda );
+
lapack_int LAPACKE_slaset( int matrix_layout, char uplo, lapack_int m,
lapack_int n, float alpha, float beta, float* a,
lapack_int lda );
lapack_int LAPACKE_zlarnv_work( lapack_int idist, lapack_int* iseed,
lapack_int n, lapack_complex_double* x );
+
+lapack_int LAPACKE_slascl_work( int matrix_layout, char type, lapack_int kl,
+ lapack_int ku, float cfrom, float cto,
+ lapack_int m, lapack_int n, float* a,
+ lapack_int lda );
+lapack_int LAPACKE_dlascl_work( int matrix_layout, char type, lapack_int kl,
+ lapack_int ku, double cfrom, double cto,
+ lapack_int m, lapack_int n, double* a,
+ lapack_int lda );
+lapack_int LAPACKE_clascl_work( int matrix_layout, char type, lapack_int kl,
+ lapack_int ku, float cfrom, float cto,
+ lapack_int m, lapack_int n, lapack_complex_float* a,
+ lapack_int lda );
+lapack_int LAPACKE_zlascl_work( int matrix_layout, char type, lapack_int kl,
+ lapack_int ku, double cfrom, double cto,
+ lapack_int m, lapack_int n, lapack_complex_double* a,
+ lapack_int lda );
+
lapack_int LAPACKE_slaset_work( int matrix_layout, char uplo, lapack_int m,
lapack_int n, float alpha, float beta, float* a,
lapack_int lda );
#define LAPACK_dlagge LAPACK_GLOBAL(dlagge,DLAGGE)
#define LAPACK_clagge LAPACK_GLOBAL(clagge,CLAGGE)
#define LAPACK_zlagge LAPACK_GLOBAL(zlagge,ZLAGGE)
+#define LAPACK_slascl LAPACK_GLOBAL(slascl,SLASCL)
+#define LAPACK_dlascl LAPACK_GLOBAL(dlascl,DLASCL)
+#define LAPACK_clascl LAPACK_GLOBAL(clascl,CLASCL)
+#define LAPACK_zlascl LAPACK_GLOBAL(zlascl,ZLASCL)
#define LAPACK_slaset LAPACK_GLOBAL(slaset,SLASET)
#define LAPACK_dlaset LAPACK_GLOBAL(dlaset,DLASET)
#define LAPACK_claset LAPACK_GLOBAL(claset,CLASET)
lapack_int* ku, const double* d, lapack_complex_double* a,
lapack_int* lda, lapack_int* iseed,
lapack_complex_double* work, lapack_int *info );
+void LAPACK_slascl( char* type, lapack_int* kl, lapack_int* ku, float* cfrom,
+ float* cto, lapack_int* m, lapack_int* n, float* a,
+ lapack_int* lda, lapack_int *info );
+void LAPACK_dlascl( char* type, lapack_int* kl, lapack_int* ku, double* cfrom,
+ double* cto, lapack_int* m, lapack_int* n, double* a,
+ lapack_int* lda, lapack_int *info );
+void LAPACK_clascl( char* type, lapack_int* kl, lapack_int* ku, float* cfrom,
+ float* cto, lapack_int* m, lapack_int* n, lapack_complex_float* a,
+ lapack_int* lda, lapack_int *info );
+void LAPACK_zlascl( char* type, lapack_int* kl, lapack_int* ku, double* cfrom,
+ double* cto, lapack_int* m, lapack_int* n, lapack_complex_double* a,
+ lapack_int* lda, lapack_int *info );
void LAPACK_slaset( char* uplo, lapack_int* m, lapack_int* n, float* alpha,
float* beta, float* a, lapack_int* lda );
void LAPACK_dlaset( char* uplo, lapack_int* m, lapack_int* n, double* alpha,
lapacke_clarfx_work.c
lapacke_clarnv.c
lapacke_clarnv_work.c
+lapacke_clascl.c
+lapacke_clascl_work.c
lapacke_claset.c
lapacke_claset_work.c
lapacke_claswp.c
lapacke_dlartgp_work.c
lapacke_dlartgs.c
lapacke_dlartgs_work.c
+lapacke_dlascl.c
+lapacke_dlascl_work.c
lapacke_dlaset.c
lapacke_dlaset_work.c
lapacke_dlasrt.c
lapacke_slartgp_work.c
lapacke_slartgs.c
lapacke_slartgs_work.c
+lapacke_slascl.c
+lapacke_slascl_work.c
lapacke_slaset.c
lapacke_slaset_work.c
lapacke_slasrt.c
lapacke_zlarfx_work.c
lapacke_zlarnv.c
lapacke_zlarnv_work.c
+lapacke_zlascl.c
+lapacke_zlascl_work.c
lapacke_zlaset.c
lapacke_zlaset_work.c
lapacke_zlaswp.c
lapacke_clarfx_work.o \
lapacke_clarnv.o \
lapacke_clarnv_work.o \
+lapacke_clascl.o \
+lapacke_clascl_work.o \
lapacke_claset.o \
lapacke_claset_work.o \
lapacke_claswp.o \
lapacke_dlartgp_work.o \
lapacke_dlartgs.o \
lapacke_dlartgs_work.o \
+lapacke_dlascl.o \
+lapacke_dlascl_work.o \
lapacke_dlaset.o \
lapacke_dlaset_work.o \
lapacke_dlasrt.o \
lapacke_slartgp_work.o \
lapacke_slartgs.o \
lapacke_slartgs_work.o \
+lapacke_slascl.o \
+lapacke_slascl_work.o \
lapacke_slaset.o \
lapacke_slaset_work.o \
lapacke_slasrt.o \
lapacke_zlarfx_work.o \
lapacke_zlarnv.o \
lapacke_zlarnv_work.o \
+lapacke_zlascl.o \
+lapacke_zlascl_work.o \
lapacke_zlaset.o \
lapacke_zlaset_work.o \
lapacke_zlaswp.o \
--- /dev/null
+/*****************************************************************************
+ Copyright (c) 2014, Intel Corp.
+ All rights reserved.
+
+ Redistribution and use in source and binary forms, with or without
+ modification, are permitted provided that the following conditions are met:
+
+ * Redistributions of source code must retain the above copyright notice,
+ this list of conditions and the following disclaimer.
+ * Redistributions in binary form must reproduce the above copyright
+ notice, this list of conditions and the following disclaimer in the
+ documentation and/or other materials provided with the distribution.
+ * Neither the name of Intel Corporation nor the names of its contributors
+ may be used to endorse or promote products derived from this software
+ without specific prior written permission.
+
+ THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
+ AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
+ IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
+ ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE
+ LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
+ CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
+ SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
+ INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
+ CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
+ ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF
+ THE POSSIBILITY OF SUCH DAMAGE.
+*****************************************************************************
+* Contents: Native high-level C interface to LAPACK function slaswp
+* Author: Intel Corporation
+* Generated November, 2011
+*****************************************************************************/
+
+#include "lapacke_utils.h"
+
+lapack_int LAPACKE_clascl( int matrix_layout, char type, lapack_int kl,
+ lapack_int ku, float cfrom, float cto,
+ lapack_int m, lapack_int n, lapack_complex_float* a,
+ lapack_int lda )
+{
+ if( matrix_layout != LAPACK_COL_MAJOR && matrix_layout != LAPACK_ROW_MAJOR ) {
+ LAPACKE_xerbla( "LAPACKE_clascl", -1 );
+ return -1;
+ }
+#ifndef LAPACK_DISABLE_NAN_CHECK
+ /* Optionally check input matrices for NaNs */
+ switch (type) {
+ case 'G':
+ if( LAPACKE_cge_nancheck( matrix_layout, lda, n, a, lda ) ) {
+ return -9;
+ }
+ break;
+ case 'L':
+ // TYPE = 'L' - lower triangular matrix.
+ if( LAPACKE_ctr_nancheck( matrix_layout, 'L', 'N', n, a, lda ) ) {
+ return -9;
+ }
+ break;
+ case 'U':
+ // TYPE = 'U' - upper triangular matrix
+ if( LAPACKE_ctr_nancheck( matrix_layout, 'U', 'N', n, a, lda ) ) {
+ return -9;
+ }
+ break;
+ case 'H':
+ // TYPE = 'H' - upper Hessenberg matrix
+ if( LAPACKE_chs_nancheck( matrix_layout, n, a, lda ) ) {
+ return -9;
+ }
+ break;
+ case 'B':
+ // TYPE = 'B' - A is a symmetric band matrix with lower bandwidth KL
+ // and upper bandwidth KU and with the only the lower
+ // half stored.
+ if( LAPACKE_csb_nancheck( matrix_layout, 'L', n, kl, a, lda ) ) {
+ return -9;
+ }
+ break;
+ case 'Q':
+ // TYPE = 'Q' - A is a symmetric band matrix with lower bandwidth KL
+ // and upper bandwidth KU and with the only the upper
+ // half stored.
+ if( LAPACKE_csb_nancheck( matrix_layout, 'U', n, ku, a, lda ) ) {
+ return -9;
+ }
+ break;
+ case 'Z':
+ // TYPE = 'Z' - A is a band matrix with lower bandwidth KL and upper
+ // bandwidth KU. See DGBTRF for storage details.
+ if( LAPACKE_cgb_nancheck( matrix_layout, n, n, kl, kl+ku, a, lda ) ) {
+ return -6;
+ }
+ break;
+ }
+#endif
+ return LAPACKE_clascl_work( matrix_layout, type, kl, ku, cfrom, cto, m, n, a, lda );
+}
--- /dev/null
+/*****************************************************************************
+ Copyright (c) 2014, Intel Corp.
+ All rights reserved.
+
+ Redistribution and use in source and binary forms, with or without
+ modification, are permitted provided that the following conditions are met:
+
+ * Redistributions of source code must retain the above copyright notice,
+ this list of conditions and the following disclaimer.
+ * Redistributions in binary form must reproduce the above copyright
+ notice, this list of conditions and the following disclaimer in the
+ documentation and/or other materials provided with the distribution.
+ * Neither the name of Intel Corporation nor the names of its contributors
+ may be used to endorse or promote products derived from this software
+ without specific prior written permission.
+
+ THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
+ AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
+ IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
+ ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE
+ LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
+ CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
+ SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
+ INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
+ CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
+ ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF
+ THE POSSIBILITY OF SUCH DAMAGE.
+*****************************************************************************
+* Contents: Native middle-level C interface to LAPACK function slaswp
+* Author: Intel Corporation
+* Generated November, 2011
+*****************************************************************************/
+
+#include "lapacke_utils.h"
+
+lapack_int LAPACKE_clascl_work( int matrix_layout, char type, lapack_int kl,
+ lapack_int ku, float cfrom, float cto,
+ lapack_int m, lapack_int n, lapack_complex_float* a,
+ lapack_int lda )
+{
+ lapack_int info = 0;
+ if( matrix_layout == LAPACK_COL_MAJOR ) {
+ /* Call LAPACK function and adjust info */
+ LAPACK_clascl( &type, &kl, &ku, &cfrom, &cto, &m, &n, a, &lda, &info);
+ if( info < 0 ) {
+ info = info - 1;
+ }
+ } else if( matrix_layout == LAPACK_ROW_MAJOR ) {
+ lapack_int lda_t = MAX(1,lda);
+ lapack_complex_float* a_t = NULL;
+ /* Check leading dimension(s) */
+ if( lda < n ) {
+ info = -9;
+ LAPACKE_xerbla( "LAPACKE_clascl_work", info );
+ return info;
+ }
+ /* Allocate memory for temporary array(s) */
+ a_t = (lapack_complex_float*)
+ LAPACKE_malloc( sizeof(lapack_complex_float) * lda_t * MAX(1,n) );
+ if( a_t == NULL ) {
+ info = LAPACK_TRANSPOSE_MEMORY_ERROR;
+ goto exit_level_0;
+ }
+ /* Transpose input matrices */
+ LAPACKE_cge_trans( matrix_layout, lda, n, a, lda, a_t, lda_t );
+ /* Call LAPACK function and adjust info */
+ LAPACK_clascl( &type, &kl, &ku, &cfrom, &cto, &m, &n, a_t, &lda_t, &info);
+ info = 0; /* LAPACK call is ok! */
+ /* Transpose output matrices */
+ LAPACKE_cge_trans( LAPACK_COL_MAJOR, lda, n, a_t, lda_t, a, lda );
+ /* Release memory and exit */
+ LAPACKE_free( a_t );
+exit_level_0:
+ if( info == LAPACK_TRANSPOSE_MEMORY_ERROR ) {
+ LAPACKE_xerbla( "LAPACKE_clascl_work", info );
+ }
+ } else {
+ info = -1;
+ LAPACKE_xerbla( "LAPACKE_clascl_work", info );
+ }
+ return info;
+}
--- /dev/null
+/*****************************************************************************
+ Copyright (c) 2014, Intel Corp.
+ All rights reserved.
+
+ Redistribution and use in source and binary forms, with or without
+ modification, are permitted provided that the following conditions are met:
+
+ * Redistributions of source code must retain the above copyright notice,
+ this list of conditions and the following disclaimer.
+ * Redistributions in binary form must reproduce the above copyright
+ notice, this list of conditions and the following disclaimer in the
+ documentation and/or other materials provided with the distribution.
+ * Neither the name of Intel Corporation nor the names of its contributors
+ may be used to endorse or promote products derived from this software
+ without specific prior written permission.
+
+ THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
+ AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
+ IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
+ ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE
+ LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
+ CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
+ SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
+ INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
+ CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
+ ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF
+ THE POSSIBILITY OF SUCH DAMAGE.
+*****************************************************************************
+* Contents: Native high-level C interface to LAPACK function dlaswp
+* Author: Intel Corporation
+* Generated November, 2011
+*****************************************************************************/
+
+#include "lapacke_utils.h"
+
+lapack_int LAPACKE_dlascl( int matrix_layout, char type, lapack_int kl,
+ lapack_int ku, double cfrom, double cto,
+ lapack_int m, lapack_int n, double* a,
+ lapack_int lda )
+{
+ if( matrix_layout != LAPACK_COL_MAJOR && matrix_layout != LAPACK_ROW_MAJOR ) {
+ LAPACKE_xerbla( "LAPACKE_dlascl", -1 );
+ return -1;
+ }
+#ifndef LAPACK_DISABLE_NAN_CHECK
+ /* Optionally check input matrices for NaNs */
+ switch (type) {
+ case 'G':
+ if( LAPACKE_dge_nancheck( matrix_layout, lda, n, a, lda ) ) {
+ return -9;
+ }
+ break;
+ case 'L':
+ // TYPE = 'L' - lower triangular matrix.
+ if( LAPACKE_dtr_nancheck( matrix_layout, 'L', 'N', n, a, lda ) ) {
+ return -9;
+ }
+ break;
+ case 'U':
+ // TYPE = 'U' - upper triangular matrix
+ if( LAPACKE_dtr_nancheck( matrix_layout, 'U', 'N', n, a, lda ) ) {
+ return -9;
+ }
+ break;
+ case 'H':
+ // TYPE = 'H' - upper Hessenberg matrix
+ if( LAPACKE_dhs_nancheck( matrix_layout, n, a, lda ) ) {
+ return -9;
+ }
+ break;
+ case 'B':
+ // TYPE = 'B' - A is a symmetric band matrix with lower bandwidth KL
+ // and upper bandwidth KU and with the only the lower
+ // half stored.
+ if( LAPACKE_dsb_nancheck( matrix_layout, 'L', n, kl, a, lda ) ) {
+ return -9;
+ }
+ break;
+ case 'Q':
+ // TYPE = 'Q' - A is a symmetric band matrix with lower bandwidth KL
+ // and upper bandwidth KU and with the only the upper
+ // half stored.
+ if( LAPACKE_dsb_nancheck( matrix_layout, 'U', n, ku, a, lda ) ) {
+ return -9;
+ }
+ break;
+ case 'Z':
+ // TYPE = 'Z' - A is a band matrix with lower bandwidth KL and upper
+ // bandwidth KU. See DGBTRF for storage details.
+ if( LAPACKE_dgb_nancheck( matrix_layout, n, n, kl, kl+ku, a, lda ) ) {
+ return -6;
+ }
+ break;
+ }
+#endif
+ return LAPACKE_dlascl_work( matrix_layout, type, kl, ku, cfrom, cto, m, n, a, lda );
+}
--- /dev/null
+/*****************************************************************************
+ Copyright (c) 2014, Intel Corp.
+ All rights reserved.
+
+ Redistribution and use in source and binary forms, with or without
+ modification, are permitted provided that the following conditions are met:
+
+ * Redistributions of source code must retain the above copyright notice,
+ this list of conditions and the following disclaimer.
+ * Redistributions in binary form must reproduce the above copyright
+ notice, this list of conditions and the following disclaimer in the
+ documentation and/or other materials provided with the distribution.
+ * Neither the name of Intel Corporation nor the names of its contributors
+ may be used to endorse or promote products derived from this software
+ without specific prior written permission.
+
+ THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
+ AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
+ IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
+ ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE
+ LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
+ CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
+ SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
+ INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
+ CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
+ ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF
+ THE POSSIBILITY OF SUCH DAMAGE.
+*****************************************************************************
+* Contents: Native middle-level C interface to LAPACK function dlaswp
+* Author: Intel Corporation
+* Generated November, 2011
+*****************************************************************************/
+
+#include "lapacke_utils.h"
+
+lapack_int LAPACKE_dlascl_work( int matrix_layout, char type, lapack_int kl,
+ lapack_int ku, double cfrom, double cto,
+ lapack_int m, lapack_int n, double* a,
+ lapack_int lda )
+{
+ lapack_int info = 0;
+ if( matrix_layout == LAPACK_COL_MAJOR ) {
+ /* Call LAPACK function and adjust info */
+ LAPACK_dlascl( &type, &kl, &ku, &cfrom, &cto, &m, &n, a, &lda, &info);
+ if( info < 0 ) {
+ info = info - 1;
+ }
+ } else if( matrix_layout == LAPACK_ROW_MAJOR ) {
+ lapack_int lda_t = MAX(1,lda);
+ double* a_t = NULL;
+ /* Check leading dimension(s) */
+ if( lda < n ) {
+ info = -9;
+ LAPACKE_xerbla( "LAPACKE_dlascl_work", info );
+ return info;
+ }
+ /* Allocate memory for temporary array(s) */
+ a_t = (double*)LAPACKE_malloc( sizeof(double) * lda_t * MAX(1,n) );
+ if( a_t == NULL ) {
+ info = LAPACK_TRANSPOSE_MEMORY_ERROR;
+ goto exit_level_0;
+ }
+ /* Transpose input matrices */
+ LAPACKE_dge_trans( matrix_layout, lda, n, a, lda, a_t, lda_t );
+ /* Call LAPACK function and adjust info */
+ LAPACK_dlascl( &type, &kl, &ku, &cfrom, &cto, &m, &n, a_t, &lda_t, &info);
+ info = 0; /* LAPACK call is ok! */
+ /* Transpose output matrices */
+ LAPACKE_dge_trans( LAPACK_COL_MAJOR, lda, n, a_t, lda_t, a, lda );
+ /* Release memory and exit */
+ LAPACKE_free( a_t );
+exit_level_0:
+ if( info == LAPACK_TRANSPOSE_MEMORY_ERROR ) {
+ LAPACKE_xerbla( "LAPACKE_dlascl_work", info );
+ }
+ } else {
+ info = -1;
+ LAPACKE_xerbla( "LAPACKE_dlascl_work", info );
+ }
+ return info;
+}
--- /dev/null
+/*****************************************************************************
+ Copyright (c) 2014, Intel Corp.
+ All rights reserved.
+
+ Redistribution and use in source and binary forms, with or without
+ modification, are permitted provided that the following conditions are met:
+
+ * Redistributions of source code must retain the above copyright notice,
+ this list of conditions and the following disclaimer.
+ * Redistributions in binary form must reproduce the above copyright
+ notice, this list of conditions and the following disclaimer in the
+ documentation and/or other materials provided with the distribution.
+ * Neither the name of Intel Corporation nor the names of its contributors
+ may be used to endorse or promote products derived from this software
+ without specific prior written permission.
+
+ THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
+ AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
+ IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
+ ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE
+ LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
+ CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
+ SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
+ INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
+ CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
+ ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF
+ THE POSSIBILITY OF SUCH DAMAGE.
+*****************************************************************************
+* Contents: Native high-level C interface to LAPACK function slaswp
+* Author: Intel Corporation
+* Generated November, 2011
+*****************************************************************************/
+
+#include "lapacke_utils.h"
+
+lapack_int LAPACKE_slascl( int matrix_layout, char type, lapack_int kl,
+ lapack_int ku, float cfrom, float cto,
+ lapack_int m, lapack_int n, float* a,
+ lapack_int lda )
+{
+ if( matrix_layout != LAPACK_COL_MAJOR && matrix_layout != LAPACK_ROW_MAJOR ) {
+ LAPACKE_xerbla( "LAPACKE_slascl", -1 );
+ return -1;
+ }
+#ifndef LAPACK_DISABLE_NAN_CHECK
+ /* Optionally check input matrices for NaNs */
+ switch (type) {
+ case 'G':
+ if( LAPACKE_sge_nancheck( matrix_layout, lda, n, a, lda ) ) {
+ return -9;
+ }
+ break;
+ case 'L':
+ // TYPE = 'L' - lower triangular matrix.
+ if( LAPACKE_str_nancheck( matrix_layout, 'L', 'N', n, a, lda ) ) {
+ return -9;
+ }
+ break;
+ case 'U':
+ // TYPE = 'U' - upper triangular matrix
+ if( LAPACKE_str_nancheck( matrix_layout, 'U', 'N', n, a, lda ) ) {
+ return -9;
+ }
+ break;
+ case 'H':
+ // TYPE = 'H' - upper Hessenberg matrix
+ if( LAPACKE_shs_nancheck( matrix_layout, n, a, lda ) ) {
+ return -9;
+ }
+ break;
+ case 'B':
+ // TYPE = 'B' - A is a symmetric band matrix with lower bandwidth KL
+ // and upper bandwidth KU and with the only the lower
+ // half stored.
+ if( LAPACKE_ssb_nancheck( matrix_layout, 'L', n, kl, a, lda ) ) {
+ return -9;
+ }
+ break;
+ case 'Q':
+ // TYPE = 'Q' - A is a symmetric band matrix with lower bandwidth KL
+ // and upper bandwidth KU and with the only the upper
+ // half stored.
+ if( LAPACKE_ssb_nancheck( matrix_layout, 'U', n, ku, a, lda ) ) {
+ return -9;
+ }
+ break;
+ case 'Z':
+ // TYPE = 'Z' - A is a band matrix with lower bandwidth KL and upper
+ // bandwidth KU. See DGBTRF for storage details.
+ if( LAPACKE_sgb_nancheck( matrix_layout, n, n, kl, kl+ku, a, lda ) ) {
+ return -6;
+ }
+ break;
+ }
+#endif
+ return LAPACKE_slascl_work( matrix_layout, type, kl, ku, cfrom, cto, m, n, a, lda );
+}
--- /dev/null
+/*****************************************************************************
+ Copyright (c) 2014, Intel Corp.
+ All rights reserved.
+
+ Redistribution and use in source and binary forms, with or without
+ modification, are permitted provided that the following conditions are met:
+
+ * Redistributions of source code must retain the above copyright notice,
+ this list of conditions and the following disclaimer.
+ * Redistributions in binary form must reproduce the above copyright
+ notice, this list of conditions and the following disclaimer in the
+ documentation and/or other materials provided with the distribution.
+ * Neither the name of Intel Corporation nor the names of its contributors
+ may be used to endorse or promote products derived from this software
+ without specific prior written permission.
+
+ THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
+ AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
+ IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
+ ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE
+ LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
+ CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
+ SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
+ INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
+ CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
+ ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF
+ THE POSSIBILITY OF SUCH DAMAGE.
+*****************************************************************************
+* Contents: Native middle-level C interface to LAPACK function slaswp
+* Author: Intel Corporation
+* Generated November, 2011
+*****************************************************************************/
+
+#include "lapacke_utils.h"
+
+lapack_int LAPACKE_slascl_work( int matrix_layout, char type, lapack_int kl,
+ lapack_int ku, float cfrom, float cto,
+ lapack_int m, lapack_int n, float* a,
+ lapack_int lda )
+{
+ lapack_int info = 0;
+ if( matrix_layout == LAPACK_COL_MAJOR ) {
+ /* Call LAPACK function and adjust info */
+ LAPACK_slascl( &type, &kl, &ku, &cfrom, &cto, &m, &n, a, &lda, &info);
+ if( info < 0 ) {
+ info = info - 1;
+ }
+ } else if( matrix_layout == LAPACK_ROW_MAJOR ) {
+ lapack_int lda_t = MAX(1,lda);
+ float* a_t = NULL;
+ /* Check leading dimension(s) */
+ if( lda < n ) {
+ info = -9;
+ LAPACKE_xerbla( "LAPACKE_slascl_work", info );
+ return info;
+ }
+ /* Allocate memory for temporary array(s) */
+ a_t = (float*)LAPACKE_malloc( sizeof(float) * lda_t * MAX(1,n) );
+ if( a_t == NULL ) {
+ info = LAPACK_TRANSPOSE_MEMORY_ERROR;
+ goto exit_level_0;
+ }
+ /* Transpose input matrices */
+ LAPACKE_sge_trans( matrix_layout, lda, n, a, lda, a_t, lda_t );
+ /* Call LAPACK function and adjust info */
+ LAPACK_slascl( &type, &kl, &ku, &cfrom, &cto, &m, &n, a_t, &lda_t, &info);
+ info = 0; /* LAPACK call is ok! */
+ /* Transpose output matrices */
+ LAPACKE_sge_trans( LAPACK_COL_MAJOR, lda, n, a_t, lda_t, a, lda );
+ /* Release memory and exit */
+ LAPACKE_free( a_t );
+exit_level_0:
+ if( info == LAPACK_TRANSPOSE_MEMORY_ERROR ) {
+ LAPACKE_xerbla( "LAPACKE_slascl_work", info );
+ }
+ } else {
+ info = -1;
+ LAPACKE_xerbla( "LAPACKE_slascl_work", info );
+ }
+ return info;
+}
--- /dev/null
+/*****************************************************************************
+ Copyright (c) 2014, Intel Corp.
+ All rights reserved.
+
+ Redistribution and use in source and binary forms, with or without
+ modification, are permitted provided that the following conditions are met:
+
+ * Redistributions of source code must retain the above copyright notice,
+ this list of conditions and the following disclaimer.
+ * Redistributions in binary form must reproduce the above copyright
+ notice, this list of conditions and the following disclaimer in the
+ documentation and/or other materials provided with the distribution.
+ * Neither the name of Intel Corporation nor the names of its contributors
+ may be used to endorse or promote products derived from this software
+ without specific prior written permission.
+
+ THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
+ AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
+ IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
+ ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE
+ LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
+ CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
+ SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
+ INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
+ CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
+ ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF
+ THE POSSIBILITY OF SUCH DAMAGE.
+*****************************************************************************
+* Contents: Native high-level C interface to LAPACK function dlaswp
+* Author: Intel Corporation
+* Generated November, 2011
+*****************************************************************************/
+
+#include "lapacke_utils.h"
+
+lapack_int LAPACKE_zlascl( int matrix_layout, char type, lapack_int kl,
+ lapack_int ku, double cfrom, double cto,
+ lapack_int m, lapack_int n, lapack_complex_double* a,
+ lapack_int lda )
+{
+ if( matrix_layout != LAPACK_COL_MAJOR && matrix_layout != LAPACK_ROW_MAJOR ) {
+ LAPACKE_xerbla( "LAPACKE_zlascl", -1 );
+ return -1;
+ }
+#ifndef LAPACK_zISABLE_NAN_CHECK
+ /* Optionally check input matrices for NaNs */
+ switch (type) {
+ case 'G':
+ if( LAPACKE_zge_nancheck( matrix_layout, lda, n, a, lda ) ) {
+ return -9;
+ }
+ break;
+ case 'L':
+ // TYPE = 'L' - lower triangular matrix.
+ if( LAPACKE_ztr_nancheck( matrix_layout, 'L', 'N', n, a, lda ) ) {
+ return -9;
+ }
+ break;
+ case 'U':
+ // TYPE = 'U' - upper triangular matrix
+ if( LAPACKE_ztr_nancheck( matrix_layout, 'U', 'N', n, a, lda ) ) {
+ return -9;
+ }
+ break;
+ case 'H':
+ // TYPE = 'H' - upper Hessenberg matrix
+ if( LAPACKE_zhs_nancheck( matrix_layout, n, a, lda ) ) {
+ return -9;
+ }
+ break;
+ case 'B':
+ // TYPE = 'B' - A is a symmetric band matrix with lower bandwidth KL
+ // and upper bandwidth KU and with the only the lower
+ // half stored.
+ if( LAPACKE_zsb_nancheck( matrix_layout, 'L', n, kl, a, lda ) ) {
+ return -9;
+ }
+ break;
+ case 'Q':
+ // TYPE = 'Q' - A is a symmetric band matrix with lower bandwidth KL
+ // and upper bandwidth KU and with the only the upper
+ // half stored.
+ if( LAPACKE_zsb_nancheck( matrix_layout, 'U', n, ku, a, lda ) ) {
+ return -9;
+ }
+ break;
+ case 'Z':
+ // TYPE = 'Z' - A is a band matrix with lower bandwidth KL and upper
+ // bandwidth KU. See DGBTRF for storage details.
+ if( LAPACKE_zgb_nancheck( matrix_layout, n, n, kl, kl+ku, a, lda ) ) {
+ return -6;
+ }
+ break;
+ }
+#endif
+ return LAPACKE_zlascl_work( matrix_layout, type, kl, ku, cfrom, cto, m, n, a, lda );
+}
--- /dev/null
+/*****************************************************************************
+ Copyright (c) 2014, Intel Corp.
+ All rights reserved.
+
+ Redistribution and use in source and binary forms, with or without
+ modification, are permitted provided that the following conditions are met:
+
+ * Redistributions of source code must retain the above copyright notice,
+ this list of conditions and the following disclaimer.
+ * Redistributions in binary form must reproduce the above copyright
+ notice, this list of conditions and the following disclaimer in the
+ documentation and/or other materials provided with the distribution.
+ * Neither the name of Intel Corporation nor the names of its contributors
+ may be used to endorse or promote products derived from this software
+ without specific prior written permission.
+
+ THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
+ AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
+ IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
+ ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE
+ LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
+ CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
+ SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
+ INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
+ CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
+ ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF
+ THE POSSIBILITY OF SUCH DAMAGE.
+*****************************************************************************
+* Contents: Native middle-level C interface to LAPACK function dlaswp
+* Author: Intel Corporation
+* Generated November, 2011
+*****************************************************************************/
+
+#include "lapacke_utils.h"
+
+lapack_int LAPACKE_zlascl_work( int matrix_layout, char type, lapack_int kl,
+ lapack_int ku, double cfrom, double cto,
+ lapack_int m, lapack_int n, lapack_complex_double* a,
+ lapack_int lda )
+{
+ lapack_int info = 0;
+ if( matrix_layout == LAPACK_COL_MAJOR ) {
+ /* Call LAPACK function and adjust info */
+ LAPACK_zlascl( &type, &kl, &ku, &cfrom, &cto, &m, &n, a, &lda, &info);
+ if( info < 0 ) {
+ info = info - 1;
+ }
+ } else if( matrix_layout == LAPACK_ROW_MAJOR ) {
+ lapack_int lda_t = MAX(1,lda);
+ lapack_complex_double* a_t = NULL;
+ /* Check leading dimension(s) */
+ if( lda < n ) {
+ info = -9;
+ LAPACKE_xerbla( "LAPACKE_zlascl_work", info );
+ return info;
+ }
+ /* Allocate memory for temporary array(s) */
+ a_t = (lapack_complex_double*)
+ LAPACKE_malloc( sizeof(lapack_complex_double) * lda_t * MAX(1,n) );
+ if( a_t == NULL ) {
+ info = LAPACK_TRANSPOSE_MEMORY_ERROR;
+ goto exit_level_0;
+ }
+ /* Transpose input matrices */
+ LAPACKE_zge_trans( matrix_layout, lda, n, a, lda, a_t, lda_t );
+ /* Call LAPACK function and adjust info */
+ LAPACK_zlascl( &type, &kl, &ku, &cfrom, &cto, &m, &n, a_t, &lda_t, &info);
+ info = 0; /* LAPACK call is ok! */
+ /* Transpose output matrices */
+ LAPACKE_zge_trans( LAPACK_COL_MAJOR, lda, n, a_t, lda_t, a, lda );
+ /* Release memory and exit */
+ LAPACKE_free( a_t );
+exit_level_0:
+ if( info == LAPACK_TRANSPOSE_MEMORY_ERROR ) {
+ LAPACKE_xerbla( "LAPACKE_zlascl_work", info );
+ }
+ } else {
+ info = -1;
+ LAPACKE_xerbla( "LAPACKE_zlascl_work", info );
+ }
+ return info;
+}