bl1_symv.c File Reference

(r12600)

Functions
void	bl1_ssymv (uplo1_t uplo, int m, float *alpha, float *a, int a_rs, int a_cs, float *x, int incx, float *beta, float *y, int incy)
void	bl1_dsymv (uplo1_t uplo, int m, double *alpha, double *a, int a_rs, int a_cs, double *x, int incx, double *beta, double *y, int incy)
void	bl1_csymv (uplo1_t uplo, int m, scomplex *alpha, scomplex *a, int a_rs, int a_cs, scomplex *x, int incx, scomplex *beta, scomplex *y, int incy)
void	bl1_zsymv (uplo1_t uplo, int m, dcomplex *alpha, dcomplex *a, int a_rs, int a_cs, dcomplex *x, int incx, dcomplex *beta, dcomplex *y, int incy)
void	bl1_ssymv_blas (uplo1_t uplo, int m, float *alpha, float *a, int lda, float *x, int incx, float *beta, float *y, int incy)
void	bl1_dsymv_blas (uplo1_t uplo, int m, double *alpha, double *a, int lda, double *x, int incx, double *beta, double *y, int incy)
void	bl1_csymv_blas (uplo1_t uplo, int m, scomplex *alpha, scomplex *a, int lda, scomplex *x, int incx, scomplex *beta, scomplex *y, int incy)
void	bl1_zsymv_blas (uplo1_t uplo, int m, dcomplex *alpha, dcomplex *a, int lda, dcomplex *x, int incx, dcomplex *beta, dcomplex *y, int incy)

---

Function Documentation

void bl1_csymv	(	uplo1_t	uplo,
		int	m,
		scomplex *	alpha,
		scomplex *	a,
		int	a_rs,
		int	a_cs,
		scomplex *	x,
		int	incx,
		scomplex *	beta,
		scomplex *	y,
		int	incy
	)

References bl1_ccreate_contigmr(), bl1_cfree_contigm(), bl1_csymv_blas(), bl1_is_row_storage(), and bl1_zero_dim1().

Referenced by FLA_Symv_external().

{
    scomplex* a_save    = a;
    int       a_rs_save = a_rs;
    int       a_cs_save = a_cs;
    int       lda, inca;

    // Return early if possible.
    if ( bl1_zero_dim1( m ) ) return;

    // If necessary, allocate, initialize, and use a temporary contiguous
    // copy of the matrix rather than the original matrix.
    bl1_ccreate_contigmr( uplo,
                          m,
                          m,
                          a_save, a_rs_save, a_cs_save,
                          &a,     &a_rs,     &a_cs );

    // Initialize with values assuming column-major storage.
    lda  = a_cs;
    inca = a_rs;

    // If A is a row-major matrix, then we can use the underlying column-major
    // BLAS implementation by fiddling with the parameters.
    if ( bl1_is_row_storage( a_rs, a_cs ) )
    {
        bl1_swap_ints( lda, inca );
        bl1_toggle_uplo( uplo );
    }

    bl1_csymv_blas( uplo,
                    m,
                    alpha,
                    a, lda,
                    x, incx,
                    beta,
                    y, incy );

    // Free the temporary contiguous matrix.
    bl1_cfree_contigm( a_save, a_rs_save, a_cs_save,
                       &a,     &a_rs,     &a_cs );
}

void bl1_csymv_blas	(	uplo1_t	uplo,
		int	m,
		scomplex *	alpha,
		scomplex *	a,
		int	lda,
		scomplex *	x,
		int	incx,
		scomplex *	beta,
		scomplex *	y,
		int	incy
	)

References bl1_callocv(), bl1_ccopyv(), bl1_cfree(), bl1_param_map_to_netlib_side(), bl1_param_map_to_netlib_uplo(), BLIS1_LEFT, BLIS1_NO_CONJUGATE, cblas_csymm(), CblasColMajor, and F77_csymm().

Referenced by bl1_csymv().

