Znn3bq.jpeg
²é¿´: 883  |  »Ø¸´: 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 µÄÖ÷Ìâ¸üÐÂ
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] 295·ÖÇóµ÷¼Á +5 ?ÒªÉϰ¶? 2026-04-17 5/250 2026-04-17 16:51 by fenglj492
[¿¼ÑÐ] 294Çóµ÷¼Á +7 µ­È»654321 2026-04-17 8/400 2026-04-17 16:36 by wutongshun
[¿¼ÑÐ] 335Çóµ÷¼Á +20 ÏëÉϰ¶Ñ½£¡£¡ 2026-04-12 23/1150 2026-04-17 10:50 by cuisz
[ÂÛÎÄͶ¸å] ÓÐûÓнÓÊձȽϿìµÄsciÆÚ¿¯Ñ½£¬×îºÃÔÚÒ»¸öÔÂÖ®Äڵģ¬ÑÐÈýº¢×ÓÇó±ÏÒµ 20+4 Ö®»¤×Å 2026-04-16 5/250 2026-04-17 10:02 by bobvan
[¿¼ÑÐ] 297£¬¹¤¿Æµ÷¼Á? +4 ºÓÄÏũҵ´óѧ-ÄÜ 2026-04-14 4/200 2026-04-16 22:52 by wulijun2012
[¿¼ÑÐ] »¯Ñ§070300 Çóµ÷¼Á +28 ¹þ¹þ¹þ^_^ 2026-04-12 28/1400 2026-04-16 21:36 by ´óÁ¦Ë®ÊÖÁ¦´óÎÞÇ
[¿¼ÑÐ] 297£¬¹¤¿Æµ÷¼Á?ºÓÄÏũҵ´óѧ±¾¿Æ +14 ºÓÄÏũҵ´óѧ-ÄÜ 2026-04-14 14/700 2026-04-16 14:41 by dingyanbo1
[¿¼ÑÐ] Çóµ÷¼Á +11 СÄô°®Ñ§Ï° 2026-04-11 15/750 2026-04-15 21:57 by noqvsozv
[¿¼ÑÐ] Ò»Ö¾Ô¸AÇø211£¬22408 321Çóµ÷¼Á +6 ËæÐÄËùÓû¡î 2026-04-15 7/350 2026-04-15 21:45 by lbsjt
[¿¼ÑÐ] ͨÐŹ¤³ÌÇóµ÷¼Á£¡£¡£¡ +6 zlb770521 2026-04-14 6/300 2026-04-15 20:00 by ѧԱJpLReM
[¿¼ÑÐ] 297¹¤¿Æµ÷¼Á? +14 ºÓÄÏũҵ´óѧ-ÄÜ 2026-04-13 15/750 2026-04-15 13:25 by ºÚ¿Æ¼¼¿óÒµ
[¿¼ÑÐ] 211±¾¿Æ²ÄÁÏ»¯¹¤Çóµ÷¼Á +19 YHLAH 2026-04-11 23/1150 2026-04-14 22:25 by fenglj492
[½Ìʦ֮¼Ò] ת³¤Æ¸ÁË +7 ¼òµ¥»¯xn 2026-04-13 7/350 2026-04-14 14:50 by xindong
[¿¼ÑÐ] Çóµ÷¼Á +12 ºÎÆøÕý 2026-04-13 13/650 2026-04-14 14:47 by zs92450
[¿¼ÑÐ] 105500ҩѧÇóµ÷¼Á +4 x_skys 2026-04-12 4/200 2026-04-14 13:37 by rndfc
[¿¼ÑÐ] 085600²ÄÁÏÓ뻯¹¤349·ÖÇóµ÷¼Á +16 Àîľ×Ó°¡¹þ¹þ 2026-04-12 17/850 2026-04-14 09:11 by fenglj492
[¿¼ÑÐ] ¿¼ÑÐÓ¢Ò»ÊýÒ»338·Ö +9 ³¤½­´óѧ¶«Ð£Çø 2026-04-13 10/500 2026-04-14 00:41 by Íõ¬Bè±
[¿¼ÑÐ] 339Çóµ÷¼Á +4 hanwudada 2026-04-12 4/200 2026-04-13 12:03 by À¶ÔÆË¼Óê
[¿¼ÑÐ] ¡£ +10 Àî¶àÃ×lee. 2026-04-12 11/550 2026-04-12 22:58 by yuyin1233
[¿¼ÑÐ] Ò»Ö¾Ô¸Ïôó0856£¬306Çóµ÷¼Á +15 Bblinging 2026-04-11 15/750 2026-04-11 22:53 by 314126402
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û