h***o 发帖数: 539 | 1 一开始setup一堆网格点,比如1x1的正方形分成10x10个格子,每个格子0.1x0.1
但程序跑着跑着,发现有些地方0.1x0.1的resolution不够,需要调整这些地方
的grid size...这个adaptive grid size的问题在fortran里用啥数据解构来实现呢?
谁有这个方面的经验吗?谢了先 |
a******s 发帖数: 232 | 2 这个我也一直想学来着,你要是都是structure grid的话,这个应该叫
adaptive mesh refinement(AMR).你可以看看这个
http://seesar.lbl.gov/anag/staff/martin/AMRPoisson.html
你如果弄出来了一定教教我阿:)
【在 h***o 的大作中提到】 : 一开始setup一堆网格点,比如1x1的正方形分成10x10个格子,每个格子0.1x0.1 : 但程序跑着跑着,发现有些地方0.1x0.1的resolution不够,需要调整这些地方 : 的grid size...这个adaptive grid size的问题在fortran里用啥数据解构来实现呢? : 谁有这个方面的经验吗?谢了先
|
h***o 发帖数: 539 | 3 我倒是有办法,用链表来实现。哪个grid 附近需要增加网格,我就给
它加一串链表。不过f77里既没有动态分配也没有指针,所以实现
起来比较别扭。
【在 a******s 的大作中提到】 : 这个我也一直想学来着,你要是都是structure grid的话,这个应该叫 : adaptive mesh refinement(AMR).你可以看看这个 : http://seesar.lbl.gov/anag/staff/martin/AMRPoisson.html : 你如果弄出来了一定教教我阿:)
|
a******s 发帖数: 232 | 4 这个不那么容易的,你解的是什么方程?poisson?
【在 h***o 的大作中提到】 : 我倒是有办法,用链表来实现。哪个grid 附近需要增加网格,我就给 : 它加一串链表。不过f77里既没有动态分配也没有指针,所以实现 : 起来比较别扭。
|
h***o 发帖数: 539 | 5 没那么复杂吧。应该跟boltzman equations差不多
【在 a******s 的大作中提到】 : 这个不那么容易的,你解的是什么方程?poisson?
|
c******m 发帖数: 599 | 6 ft, poisson equation is almost the simplest model problem
【在 h***o 的大作中提到】 : 没那么复杂吧。应该跟boltzman equations差不多
|
h***o 发帖数: 539 | 7 我也觉得比poisson的表达式复杂些。hoho
【在 c******m 的大作中提到】 : ft, poisson equation is almost the simplest model problem
|
t*****o 发帖数: 74 | 8 又是boltzman equation
天天听到这个名词boltzman, 纳米, MEMS, Fuel Cell 。。。
F77搞adaptive grid比较困难把
要不混着, C负责网格, f77负责算? HOHO
【在 h***o 的大作中提到】 : 没那么复杂吧。应该跟boltzman equations差不多
|
s*****l 发帖数: 167 | 9 F90 supports pointers
【在 t*****o 的大作中提到】 : 又是boltzman equation : 天天听到这个名词boltzman, 纳米, MEMS, Fuel Cell 。。。 : F77搞adaptive grid比较困难把 : 要不混着, C负责网格, f77负责算? HOHO
|
s*****l 发帖数: 167 | 10 Marsha Berger and Randall Leveque have some papers on this.
【在 h***o 的大作中提到】 : 一开始setup一堆网格点,比如1x1的正方形分成10x10个格子,每个格子0.1x0.1 : 但程序跑着跑着,发现有些地方0.1x0.1的resolution不够,需要调整这些地方 : 的grid size...这个adaptive grid size的问题在fortran里用啥数据解构来实现呢? : 谁有这个方面的经验吗?谢了先
|