Exact computation of the 3-f and 6-f symbols - iucaa

More documents

Recommendations

Info

S. -T Lai~Y.-N. C’hiu / Exact computation of the 3-f and 6-f symbols 359 455 QP.i.OQO L1—INT)Si+32+33+1) WRITE)2i,429) N9—.5 429 FORMAT(11, ‘APPROXIMATE VALUE OF 3—) SYMBOL’/) CALL DECOM(Li.QR.R9) 00430 II.i,ui0 Li—INT)34535—33) IF)ABS)PP)II)).LT.i.OE—5) PP)II).0.S RB—S IF)ABS(PP)II)).LY.i.SE—5) GOYO 430 CALL DECOM)Li,QK,R9) QPQP’Qi(II)’’PP(II) Li—IIOT)Ji+32—33) 430 CONTINUE N9—. Q3)—CC’QP’QABS(DIDS) CALL DECOM)L1,QK,R9) WRIYE)2i,450) Q3) Li.INY(34—35*33) 450 ?ORMAT(iI, ‘TUB)—’ ,F37.34) R9—.5 0010 600 CALL DECOM)L1.Qk.B9) 550 TBR)—S.0 Li.INT)—Ji+32+33) IP)TNR).EQ.S.0) WRITE(21,550) R9—.5 550 PORMAT(iI.’YHR)—O.S’) CALL DECOM)Li.QE.R9) 600 CONTINUE Li.INT)—34.35+33) RETURN R9.S END CALL DECOM)Li,QK,R9) C L1ENY(24+35+33+i) SUBROUTINE A6)SY R9—.S IMPLICIT REAL’9(A—H,O—Z) CALL DECOM)Li,QK,R9) REAL’S 31,32.33,34,35,36,) Li.INT)34s.JS—36) REAL’i6 DIAi,DIA2,01A3,DIA4,DIB1,DIB2,DIB3,DIB4 R9.S REAL’i6 DIBS,DIB6 ,DIB7,DIBB ,DEFAC.0Y3 ,DIC1 ,DIC2 CALL DECOM(Li,QE.R9) REAL*i6 DIC3,DIC4,DICC,DIDS,01A,KQ,QE,QP,Q3) Li—INT)34+)2+36+i) INTEOER’2 Qi(uio) RB.5 DIMENSION PP(ii0) CALL DECOM)Li,QK,R9) COMMON/PRI/Q1,PP LiINT))4—32+36) 5 DO 10 lu—illS PP)IU)0.O CALL DECOM)L1.QE,R9) 10 CONTINUE Li—INT)—34+)S’36) QK.i.OQO R9.5 12 WRITE)2i,iS) CALL DECOM)Li,QK,R9) 15 FORMAT(iX.’)i—?’,’32.?’,’33—P’.’)4.?’.’)5.?,’36—?’/) L5—INT)Ji+3S—36) REAO)S2,i9) Ji,32,33,34,3S,36 R9..S 19 FORMAT(6CPiO.i)) CALL DECOM)L1,QE,R9) WRIYE)2i .20)31 .32,33 L1INT)ii+)S+36+i) SO FORMAT)iX,3)F9.i)/) R9—.S BORIYE)2i.2S)34,JS,36 CALL DECOM)Li.QK,R9) 22 FORMAT)SX,3)FB.i)/) L1INT)J1—)5S)6) C R9—.5 iMi.Ji+JS+33 CALL DECOM)Li,QK,R9) 3MS—)i*35+36 Li.INT(—31,3S+)6) 3N3.i3*34÷)6 RB— .5 3M4)3+)4.35 CALL DECOM)L1.QE,R9) )Xi)i+3S+)4+35 IP)Ji.GY.S0.0) OOTO 400 3X2—)5+33+)5+36 IF)32.GT.20.0) 0010 400 3X3.)i+)3+34+.76 IP))3.GT.20.0) 0010 400 C IF)34.GT.20.0) OOTO 400 )M)Mi IF)35.OT.S0.0) 0010 400 IP)3M2.G1.JM)JM—3M2 IF)36.GY.S0.S) OOYO 400 IF))M3.GT.JM)JM.3M3 WBITE)2i,4i9) IP(7M4.GT.3M)JM.)M4 419 PORMAT)1X, ‘EXACT VALUE OP 6—) SYMBOL’/) )XJXi VRITE(Si,4S0) CC IP))XS.LT.3X))X—3X2 420 FORMAT)iX,’SIXJ.’,F4.i) IP))X3.LT.)X))X.)X3 0042511—1,110 KM.IMT))M) IF)ABS(PP)II)).LY.i.OE—5) 8010 425 KX.INT))X) WRITE)Si,422) Qi(II),PP)II) IF(EM.GY.RX) GOYO 500 422 FORMAT)iOX,’’,I4,’’’’,F7.i) DIDS—0.008 425 CONTINUE DO 250 LKM,KX WRITE)21,420) PL*QFLOAT)L) 428 FORMAT(1X,’MULTIPLY) DIA1—DIPAC)L.1) KQ..QINY)QK+.500) DIA2.DIFAC(INT)Ji,J2—33)) IF)KQ.LT.601.OQO) tCQ—i.OQO DEA3—DIPAC)INT)34—)S+36)) WRITE)2i,’)EQ DIA4—DIFAC)INT)—Ji+JS+JS)) GOTO 650 DEN1DIPAC)INT)PL—.21—)2—33)) 400 QP1.OQO 0182.DEFAC)INT))4*.32+36+i,)) WRITE)21,429) DIB3—DIFAC)INT)PL—J4—35—33)) 429 FORMAT)1X,’APPROXIMATE VALUE OF 6—) SYMBOL’!) DIB4—DIPAC)INT)JS+3i,)4+35—PL)) DO 438 1I1,iiO DIBS.DIFAC)INT)PL—)4—32—)6)) IF)ABS(PP)II)).LT.i.SE—5) PP)II)—S.S DIB6.DIFAC)INT)33+)6+)1+34—pL)) IF)ABS(PP)II)).LT.i.OE—5) 0010 430 DIB7—DIPAC)INT)PL—31—)5—)6)) QPQP’Q1)II)”PP(II) DIBO—DIFAC)INT)35+334)S+.76—pL)) 430 CONTINUE DIC1DIA1/)DIBSDIB2) Q3).CCQP*QABS(DIOS) DICSDIAS/)DIB3’DIB4) WRITE)2i,450) Q33 DIC3.DIA3/)DIBiDIB6) 450 FORMAT)iX,’SII)*’,F37.34) DIC4.DIA4/)DLB7’DIB8) GOTO 650 DICC—OICi’DIC2DIC3’01C4 500 SIX)—S.0 DIDSDIDS+)—i)”L’DICC IF(SII).EQ.0 .0) WRITE(2i,600) 250 CONTINUE 600 PORMAT(iI,’SIX)=0.3) IF)DIDS.EQ.0.500) 0010 580 650 CONTINUE IF(DIDS.LT.S.OQO) CC.—i.0 RETURN IF(DIDS.GT.0.SQO) CC.i.0 END DIAQABS)D105) C IF))i.GT.2S.0) GOTO 255 C ‘CALCULATION OF Nt IF))S.OT.S0.0) OOTO 255 C IF)33.GT.20.O) GOTO 2SS FUNCTION DIPAC(N) IP)34.0Y.20.0) OOTO 255 REAL16 DIFAC,Pi IF)35.GT.2B.0) GOYO 255 INTEOER’2 K,N IF)36.GY.SO,O) 0010 255 DIFAC.i.OQQ NO.1.0 SF (N.LE.i) GOTO 20 CALL PRIME(DIA,OE,R9) 00 10 K1,N 255 L1.INT)31—32+33) Pi.QFLOAT(K) R9.5 OIFACOIFACP1 CALL D!COM)Li.QX,R9) 10 CONTINUE
360 S.-T Lai, Y -N. Chiu / Exact computation of the 3-f and 6-f symbols 20 RETURN 15 CONTINUE END RETURN C END SUBROUTINE PBIMi(Qit C IMPLICIT REALB)A—U,O—I) SUBROUTINE PRIME)DIA,QK,BB) INTEOER*2 Qi(11S) IMPLICIT REALO(A—N,O—Z) LOGICAL KR INTEGER’S Qi)liS),LY Qi(i)—2 REALi6 DIA.DY3,QE NI2 LOGICAL LL,LQ 00101—3,601 DIMENSION PP)ii0) DO 25 3—2,1—i COMMON/PRI/Q1,PP KK—ABS(FL0AT(I/))—FLOAT(I)/FLOAT~fl).LT.i.0E—5 LY—0 IF(EE) GOTO iS LM1 20 CONTINUE DIA—QINT)DEA+.5Q5) Qi(NI)—I 15 IF)LM.OT.LY) LYLM NI—N1+i IF)DIA.EQ.i.558.OR.DIA.EQ.0.SQO) 001040 iO CONTINUE IF)01)LM).EQ.0) QKDIA RETURN IF)Q1(LM).EQ.0) 001040 END SO DT3DEA/QFLOAT)Qi)LM)) C LQ.QABS(DT3—i.000).LT.i.OQ—5 SUBROUTINE DECOM)Li,QE.R9) IF)LQ) PP)LM)—PP)LM)+R9 IMPLICIT REALB)A—U.O—Z) IF)LQ) 001040 REAL’16 DIA,DT3.QE LL.QABS(D13—QINT)DY3*.SQO)).LT. 1.SQ—iO INTEOERS Qi)iio) IF)LL) PP)LM)—PP)LM)+R9 DIMENSION PP)uio) IF)LL) DIA:DY3 COMMON/PRI/Ql,PP IF)LL) 0010 50 IF(Li.EQ.0.OR.Li .EQ.i) RETURN LMLM+i DO 10 LEi.Li 0010 10 DIA—QFLOAT)LK) 40 RETURN CALL PRIME)DIA.QK,R9) END References [11 M. Rotenberg, R. Bivins, N. Metropolis and J.K. Wooten Jr.. The 3-f and 6-f Symbols (The Technology Press. MIT, Cambridge, MA, ~959). [2] K. Schulten and R.G. Gordon, J. Math. Phys. 16 (1975) 1961; Comput. Phys. Commun. 11(1976) 269. [3] RN. Zare, Angular Momentum (Wiley, New York, 1988).
Page 1 and 2: 350 Computer Physics Communications
Page 3 and 4: 352 S. -T. Lai, Y. -N. Chiu / Exact
Page 5 and 6: 354 5. -T Lai, 1’. -N. Chiu / Exa
Page 7 and 8: 356 5. -T Lai, Y. -N. Chiu / Exact
Page 9: 358 S.-T Liii, Y. -N. Chiu / Exact

Exact computation of the 3-f and 6-f symbols - iucaa

You also want an ePaper? Increase the reach of your titles

Delete template?

Save as template?