Compare
Compare
按指定的精度比较两个矩阵/向量的元素。
ulong vector::Compare(
const vector& vec, // 欲比较的向量
const double epsilon // 精度
);
ulong matrix::Compare(
const matrix& mat, // 欲比较的矩阵
const double epsilon // 精度
);参数
- vector_b
[输入] 欲比较的向量。
- epsilon
[输入] 精度。
返回值
所比较的矩阵或向量的不匹配元素数量:如果矩阵相等,则为 0,否则大于 0。
注意
比较运算符 == 或 != 执行精确的元素级比较。 众所周知,实数的精确比较用途有限,因此增加了 epsilon 比较方法。 一个矩阵也许包含某个范围内的元素,例如从 1e-20 到 1e+20。 这些矩阵可逐元素比较来进行处理,直到有效数字。
对于复数矩阵/向量,该比较包括估计复数之间的距离。距离的计算公式为sqrt(pow(r1-r2, 2) + pow(i1-i2, 2),并且是一个可以与epsilon进行比较的实数。
举例
matrix matrix_a={{10,3,2},{1,8,12},{6,5,4}};
matrix matrix_i=matrix::Identity(3,3);
matrix matrix_c=matrix_a.Inv();
matrix matrix_check=matrix_a.MatMul(matrix_c);
Print("matrix_check\n",matrix_check);
ulong errors=matrix_check.Compare(matrix::Identity(3,3),1e-15);
Print("errors=",errors);
/*
matrix_check
[[1,0,0]
[4.440892098500626e-16,1,8.881784197001252e-16]
[4.440892098500626e-16,2.220446049250313e-16,0.9999999999999996]]
errors=0
*/