应使用 vec3 结构体或 eigen 库实现点积与叉积：vec3 明确三维、重载操作符区分点/叉积、默认零初始化；eigen 提供编译期维度检查和高效接口，避免 std::vector 或 std::valarray 的越界、歧义与精度问题。
用 std::valarray 或手写结构体做点积，别碰 std::vector 直接运算
标准库的 std::vector 没有内置点积或叉积接口，强行遍历累加容易写错索引或漏掉维度检查。用 std::valarray 能靠 operator* 和 sum() 快速实现点积，但要注意它只支持一维且无空间维度语义；更稳妥的是自己封装一个 Vec3 结构体——明确限定为三维、重载 operator* 区分点积与叉积、默认初始化为零向量。
常见错误是把 std::vector<double></double> 当作数学向量传给多个函数反复计算，结果因缺少长度校验导致越界或静默错误。比如两个 std::vector 长度不等时点积应报错，但裸循环可能只算到短的那个长度。
- 点积必须要求两向量维度一致，
Vec3构造时可强制固定为 3，避免运行时判断 -
std::valarray的sum()是double类型，注意精度丢失风险（尤其大量小数累加） - 别用
std::vector::data()+std::inner_product—— 它不检查长度，且需手动传迭代器范围，易出错
叉积只对三维有效，Vec3::cross() 必须硬编码 x/y/z 分量公式
二维向量没有叉积（只有标量伪叉积），四维及以上没有唯一垂直向量定义。C++ 里所有“通用 N 维叉积”实现都是错的或自欺欺人的。真正可用的只有三维叉积，公式固定：v1.y * v2.z - v1.z * v2.y 等三项。必须手写，不能靠模板推导维度。
有人用 std::array<double></double> 加 constexpr 函数封装，可行但调用略啰嗦；struct Vec3 成员函数更直观，也方便后续加归一化、长度等方法。
立即学习“C++免费学习笔记（深入）”；
- 叉积结果方向依赖于坐标系手性（右手系下
i×j=k），代码里不注释清楚会导致物理仿真翻车 - 输入向量共线时结果为零向量，这是正确行为，不是 bug
- 不要尝试用行列式形式（如嵌套
std::array模拟矩阵）——编译慢、无实际收益、还掩盖了三维特异性
Eigen 库的 dot() 和 cross() 是最省心的选择，但得接受它的类型体系
如果项目允许第三方依赖，Eigen::Vector3d 的 dot() 和 cross() 不仅正确、高效，还带编译期维度检查。调用 v1.cross(v2) 时若 v1 不是 Vector3d，直接编译失败，比运行时断言强得多。
代价是引入整个 Eigen，且它的向量类型不能和 std::vector 混用——你得把原始数据拷进 Eigen::Map 或显式构造，中间多一次内存访问。不过对数学密集型逻辑来说，这点开销远小于手写 bug 的调试成本。
-
Eigen::Vector3d::dot()返回double，cross()返回Vector3d，类型清晰，不会误当指针用 - 不支持
std::vector直接转Eigen::VectorX—— 必须用Eigen::Map映射已有内存，否则拷贝 - 头文件模式集成简单，但注意
EIGEN_DONT_VECTORIZE等宏会影响性能，调试时默认关闭向量化
自定义 Vec3 时，operator* 重载必须拆成两个：点积用 double operator*(const Vec3&, const Vec3&)，叉积用 Vec3 Vec3::cross(const Vec3&)
用同一个操作符 * 同时表示点积和叉积是典型歧义设计。C++ 不支持按返回类型重载，所以 v1 * v2 只能是一种含义。要么全用成员函数（推荐），要么点积用全局 operator*，叉积坚决用命名函数——这是避免调用方困惑的底线。
另外，Vec3 的构造函数建议带默认参数：Vec3(double x = 0, double y = 0, double z = 0)，这样 Vec3 v; 就是零向量，不用写 Vec3{0,0,0}。
- 别把叉积写成
operator%或operator^—— 这些不是常规用法，别人读代码第一反应是位运算 - 如果需要支持 float/double 模板，用
template<typename t> struct Vec3</typename>，但先确保 double 版本跑通再泛化 - 所有运算函数标记
constexpr和noexcept，方便编译器优化，也表明无副作用
真正麻烦的从来不是写那几行公式，而是统一团队对“向量”的理解：它到底是一个动态数组，还是一个有明确维度和代数规则的数学对象。一旦混用，调试时连该查内存布局还是查手性都会犹豫。
