bl1_invertv.c File Reference

(r12600)

Functions
void	bl1_sinvertv (conj1_t conj, int n, float *x, int incx)
void	bl1_dinvertv (conj1_t conj, int n, double *x, int incx)
void	bl1_cinvertv (conj1_t conj, int n, scomplex *x, int incx)
void	bl1_zinvertv (conj1_t conj, int n, dcomplex *x, int incx)

---

Function Documentation

void bl1_cinvertv	(	conj1_t	conj,
		int	n,
		scomplex *	x,
		int	incx
	)

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

Referenced by FLA_Invert().

{
    float     one = 1.0F;
    float     temp;
    float     s, xr_s, xi_s;
    float     conjsign;
    scomplex* chi;
    int       i;

    if ( bl1_is_conj( conj ) ) conjsign =  one;
    else                       conjsign = -one;

    for ( i = 0; i < n; ++i )
    {
        chi = x + i*incx;

        s         = bl1_fmaxabs( chi->real, chi->imag ); \
        xr_s      = chi->real / s;
        xi_s      = chi->imag / s;
        temp      = xr_s * chi->real + xi_s * chi->imag;

        chi->real =            xr_s / temp;
        chi->imag = conjsign * xi_s / temp;
    }
}

void bl1_dinvertv	(	conj1_t	conj,
		int	n,
		double *	x,
		int	incx
	)

Referenced by FLA_Invert().

{
    double  one = 1.0;
    double* chi;
    int     i;

    for ( i = 0; i < n; ++i )
    {
        chi = x + i*incx;

        *chi = one / *chi;
    }
}

void bl1_sinvertv	(	conj1_t	conj,
		int	n,
		float *	x,
		int	incx
	)

Referenced by FLA_Invert().

{
    float  one = 1.0F;
    float* chi;
    int    i;

    for ( i = 0; i < n; ++i )
    {
        chi = x + i*incx;

        *chi = one / *chi;
    }
}

void bl1_zinvertv	(	conj1_t	conj,
		int	n,
		dcomplex *	x,
		int	incx
	)

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

Referenced by FLA_Invert().

{
    double    one = 1.0;
    double    temp;
    double    s, xr_s, xi_s;
    double    conjsign;
    dcomplex* chi;
    int       i;

    if ( bl1_is_conj( conj ) ) conjsign =  one;
    else                       conjsign = -one;

    for ( i = 0; i < n; ++i )
    {
        chi = x + i*incx;

        s         = bl1_fmaxabs( chi->real, chi->imag ); \
        xr_s      = chi->real / s;
        xi_s      = chi->imag / s;
        temp      = xr_s * chi->real + xi_s * chi->imag;

        chi->real =            xr_s / temp;
        chi->imag = conjsign * xi_s / temp;
    }
}