²é¿´: 2113  |  »Ø¸´: 9
µ±Ç°Ö»ÏÔʾÂú×ãÖ¸¶¨Ìõ¼þµÄ»ØÌû£¬µã»÷ÕâÀï²é¿´±¾»°ÌâµÄËùÓлØÌû

yjcmwgk

½û³æ (ÎÄ̳¾«Ó¢)

Ãܶȷºº¯¡¤Ð¡×ä

ÓÅÐã°æÖ÷

[ÇóÖú] ÓÃPython¶ÁGaussianµÄcubeÎļþ£¬ËٶȺÜÂý£¬ÔõÃ´ÆÆ£¿Ð»Ð»£¡£¡£¡ ÒÑÓÐ4È˲ÎÓë

Ïò¸÷λpython´óÏÀÇó½ÌÒ»¸öÎÊÌ⣬лл¸÷λ£¡
ÎÒÕýÔÚÓÃpython¶ÁÒ»¸öÓÉGaussianÉú³ÉµÄcubeÎļþ¡£µ«ÊǶÁÎļþ¶ÁµÃ·Ç³£Âý£¬¶ÁÒ»¸ö6M¶àµÄcubeÎļþ£¬¾ÓÈ»ÓÃÁË3·ÖÖÓ£¡

¾ßÌåÎÊÌâÈçÏ£º
ÓÃcubefileÀ´´ò¿ªÄ³¸öcubeÎļþ¡£È»ºó¼ÌÐø
CODE:
    tmpcubedata = [] !ÏÈÉèÒ»¸ö¿ÕÁÐ±í£¬´ý»á¶ùÓÃÀ´´æ´¢Êý¾Ý
    tmpcubedata = tmpcubedata + tmpline.split() !tmplineÊÇÓÃÀ´´æ´¢Ä³Ò»ÐеÄÊý¾ÝµÄ¡£ÓÉÓÚÏÈǰÒѾ­¾­¹ýÁ˶à´Îtmpline.split()Ö®ºó£¬Êµ¼ÊÉÏÏÖÔÚµÄtmpline´æ´¢µÄ¾ÍÊÇÎÒÒѾ­ÏëÒªµÄÄÇЩÊý¾ÝµÄµÚÒ»ÐÐÁË¡£ËùÒÔ£¬ÏȶÁһϲ¢´æÈëtmpcubedataÖС£
    i = 1 !ÕâÊǸö¼ÆÊýÆ÷¡£Ã¿¶ÁÒ»ÐÐÊý¾Ý¾Í¼Ó1
    tmplines=cubefile.readlines() !°ÑʣϵÄÐÐÈ«¶Á½øÀ´
    for i in range(1,linecount): !linecount¾ÍÊÇÓ¦¸Ã¶ÁÈëµÄcubeÎļþÖÐÊý¾Ý²¿·ÖµÄ×ÜÐÐÊý£¬Ò»°ãÀ´Ëµ£¬GaussianÉú³ÉµÄcubeÎļþÖУ¬Èç¹ûÉú³É¹ý³ÌȫĬÈÏ£¬ÄÇôcubeÎļþµÄÊý¾Ý²¿·ÖÓ¦¸ÃÊÇ8-9ÍòÐеÄÑù×Ó¡£´ó²¿·ÖÐУ¬Ã¿ÐÐÓÐ6¸öÊý×Ö£¨´ó¸Å95%µÄÐÐÊÇÕâÑùµÄ£©£¬µ«Ò²²»ÅųýÖ»ÓÐ1-5¸öÊý×ÖµÄÇé¿ö£¨´ó¸Å5%µÄÐÐÊÇÕâÑùµÄ£©¡£
        tmpline = tmplines[i-1] !´ÓtmplinesÖлñµÃijһÐÐ
        if tmpline == 'EOF' : break
        if tmpline == [] : break
        if tmpline == None : break
        if i == linecount: break !Èç¹û¼ÆÊýÆ÷i´ïµ½ÁË×ÜÐÐÊý£¬¾Í²»¶ÁÁË
        if i%10000 == 0: print "Converting", cubefilename, "......  ", int(100.0*float(i)/float(linecount)), "%" !ÿ¶ÁÒ»ÍòÐÐÊý¾Ý£¬¾ÍÏÔʾһ¸ö½ø¶È°Ù·Ö±È
        if i == linecount-1 : print "Converting", cubefilename, "...... ", "100 %" !¶Áµ½×îºóÒ»ÐÐÊý¾ÝµÄʱºò£¬¾ÍÏÔʾһ¸ö°Ù·ÖÖ®°Ù
        i = i+1
        tmpcubedata = tmpcubedata + tmpline.split() !½«¶Á½øÀ´µÄÊý¾Ý´æÈëtmpcubedataÖÐ

