40 float *buff_a = (
float * ) FLA_FLOAT_PTR( A );
41 float *buff_b = (
float * ) FLA_FLOAT_PTR( B );
43 for( j = 0; j < n_A; j++ )
45 for( i = 0; i < m_A; i++ )
47 diff = ( double ) ( buff_a[ j*cs_A + i*rs_A ] - buff_b[ j*cs_B + i*rs_B ] );
49 if( fabs(diff) > d_max )
59 double *buff_a = (
double * ) FLA_DOUBLE_PTR( A );
60 double *buff_b = (
double * ) FLA_DOUBLE_PTR( B );
62 for( j = 0; j < n_A; j++ )
64 for( i = 0; i < m_A; i++ )
66 diff = ( double ) ( buff_a[ j*cs_A + i*rs_A ] - buff_b[ j*cs_B + i*rs_B ] );
68 if( fabs(diff) > d_max )
81 for( j = 0; j < n_A; j++ )
83 for( i = 0; i < m_A; i++ )
85 diff = ( double ) ( buff_a[ j*cs_A + i*rs_A ].
real - buff_b[ j*cs_B + i*rs_B ].
real );
87 if( fabs(diff) > d_max )
90 diff = ( double ) ( buff_a[ j*cs_A + i*rs_A ].
imag - buff_b[ j*cs_B + i*rs_B ].
imag );
92 if( fabs(diff) > d_max )
100 case FLA_DOUBLE_COMPLEX:
105 for( j = 0; j < n_A; j++ )
107 for( i = 0; i < m_A; i++ )
109 diff = ( double ) ( buff_a[ j*cs_A + i*rs_A ].real - buff_b[ j*cs_B + i*rs_B ].real );
111 if( fabs(diff) > d_max )
114 diff = ( double ) ( buff_a[ j*cs_A + i*rs_A ].imag - buff_b[ j*cs_B + i*rs_B ].imag );
116 if( fabs(diff) > d_max )
unsigned long dim_t
Definition: FLA_type_defs.h:71
dim_t FLA_Obj_row_stride(FLA_Obj obj)
Definition: FLA_Query.c:167
float real
Definition: FLA_f2c.h:30
FLA_Error FLA_Max_elemwise_diff_check(FLA_Obj A, FLA_Obj B)
Definition: FLA_Max_elemwise_diff_check.c:13
FLA_Datatype FLA_Obj_datatype(FLA_Obj obj)
Definition: FLA_Query.c:13
rho_c imag
Definition: bl1_axpyv2bdotaxpy.c:483
dim_t FLA_Obj_width(FLA_Obj obj)
Definition: FLA_Query.c:123
Definition: blis_type_defs.h:132
unsigned int FLA_Check_error_level(void)
Definition: FLA_Check.c:18
int FLA_Datatype
Definition: FLA_type_defs.h:49
dim_t FLA_Obj_col_stride(FLA_Obj obj)
Definition: FLA_Query.c:174
int i
Definition: bl1_axmyv2.c:145
dim_t FLA_Obj_length(FLA_Obj obj)
Definition: FLA_Query.c:116
Definition: blis_type_defs.h:137