bl1_dot2s.c File Reference

(r12600)

Functions
void	bl1_sdot2s (conj1_t conj, int n, float *alpha, float *x, int incx, float *y, int incy, float *beta, float *rho)
void	bl1_ddot2s (conj1_t conj, int n, double *alpha, double *x, int incx, double *y, int incy, double *beta, double *rho)
void	bl1_cdot2s (conj1_t conj, int n, scomplex *alpha, scomplex *x, int incx, scomplex *y, int incy, scomplex *beta, scomplex *rho)
void	bl1_zdot2s (conj1_t conj, int n, dcomplex *alpha, dcomplex *x, int incx, dcomplex *y, int incy, dcomplex *beta, dcomplex *rho)

---

Function Documentation

void bl1_cdot2s	(	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_Dot2cs_external(), FLA_Dot2s_external(), FLA_Eig_gest_il_opc_var1(), FLA_Eig_gest_il_opc_var2(), FLA_Eig_gest_il_opc_var3(), FLA_Eig_gest_iu_opc_var1(), FLA_Eig_gest_iu_opc_var2(), FLA_Eig_gest_iu_opc_var3(), FLA_Eig_gest_nl_opc_var1(), FLA_Eig_gest_nl_opc_var2(), FLA_Eig_gest_nu_opc_var1(), FLA_Eig_gest_nu_opc_var2(), FLA_Lyap_h_opc_var1(), FLA_Lyap_h_opc_var2(), FLA_Lyap_h_opc_var3(), FLA_Lyap_n_opc_var1(), FLA_Lyap_n_opc_var2(), and FLA_Lyap_n_opc_var3().

{
    scomplex dotxy;
    scomplex dotyx;
    scomplex alpha_d    = *alpha;
    scomplex alphac_d   = *alpha;
    scomplex beta_d     = *beta;
    scomplex rho_d      = *rho;

    alphac_d.imag *= -1.0F;

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

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

    rho->real = beta_d.real   * rho_d.real - beta_d.imag   * rho_d.imag +
                alpha_d.real  * dotxy.real - alpha_d.imag  * dotxy.imag +
                alphac_d.real * dotyx.real - alphac_d.imag * dotyx.imag; 
    rho->imag = beta_d.real   * rho_d.imag + beta_d.imag   * rho_d.real +
                alpha_d.real  * dotxy.imag + alpha_d.imag  * dotxy.real +
                alphac_d.real * dotyx.imag + alphac_d.imag * dotyx.real; 
}

void bl1_ddot2s	(	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_Dot2cs_external(), FLA_Dot2s_external(), FLA_Eig_gest_il_opd_var1(), FLA_Eig_gest_il_opd_var2(), FLA_Eig_gest_il_opd_var3(), FLA_Eig_gest_iu_opd_var1(), FLA_Eig_gest_iu_opd_var2(), FLA_Eig_gest_iu_opd_var3(), FLA_Eig_gest_nl_opd_var1(), FLA_Eig_gest_nl_opd_var2(), FLA_Eig_gest_nu_opd_var1(), FLA_Eig_gest_nu_opd_var2(), FLA_Lyap_h_opd_var1(), FLA_Lyap_h_opd_var2(), FLA_Lyap_h_opd_var3(), FLA_Lyap_n_opd_var1(), FLA_Lyap_n_opd_var2(), and FLA_Lyap_n_opd_var3().

{
    double dot;

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

    *rho = (*beta) * (*rho) + 2.0 * (*alpha) * dot;
}

void bl1_sdot2s	(	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_Dot2cs_external(), FLA_Dot2s_external(), FLA_Eig_gest_il_ops_var1(), FLA_Eig_gest_il_ops_var2(), FLA_Eig_gest_il_ops_var3(), FLA_Eig_gest_iu_ops_var1(), FLA_Eig_gest_iu_ops_var2(), FLA_Eig_gest_iu_ops_var3(), FLA_Eig_gest_nl_ops_var1(), FLA_Eig_gest_nl_ops_var2(), FLA_Eig_gest_nu_ops_var1(), FLA_Eig_gest_nu_ops_var2(), FLA_Lyap_h_ops_var1(), FLA_Lyap_h_ops_var2(), FLA_Lyap_h_ops_var3(), FLA_Lyap_n_ops_var1(), FLA_Lyap_n_ops_var2(), and FLA_Lyap_n_ops_var3().

{
    float dot;

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

    *rho = (*beta) * (*rho) + 2.0F * (*alpha) * dot;
}

void bl1_zdot2s	(	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_Dot2cs_external(), FLA_Dot2s_external(), FLA_Eig_gest_il_opz_var1(), FLA_Eig_gest_il_opz_var2(), FLA_Eig_gest_il_opz_var3(), FLA_Eig_gest_iu_opz_var1(), FLA_Eig_gest_iu_opz_var2(), FLA_Eig_gest_iu_opz_var3(), FLA_Eig_gest_nl_opz_var1(), FLA_Eig_gest_nl_opz_var2(), FLA_Eig_gest_nu_opz_var1(), FLA_Eig_gest_nu_opz_var2(), FLA_Lyap_h_opz_var1(), FLA_Lyap_h_opz_var2(), FLA_Lyap_h_opz_var3(), FLA_Lyap_n_opz_var1(), FLA_Lyap_n_opz_var2(), and FLA_Lyap_n_opz_var3().

{
    dcomplex dotxy;
    dcomplex dotyx;
    dcomplex alpha_d    = *alpha;
    dcomplex alphac_d   = *alpha;
    dcomplex beta_d     = *beta;
    dcomplex rho_d      = *rho;

    alphac_d.imag *= -1.0;

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

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

    rho->real = beta_d.real   * rho_d.real - beta_d.imag   * rho_d.imag +
                alpha_d.real  * dotxy.real - alpha_d.imag  * dotxy.imag +
                alphac_d.real * dotyx.real - alphac_d.imag * dotyx.imag; 
    rho->imag = beta_d.real   * rho_d.imag + beta_d.imag   * rho_d.real +
                alpha_d.real  * dotxy.imag + alpha_d.imag  * dotxy.real +
                alphac_d.real * dotyx.imag + alphac_d.imag * dotyx.real; 
}