求解线性方程组

aliouying

如果单纯给代码，咱们没法给您看，没有任何说明，没有任何算法原理

最好给出输入文件，别人才好测试你的程序是否正确

aliouying

麦田守望者 发表于 2014-4-17 19:55
恩，第一次发帖问问题，见谅哦，整体的代码有几百行，发上来也没人愿意看吧，这个应该是一个单独的方程组 ...

我说的输入文件的意思是指矩阵，不是其他，不要拿专业的东西，能碰到同一个专业且热心帮你解决的可能性很小，应该从纯数学或者fortran程序的角度上描述问题

aliouying

麦田守望者 发表于 2014-4-18 23:04
嗯，还是谢谢你了，那样的话我还得考虑下怎么描述

如果这个程序是以前别人写的，而且是对的，那么肯定是你的矩阵组装的时候或者单元矩阵有问题
一般解方程组的这些代码是公用的，经过测试了的

如果非要检查，可以把矩阵的行号、列号、非零值 输出来，然后用matlab或者其他的求解器来对比，看求解方程组是否有问题，仅仅给定这个代码，没有算法，那是真心不好检查，只能自己猜测的事情没人愿意干。