标签： performance fortran intel-mkl matrix-inverse
我在Fortran 90中有一段代码，我必须解决非线性（我必须反转雅可比矩阵）和线性方程组。当我说小时我的意思是 n 两个操作的未知数， n <= 4 。遗憾的是， n 未知 a priori 。你认为最好的选择是什么？
我想为 n = 1,2 并使用 n = 3,4 的其他方法编写显式公式（例如英特尔MKL库的某些函数） ），为了表现。这是明智的，还是应该为 n = 3,4 的逆矩阵编写显式公式？