bl1_ger.c File Reference

(r12600)

Functions
void	bl1_sger (conj1_t conjx, conj1_t conjy, int m, int n, float *alpha, float *x, int incx, float *y, int incy, float *a, int a_rs, int a_cs)
void	bl1_dger (conj1_t conjx, conj1_t conjy, int m, int n, double *alpha, double *x, int incx, double *y, int incy, double *a, int a_rs, int a_cs)
void	bl1_cger (conj1_t conjx, conj1_t conjy, int m, int n, scomplex *alpha, scomplex *x, int incx, scomplex *y, int incy, scomplex *a, int a_rs, int a_cs)
void	bl1_zger (conj1_t conjx, conj1_t conjy, int m, int n, dcomplex *alpha, dcomplex *x, int incx, dcomplex *y, int incy, dcomplex *a, int a_rs, int a_cs)
void	bl1_sger_blas (int m, int n, float *alpha, float *x, int incx, float *y, int incy, float *a, int lda)
void	bl1_dger_blas (int m, int n, double *alpha, double *x, int incx, double *y, int incy, double *a, int lda)
void	bl1_cgerc_blas (int m, int n, scomplex *alpha, scomplex *x, int incx, scomplex *y, int incy, scomplex *a, int lda)
void	bl1_cgeru_blas (int m, int n, scomplex *alpha, scomplex *x, int incx, scomplex *y, int incy, scomplex *a, int lda)
void	bl1_zgerc_blas (int m, int n, dcomplex *alpha, dcomplex *x, int incx, dcomplex *y, int incy, dcomplex *a, int lda)
void	bl1_zgeru_blas (int m, int n, dcomplex *alpha, dcomplex *x, int incx, dcomplex *y, int incy, dcomplex *a, int lda)

---

Function Documentation

void bl1_cger	(	conj1_t	conjx,
		conj1_t	conjy,
		int	m,
		int	n,
		scomplex *	alpha,
		scomplex *	x,
		int	incx,
		scomplex *	y,
		int	incy,
		scomplex *	a,
		int	a_rs,
		int	a_cs
	)

References bl1_callocv(), bl1_ccopyv(), bl1_ccreate_contigm(), bl1_cfree(), bl1_cfree_saved_contigm(), bl1_cgerc_blas(), bl1_cgeru_blas(), bl1_is_conj(), bl1_is_row_storage(), bl1_zero_dim2(), and BLIS1_CONJUGATE.

Referenced by FLA_Apply_H2_UT_l_opc_var1(), FLA_Apply_H2_UT_r_opc_var1(), FLA_Apply_HUD_UT_l_opc_var1(), FLA_Bidiag_UT_u_step_ofc_var3(), FLA_Bidiag_UT_u_step_opc_var2(), FLA_Bidiag_UT_u_step_opc_var3(), FLA_Eig_gest_il_opc_var3(), FLA_Eig_gest_il_opc_var4(), FLA_Eig_gest_iu_opc_var3(), FLA_Eig_gest_iu_opc_var4(), FLA_Eig_gest_nl_opc_var4(), FLA_Eig_gest_nu_opc_var4(), FLA_Ger_external(), FLA_Gerc_external(), FLA_Hess_UT_step_ofc_var2(), FLA_Hess_UT_step_ofc_var3(), FLA_Hess_UT_step_ofc_var4(), FLA_Hess_UT_step_opc_var2(), FLA_Hess_UT_step_opc_var3(), FLA_Hess_UT_step_opc_var4(), FLA_LU_nopiv_opc_var5(), FLA_LU_piv_opc_var5(), FLA_Lyap_h_opc_var2(), FLA_Lyap_n_opc_var2(), FLA_SA_LU_unb(), FLA_Trinv_ln_opc_var3(), FLA_Trinv_ln_opc_var4(), FLA_Trinv_lu_opc_var3(), FLA_Trinv_lu_opc_var4(), FLA_Trinv_un_opc_var3(), FLA_Trinv_un_opc_var4(), FLA_Trinv_uu_opc_var3(), and FLA_Trinv_uu_opc_var4().

