V700: Volume of Intersection of a Circular Cylinder with a Sphere

Author(s): K.S. Kölbig, F. Lamarche, C. Leroy	Library: MATHLIB
Submitter:	Submitted: 07.06.1992
Language: Fortran	Revised:

Function subprograms RVNSPC and DVNSPC calculate the volume of intersection of a circular cylinder of radius with a sphere of radius , the distance from the center of the sphere to the axis of the cylinder being .

This volume is given by

where the integration is performed over the intersection, if any, of the two circular disks and . If this is equal to

Otherwise .

On CDC and Cray computers, the double-precision version DVNSPC is not provided.

Structure:

FUNCTION subprograms
User Entry Names: RVNSPCC347, DVNSPCC347
External References: DELI3C, DELIKC, DELIEC

Usage:

In any arithmetic expression,

RVNSPC(R,RHO,D) or DVNSPC(R,RHO,D) has the value .

RVNSPC is if type REAL, DVNSPC is of type DOUBLE PRECISION, and R, RHO and D are of the same type as the function name.

Method:

The integral given above can be expressed in closed form in terms of complete elliptic integrals of the first, second, and third kind. For details see Ref. 1.

Notes:

Any negative sign in the parameters is ignored.

In the single-precision version RVNSPC on machines other than CDC or Cray, the complete elliptic integrals are calculated inside the subprogram. This version, faster than DVNSPC, is intended mainly for applications in experimental physics, where its limited accuracy of about 6 digits can be tolerated.

References:

1. F. Lamarche and C. Leroy, Evaluation of the volume of intersection of a sphere with a cylinder by elliptic integrals, Computer Phys. Comm. 59 (1990) 359-369.