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。这些矩阵可逐元素比较来进行处理，直到有效数字。

举例

   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



  */