void libblis_test_subv_check( obj_t* alpha,
obj_t* beta,
obj_t* x,
obj_t* y,
double* resid )
{
num_t dt = bli_obj_datatype( *x );
num_t dt_real = bli_obj_datatype_proj_to_real( *x );
dim_t m = bli_obj_vector_dim( *x );
conj_t conjx = bli_obj_conj_status( *x );
obj_t aminusb;
obj_t alpha_conj;
obj_t norm_r, m_r, temp_r;
double junk;
//
// Pre-conditions:
// - x is set to alpha.
// - y_orig is set to beta.
// Note:
// - alpha and beta should have non-zero imaginary components in the
// complex cases in order to more fully exercise the implementation.
//
// Under these conditions, we assume that the implementation for
//
// y := y_orig - conjx(x)
//
// is functioning correctly if
//
// normfv(y) - sqrt( absqsc( beta - conjx(alpha) ) * m )
//
// is negligible.
//
bli_obj_scalar_init_detached( dt, &aminusb );
bli_obj_scalar_init_detached( dt_real, &temp_r );
bli_obj_scalar_init_detached( dt_real, &norm_r );
bli_obj_scalar_init_detached( dt_real, &m_r );
bli_obj_scalar_init_detached_copy_of( dt, conjx, alpha, &alpha_conj );
bli_normfv( y, &norm_r );
bli_copysc( beta, &aminusb );
bli_subsc( &alpha_conj, &aminusb );
bli_setsc( ( double )m, 0.0, &m_r );
bli_absqsc( &aminusb, &temp_r );
bli_mulsc( &m_r, &temp_r );
bli_sqrtsc( &temp_r, &temp_r );
bli_subsc( &temp_r, &norm_r );
bli_getsc( &norm_r, resid, &junk );
}
void libblis_test_normfm_check
(
test_params_t* params,
obj_t* beta,
obj_t* x,
obj_t* norm,
double* resid
)
{
num_t dt_real = bli_obj_datatype_proj_to_real( *x );
dim_t m = bli_obj_length( *x );
dim_t n = bli_obj_width( *x );
obj_t m_r, n_r, temp_r;
double junk;
//
// Pre-conditions:
// - x is set to beta.
// Note:
// - beta should have a non-zero imaginary component in the complex
// cases in order to more fully exercise the implementation.
//
// Under these conditions, we assume that the implementation for
//
// norm := normf( x )
//
// is functioning correctly if
//
// norm = sqrt( absqsc( beta ) * m * n )
//
// where m and n are the dimensions of x.
//
bli_obj_scalar_init_detached( dt_real, &temp_r );
bli_obj_scalar_init_detached( dt_real, &m_r );
bli_obj_scalar_init_detached( dt_real, &n_r );
bli_setsc( ( double )m, 0.0, &m_r );
bli_setsc( ( double )n, 0.0, &n_r );
bli_absqsc( beta, &temp_r );
bli_mulsc( &m_r, &temp_r );
bli_mulsc( &n_r, &temp_r );
bli_sqrtsc( &temp_r, &temp_r );
bli_subsc( &temp_r, norm );
bli_getsc( norm, resid, &junk );
}
void libblis_test_fnormv_check( obj_t* beta,
obj_t* x,
obj_t* norm,
double* resid )
{
num_t dt_real = bli_obj_datatype_proj_to_real( *x );
dim_t m = bli_obj_vector_dim( *x );
obj_t m_r, temp_r;
double junk;
//
// Pre-conditions:
// - x is set to beta.
// Note:
// - beta should have a non-zero imaginary component in the complex
// cases in order to more fully exercise the implementation.
//
// Under these conditions, we assume that the implementation for
//
// norm := fnorm( x )
//
// is functioning correctly if
//
// norm = sqrt( absqsc( beta ) * m )
//
// where m is the length of x.
//
bli_obj_init_scalar( dt_real, &temp_r );
bli_obj_init_scalar( dt_real, &m_r );
bli_setsc( ( double )m, 0.0, &m_r );
bli_absqsc( beta, &temp_r );
bli_mulsc( &m_r, &temp_r );
bli_sqrtsc( &temp_r, &temp_r );
bli_subsc( &temp_r, norm );
bli_getsc( norm, resid, &junk );
}
void libblis_test_dotaxpyv_check( obj_t* alpha,
obj_t* xt,
obj_t* x,
obj_t* y,
obj_t* rho,
obj_t* z,
obj_t* z_orig,
double* resid )
{
num_t dt = bli_obj_datatype( *z );
num_t dt_real = bli_obj_datatype_proj_to_real( *z );
dim_t m = bli_obj_vector_dim( *z );
obj_t rho_temp;
obj_t z_temp;
obj_t norm_z;
double resid1, resid2;
double junk;
//
// Pre-conditions:
// - x is randomized.
// - y is randomized.
// - z_orig is randomized.
// - xt is an alias to x.
// Note:
// - alpha should have a non-zero imaginary component in the complex
// cases in order to more fully exercise the implementation.
//
// Under these conditions, we assume that the implementation for
//
// rho := conjxt(x^T) conjy(y)
// z := z_orig + alpha * conjx(x)
//
// is functioning correctly if
//
// ( rho - rho_temp )
//
// and
//
// normf( z - z_temp )
//
// are negligible, where rho_temp and z_temp contain rho and z as
// computed by dotv and axpyv, respectively.
//
bli_obj_scalar_init_detached( dt, &rho_temp );
bli_obj_scalar_init_detached( dt_real, &norm_z );
bli_obj_create( dt, m, 1, 0, 0, &z_temp );
bli_copyv( z_orig, &z_temp );
bli_dotv( xt, y, &rho_temp );
bli_axpyv( alpha, x, &z_temp );
bli_subsc( rho, &rho_temp );
bli_getsc( &rho_temp, &resid1, &junk );
bli_subv( &z_temp, z );
bli_normfv( z, &norm_z );
bli_getsc( &norm_z, &resid2, &junk );
*resid = bli_fmaxabs( resid1, resid2 );
bli_obj_free( &z_temp );
}
