24小时热门版块排行榜    

查看: 1451  |  回复: 2

feynman1965

金虫 (正式写手)

[求助] doscar数据处理

A16B15(A,B分别代表一种原子)的体系,用 doscar文件除了有总体系的态密度,还有16个A原子每个的态密度,15原子个B每个的态密度。现在我已经把每个原子的态密度数据分割出来,总共16+15=31个态密度文件。现在我想画16个A原子的总态密度和15个B原子总的态密度,可以面临着两个难题:(1)每个原子的态密度数据格式为:ene  S+  S-  Py+  Py- Pz+  Pz-  Px+  Px-  Dxy+  Dxy-  Dyz+  Dyz-  Dz2+   Dz2-  Dxz+  Dxz-  Dx2  Dx2(每行19个数据),然后对每个原子态密度数据进行对应分量列相加得到:ene  S+  S-  P+  P-  d+  d-。这个工作量很大,需要对31个原子的态密度数据都进行这样的处理;(2)通过上步得到的31个原子态密度数据:ene  S+  S-  P+  P-  d+  d-,需要对16个A原子和15个B原子态密度数据对应每个分量列相加,得到总的A、B,这个操作也很复杂。请问各位大侠,有没有简单的方法实现以上两个步骤。最好是有程序!
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hn19870519

木虫 (正式写手)

就像你说的 写个程序啊。你莫不是想手动来解决问题。。。
2楼2013-01-05 21:57:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

了了工作室

新虫 (正式写手)

楼主解决了吗?
3楼2018-06-27 19:53:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 feynman1965 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] FileCode能看出啥? +8 要乐观耀哥 2026-08-10 26/1300 2026-08-13 12:01 by Tide man
[基金申请] 关于Filecode分析方法 +9 majunge000 2026-08-10 11/550 2026-08-13 10:40 by ZH远山
[基金申请] 重要来源:本周末出结果 +9 瞬息宇宙 2026-08-12 9/450 2026-08-13 10:12 by Vivilian
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +9 家与远方 2026-08-10 14/700 2026-08-13 09:41 by 臭臭不臭01
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 好奇怪的filecode +5 布布和一二 2026-08-08 6/300 2026-08-12 16:11 by 云上清扬
[基金申请] filecode +13 documentary 2026-08-10 15/750 2026-08-12 15:56 by 云上清扬
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +7 布布和一二 2026-08-10 10/500 2026-08-12 09:12 by 布布和一二
[基金申请] 是这周出结果还是下周出结果? +3 yuleib84 2026-08-11 3/150 2026-08-11 21:53 by jnhyjjm
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +3 majunge000 2026-08-11 4/200 2026-08-11 20:13 by lch2012
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[硕博家园] 读博的好处 +3 lnee 2026-08-11 3/150 2026-08-11 18:10 by 希望我好好的
[基金申请] 基金中了 +15 laoda193707 2026-08-06 15/750 2026-08-11 00:11 by jiafei2190
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] 这样的filecode谁见过 +11 布布和一二 2026-08-08 22/1100 2026-08-10 11:10 by wmfsnow
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
[基金申请] filecode与中标关系的预测 +5 布布和一二 2026-08-07 5/250 2026-08-09 16:15 by 袁向阳007
[基金申请] 固定端突然变了,今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
信息提示
请填处理意见