bl1_invert2s.c File Reference

(r)

Functions
void	bl1_sinvert2s (conj1_t conj, float *alpha, float *beta)
void	bl1_dinvert2s (conj1_t conj, double *alpha, double *beta)
void	bl1_cinvert2s (conj1_t conj, scomplex *alpha, scomplex *beta)
void	bl1_zinvert2s (conj1_t conj, dcomplex *alpha, dcomplex *beta)

Function Documentation

bl1_cinvert2s()

void bl1_cinvert2s	(	conj1_t	conj,
		scomplex *	alpha,
		scomplex *	beta )

{
    float  temp;
    float  s, xr_s, xi_s;
 
    s           = bl1_fmaxabs( alpha->real, alpha->imag ); \
    xr_s        = alpha->real / s;
    xi_s        = alpha->imag / s;
    temp        = xr_s * alpha->real + xi_s * alpha->imag;
 
    beta->real =  xr_s / temp;
    beta->imag = -xi_s / temp;
 
    if ( bl1_is_conj( conj ) )
        bl1_cconjs( beta );
}

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

Referenced by bl1_cinvscalm(), and bl1_cinvscalv().

bl1_dinvert2s()

void bl1_dinvert2s	(	conj1_t	conj,
		double *	alpha,
		double *	beta )

{
    double one = 1.0;
 
    *beta = one / *alpha;
}

Referenced by bl1_dinvscalm(), and bl1_zdinvscalm().

bl1_sinvert2s()

void bl1_sinvert2s	(	conj1_t	conj,
		float *	alpha,
		float *	beta )

{
    float  one = 1.0F;
 
    *beta = one / *alpha;
}

Referenced by bl1_csinvscalm(), and bl1_sinvscalm().

bl1_zinvert2s()

void bl1_zinvert2s	(	conj1_t	conj,
		dcomplex *	alpha,
		dcomplex *	beta )

{
    double temp;
    double s, xr_s, xi_s;
 
    s           = bl1_fmaxabs( alpha->real, alpha->imag ); \
    xr_s        = alpha->real / s;
    xi_s        = alpha->imag / s;
    temp        = xr_s * alpha->real + xi_s * alpha->imag;
 
    beta->real =  xr_s / temp;
    beta->imag = -xi_s / temp;
 
    if ( bl1_is_conj( conj ) )
        bl1_zconjs( beta );
}

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

Referenced by bl1_zinvscalm(), and bl1_zinvscalv().