单元细分源代码和线性方程组求解器
本帖最后由 Esimulate 于 2011-3-13 17:10 编辑Fortran的程序，完全原创，2D和3D的单元，生成Tecplot格式的输出文件用于检查细分是否正确；主要用于扩展有限元编程，需要通过裂纹面把单元细分。代码开源，期待能够共同进步。
刚来论坛，希望这第一帖对大家有用。联系邮箱：Esimulate@gmail.com 在网上能找到Suku的类似代码，不过他是用C++写的。 开源第二贴
最近写程序，需要一个高效简单的线性方程组求解器，就在别人的基础上改造了一套代码。
这套有两个版本:fortran/C++，主要用于求解稀疏对称矩阵；
本代码完全开源，如果您觉得有用，并要加以利用，请告诉大家该代码来源于SimWe仿真论坛的Esimulate；如果发现有bug，或者愿意共享您的改进版本，或者有其他的需求都可与我们联系； C++ version 来源于网络。 fortran version 由Esimulate改写。 fortran version 的输入文件为C++ version\LDL\Matrix 中任意文件，名字固定为A01。
联系邮箱：Esimulate@gmail.com 谢谢！下来试用一下 真是杰出的工作，非常感谢！下来试用一下！ 开源第二贴
最近写程序，需要一个高效简单的线性方程组求解器，就在别人的基础上改造了一套代码。
这套有两个版本:fortran/C++，主要用于求解稀疏对称矩阵；
本代码完全开源，如果您觉得有用，并要加以利用，请告 ...
Esimulate 发表于 2011-3-11 21:13 http://forum.simwe.com/images/common/back.gif
请问lz，Fortran版的稀疏实对称矩阵求解器采用了什么算法呢？
效率方面和pardiso、UMFPACK、MA47、MA57等相比，优势明显吗？ 6# pasuka 6# pasuka 算法还是LDL的三角分解 但是和一般书上的LDL相比，其最大的优势是节省了大量的存储，你可以发现的他的稀疏矩阵存储格式是CSC的。当然我也没有和那些大的开源求解器做过对比，不过这个求解器真的很精致，更方便大家利用和再开发 6# pasuka 算法还是LDL的三角分解 但是和一般书上的LDL相比，其最大的优势是节省了大量的存储，你可以发现的他的稀疏矩阵存储格式是CSC的。当然我也没有和那些大的开源求解器做过对比，不过这个求解器真的很精致，更 ...
Esimulate 发表于 2011-3-15 17:28 http://forum.simwe.com/images/common/back.gif
呵呵，这个见仁见智吧
即便是矩阵三角分解，也有若干种不同方法 自己顶，希望大家评价 希望大家多多给意见 建议试用一下MUMPS parallel sparse direct solver,with GotoBLAS library...
请问lz，Fortran版的稀疏实对称矩阵求解器采用了什么算法呢？
效率方面和pardiso、UMFPACK、MA47、MA57等相比，优势明显吗？
pasuka 发表于 2011-3-14 16:19 http://forum.simwe.com/images/common/back.gif MUMPS用过。太大了，后来不想用了。而且windows下编不出来。
看到这个程序比较精致，就尝试了一下。
12# iomega 建议试用一下MUMPS parallel sparse direct solver,with GotoBLAS library...
iomega 发表于 2011-4-9 14:40 http://forum.simwe.com/images/common/back.gif
用下来感觉还是pardiso使用起来方便，装ivf专业版就自带了
MA47、57、87系列的话，不给带注释源码版本，只能是一个灰箱，不过小巧、方便
国产的GSS没有使用过，据称性能超过pardiso I don't have problem of compiling MUMPS under windows. I was using microsoft C compiler cl and g95.
13# Esimulate 好贴，热烈欢迎 求解器的效率应该对比下看看到底怎么样。 谢谢分享
页:
[1]