SCULIB_BESSEL_WTINIT - Generate a weighting function for rebinning

Description:

This is a FORTRAN version of the C code written for the transputer rebinning. Here is the C description: Initialise the BESSEL weighting function for Bessel interpolation. The 1-D function 2.0$*$J1(x)/x is initialised in the declaration, and the weighting function lookup table is calculated from this. The tabulation is at intervals of pi/RES for the range 0.0 to RADIUS$*$pi, resulting in (RES$*$RADIUS)$*$$*$2 values. The weight function is also multiplied by a cosine over the oute third to reduce edge effects.

The look up table is initialised such that one can access the weight directly given the square of the distance between the input and output pixel. We use a fixed resolution out to 10PI.

Invocation:

CALL SCULIB_BESSEL_WTINIT( WTFN, RADIUES, RES, STATUS)

Arguments:

WTFN ( RADIUS $*$ RADIUS $*$ RES $*$ RES $+$ 1) = REAL (Returned)

The weighting function generated by this routine. The index corresponds to the square of the distance from the centre in scale units.

RADIUS = INTEGER (Given)

Size of the weighting function in scale units.

RES = INTEGER (Given)

Number of points per scale lenghth.

STATUS = INTEGER (Given & Returned)

Global status

Copyright

Copyright (C) 1995,1996,1997,1998,1999 Particle Physics and Astronomy Research Council. All Rights Reserved.