{
    scomplex* x_copy;
    scomplex* y_copy;
    int       n   = 1;
    int       ldx = m;
    int       ldy = m;

#ifdef BLIS1_ENABLE_CBLAS_INTERFACES
    enum CBLAS_ORDER cblas_order = CblasColMajor;
    enum CBLAS_SIDE  cblas_side;
    enum CBLAS_UPLO  cblas_uplo;

    bl1_param_map_to_netlib_side( BLIS1_LEFT, &cblas_side );
    bl1_param_map_to_netlib_uplo( uplo, &cblas_uplo );

    x_copy = bl1_callocv( m );
    y_copy = bl1_callocv( m );

    bl1_ccopyv( BLIS1_NO_CONJUGATE,
                m,
                x,      incx,
                x_copy, 1 );

    bl1_ccopyv( BLIS1_NO_CONJUGATE,
                m,
                y,      incy,
                y_copy, 1 );

    cblas_csymm( cblas_order,
                 cblas_side,
                 cblas_uplo,
                 m,
                 n,
                 alpha,
                 a,      lda,
                 x_copy, ldx,
                 beta,
                 y_copy, ldy );

    bl1_ccopyv( BLIS1_NO_CONJUGATE,
                m,
                y_copy, 1,
                y,      incy );

    bl1_cfree( x_copy );
    bl1_cfree( y_copy );

#else
    char blas_side;
    char blas_uplo;

    bl1_param_map_to_netlib_side( BLIS1_LEFT, &blas_side );
    bl1_param_map_to_netlib_uplo( uplo, &blas_uplo );

    x_copy = bl1_callocv( m );
    y_copy = bl1_callocv( m );

    bl1_ccopyv( BLIS1_NO_CONJUGATE,
                m,
                x,      incx,
                x_copy, 1 );

    bl1_ccopyv( BLIS1_NO_CONJUGATE,
                m,
                y,      incy,
                y_copy, 1 );

    F77_csymm ( &blas_side,
                &blas_uplo,
                &m,
                &n,
                alpha,
                a,      &lda,
                x_copy, &ldx,
                beta,
                y_copy, &ldy );

    bl1_ccopyv( BLIS1_NO_CONJUGATE,
                m,
                y_copy, 1,
                y,      incy );

    bl1_cfree( x_copy );
    bl1_cfree( y_copy );
#endif
}

void bl1_dsymv	(	uplo1_t	uplo,
		int	m,
		double *	alpha,
		double *	a,
		int	a_rs,
		int	a_cs,
		double *	x,
		int	incx,
		double *	beta,
		double *	y,
		int	incy
	)

References bl1_dcreate_contigmr(), bl1_dfree_contigm(), bl1_dsymv_blas(), bl1_is_row_storage(), and bl1_zero_dim1().

Referenced by bl1_dhemv(), FLA_Hemv_external(), FLA_Hemvc_external(), FLA_Symv_external(), FLA_Tridiag_UT_l_step_ofd_var2(), FLA_Tridiag_UT_l_step_ofd_var3(), FLA_Tridiag_UT_l_step_opd_var1(), FLA_Tridiag_UT_l_step_opd_var2(), and FLA_Tridiag_UT_l_step_opd_var3().

{
    double*   a_save    = a;
    int       a_rs_save = a_rs;
    int       a_cs_save = a_cs;
    int       lda, inca;

    // Return early if possible.
    if ( bl1_zero_dim1( m ) ) return;

    // If necessary, allocate, initialize, and use a temporary contiguous
    // copy of the matrix rather than the original matrix.
    bl1_dcreate_contigmr( uplo,
                          m,
                          m,
                          a_save, a_rs_save, a_cs_save,
                          &a,     &a_rs,     &a_cs );

    // Initialize with values assuming column-major storage.
    lda  = a_cs;
    inca = a_rs;

    // If A is a row-major matrix, then we can use the underlying column-major
    // BLAS implementation by fiddling with the parameters.
    if ( bl1_is_row_storage( a_rs, a_cs ) )
    {
        bl1_swap_ints( lda, inca );
        bl1_toggle_uplo( uplo );
    }

    bl1_dsymv_blas( uplo,
                    m,
                    alpha,
                    a, lda,
                    x, incx,
                    beta,
                    y, incy );

    // Free the temporary contiguous matrix.
    bl1_dfree_contigm( a_save, a_rs_save, a_cs_save,
                       &a,     &a_rs,     &a_cs );
}

void bl1_dsymv_blas	(	uplo1_t	uplo,
		int	m,
		double *	alpha,
		double *	a,
		int	lda,
		double *	x,
		int	incx,
		double *	beta,
		double *	y,
		int	incy
	)

References bl1_param_map_to_netlib_uplo(), cblas_dsymv(), CblasColMajor, and F77_dsymv().

Referenced by bl1_dsymv().

