Znn3bq.jpeg
²é¿´: 879  |  »Ø¸´: 3

voleyes

Òø³æ (ÕýʽдÊÖ)

[ÇóÖú] ÓÐûÓÐÄܹ»¸Ä±äµ°°×ÖÊ·Ö×Ó×ÜÌå¹¹Ð͵ÄÈí¼þ ÒÑÓÐ1È˲ÎÓë

±ÈÈçÎÒÓÐÒ»¸öÖ±Á´Ðεĵ°°×ÖÊ£¬Ïë°ÑËüÔÚÊÓ¾õЧ¹ûÉÏŪÍ䣬Õâ¸öÓÐÈí¼þÄܹ»²Ù×÷Âð£¿
Èç¹ûÒ»µãÒ»µãµÄ¸Ä¹¹ÐÍ»á±È½ÏÂé·³£¬ÒòΪµ°°×Öʹ¹Ðͺܴ󣬲»ºÃ²Ù×÷
»Ø¸´´ËÂ¥

» ²ÂÄãϲ»¶

» ±¾Ö÷ÌâÏà¹ØÉ̼ÒÍÆ¼ö: (ÎÒÒ²ÒªÔÚÕâÀïÍÆ¹ã)

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

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

polypro

ľ³æ (ÕýʽдÊÖ)

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

¡ï ¡ï ¡ï ¡ï ¡ï
¸Ðл²ÎÓ룬ӦÖúÖ¸Êý +1
voleyes: ½ð±Ò+5, ¡ï¡ï¡ïºÜÓаïÖú 2014-09-26 11:46:58
Õâ¸ö¼òµ¥£¬ÔÚÖ÷Á´Ä³¸öλÖÃתһ¸ö¶þÃæ½Ç(PHI»òPSI),ºóÃæµÄËùÓÐÔ­×Ó¶¼×ªÒ»Ï£¬ÐèÒªµÄ»°ÎÒ¸øÄãÒ»¸ö¶þÃæ½ÇתµÄ×Ó³ÌÐò£¬Íâ¼ÓÒ»¸öÑ­»·¾ÍÐÐÁË¡£
ȪºÔ,ÓãÏàÓë´¦ÓÚ½,Ïà…éÒÔʪ,Ïàå¦ÒÔÄ­,²»ÈçÏàÍüÓÚ½­ºþ¡£
2Â¥2014-09-24 20:42:16
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

voleyes

Òø³æ (ÕýʽдÊÖ)

ÒýÓûØÌû:
2Â¥: Originally posted by polypro at 2014-09-24 20:42:16
Õâ¸ö¼òµ¥£¬ÔÚÖ÷Á´Ä³¸öλÖÃתһ¸ö¶þÃæ½Ç(PHI»òPSI),ºóÃæµÄËùÓÐÔ­×Ó¶¼×ªÒ»Ï£¬ÐèÒªµÄ»°ÎÒ¸øÄãÒ»¸ö¶þÃæ½ÇתµÄ×Ó³ÌÐò£¬Íâ¼ÓÒ»¸öÑ­»·¾ÍÐÐÁË¡£

лл£¬³ÌÐòÄÜ·ñ·¢¸øÎÒ  okakasdu@126.com
3Â¥2014-09-26 11:49:27
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

polypro

ľ³æ (ÕýʽдÊÖ)

ÒýÓûØÌû:
3Â¥: Originally posted by voleyes at 2014-09-26 11:49:27
лл£¬³ÌÐòÄÜ·ñ·¢¸øÎÒ  okakasdu@126.com...

!!rotate atoms from iatom0 to iatom1     
do i=iatom0,iatom1   
        x=xx(i)
        y=yy(i)
        z=zz(i)
        call rotate_matrix2(x,y,z,ax,ay,az,xi,yi,zi,angg)  !C
        xx(i)=x
        yy(i)=y
        zz(i)=z        
     enddo
  subroutine rotate_matrix2(x,y,z,ax,ay,az,xi,yi,zi,an) !rotate (x,y,z) with origin (ax,ay,az) around axis (xi,yi,zi) with angle an
    implicit none
    real x,y,z,ax,ay,az,xi,yi,zi,an
    real r,r12,r23,r13,rcor
    real xr,yr,zr
    real x0,y0,z0    !!only use in test
    integer i,j,k
r=sqrt(xi*xi+yi*yi+zi*zi)
if(r<0.1)goto 101
x0=x
y0=y
x0=z
xi=xi/r
yi=yi/r
zi=zi/r
r12=xi*xi+yi*yi
r23=yi*yi+zi*zi
r13=xi*xi+zi*zi
rcor=xi*x+yi*y+zi*z
xr=ax*r23+xi*(-ay*yi-az*zi+rcor)+((x-ax)*r23+xi*(ay*yi+az*zi-yi*y-zi*z))*cos(an)+(ay*zi-az*yi-zi*y+yi*z)*sin(an)
yr=ay*r13+yi*(-ax*xi-az*zi+rcor)+((y-ay)*r13+yi*(ax*xi+az*zi-xi*x-zi*z))*cos(an)+(-ax*zi+az*xi+zi*x-xi*z)*sin(an)
zr=az*r12+zi*(-ax*xi-ay*yi+rcor)+((z-az)*r12+zi*(ax*xi+ay*yi-xi*x-yi*y))*cos(an)+(ax*yi-ay*xi-yi*x+xi*y)*sin(an)
x=xr
y=yr
z=zr
!if(abs(x)+abs(y)+abs(z)<0.001)then
!   print*,"warning",xi,yi,zi,r,x,y,z,an,x0,y0,z0
!endif
101 continue
  end subroutine rotate_matrix2
