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
  */

Compare

Flat