{
#ifdef BLIS1_ENABLE_CBLAS_INTERFACES
    enum CBLAS_ORDER cblas_order = CblasColMajor;
    enum CBLAS_UPLO  cblas_uplo;

    bl1_param_map_to_netlib_uplo( uplo, &cblas_uplo );

    cblas_dsymv( cblas_order,
                 cblas_uplo,
                 m,
                 *alpha,
                 a, lda,
                 x, incx,
                 *beta,
                 y, incy );
#else
    char blas_uplo;

    bl1_param_map_to_netlib_uplo( uplo, &blas_uplo );

    F77_dsymv( &blas_uplo,
               &m,
               alpha,
               a, &lda,
               x, &incx,
               beta,
               y, &incy );
#endif
}

void bl1_ssymv	(	uplo1_t	uplo,
		int	m,
		float *	alpha,
		float *	a,
		int	a_rs,
		int	a_cs,
		float *	x,
		int	incx,
		float *	beta,
		float *	y,
		int	incy
	)

References bl1_is_row_storage(), bl1_screate_contigmr(), bl1_sfree_contigm(), bl1_ssymv_blas(), and bl1_zero_dim1().

Referenced by bl1_shemv(), FLA_Hemv_external(), FLA_Hemvc_external(), FLA_Symv_external(), FLA_Tridiag_UT_l_step_ofs_var2(), FLA_Tridiag_UT_l_step_ofs_var3(), FLA_Tridiag_UT_l_step_ops_var1(), FLA_Tridiag_UT_l_step_ops_var2(), and FLA_Tridiag_UT_l_step_ops_var3().

{
    float*    a_save    = a;
    int       a_rs_save = a_rs;
    int       a_cs_save = a_cs;
    int       lda, inca;

    // Return early if possible.
    if ( bl1_zero_dim1( m ) ) return;

    // If necessary, allocate, initialize, and use a temporary contiguous
    // copy of the matrix rather than the original matrix.
    bl1_screate_contigmr( uplo,
                          m,
                          m,
                          a_save, a_rs_save, a_cs_save,
                          &a,     &a_rs,     &a_cs );

    // Initialize with values assuming column-major storage.
    lda  = a_cs;
    inca = a_rs;

    // If A is a row-major matrix, then we can use the underlying column-major
    // BLAS implementation by fiddling with the parameters.
    if ( bl1_is_row_storage( a_rs, a_cs ) )
    {
        bl1_swap_ints( lda, inca );
        bl1_toggle_uplo( uplo );
    }

    bl1_ssymv_blas( uplo,
                    m,
                    alpha,
                    a, lda,
                    x, incx,
                    beta,
                    y, incy );

    // Free the temporary contiguous matrix.
    bl1_sfree_contigm( a_save, a_rs_save, a_cs_save,
                       &a,     &a_rs,     &a_cs );
}

void bl1_ssymv_blas	(	uplo1_t	uplo,
		int	m,
		float *	alpha,
		float *	a,
		int	lda,
		float *	x,
		int	incx,
		float *	beta,
		float *	y,
		int	incy
	)

References bl1_param_map_to_netlib_uplo(), cblas_ssymv(), CblasColMajor, and F77_ssymv().

Referenced by bl1_ssymv().

{
#ifdef BLIS1_ENABLE_CBLAS_INTERFACES
    enum CBLAS_ORDER cblas_order = CblasColMajor;
    enum CBLAS_UPLO  cblas_uplo;

    bl1_param_map_to_netlib_uplo( uplo, &cblas_uplo );

    cblas_ssymv( cblas_order,
                 cblas_uplo,
                 m,
                 *alpha,
                 a, lda,
                 x, incx,
                 *beta,
                 y, incy );
#else
    char blas_uplo;

    bl1_param_map_to_netlib_uplo( uplo, &blas_uplo );

    F77_ssymv( &blas_uplo,
               &m,
               alpha,
               a, &lda,
               x, &incx,
               beta,
               y, &incy );
#endif
}

void bl1_zsymv	(	uplo1_t	uplo,
		int	m,
		dcomplex *	alpha,
		dcomplex *	a,
		int	a_rs,
		int	a_cs,
		dcomplex *	x,
		int	incx,
		dcomplex *	beta,
		dcomplex *	y,
		int	incy
	)

