2 ! { dg
-additional
-options
"-Ofast" }
3 SUBROUTINE GAUBON
(NV
,PTS
,PP
)
4 IMPLICIT DOUBLE PRECISION(A
-H
,O
-Z
)
6 DIMENSION PTS
(3,10),PP
(3)
7 COMMON /PCMPLY
/ XE
(MXSP
),YE
(MXSP
),ZE
(MXSP
)
8 DATA PI
/3.141592653589793D
+00/
11 PP
(1) = PP
(1) + (PTS
(1,I
)-XE
(NS
))
12 PP
(2) = PP
(2) + (PTS
(2,I
)-YE
(NS
))
13 PP
(3) = PP
(3) + (PTS
(3,I
)-ZE
(NS
))