Fortran Coder

查看: 3718|回复: 9
打印 上一主题 下一主题

[非线性] Fortran程序中求解非线性方程组

[复制链接]

6

帖子

1

主题

0

精华

新人

F 币
28 元
贡献
13 点
跳转到指定楼层
楼主
发表于 2022-6-10 22:46:51 | 只看该作者 回帖奖励 |倒序浏览 |阅读模式
请教各位大佬,Fortran程序中有没有现成的求解非线性方程组的子程序,刚开始接触,需要联立求解两个体积守恒方程,因此需要用到
分享到:  微信微信
收藏收藏 点赞点赞 点踩点踩

44

帖子

4

主题

0

精华

熟手

F 币
179 元
贡献
90 点
沙发
发表于 2022-6-11 12:44:06 | 只看该作者
PETSC 可以并行的开源求解器,线性非线性都有,你说的体积守恒方程也太不清楚了。不知道你到底要什么样的求解器。

44

帖子

4

主题

0

精华

熟手

F 币
179 元
贡献
90 点
板凳
发表于 2022-6-11 12:49:28 | 只看该作者
intel 的MKL应该也都够用了吧,

6

帖子

1

主题

0

精华

新人

F 币
28 元
贡献
13 点
地板
 楼主| 发表于 2022-6-13 16:45:19 | 只看该作者
zjk0112 发表于 2022-6-11 12:49
intel 的MKL应该也都够用了吧,

用过这个,但一方面他经常报错,或者给不出结果(可能我的公式比较复杂吧,方程本身是个积分),另一方面,有时候需要在linux上运行,所以想的是不用函数库,这些功能尽量用程序本身的子程序或者函数实现

6

帖子

1

主题

0

精华

新人

F 币
28 元
贡献
13 点
5#
 楼主| 发表于 2022-6-13 16:55:36 | 只看该作者
zjk0112 发表于 2022-6-11 12:44
PETSC 可以并行的开源求解器,线性非线性都有,你说的体积守恒方程也太不清楚了。不知道你到底要什么样的求 ...

谢谢您的推荐,倒也不是求解器,就是需要个子程序,至于体积守恒方程就是,体积已知,它等于个积分,未知数存在于积分中

59

帖子

2

主题

0

精华

大师

F 币
810 元
贡献
476 点
6#
发表于 2022-6-15 06:29:45 | 只看该作者
你是积分方程组,还是非线性方程组,你的积分是固定区间积分还是积分限是变量,你这种笼统问,得不到什么有用信息

44

帖子

4

主题

0

精华

熟手

F 币
179 元
贡献
90 点
7#
发表于 2022-6-15 20:19:06 | 只看该作者
hanch 发表于 2022-6-13 16:55
谢谢您的推荐,倒也不是求解器,就是需要个子程序,至于体积守恒方程就是,体积已知,它等于个积分,未知 ...

感觉你人有点糊里糊涂的。你到底想要啥啊

6

帖子

1

主题

0

精华

新人

F 币
28 元
贡献
13 点
8#
 楼主| 发表于 2022-9-28 16:09:22 | 只看该作者
kerb 发表于 2022-6-15 06:29
你是积分方程组,还是非线性方程组,你的积分是固定区间积分还是积分限是变量,你这种笼统问,得不到什么有 ...

两个方程,都是先对一个变量积分,积分后的结果都是关于另外两个变量的方程,联立解方程组。

6

帖子

1

主题

0

精华

新人

F 币
28 元
贡献
13 点
9#
 楼主| 发表于 2022-9-28 16:10:56 | 只看该作者
hanch 发表于 2022-9-28 16:09
两个方程,都是先对一个变量积分,积分后的结果都是关于另外两个变量的方程,联立解方程组。 ...

积分限和需要积分的表达式都有所求变量

1

帖子

0

主题

0

精华

新人

F 币
13 元
贡献
2 点
10#
发表于 2023-4-11 14:56:41 | 只看该作者
我也遇到此类非线性的积分方程,用MKL库求解总是不如人意,没能正常迭代,最后一气之下自己写了求解程序,虽然重造轮子,但自己写的用起来比较得心应手
您需要登录后才可以回帖 登录 | 极速注册

本版积分规则

捐赠本站|Archiver|关于我们 About Us|小黑屋|Fcode ( 京ICP备18005632-2号 )

GMT+8, 2024-12-23 21:49

Powered by Tencent X3.4

© 2013-2024 Tencent

快速回复 返回顶部 返回列表