bl1_dots.c File Reference

(r12600)

Functions
void	bl1_sdots (conj1_t conj, int n, float *alpha, float *x, int incx, float *y, int incy, float *beta, float *rho)
void	bl1_ddots (conj1_t conj, int n, double *alpha, double *x, int incx, double *y, int incy, double *beta, double *rho)
void	bl1_cdots (conj1_t conj, int n, scomplex *alpha, scomplex *x, int incx, scomplex *y, int incy, scomplex *beta, scomplex *rho)
void	bl1_zdots (conj1_t conj, int n, dcomplex *alpha, dcomplex *x, int incx, dcomplex *y, int incy, dcomplex *beta, dcomplex *rho)

---

Function Documentation

void bl1_cdots	(	conj1_t	conj,
		int	n,
		scomplex *	alpha,
		scomplex *	x,
		int	incx,
		scomplex *	y,
		int	incy,
		scomplex *	beta,
		scomplex *	rho
	)

References bl1_cdot(), scomplex::imag, and scomplex::real.

Referenced by FLA_Chol_l_opc_var1(), FLA_Chol_l_opc_var2(), FLA_Chol_u_opc_var1(), FLA_Chol_u_opc_var2(), FLA_Dotcs_external(), FLA_Dots_external(), FLA_Fused_Ahx_Axpy_Ax_opc_var1(), FLA_Hess_UT_step_opc_var5(), FLA_LU_nopiv_opc_var1(), FLA_LU_nopiv_opc_var2(), FLA_LU_nopiv_opc_var3(), FLA_LU_nopiv_opc_var4(), FLA_LU_piv_opc_var3(), FLA_LU_piv_opc_var4(), FLA_Ttmm_l_opc_var2(), FLA_Ttmm_l_opc_var3(), FLA_Ttmm_u_opc_var2(), and FLA_Ttmm_u_opc_var3().

{
    scomplex rho_orig = *rho;
    scomplex dot_prod;

    bl1_cdot( conj,
              n,
              x, incx,
              y, incy,
              &dot_prod );

    rho->real = beta->real  * rho_orig.real - beta->imag  * rho_orig.imag +
                alpha->real * dot_prod.real - alpha->imag * dot_prod.imag;
    rho->imag = beta->real  * rho_orig.imag + beta->imag  * rho_orig.real +
                alpha->real * dot_prod.imag + alpha->imag * dot_prod.real;
}

void bl1_ddots	(	conj1_t	conj,
		int	n,
		double *	alpha,
		double *	x,
		int	incx,
		double *	y,
		int	incy,
		double *	beta,
		double *	rho
	)

References bl1_ddot().

Referenced by FLA_Chol_l_opd_var1(), FLA_Chol_l_opd_var2(), FLA_Chol_u_opd_var1(), FLA_Chol_u_opd_var2(), FLA_Dotcs_external(), FLA_Dots_external(), FLA_Hess_UT_step_opd_var5(), FLA_LU_nopiv_opd_var1(), FLA_LU_nopiv_opd_var2(), FLA_LU_nopiv_opd_var3(), FLA_LU_nopiv_opd_var4(), FLA_LU_piv_opd_var3(), FLA_LU_piv_opd_var4(), FLA_Ttmm_l_opd_var2(), FLA_Ttmm_l_opd_var3(), FLA_Ttmm_u_opd_var2(), and FLA_Ttmm_u_opd_var3().

{
    double dot_prod;

    bl1_ddot( conj,
              n,
              x, incx,
              y, incy,
              &dot_prod );

    *rho = (*beta) * (*rho) + (*alpha) * dot_prod;
}

void bl1_sdots	(	conj1_t	conj,
		int	n,
		float *	alpha,
		float *	x,
		int	incx,
		float *	y,
		int	incy,
		float *	beta,
		float *	rho
	)

References bl1_sdot().

Referenced by FLA_Chol_l_ops_var1(), FLA_Chol_l_ops_var2(), FLA_Chol_u_ops_var1(), FLA_Chol_u_ops_var2(), FLA_Dotcs_external(), FLA_Dots_external(), FLA_Fused_Ahx_Axpy_Ax_ops_var1(), FLA_Hess_UT_step_ops_var5(), FLA_LU_nopiv_ops_var1(), FLA_LU_nopiv_ops_var2(), FLA_LU_nopiv_ops_var3(), FLA_LU_nopiv_ops_var4(), FLA_LU_piv_ops_var3(), FLA_LU_piv_ops_var4(), FLA_Ttmm_l_ops_var2(), FLA_Ttmm_l_ops_var3(), FLA_Ttmm_u_ops_var2(), and FLA_Ttmm_u_ops_var3().

{
    float dot_prod;

    bl1_sdot( conj,
              n,
              x, incx,
              y, incy,
              &dot_prod );

    *rho = (*beta) * (*rho) + (*alpha) * dot_prod;
}

void bl1_zdots	(	conj1_t	conj,
		int	n,
		dcomplex *	alpha,
		dcomplex *	x,
		int	incx,
		dcomplex *	y,
		int	incy,
		dcomplex *	beta,
		dcomplex *	rho
	)

References bl1_zdot(), dcomplex::imag, and dcomplex::real.

Referenced by FLA_Chol_l_opz_var1(), FLA_Chol_l_opz_var2(), FLA_Chol_u_opz_var1(), FLA_Chol_u_opz_var2(), FLA_Dotcs_external(), FLA_Dots_external(), FLA_Hess_UT_step_opz_var5(), FLA_LU_nopiv_opz_var1(), FLA_LU_nopiv_opz_var2(), FLA_LU_nopiv_opz_var3(), FLA_LU_nopiv_opz_var4(), FLA_LU_piv_opz_var3(), FLA_LU_piv_opz_var4(), FLA_Ttmm_l_opz_var2(), FLA_Ttmm_l_opz_var3(), FLA_Ttmm_u_opz_var2(), and FLA_Ttmm_u_opz_var3().

{
    dcomplex rho_orig = *rho;
    dcomplex dot_prod;

    bl1_zdot( conj,
              n,
              x, incx,
              y, incy,
              &dot_prod );

    rho->real = beta->real  * rho_orig.real - beta->imag  * rho_orig.imag +
                alpha->real * dot_prod.real - alpha->imag * dot_prod.imag;
    rho->imag = beta->real  * rho_orig.imag + beta->imag  * rho_orig.real +
                alpha->real * dot_prod.imag + alpha->imag * dot_prod.real;
}