bl1_invertv.c File Reference

(r)

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

bl1_cinvertv()

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

{
    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;
    }
}

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

Referenced by FLA_Invert().

bl1_dinvertv()

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

{
    double  one = 1.0;
    double* chi;
    int     i;
 
    for ( i = 0; i < n; ++i )
    {
        chi = x + i*incx;
 
        *chi = one / *chi;
    }
}

References i.

Referenced by FLA_Invert().

bl1_sinvertv()

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

{
    float  one = 1.0F;
    float* chi;
    int    i;
 
    for ( i = 0; i < n; ++i )
    {
        chi = x + i*incx;
 
        *chi = one / *chi;
    }
}

References i.

Referenced by FLA_Invert().

bl1_zinvertv()

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

{
    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;
    }
}

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

Referenced by FLA_Invert().