nalgebra中如何正确计算向量内积得到f64标量结果？

nalgebra的矩阵运算遵循严格的维度类型规则：行向量（转置后的列向量本质是1×N的矩阵）和列向量（N×1的矩阵）做通用矩阵乘法时，返回结果维度为1×1，类型保持为OMatrix，不会自动隐式转换为f64这类原生标量类型，因此直接用转置相乘的结果参与标量和矩阵的乘法运算时，会触发类型不匹配的编译错误。

优先使用库内置的专用内积API，其次可以对1×1矩阵手动提取标量，两种方式都能满足需求：

• 方案1：使用向量内置的.dot()方法（最推荐，语义清晰且性能最优）
  .dot()是nalgebra专门为向量内积（点积）设计的方法，直接返回对应数值类型的标量结果，不会产生1×1矩阵，内部也有针对点积的专属性能优化。
  示例代码：

  use nalgebra::{OVector, OMatrix, U2};
  
  let m = OMatrix::<f64, U2, U2>::from_element(1.0);
  let v = OVector::<f64, U2>::from_column_slice(&[3_f64, 4_f64]);
  
  let inner_product = v.dot(&v); // 直接得到f64类型的25.0
  let y = inner_product * m; // 编译通过，运算结果和25.0 * m完全一致
  

• 方案2：从1×1矩阵结果中手动提取标量
  如果已经通过矩阵乘法得到了1×1的结果，可以通过索引或者专属快捷属性取出内部存储的标量值：

  use nalgebra::{OVector, OMatrix, U2};
  
  let v = OVector::<f64, U2>::from_column_slice(&[3_f64, 4_f64]);
  let x_1x1_mat = &v.transpose() * &v;
  let scalar_x = x_1x1_mat[(0, 0)]; // 也可以用x_1x1_mat.x，1×1矩阵专属的取值快捷属性
  let y = scalar_x * OMatrix::<f64, U2, U2>::from_element(1.0); // 正常编译运行
  

注意：不建议在向量内积场景下使用转置相乘的写法，除了类型不符合预期外，通用矩阵乘法的执行路径性能也略低于专用的.dot()方法。

内容的提问来源于stack exchange，提问作者xc wang