JAW,WA functions in Passarino-Veltman reduction

We continue to study a special class of Passarino-Veltman functions J arising at the reduction of infrared divergent box diagrams. We describe a procedure of separation of two types of singularities, infrared and mass singularities, which are absorbed in simple C0 functions. The infrared divergences of C0's can be regularized then by any method: photon mass, dimensionally or by the width of an unstable particle. Functions D0 are represented as certain linear combinations of the standard C0 Passarino-Veltman functions and infrared finite functions J. Then mass singularities are extracted from J to other combinations of C0. The rests are free of both types of singularities and are expressed as explicit and compact linear combination of logarithms and dilogarithm functions. The extensive comparison of numerical results with those obtained with the aid of the Loop Tools package is presented.