File indexing completed on 2025-08-05 08:21:15
0001
0002
0003
0004
0005
0006
0007 FUNCTION PYRVG4(Y)
0008
0009 IMPLICIT NONE
0010 COMMON/PYRVPM/RM(0:3),A(2),B(2),RESM(2),RESW(2),MFLAG
0011 COMMON/PYG2DX/X
0012 DOUBLE PRECISION X, Y, PYRVG4, RM, A, B, RESM, RESW, RVS, PYRVS
0013 LOGICAL MFLAG
0014 SAVE /PYRVPM/,/PYG2DX/
0015 PYRVG4=0D0
0016 RVS=PYRVS(X,Y,RESM(1),RESW(1),RESM(2),RESW(2))
0017 IF (.NOT.MFLAG) THEN
0018 PYRVG4 = RVS*B(1)*B(2)*X*Y
0019 ELSE
0020 PYRVG4 = RVS*(RM(1)*RM(3)*A(1)*A(2)*(X+Y-RM(1)**2-RM(3)**2)
0021 & + RM(1)*RM(0)*B(1)*A(2)*(Y-RM(2)**2-RM(3)**2)
0022 & + RM(3)*RM(0)*A(1)*B(2)*(X-RM(1)**2-RM(2)**2)
0023 & + B(1)*B(2)*(X*Y-(RM(1)*RM(3))**2-(RM(0)*RM(2))**2))
0024 ENDIF
0025 RETURN
0026 END