ÎÒ»³ÒÉtmpline = tmplines[i-1]Õâ¾ä»°ÍÏÀÛµÄËÙ¶È£¿Ó¦¸ÃÈçºÎÐ޸ģ¿·Ç³£¸Ðл¸ßÊÖ°ïÖú£¡£¡
¸Ðл£¡£¡
»Ø¸´´ËÂ¥

» ÊÕ¼±¾ÌûµÄÌÔÌûר¼­ÍƼö

SOB ¼¯½õ ¾­Ñé ÔÓÊé

» ²ÂÄãϲ»¶

» ±¾Ö÷ÌâÏà¹Ø¼ÛÖµÌùÍÆ¼ö£¬¶ÔÄúͬÑùÓаïÖú:

ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

frk2008

ľ³æ (ÕýʽдÊÖ)

¡¾´ð°¸¡¿Ó¦Öú»ØÌû

¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï
¸Ðл²ÎÓ룬ӦÖúÖ¸Êý +1
yjcmwgk: ½ð±Ò+20, ¡ïÓаïÖú, ·Ç³£¸Ðл£¬ÎÒºóÀ´ÓÃtmpp = ((' '.join(cubefile.readlines())).replace("/n"," ")).split()½â¾öÁË 2014-10-05 16:10:10
gmy1990: ½ð±Ò+3 2014-10-05 18:49:03
ÊÔÊÔÕâ¸ö£¬¿´¿´Ð§¹ûÈçºÎ£º

[ code]

tmplines = cubefile.readlines()  
tempstr = "".join(tmplines)             #°Ñlistת»¯Îªstr£¬ÐÎʽӦ¸ÃÊÇ'num  num  ...num\nnum  ...'
tempstr = tempstr.replace('\n','  ')    #°ÑstrÖеÄ'\n'Ìæ»»Îª'  '
tmpcubedata = tempstr.split('  ')       #°ÑstrתΪlist
n=0
numcount=len(tmpcubedata)               #listµÄÊý¾Ý×ÜÊý
for i in tmpcubedata:
    n+=1
    if i%60000==0:
        print "Converting", cubefilename, "......  ", int(100.0*float(i)/float(numcount)), "%"
    #if i=="":
        #tmpcubedata.remove(i)          #Èç¹ûcubefile ÎļþÖÐÖмäÓпÕÐУ¬ÔòÐèÒªÕâ¸ö¹¦ÄÜ
if tmpcubedata[-1]=="":
    del tmpcubedata[-1]                 #Èç¹ûcubefile ÎļþÒÔ¡®num\n¡®½áβ
print "Converting", cubefilename, "...... ", "100 %"

[ /code]
3Â¥2014-10-04 15:23:36
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
²é¿´È«²¿ 10 ¸ö»Ø´ð

yjcmwgk

½û³æ (ÎÄ̳¾«Ó¢)

Ãܶȷºº¯¡¤Ð¡×ä

ÓÅÐã°æÖ÷