References bl1_is_row_storage(), bl1_zcreate_contigmr(), bl1_zero_dim1(), bl1_zfree_contigm(), and bl1_zsymv_blas().

Referenced by FLA_Symv_external().

{
    dcomplex* a_save    = a;
    int       a_rs_save = a_rs;
    int       a_cs_save = a_cs;
    int       lda, inca;

    // Return early if possible.
    if ( bl1_zero_dim1( m ) ) return;

    // If necessary, allocate, initialize, and use a temporary contiguous
    // copy of the matrix rather than the original matrix.
    bl1_zcreate_contigmr( uplo,
                          m,
                          m,
                          a_save, a_rs_save, a_cs_save,
                          &a,     &a_rs,     &a_cs );

    // Initialize with values assuming column-major storage.
    lda  = a_cs;
    inca = a_rs;

    // If A is a row-major matrix, then we can use the underlying column-major
    // BLAS implementation by fiddling with the parameters.
    if ( bl1_is_row_storage( a_rs, a_cs ) )
    {
        bl1_swap_ints( lda, inca );
        bl1_toggle_uplo( uplo );
    }

    bl1_zsymv_blas( uplo,
                    m,
                    alpha,
                    a, lda,
                    x, incx,
                    beta,
                    y, incy );

    // Free the temporary contiguous matrix.
    bl1_zfree_contigm( a_save, a_rs_save, a_cs_save,
                       &a,     &a_rs,     &a_cs );
}

void bl1_zsymv_blas	(	uplo1_t	uplo,
		int	m,
		dcomplex *	alpha,
		dcomplex *	a,
		int	lda,
		dcomplex *	x,
		int	incx,
		dcomplex *	beta,
		dcomplex *	y,
		int	incy
	)

References bl1_param_map_to_netlib_side(), bl1_param_map_to_netlib_uplo(), bl1_zallocv(), bl1_zcopyv(), bl1_zfree(), BLIS1_LEFT, BLIS1_NO_CONJUGATE, cblas_zsymm(), CblasColMajor, and F77_zsymm().

Referenced by bl1_zsymv().

{
    dcomplex*        x_copy;
    dcomplex*        y_copy;
    int              n   = 1;
    int              ldx = m;
    int              ldy = m;

#ifdef BLIS1_ENABLE_CBLAS_INTERFACES
    enum CBLAS_ORDER cblas_order = CblasColMajor;
    enum CBLAS_SIDE  cblas_side;
    enum CBLAS_UPLO  cblas_uplo;

    bl1_param_map_to_netlib_side( BLIS1_LEFT, &cblas_side );
    bl1_param_map_to_netlib_uplo( uplo, &cblas_uplo );

    x_copy = bl1_zallocv( m );
    y_copy = bl1_zallocv( m );

    bl1_zcopyv( BLIS1_NO_CONJUGATE,
                m,
                x,      incx,
                x_copy, 1 );

    bl1_zcopyv( BLIS1_NO_CONJUGATE,
                m,
                y,      incy,
                y_copy, 1 );

    cblas_zsymm( cblas_order,
                 cblas_side,
                 cblas_uplo,
                 m,
                 n,
                 alpha,
                 a,      lda,
                 x_copy, ldx,
                 beta,
                 y_copy, ldy );

    bl1_zcopyv( BLIS1_NO_CONJUGATE,
                m,
                y_copy, 1,
                y,      incy );

    bl1_zfree( x_copy );
    bl1_zfree( y_copy );

#else
    char blas_side;
    char blas_uplo;

    bl1_param_map_to_netlib_side( BLIS1_LEFT, &blas_side );
    bl1_param_map_to_netlib_uplo( uplo, &blas_uplo );

    x_copy = bl1_zallocv( m );
    y_copy = bl1_zallocv( m );

    bl1_zcopyv( BLIS1_NO_CONJUGATE,
                m,
                x,      incx,
                x_copy, 1 );

    bl1_zcopyv( BLIS1_NO_CONJUGATE,
                m,
                y,      incy,
                y_copy, 1 );

    F77_zsymm ( &blas_side,
                &blas_uplo,
                &m,
                &n,
                alpha,
                a,      &lda,
                x_copy, &ldx,
                beta,
                y_copy, &ldy );

    bl1_zcopyv( BLIS1_NO_CONJUGATE,
                m,
                y_copy, 1,
                y,      incy );

    bl1_zfree( x_copy );
    bl1_zfree( y_copy );
#endif
}