单元细分源代码和线性方程组求解器

Esimulate · 发表于 2011-3-8 16:58:50

本帖最后由 Esimulate 于 2011-3-13 17:10 编辑

Fortran的程序，完全原创，2D和3D的单元，生成Tecplot格式的输出文件用于检查细分是否正确；主要用于扩展有限元编程，需要通过裂纹面把单元细分。代码开源，期待能够共同进步。
刚来论坛，希望这第一帖对大家有用。联系邮箱：Esimulate@gmail.com

bbssbb · 发表于 2011-3-8 19:33:02

在网上能找到Suku的类似代码，不过他是用C++写的。

Esimulate · 发表于 2011-3-11 21:13:57

开源第二贴
最近写程序，需要一个高效简单的线性方程组求解器，就在别人的基础上改造了一套代码。
这套有两个版本:fortran/C++，主要用于求解稀疏对称矩阵；
本代码完全开源，如果您觉得有用，并要加以利用，请告诉大家该代码来源于SimWe仿真论坛的Esimulate；如果发现有bug，或者愿意共享您的改进版本，或者有其他的需求都可与我们联系； C++ version 来源于网络。 fortran version 由Esimulate改写。 fortran version 的输入文件为C++ version\LDL\Matrix 中任意文件，名字固定为A01。
联系邮箱：Esimulate@gmail.com

fanxiaochen · 发表于 2011-3-14 21:05:20

谢谢！下来试用一下

colinmi · 发表于 2011-3-14 23:04:59

真是杰出的工作，非常感谢！下来试用一下！

pasuka · 发表于 2011-3-15 08:19:53

开源第二贴
最近写程序，需要一个高效简单的线性方程组求解器，就在别人的基础上改造了一套代码。
这套有两个版本:fortran/C++，主要用于求解稀疏对称矩阵；
本代码完全开源，如果您觉得有用，并要加以利用，请告 ...
Esimulate 发表于 2011-3-11 21:13

请问lz，Fortran版的稀疏实对称矩阵求解器采用了什么算法呢？
效率方面和pardiso、UMFPACK、MA47、MA57等相比，优势明显吗？

Esimulate · 发表于 2011-3-15 17:25:17

6# pasuka

Esimulate · 发表于 2011-3-15 17:28:51

6# pasuka 算法还是LDL的三角分解但是和一般书上的LDL相比，其最大的优势是节省了大量的存储，你可以发现的他的稀疏矩阵存储格式是CSC的。当然我也没有和那些大的开源求解器做过对比，不过这个求解器真的很精致，更方便大家利用和再开发

pasuka · 发表于 2011-3-16 07:54:06

6# pasuka 算法还是LDL的三角分解但是和一般书上的LDL相比，其最大的优势是节省了大量的存储，你可以发现的他的稀疏矩阵存储格式是CSC的。当然我也没有和那些大的开源求解器做过对比，不过这个求解器真的很精致，更 ...
Esimulate 发表于 2011-3-15 17:28

呵呵，这个见仁见智吧
即便是矩阵三角分解，也有若干种不同方法

Esimulate · 发表于 2011-3-23 22:28:15

自己顶，希望大家评价

Esimulate · 发表于 2011-4-5 22:43:25

希望大家多多给意见

iomega · 发表于 2011-4-9 14:40:07

建议试用一下MUMPS parallel sparse direct solver,with GotoBLAS library...

请问lz，Fortran版的稀疏实对称矩阵求解器采用了什么算法呢？
效率方面和pardiso、UMFPACK、MA47、MA57等相比，优势明显吗？
pasuka 发表于 2011-3-14 16:19

Esimulate · 发表于 2011-4-9 19:46:44

MUMPS用过。太大了，后来不想用了。而且windows下编不出来。
看到这个程序比较精致，就尝试了一下。

12# iomega

pasuka · 发表于 2011-4-10 09:02:59

建议试用一下MUMPS parallel sparse direct solver,with GotoBLAS library...

iomega 发表于 2011-4-9 14:40

用下来感觉还是pardiso使用起来方便，装ivf专业版就自带了
MA47、57、87系列的话，不给带注释源码版本，只能是一个灰箱，不过小巧、方便
国产的GSS没有使用过，据称性能超过pardiso

iomega · 发表于 2011-4-10 14:06:22

I don't have problem of compiling MUMPS under windows. I was using microsoft C compiler cl and g95.

13# Esimulate

z03063023 · 发表于 2011-4-10 14:50:35

好贴，热烈欢迎

zgj-chaoyue · 发表于 2011-12-15 12:13:00

求解器的效率应该对比下看看到底怎么样。

S170200232 · 发表于 2017-11-9 18:11:02

谢谢分享

账号		自动登录	找回密码
密码			注册

账号		自动登录	找回密码
密码			立即注册

[H. 有限元编程] 单元细分源代码和线性方程组求解器

本帖子中包含更多资源

评分

相关帖子

本帖子中包含更多资源

评分