ȪºÔ,ÓãÏàÓë´¦ÓÚ½,Ïà…éÒÔʪ,Ïàå¦ÒÔÄ­,²»ÈçÏàÍüÓÚ½­ºþ¡£
4Â¥2014-09-27 21:07:41
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
Ïà¹Ø°æ¿éÌø×ª ÎÒÒª¶©ÔÄÂ¥Ö÷ voleyes µÄÖ÷Ìâ¸üÐÂ
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] ͨÐŹ¤³ÌÇóµ÷¼Á£¡£¡£¡ +5 zlb770521 2026-04-14 5/250 2026-04-15 19:11 by ºÚ¿Æ¼¼»¨¸ÚÑÒ
[¿¼ÑÐ] Ò»Ö¾Ô¸»ªÖÐũҵ071010£¬320Çóµ÷¼Á +11 À§À§À§À§À¤À¤ 2026-04-14 12/600 2026-04-15 18:58 by zhuwenxu
[¿¼ÑÐ] 296Çóµ÷¼Á +11 Íô£¡£¿£¡ 2026-04-09 12/600 2026-04-15 18:32 by ANÁ÷800
[¿¼ÑÐ] 291Çóµ÷¼Á +4 ¹ØÒä±±. 2026-04-14 4/200 2026-04-15 10:49 by shenshen207
[¿¼ÑÐ] һ־Ը³¶«´óѧ071000ÉúÎïѧѧ˶³õÊÔ·ÖÊý276Çóµ÷¼Á +26 Ľ¾øcc 2026-04-09 30/1500 2026-04-14 18:50 by ²ÌËÕÑô
[¿¼ÑÐ] Ò»Ö¾Ô¸Öйú¿ÆÑ§ÔºÉϺ£ÓлúËù£¬Óлú»¯Ñ§356·ÖÕÒµ÷¼Á +12 Nadiums 2026-04-09 13/650 2026-04-14 17:54 by lhj2009
[¿¼ÑÐ] 085600²ÄÁÏÓ뻯¹¤329·ÖÇóµ÷¼Á +24 Ò¶zilin 2026-04-13 25/1250 2026-04-14 09:20 by ÊÔ¹ÜÆÆÁÑ
[¿¼ÑÐ] 293Çóµ÷¼Á +16 ÎÒ°®¸ßÊý¸ßÊý°®Î 2026-04-12 18/900 2026-04-13 21:47 by ѧԱJpLReM
[¿¼ÑÐ] 302Çóµ÷¼Á +10 Ò×£¡? 2026-04-13 10/500 2026-04-13 19:04 by lbsjt
[¿¼ÑÐ] 290Çóµ÷¼Á +18 ¿Â»´È» 2026-04-12 20/1000 2026-04-13 12:56 by cyh¡ª315
[½Ìʦ֮¼Ò] ɽ¶«Ë«·ÇԺУ¿¼ºË³¬¼¶ÎÞµ×Ïߣ¬Áìµ¼ÐÒÔÖÀÖ»ö£¬½ÌʦÔâÑê¿Ö +3 qut2026 2026-04-11 7/350 2026-04-12 20:24 by qut2026
[¿¼ÑÐ] ҩѧר˶µ÷¼Á +8 ? һ·Éú?»¨? 2026-04-10 10/500 2026-04-11 21:21 by zhouxiaoyu
[¿¼ÑÐ] 359Çóµ÷¼Á +5 θ¾·ÂÎÀÛÁË 2026-04-11 5/250 2026-04-11 19:55 by lbsjt
[¿¼ÑÐ] ¹ã¶«Ê¡ 085601 329·ÖÇóµ÷¼Á +14 Eddieddd 2026-04-10 14/700 2026-04-11 09:58 by bljnqdcc
[¿¼ÑÐ] 311Çóµ÷¼Á +13 xypÏë¶ÁÊé 2026-04-10 14/700 2026-04-11 09:41 by Öí»á·É
[¿¼ÑÐ] 22408µ÷¼ÁÇóÖú +7 ì±12 2026-04-09 9/450 2026-04-11 09:23 by ŶŶ123
[¿¼ÑÐ] 284Çóµ÷¼Á +12 archer.. 2026-04-10 13/650 2026-04-11 08:44 by zhq0425
[¿¼ÑÐ] ±¾9 Ò»Ö¾Ô¸Î÷¹¤´ó085601 324Çóµ÷¼Á +5 wysyjs25 2026-04-10 5/250 2026-04-10 16:57 by luoyongfeng
[¿¼ÑÐ] 298Çóµ÷¼Á +13 ¶¤¶£ß˶¬¹Ï 2026-04-09 13/650 2026-04-10 15:49 by jiajinhpu
[¿¼ÑÐ] Ò»Ö¾Ô¸»¦9£¬326ÉúÎïѧÇóÏà¹Ø×¨Òµµ÷¼Á +4 Áõīī 2026-04-09 4/200 2026-04-10 12:07 by pengliang8036
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û