CompareByDigits
CompareByDigits
按有效位数精度比较两个矩阵/向量的元素。
ulong vector::CompareByDigits(
const vector& vec, // 欲比较的向量
const int digits // 有效位数
);
ulong matrix::CompareByDigits(
const matrix& mat, // 欲比较的矩阵
const int digits // 有效位数
);参数
- vector_b
[输入] 欲比较的向量。
- digits
[输入] 进行比较采用的有效位数。
- epsilon
[输入] 比较精度。 如果两个值的绝对值相差小于指定的精度,则认为它们相等。
返回值
所比较的矩阵或向量的不匹配元素数量:如果矩阵相等,则为 0,否则大于 0。
注意
比较运算符 == 或 != 执行精确的元素级比较。 众所周知,实数的精确比较用途有限,因此增加了 epsilon 比较方法。 一个矩阵也许包含某个范围内的元素,例如从 1e-20 到 1e+20。 此类矩阵可用元素级比较来处理,直至有效位数。
举例
int size_m=128;
int size_k=256;
matrix matrix_a(size_m,size_k);
//--- 填充矩阵
double value=0.0;
for(int i=0; i<size_m; i++)
{
for(int j=0; j<size_k; j++)
{
if(i==j)
matrix_a[i][j]=1.0+i;
else
{
value+=1.0;
matrix_a[i][j]=value/1e+20;
}
}
}
//--- 获取其它矩阵
matrix matrix_c = matrix_a * -1;
ulong errors_epsilon=matrix_a.Compare(matrix_c,1e-15);
ulong errors_digits=matrix_a.CompareByDigits(matrix_c,15);
printf("Compare matrix %d x %d errors_epsilon=%I64u errors_digits=%I64u",size_m,size_k,errors_epsilon,errors_digits);
/*
比较 矩阵 128 x 256 errors_epsilon=128 errors_digits=32768
*/