cubeÎļþµÄ»ù±¾½á¹¹¾ÙÀýÈçÏÂ
CODE:
Title Card Required mo=1
MO coefficients
   -1   -2.916225   -2.916225   -2.916225
   80    0.073828    0.000000    0.000000
   80    0.000000    0.073828    0.000000
   80    0.000000    0.000000    0.073828
    2    2.000000    0.000000    0.000000    0.000000
    1    1
  7.36784E-05  8.36290E-05  9.46155E-05  1.06698E-04  1.19934E-04  1.34374E-04
  1.50064E-04  1.67044E-04  1.85341E-04  2.04976E-04  2.25955E-04  2.48274E-04
  2.71913E-04  2.96838E-04  3.22995E-04  3.50318E-04  3.78721E-04  4.08098E-04
  4.38328E-04  4.69271E-04  5.00768E-04  5.32647E-04  5.64718E-04  5.96779E-04
  6.28614E-04  6.60000E-04  6.90707E-04  7.20497E-04  7.49136E-04  7.76386E-04
  8.02019E-04  8.25811E-04  8.47551E-04  8.67043E-04  8.84107E-04  8.98583E-04
  9.10334E-04  9.19248E-04  9.25239E-04  9.28250E-04  9.28250E-04  9.25239E-04
  9.19248E-04  9.10334E-04  8.98583E-04  8.84107E-04  8.67043E-04  8.47551E-04
  8.25811E-04  8.02019E-04  7.76386E-04  7.49136E-04  7.20497E-04  6.90707E-04
  6.60000E-04  6.28614E-04  5.96779E-04  5.64718E-04  5.32647E-04  5.00768E-04
  4.69271E-04  4.38328E-04  4.08098E-04  3.78721E-04  3.50318E-04  3.22995E-04
  2.96838E-04  2.71913E-04  2.48274E-04  2.25955E-04  2.04976E-04  1.85341E-04
  1.67044E-04  1.50064E-04  1.34374E-04  1.19934E-04  1.06698E-04  9.46155E-05
  8.36290E-05  7.36784E-05
  8.36290E-05  9.49233E-05  1.07394E-04  1.21108E-04  1.36131E-04  1.52521E-04
  1.70331E-04  1.89603E-04  2.10372E-04  2.32658E-04  2.56471E-04  2.81805E-04
  3.08636E-04  3.36926E-04  3.66617E-04  3.97630E-04  4.29868E-04  4.63213E-04
ÏÂÃæ»¹ÓаËÍò¶àÐУ¬ÂÔÁË

¶ÔÓÚÕâ¸öcubeÎļþÀ´Ëµ£¬´ÓµÚ9ÐпªÊ¼ÊǵÚÒ»¸öÊý¾ÝÐУ¬Ö±µ½Îļþ½áÊø¡£ÖµµÃ×¢ÒâµÄÊÇ£¬¾ø´ó¶àÊýÊý¾ÝÐУ¬Ã¿ÐÐÓÐ6¸öÊý×Ö£¬µ«µÚ22ÐоÍÖ»ÓÐ2¸öÊý×ÖÁË¡£
2Â¥2014-10-04 11:08:58
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

rh1127

¾èÖú¹ó±ö (ÖøÃûдÊÖ)

Ò»Ãû¹¤³Ìʦ@ÉϺ£

¡¾´ð°¸¡¿Ó¦Öú»ØÌû

¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï
¸Ðл²ÎÓ룬ӦÖúÖ¸Êý +1
yjcmwgk: ½ð±Ò+20, ¡ïÓаïÖú, ·Ç³£¸Ðл£¬ÎÒºóÀ´ÓÃtmpp = ((' '.join(cubefile.readlines())).replace("/n"," ")).split()½â¾öÁË 2014-10-05 16:10:21
gmy1990: ½ð±Ò+2 2014-10-05 18:49:14
ÓÃnumpyºÍpandasÄ£¿é¶ùÀ´´¦Àí´óÐÍÊý×飬±ÈÖ»ÓÃforÑ­»·Ëٶȿ켸ʮ±¶²»Ö¹¡£python×Ô´øµÄÑ­»·Î´¾­ÓÅ»¯£¬ºÜ¼É»äÀ´´¦Àí´óÐÍÊý×飡
Ï£Íû¶ÔÄãÓаïÖú¡£
Actionnoworachievenever!
4Â¥2014-10-05 09:21:10
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

magicmonk

ÖÁ×ðľ³æ (ÖøÃûдÊÖ)

¡¾´ð°¸¡¿Ó¦Öú»ØÌû

