处理线性代数的最佳基类型

我正在用C++为一个项目编写一个小而且不充分的线性代数库(对不起).我正在使用双精度数来实现矩阵和运算.我做得对吗？我应该实现模板类吗？周围有更精确的类型吗？

我将使用模板实现类/结构.一开始,你很可能会对just感到满意double,但我发现在我没有将矩阵实现为模板的每个项目中,我后来后悔了.

此外,它为您提供了使用更有趣的元素代数的机会 - 区间算法,概率分布,复杂数学,定点匹配,子矩阵,简单数学:-)等.