FLA_Bsvd_find_max_min.c File Reference

(r)

Functions
FLA_Error	FLA_Bsvd_find_max (FLA_Obj d, FLA_Obj e, FLA_Obj smax, FLA_Obj smin)
FLA_Error	FLA_Bsvd_find_max_min_ops (int m_A, float *buff_d, int inc_d, float *buff_e, int inc_e, float *smax, float *smin)
FLA_Error	FLA_Bsvd_find_max_min_opd (int m_A, double *buff_d, int inc_d, double *buff_e, int inc_e, double *smax, double *smin)

Function Documentation

FLA_Bsvd_find_max()

FLA_Error FLA_Bsvd_find_max	(	FLA_Obj	d,
		FLA_Obj	e,
		FLA_Obj	smax,
		FLA_Obj	smin
	)

{
    FLA_Datatype datatype;
    int          m_A;
    int          inc_d;
    int          inc_e;
 
    datatype = FLA_Obj_datatype( d );
 
    m_A      = FLA_Obj_vector_dim( d );
 
    inc_d    = FLA_Obj_vector_inc( d );
    inc_e    = FLA_Obj_vector_inc( e );
 
 
    switch ( datatype )
    {
    case FLA_FLOAT:
    {
        float*    buff_d    = FLA_FLOAT_PTR( d );
        float*    buff_e    = FLA_FLOAT_PTR( e );
        float*    buff_smax = FLA_FLOAT_PTR( smax );
        float*    buff_smin = FLA_FLOAT_PTR( smin );
 
        FLA_Bsvd_find_max_min_ops( m_A,
                                   buff_d, inc_d,
                                   buff_e, inc_e,
                                   buff_smax,
                                   buff_smin );
 
        break;
    }
 
    case FLA_DOUBLE:
    {
        double*   buff_d    = FLA_DOUBLE_PTR( d );
        double*   buff_e    = FLA_DOUBLE_PTR( e );
        double*   buff_smax = FLA_DOUBLE_PTR( smax );
        double*   buff_smin = FLA_DOUBLE_PTR( smin );
 
        FLA_Bsvd_find_max_min_opd( m_A,
                                   buff_d, inc_d,
                                   buff_e, inc_e,
                                   buff_smax,
                                   buff_smin );
 
        break;
    }
    }
 
    return FLA_SUCCESS;
}

References FLA_Bsvd_find_max_min_opd(), FLA_Bsvd_find_max_min_ops(), FLA_Obj_datatype(), FLA_Obj_vector_dim(), FLA_Obj_vector_inc(), and i.

FLA_Bsvd_find_max_min_opd()

FLA_Error FLA_Bsvd_find_max_min_opd	(	int	m_A,
		double *	buff_d,
		int	inc_d,
		double *	buff_e,
		int	inc_e,
		double *	smax,
		double *	smin
	)

{
    double smax_cand;
    double smin_cand;
    int    i;
 
    smax_cand = fabs( buff_d[ (m_A-1)*inc_d ] );
    smin_cand = smax_cand;
 
    for ( i = 0; i < m_A - 1; ++i )
    {
        double abs_di = fabs( buff_d[ i*inc_d ] );
        double abs_ei = fabs( buff_e[ i*inc_e ] );
 
        // Track the minimum element.
        smin_cand = min( smin_cand, abs_di );
 
        // Track the maximum element.
        smax_cand = max( smax_cand, abs_di );
        smax_cand = max( smax_cand, abs_ei );
    }
 
    // Save the results of the search.
    *smax = smax_cand;
    *smin = smin_cand;
 
    return FLA_SUCCESS;
}

References i.

Referenced by FLA_Bsvd_find_max(), and FLA_Bsvd_sinval_v_opd_var1().

FLA_Bsvd_find_max_min_ops()

FLA_Error FLA_Bsvd_find_max_min_ops	(	int	m_A,
		float *	buff_d,
		int	inc_d,
		float *	buff_e,
		int	inc_e,
		float *	smax,
		float *	smin
	)

{
    float  smax_cand;
    float  smin_cand;
    int    i;
 
    smax_cand = fabsf( buff_d[ (m_A-1)*inc_d ] );
    smin_cand = smax_cand;
 
    for ( i = 0; i < m_A - 1; ++i )
    {
        float abs_di = fabsf( buff_d[ i*inc_d ] );
        float abs_ei = fabsf( buff_e[ i*inc_e ] );
 
        // Track the minimum element.
        smin_cand = min( smin_cand, abs_di );
 
        // Track the maximum element.
        smax_cand = max( smax_cand, abs_di );
        smax_cand = max( smax_cand, abs_ei );
    }
 
    // Save the results of the search.
    *smax = smax_cand;
    *smin = smin_cand;
 
    return FLA_SUCCESS;
}

References i.

Referenced by FLA_Bsvd_find_max(), and FLA_Bsvd_sinval_v_ops_var1().