应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » 其他 » Fortran程序，如何实现采用矩阵保存单链表表头？
71/1返回列表
查看: 1660  |  回复: 6
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

andersonlyw

新虫 (著名写手)

[求助] Fortran程序，如何实现采用矩阵保存单链表表头？

问题是这样，小弟想做三维网格的DSMC稀薄气体模拟，打算：
用N*N*N的三维矩阵存贮N*N*N个单链表表头，这样每个单链表可以用来存贮该网格内所有的气体分子信息（速度，坐标等），但是不知如何用Fortran程序实现？
请高手指教！可否提供几行程序片段供参考？

相应的物理思想见王学德，发表在2007年空气动力学学报上的论文：
http://d.g.wanfangdata.com.hk.ez ... dlxxb200701020.aspx
题名：一类二维非结构网格DSMC方法的实现策略及其应用
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)
1楼 2013-10-12 12:57:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

andersonlyw

新虫 (著名写手)
解决此问题的虫友，另有100金币可领。见：
http://muchong.com/bbs/viewthread.php?tid=6451229
赞一下 回复此楼
2楼2013-10-12 17:48:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

andersonlyw

新虫 (著名写手)
我似乎已经想到了一种方式。
如下：采用一个指针数组MP（分配三维N×N×N空间），一个单链表P，一个普通矩阵M（三维N×N×N）可完成任务。首先，单链表P记录所有计算区（三维N×N×N　Cells）内的分子信息，分子运动一步后，计算碰撞时，先记录每Cell内的分子数量于普通矩阵M内。然后根据矩阵M提供的信息，给指针数组每个点（共N×N×N个）分配（６×Ｘ的二维存贮空间，其中3行存分子速度，3行存贮分子坐标）。这类似于动态矩阵，完成后，在计算分子碰撞时，可以直接采用矩阵方式选择碰撞分子对。完成分子碰撞后，再将指针数组内保存的数据全部转存至单链表。如此循环即可。
赞一下 回复此楼
3楼2013-10-13 12:34:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

andersonlyw

新虫 (著名写手)
请教各位做稀薄气体DSMC模拟的虫友们，你们是怎么存贮分子信息的？互相交流一下，共同提高。有金币送。
赞一下 回复此楼
4楼2013-10-13 12:36:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

andersonlyw

新虫 (著名写手)
引用回帖:
3楼: Originally posted by andersonlyw at 2013-10-13 12:34:37
我似乎已经想到了一种方式。
如下：采用一个指针数组MP（分配三维N×N×N空间），一个单链表P，一个普通矩阵M（三维N×N×N）可完成任务。首先，单链表P记录所有计算区（三维N×N×N　Cells）内的分子信息，分子运 ...

最近采用此方法编写了程序，完成了分子泵的过渡流态下的压缩比计算，确认该方法可行，占用内存少，可放心使用。
赞一下 回复此楼
5楼2013-11-27 08:38:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

meigudengfei

新虫 (初入文坛)
楼主，我做的是粘性颗粒的碰撞，最近一直在纠结碰撞信息的存储，想问问运用指针线性链表怎么处理的，还有你是怎么搜索分子颗粒碰撞的，判断分子颗粒之间的距离吗？求教。。。
赞一下 回复此楼
6楼2014-03-23 20:51:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zzpjjjjj

禁虫 (初入文坛)
本帖内容被屏蔽
7楼2016-10-20 08:48:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 andersonlyw 的主题更新
71/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 上海工程技术大学 激光智能制造课题组 2027级博士研究生招生 +5 oxidpl 2026-09-20 7/350 2026-09-27 10:43 by Equinoxhua
[博后之家] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 ViqFlLSxHDGt 2026-09-26 3/150 2026-09-27 10:42 by 4IW0sJvtEqX8
[找工作] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3 ViqFlLSxHDGt 2026-09-26 3/150 2026-09-27 10:26 by 4IW0sJvtEqX8
[博后之家] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 ViqFlLSxHDGt 2026-09-26 3/150 2026-09-27 10:21 by 4IW0sJvtEqX8
[考博] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 ViqFlLSxHDGt 2026-09-26 3/150 2026-09-27 10:03 by 4IW0sJvtEqX8
[考博] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 oeyPlfmMUpOC 2026-09-26 3/150 2026-09-27 09:13 by 4IW0sJvtEqX8
[公派出国] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 GxESXJmptOtQ 2026-09-25 4/200 2026-09-27 06:44 by XgCC8uTcwILl
[考博] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 GxESXJmptOtQ 2026-09-25 4/200 2026-09-27 06:42 by XgCC8uTcwILl
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 GxESXJmptOtQ 2026-09-25 3/150 2026-09-27 06:32 by XgCC8uTcwILl
[找工作] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 vQZoDrrm7mUF 2026-09-25 3/150 2026-09-27 06:08 by XgCC8uTcwILl
[公派出国] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 vQZoDrrm7mUF 2026-09-25 3/150 2026-09-27 06:05 by XgCC8uTcwILl
[考研] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 vQZoDrrm7mUF 2026-09-25 3/150 2026-09-27 05:56 by XgCC8uTcwILl
[论文投稿] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +4 vQZoDrrm7mUF 2026-09-25 4/200 2026-09-27 05:56 by XgCC8uTcwILl
[找工作] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +5 vQZoDrrm7mUF 2026-09-25 5/250 2026-09-27 05:53 by XgCC8uTcwILl
[考研] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3 vQZoDrrm7mUF 2026-09-25 5/250 2026-09-27 05:44 by XgCC8uTcwILl
[公派出国] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +4 vQZoDrrm7mUF 2026-09-25 4/200 2026-09-27 05:44 by XgCC8uTcwILl
[考博] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 vQZoDrrm7mUF 2026-09-25 3/150 2026-09-27 05:44 by XgCC8uTcwILl
[论文投稿] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3 vQZoDrrm7mUF 2026-09-25 6/300 2026-09-27 05:42 by XgCC8uTcwILl
[找工作] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +4 PGTSwSC3F6nQ 2026-09-24 5/250 2026-09-25 11:19 by yhq9807
[基金申请] 师弟论文见刊大半年才想起来申请专利，还能抢救一下吗？ +4 13108017953 2026-09-21 4/200 2026-09-24 10:06 by kudofaye
信息提示
关闭
请填处理意见
关闭