{
    int       m_save    = m;
    int       n_save    = n;
    scomplex* a_save    = a;
    int       a_rs_save = a_rs;
    int       a_cs_save = a_cs;
    scomplex* x_conj;
    int       incx_conj;
    int       lda, inca;

    // Return early if possible.
    if ( bl1_zero_dim2( m, n ) ) return;

    // If necessary, allocate, initialize, and use a temporary contiguous
    // copy of the matrix rather than the original matrix.
    bl1_ccreate_contigm( m,
                         n,
                         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( m, n );
        bl1_swap_ints( lda, inca );
        bl1_swap_ints( incx, incy );
        bl1_swap_conj( conjx, conjy );
        bl1_cswap_pointers( x, y );
    }

    // Initialize with values assuming no conjugation of x.
    x_conj    = x;
    incx_conj = incx;

    // We need a temporary vector for the cases when x is conjugated.
    if ( bl1_is_conj( conjx ) )
    {
        x_conj    = bl1_callocv( m );
        incx_conj = 1;

        bl1_ccopyv( BLIS1_CONJUGATE,
                    m,
                    x,      incx,
                    x_conj, incx_conj );
    }

    // Conjugation of y is supported in the BLAS.
    if ( bl1_is_conj( conjy ) )
    {
        bl1_cgerc_blas( m,
                        n,
                        alpha,
                        x_conj, incx_conj,
                        y,      incy, 
                        a,      lda );
    }
    else
    {
        bl1_cgeru_blas( m,
                        n,
                        alpha,
                        x_conj, incx_conj,
                        y,      incy, 
                        a,      lda );
    }

    // Free the temporary conjugated x vector.
    if ( bl1_is_conj( conjx ) )
        bl1_cfree( x_conj );

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

void bl1_cgerc_blas	(	int	m,
		int	n,
		scomplex *	alpha,
		scomplex *	x,
		int	incx,
		scomplex *	y,
		int	incy,
		scomplex *	a,
		int	lda
	)

References cblas_cgerc(), CblasColMajor, and F77_cgerc().

Referenced by bl1_cger().

{
#ifdef BLIS1_ENABLE_CBLAS_INTERFACES
    enum CBLAS_ORDER cblas_order = CblasColMajor;

    cblas_cgerc( cblas_order,
                 m,
                 n,
                 alpha,
                 x, incx,
                 y, incy,
                 a, lda );
#else
    F77_cgerc ( &m,
                &n,
                alpha,
                x, &incx,
                y, &incy,
                a, &lda );
#endif
}

void bl1_cgeru_blas	(	int	m,
		int	n,
		scomplex *	alpha,
		scomplex *	x,
		int	incx,
		scomplex *	y,
		int	incy,
		scomplex *	a,
		int	lda
	)

References cblas_cgeru(), CblasColMajor, and F77_cgeru().

Referenced by bl1_cger().

{
#ifdef BLIS1_ENABLE_CBLAS_INTERFACES
    enum CBLAS_ORDER cblas_order = CblasColMajor;

    cblas_cgeru( cblas_order,
                 m,
                 n,
                 alpha,
                 x, incx,
                 y, incy,
                 a, lda );
#else
    F77_cgeru ( &m,
                &n,
                alpha,
                x, &incx,
                y, &incy,
                a, &lda );
#endif
}

void bl1_dger	(	conj1_t	conjx,
		conj1_t	conjy,
		int	m,
		int	n,
		double *	alpha,
		double *	x,
		int	incx,
		double *	y,
		int	incy,
		double *	a,
		int	a_rs,
		int	a_cs
	)

References bl1_dcreate_contigm(), bl1_dfree_saved_contigm(), bl1_dger_blas(), bl1_is_row_storage(), and bl1_zero_dim2().

Referenced by FLA_Apply_H2_UT_l_opd_var1(), FLA_Apply_H2_UT_r_opd_var1(), FLA_Apply_HUD_UT_l_opd_var1(), FLA_Bidiag_UT_u_step_ofd_var3(), FLA_Bidiag_UT_u_step_opd_var2(), FLA_Bidiag_UT_u_step_opd_var3(), FLA_Eig_gest_il_opd_var3(), FLA_Eig_gest_il_opd_var4(), FLA_Eig_gest_iu_opd_var3(), FLA_Eig_gest_iu_opd_var4(), FLA_Eig_gest_nl_opd_var4(), FLA_Eig_gest_nu_opd_var4(), FLA_Ger_external(), FLA_Gerc_external(), FLA_Hess_UT_step_ofd_var2(), FLA_Hess_UT_step_ofd_var3(), FLA_Hess_UT_step_ofd_var4(), FLA_Hess_UT_step_opd_var2(), FLA_Hess_UT_step_opd_var3(), FLA_Hess_UT_step_opd_var4(), FLA_LU_nopiv_opd_var5(), FLA_LU_piv_opd_var5(), FLA_Lyap_h_opd_var2(), FLA_Lyap_n_opd_var2(), FLA_SA_LU_unb(), FLA_Trinv_ln_opd_var3(), FLA_Trinv_ln_opd_var4(), FLA_Trinv_lu_opd_var3(), FLA_Trinv_lu_opd_var4(), FLA_Trinv_un_opd_var3(), FLA_Trinv_un_opd_var4(), FLA_Trinv_uu_opd_var3(), and FLA_Trinv_uu_opd_var4().

{
    int       m_save    = m;
    int       n_save    = n;
    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_dim2( m, n ) ) return;

    // If necessary, allocate, initialize, and use a temporary contiguous
    // copy of the matrix rather than the original matrix.
    bl1_dcreate_contigm( m,
                         n,
                         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( m, n );
        bl1_swap_ints( lda, inca );
        bl1_swap_ints( incx, incy );
        bl1_swap_conj( conjx, conjy );
        bl1_dswap_pointers( x, y );
    }

    // Initialize with values assuming no conjugation of x.
    bl1_dger_blas( m,
                   n,
                   alpha,
                   x, incx,
                   y, incy,
                   a, lda );

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

void bl1_dger_blas	(	int	m,
		int	n,
		double *	alpha,
		double *	x,
		int	incx,
		double *	y,
		int	incy,
		double *	a,
		int	lda
	)

References cblas_dger(), CblasColMajor, and F77_dger().

Referenced by bl1_dger().

{
#ifdef BLIS1_ENABLE_CBLAS_INTERFACES
    enum CBLAS_ORDER cblas_order = CblasColMajor;

    cblas_dger( cblas_order,
                m,
                n,
                *alpha,
                x, incx,
                y, incy,
                a, lda );
#else
    F77_dger( &m,
              &n,
              alpha,
              x, &incx,
              y, &incy,
              a, &lda );
#endif
}

void bl1_sger	(	conj1_t	conjx,
		conj1_t	conjy,
		int	m,
		int	n,
		float *	alpha,
		float *	x,
		int	incx,
		float *	y,
		int	incy,
		float *	a,
		int	a_rs,
		int	a_cs
	)

References bl1_is_row_storage(), bl1_screate_contigm(), bl1_sfree_saved_contigm(), bl1_sger_blas(), and bl1_zero_dim2().

Referenced by FLA_Apply_H2_UT_l_ops_var1(), FLA_Apply_H2_UT_r_ops_var1(), FLA_Apply_HUD_UT_l_ops_var1(), FLA_Bidiag_UT_u_step_ofs_var3(), FLA_Bidiag_UT_u_step_ops_var2(), FLA_Bidiag_UT_u_step_ops_var3(), FLA_Eig_gest_il_ops_var3(), FLA_Eig_gest_il_ops_var4(), FLA_Eig_gest_iu_ops_var3(), FLA_Eig_gest_iu_ops_var4(), FLA_Eig_gest_nl_ops_var4(), FLA_Eig_gest_nu_ops_var4(), FLA_Ger_external(), FLA_Gerc_external(), FLA_Hess_UT_step_ofs_var2(), FLA_Hess_UT_step_ofs_var3(), FLA_Hess_UT_step_ofs_var4(), FLA_Hess_UT_step_ops_var2(), FLA_Hess_UT_step_ops_var3(), FLA_Hess_UT_step_ops_var4(), FLA_LU_nopiv_ops_var5(), FLA_LU_piv_ops_var5(), FLA_Lyap_h_ops_var2(), FLA_Lyap_n_ops_var2(), FLA_SA_LU_unb(), FLA_Trinv_ln_ops_var3(), FLA_Trinv_ln_ops_var4(), FLA_Trinv_lu_ops_var3(), FLA_Trinv_lu_ops_var4(), FLA_Trinv_un_ops_var3(), FLA_Trinv_un_ops_var4(), FLA_Trinv_uu_ops_var3(), and FLA_Trinv_uu_ops_var4().

{
    int       m_save    = m;
    int       n_save    = n;
    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_dim2( m, n ) ) return;

    // If necessary, allocate, initialize, and use a temporary contiguous
    // copy of the matrix rather than the original matrix.
    bl1_screate_contigm( m,
                         n,
                         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( m, n );
        bl1_swap_ints( lda, inca );
        bl1_swap_ints( incx, incy );
        bl1_swap_conj( conjx, conjy );
        bl1_sswap_pointers( x, y );
    }

    // Initialize with values assuming no conjugation of x.
    bl1_sger_blas( m,
                   n,
                   alpha,
                   x, incx,
                   y, incy,
                   a, lda );

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

void bl1_sger_blas	(	int	m,
		int	n,
		float *	alpha,
		float *	x,
		int	incx,
		float *	y,
		int	incy,
		float *	a,
		int	lda
	)

References cblas_sger(), CblasColMajor, and F77_sger().

Referenced by bl1_sger().

{
#ifdef BLIS1_ENABLE_CBLAS_INTERFACES
    enum CBLAS_ORDER cblas_order = CblasColMajor;

    cblas_sger( cblas_order,
                m,
                n,
                *alpha,
                x, incx,
                y, incy,
                a, lda );
#else
    F77_sger( &m,
              &n,
              alpha,
              x, &incx,
              y, &incy,
              a, &lda );
#endif
}

void bl1_zger	(	conj1_t	conjx,
		conj1_t	conjy,
		int	m,
		int	n,
		dcomplex *	alpha,
		dcomplex *	x,
		int	incx,
		dcomplex *	y,
		int	incy,
		dcomplex *	a,
		int	a_rs,
		int	a_cs
	)

References bl1_is_conj(), bl1_is_row_storage(), bl1_zallocv(), bl1_zcopyv(), bl1_zcreate_contigm(), bl1_zero_dim2(), bl1_zfree(), bl1_zfree_saved_contigm(), bl1_zgerc_blas(), bl1_zgeru_blas(), and BLIS1_CONJUGATE.

Referenced by FLA_Apply_H2_UT_l_opz_var1(), FLA_Apply_H2_UT_r_opz_var1(), FLA_Apply_HUD_UT_l_opz_var1(), FLA_Bidiag_UT_u_step_ofz_var3(), FLA_Bidiag_UT_u_step_opz_var2(), FLA_Bidiag_UT_u_step_opz_var3(), FLA_Eig_gest_il_opz_var3(), FLA_Eig_gest_il_opz_var4(), FLA_Eig_gest_iu_opz_var3(), FLA_Eig_gest_iu_opz_var4(), FLA_Eig_gest_nl_opz_var4(), FLA_Eig_gest_nu_opz_var4(), FLA_Ger_external(), FLA_Gerc_external(), FLA_Hess_UT_step_ofz_var2(), FLA_Hess_UT_step_ofz_var3(), FLA_Hess_UT_step_ofz_var4(), FLA_Hess_UT_step_opz_var2(), FLA_Hess_UT_step_opz_var3(), FLA_Hess_UT_step_opz_var4(), FLA_LU_nopiv_opz_var5(), FLA_LU_piv_opz_var5(), FLA_Lyap_h_opz_var2(), FLA_Lyap_n_opz_var2(), FLA_SA_LU_unb(), FLA_Trinv_ln_opz_var3(), FLA_Trinv_ln_opz_var4(), FLA_Trinv_lu_opz_var3(), FLA_Trinv_lu_opz_var4(), FLA_Trinv_un_opz_var3(), FLA_Trinv_un_opz_var4(), FLA_Trinv_uu_opz_var3(), and FLA_Trinv_uu_opz_var4().

{
    int       m_save    = m;
    int       n_save    = n;
    dcomplex* a_save    = a;
    int       a_rs_save = a_rs;
    int       a_cs_save = a_cs;
    dcomplex* x_conj;
    int       incx_conj;
    int       lda, inca;

    // Return early if possible.
    if ( bl1_zero_dim2( m, n ) ) return;

    // If necessary, allocate, initialize, and use a temporary contiguous
    // copy of the matrix rather than the original matrix.
    bl1_zcreate_contigm( m,
                         n,
                         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( m, n );
        bl1_swap_ints( lda, inca );
        bl1_swap_ints( incx, incy );
        bl1_swap_conj( conjx, conjy );
        bl1_zswap_pointers( x, y );
    }

    // Initialize with values assuming no conjugation of x.
    x_conj    = x;
    incx_conj = incx;

    // We need a temporary vector for the cases when x is conjugated.
    if ( bl1_is_conj( conjx ) )
    {
        x_conj    = bl1_zallocv( m );
        incx_conj = 1;

        bl1_zcopyv( BLIS1_CONJUGATE,
                    m,
                    x,      incx,
                    x_conj, incx_conj );
    }

    // Conjugation of y is supported in the BLAS.
    if ( bl1_is_conj( conjy ) )
    {
        bl1_zgerc_blas( m,
                        n,
                        alpha,
                        x_conj, incx_conj,
                        y,      incy, 
                        a,      lda );
    }
    else
    {
        bl1_zgeru_blas( m,
                        n,
                        alpha,
                        x_conj, incx_conj,
                        y,      incy, 
                        a,      lda );
    }

    // Free the temporary conjugated x vector.
    if ( bl1_is_conj( conjx ) )
        bl1_zfree( x_conj );

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

void bl1_zgerc_blas	(	int	m,
		int	n,
		dcomplex *	alpha,
		dcomplex *	x,
		int	incx,
		dcomplex *	y,
		int	incy,
		dcomplex *	a,
		int	lda
	)

References cblas_zgerc(), CblasColMajor, and F77_zgerc().

Referenced by bl1_zger().

{
#ifdef BLIS1_ENABLE_CBLAS_INTERFACES
    enum CBLAS_ORDER cblas_order = CblasColMajor;

    cblas_zgerc( cblas_order,
                 m,
                 n,
                 alpha,
                 x, incx,
                 y, incy,
                 a, lda );
#else
    F77_zgerc ( &m,
                &n,
                alpha,
                x, &incx,
                y, &incy,
                a, &lda );
#endif
}

void bl1_zgeru_blas	(	int	m,
		int	n,
		dcomplex *	alpha,
		dcomplex *	x,
		int	incx,
		dcomplex *	y,
		int	incy,
		dcomplex *	a,
		int	lda
	)

References cblas_zgeru(), CblasColMajor, and F77_zgeru().

Referenced by bl1_zger().

{
#ifdef BLIS1_ENABLE_CBLAS_INTERFACES
    enum CBLAS_ORDER cblas_order = CblasColMajor;

    cblas_zgeru( cblas_order,
                 m,
                 n,
                 alpha,
                 x, incx,
                 y, incy,
                 a, lda );
#else
    F77_zgeru ( &m,
                &n,
                alpha,
                x, &incx,
                y, &incy,
                a, &lda );
#endif
}