The point-kernel (PK) technique has a long history in applied radiation shielding, originating from the early days of digital computers. The PK technique applied to gamma-ray attenuation is one of many successful applications, based on the linear superposition principle applied to distributed radiation sources. Mathematically speaking, the distributed source will produce a detector response equivalent to the numerical integration of the radiation received from an equivalent number of point sources. In this treatment, there is no interference between individual point sources, while inherent limitations of the PK method are its inability to simulate gamma scattering in shields and the usage of simple boundary conditions. The PK method generally works for gamma-ray shielding with corrective B-factor for scattering and only specifically for fast neutron attenuation in a hydrogenous medium with the definition of cross section removal. This paper presents theoretical and programming aspects of the PK program developed for a distributed source of photons (line, disc, plane, sphere, slab volume, etc.) and slab shields. The derived flux solutions go beyond classical textbooks as they include the analytical integration of Taylor B-factor, obtaining a closed form readily suitable for programming. The specific computational modules are unified with a graphical user interface (GUI), assisting users with input/output data and visualization, developed for the fast radiological characterization of simple shielding problems. Numerical results of the selected PK test cases are presented and verified with the CADIS hybrid shielding methodology of the MAVRIC/SCALE6.1.3 code package from the ORNL.