i***k 发帖数: 45 | 1 【 以下文字转载自 Chemistry 讨论区 】
发信人: icook (icook), 信区: Chemistry
标 题: 问个计算化学问题:怎么读GRID?
发信站: BBS 未名空间站 (Tue May 3 12:55:17 2011, 美东)
初入计算化学行当,现在有一个势能FUNCTION, f(a,b,c,d,e),
我们从别的组拿到了f的ab initio数据,是17万行的数据,分别是变换ABCD得来的
我老板的意思是,做个GRIDDING,就相当于做个LIBRARY,
然后从里面内插,或者外插(interpolation, extrapolation)得到我们要的f值。
可是俺比较笨,不知道该怎么做。能请各位前辈指导一下么?我们用fortran,c
也可以。
1,我是先把17万行数据都读入程序先么?这样会不会内存废掉?俺是在ibm系列
的super computer上做,不知道会不会导致内存太大
2,假设我不先读入,而是每个f(a,b,c,d,e),我都到外部文件里面去查找,
这个怎么实现呢? 非常感谢关于这个的建议! 俺编程半白痴,实在不知道怎么实现。
如果一开始都读入,俺可以想办法实现。如果外部每次都读,就不知道了。:)
谁能帮忙一下? 多谢!
要是有人有类似的gridding的,外部读入的或者内部读入的程序可以介绍看一下,
就更好了! 多谢多谢!! | g*****k 发帖数: 623 | 2 17万行的数据应该不需要你做高深的优化。
应该可以直接读入内存操作。上了billion再考虑如何分块使用buffer
【在 i***k 的大作中提到】 : 【 以下文字转载自 Chemistry 讨论区 】 : 发信人: icook (icook), 信区: Chemistry : 标 题: 问个计算化学问题:怎么读GRID? : 发信站: BBS 未名空间站 (Tue May 3 12:55:17 2011, 美东) : 初入计算化学行当,现在有一个势能FUNCTION, f(a,b,c,d,e), : 我们从别的组拿到了f的ab initio数据,是17万行的数据,分别是变换ABCD得来的 : 我老板的意思是,做个GRIDDING,就相当于做个LIBRARY, : 然后从里面内插,或者外插(interpolation, extrapolation)得到我们要的f值。 : 可是俺比较笨,不知道该怎么做。能请各位前辈指导一下么?我们用fortran,c : 也可以。
| i***k 发帖数: 45 | 3 就是说,直接当做数组读入数据先,然后直接访问数组数据来进行内插和外插?
(interpolation extrapolation)
那就简单多了! 多谢!
【在 g*****k 的大作中提到】 : 17万行的数据应该不需要你做高深的优化。 : 应该可以直接读入内存操作。上了billion再考虑如何分块使用buffer
|
|