¸Ðл²ÎÓ룬ӦÖúÖ¸Êý +1
¸öÈËÄØ½¨ÒéÓÃfortran»òMATLAB
ÎÒÃÇ×îÉîµÄ¿Ö¾å²»ÊÇÎÒÃÇÄÜÁ¦²»¹»£¬ÎÒÃÇ×îÉîµÄ¿Ö¾åÊÇÄÜÁ¦³¬Ô½Á˼«ÏÞ¡£
5Â¥2014-10-05 09:50:45
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] 317Çóµ÷¼Á +3 Éê×ÓÉêÉê 2026-03-19 6/300 2026-03-20 07:36 by Iveryant
[¿¼ÑÐ] 286·ÖÈ˹¤ÖÇÄÜרҵÇëÇóµ÷¼ÁÔ¸Òâ¿ç¿¼£¡ +3 lemonzzn 2026-03-17 3/150 2026-03-19 23:58 by 23Postgrad
[¿¼ÑÐ] ²ÄÁÏѧ˶318Çóµ÷¼Á +5 February_Feb 2026-03-19 5/250 2026-03-19 23:51 by 23Postgrad
[¿¼ÑÐ] Ò»Ö¾Ô¸ËÕÖÝ´óѧ²ÄÁÏÇóµ÷¼Á£¬×Ü·Ö315£¨Ó¢Ò»£© +3 sbdksD 2026-03-19 3/150 2026-03-19 23:21 by fmesaito
[¿¼ÑÐ] 307Çóµ÷¼Á +9 ÀäóÏ123 2026-03-17 9/450 2026-03-19 22:44 by ѧԱ8dgXkO
[¿¼ÑÐ] Çóµ÷¼Á +3 Mqqqqqq 2026-03-19 3/150 2026-03-19 14:11 by peike
[¿¼ÑÐ] ¡¾Í¬¼ÃÈí¼þ¡¿Èí¼þ£¨085405£©¿¼ÑÐÇóµ÷¼Á +3 2026eternal 2026-03-18 3/150 2026-03-18 19:09 by ²«»÷518
[¿¼ÑÐ] 085600²ÄÁÏÓ뻯¹¤ +5 °²È«Éϰ¶£¡ 2026-03-16 5/250 2026-03-18 15:33 by cmz0325
[¿¼ÑÐ] 311Çóµ÷¼Á +6 26ÑÐ0 2026-03-15 6/300 2026-03-18 14:43 by haxia
[¿¼ÑÐ] 312Çóµ÷¼Á +8 İå·Ï£ 2026-03-16 9/450 2026-03-18 12:39 by Linda Hu
[¿¼ÑÐ] 304Çóµ÷¼Á +12 СÐÜjoy 2026-03-14 13/650 2026-03-18 12:34 by Linda Hu
[¿¼ÑÐ] 0854£¬¼ÆËã»úÀàÕÐÊÕµ÷¼Á +3 ºúÀ±ÌÀ·ÅÌÇ 2026-03-15 6/300 2026-03-18 12:09 by Éϰ¶Éϰ¶¡­¡­..
[¿¼ÑÐ] 280Çóµ÷¼Á +6 ¹¾ààÏþÏþ 2026-03-18 7/350 2026-03-18 11:25 by Î޼ʵIJÝÔ­
[¿¼ÑÐ] ÉúÎïѧ071000 329·ÖÇóµ÷¼Á +3 ÎÒ°®ÉúÎïÉúÎﰮΠ2026-03-17 3/150 2026-03-18 10:12 by macy2011
[¿¼ÑÐ] 301Çóµ÷¼Á +9 yyÒªÉϰ¶Ñ½ 2026-03-17 9/450 2026-03-18 08:58 by Î޼ʵIJÝÔ­
[¿¼ÑÐ] 293Çóµ÷¼Á +11 zjlµÄºÅ 2026-03-16 16/800 2026-03-18 08:10 by zhukairuo
[¿¼ÑÐ] 301Çóµ÷¼Á +4 A_JiXing 2026-03-16 4/200 2026-03-17 17:32 by ruiyingmiao
[¿¼ÑÐ] 304Çóµ÷¼Á +5 ËØÄê¼ÀÓï 2026-03-15 5/250 2026-03-16 17:00 by ÎҵĴ¬Îҵĺ£
[¿¼ÑÐ] 070300»¯Ñ§Ñ§Ë¶Çóµ÷¼Á +6 Ì«Ïë½ø²½ÁË0608 2026-03-16 6/300 2026-03-16 16:13 by kykm678
[¿¼ÑÐ] 26¿¼ÑÐÒ»Ö¾Ô¸ÖйúʯÓÍ´óѧ(»ª¶«)305·ÖÇóµ÷¼Á +3 ¼ÎÄêÐÂ³Ì 2026-03-15 3/150 2026-03-15 13:58 by ¹þ¹þ¹þ¹þºÙºÙºÙ
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û