c09acf returns the details of the chosen three-dimensional discrete wavelet filter. For a chosen mother wavelet, discrete wavelet transform type (single-level or multi-level DWT) and end extension method, this routine returns the maximum number of levels of resolution (appropriate to a multi-level transform), the filter length, the total number of coefficients and the number of wavelet coefficients in the second and third dimensions for the single-level case. This routine must be called before any of the three-dimensional transform routines in this chapter.

2 Specification

Fortran Interface

Subroutine c09acf (

wavnam, wtrans, mode, m, n, fr, nwlmax, nf, nwct, nwcn, nwcfr, icomm, ifail)

Integer, Intent (In)	::	m, n, fr
Integer, Intent (Inout)	::	ifail
Integer, Intent (Out)	::	nwlmax, nf, nwct, nwcn, nwcfr, icomm(260)
Character (*), Intent (In)	::	wavnam
Character (1), Intent (In)	::	wtrans, mode

C Header Interface

#include <nag.h>

void

c09acf_ (const char *wavnam, const char *wtrans, const char *mode, const Integer *m, const Integer *n, const Integer *fr, Integer *nwlmax, Integer *nf, Integer *nwct, Integer *nwcn, Integer *nwcfr, Integer icomm[], Integer *ifail, const Charlen length_wavnam, const Charlen length_wtrans, const Charlen length_mode)

The routine may be called by the names c09acf or nagf_wav_dim3_init.

3 Description

Three-dimensional discrete wavelet transforms (DWT) are characterised by the mother wavelet, the end extension method and whether multiresolution analysis is to be performed. For the selected combination of choices for these three characteristics, and for given dimensions (

m \times n \times fr

) of data array

A

, c09acf returns the dimension details for the transform determined by this combination. The dimension details are:

l_{\max}

, the maximum number of levels of resolution that would be computed were a multi-level DWT applied;

n_{f}

, the filter length;

n_{ct}

the total number of wavelet coefficients (over all levels in the multi-level DWT case);

n_{cn}

, the number of coefficients in the second dimension for a single-level DWT; and

n_{cfr}

, the number of coefficients in the third dimension for a single-level DWT. These values are also stored in the communication array icomm, as are the input choices, so that they may be conveniently communicated to the three-dimensional transform routines in this chapter.

4 References

None.

5 Arguments

1: $wavnam$ – Character(*) Input

On entry: the name of the mother wavelet. See the C09 Chapter Introduction for details.

$wavnam ='HAAR' or 'DB1'$: Haar wavelet, also known as $'DB1'$ as a special case of the Daubechies wavelet.
$wavnam =' DB n'$ , where $n = 2, 3, \dots, 38$: Daubechies wavelet with $n$ vanishing moments ( $2 n$ coefficients). For example, $wavnam ='DB4'$ is the name for the Daubechies wavelet with $4$ vanishing moments ( $8$ coefficients).
$wavnam =' COIF n'$ , where $n = 1, 2, \dots, 17$: Coiflet wavelet of order $n$ .
$wavnam ='BEYL'$: Beylkin wavelet.
$wavnam ='VAID'$: Vaidyanathan wavelet.
$wavnam =' SYM n'$ , where $n = 2, 3, \dots, 20$: Symlet wavelet of order $n$ .
$wavnam =' BIOR x . y'$ , where $x . y$ can be one of 1.1, 1.3, 1.5, 2.2, 2.4, 2.6, 2.8, 3.1, 3.3, 3.5, 3.7, 3.9, 4.4, 5.5 or 6.8: Biorthogonal wavelet of order $x . y$ . For example $wavnam ='BIOR3.1'$ is the name for the biorthogonal wavelet of order $3.1$ .
$wavnam =' RBIO x . y'$ , where $x . y$ can be one of 1.1, 1.3, 1.5, 2.2, 2.4, 2.6, 2.8, 3.1, 3.3, 3.5, 3.7, 3.9, 4.4, 5.5 or 6.8: Reverse biorthogonal wavelet of order $x . y$ . For example $wavnam ='RBIO3.1'$ is the name for the reverse biorthogonal wavelet of order $3.1$ .

Constraint:

wavnam ='HAAR', ​'DB1'

'DB2'

'DB3'

'DB4'

'DB5'

'DB6'

'DB7'

'DB8'

'DB9'

'DB10'

'DB11'

'DB12'

'DB13'

'DB14'

'DB15'

'DB16'

'DB17'

'DB18'

'DB19'

'DB20'

'DB21'

'DB22'

'DB23'

'DB24'

'DB25'

'DB26'

'DB27'

'DB28'

'DB29'

'DB30'

'DB31'

'DB32'

'DB33'

'DB34'

'DB35'

'DB36'

'DB37'

'DB38'

'COIF1'

'COIF2'

'COIF3'

'COIF4'

'COIF5'

'COIF6'

'COIF7'

'COIF8'

'COIF9'

'COIF10'

'COIF11'

'COIF12'

'COIF13'

'COIF14'

'COIF15'

'COIF16'

'COIF17'

'BEYL'

'VAID'

'SYM2'

'SYM3'

'SYM4'

'SYM5'

'SYM6'

'SYM7'

'SYM8'

'SYM9'

'SYM10'

'SYM11'

'SYM12'

'SYM13'

'SYM14'

'SYM15'

'SYM16'

'SYM17'

'SYM18'

'SYM19'

'SYM20'

'BIOR1.1'

'BIOR1.3'

'BIOR1.5'

'BIOR2.2'

'BIOR2.4'

'BIOR2.6'

'BIOR2.8'

'BIOR3.1'

'BIOR3.3'

'BIOR3.5'

'BIOR3.7'

'BIOR3.9'

'BIOR4.4'

'BIOR5.5'

'BIOR6.8'

'RBIO1.1'

'RBIO1.3'

'RBIO1.5'

'RBIO2.2'

'RBIO2.4'

'RBIO2.6'

'RBIO2.8'

'RBIO3.1'

'RBIO3.3'

'RBIO3.5'

'RBIO3.7'

'RBIO3.9'

'RBIO4.4'

'RBIO5.5'

'RBIO6.8'

2: $wtrans$ – Character(1) Input

On entry: the type of discrete wavelet transform that is to be applied.

$wtrans ='S'$: Single-level decomposition or reconstruction by discrete wavelet transform.
$wtrans ='M'$: Multiresolution, by a multi-level DWT or its inverse.

Constraint:

wtrans ='S'

'M'

3: $mode$ – Character(1) Input

On entry: the end extension method.

$mode ='P'$: Periodic end extension.
$mode ='H'$: Half-point symmetric end extension.
$mode ='W'$: Whole-point symmetric end extension.
$mode ='Z'$: Zero end extension.

Constraint:

mode ='P'

'H'

'W'

'Z'

4: $m$ – Integer Input

On entry: the number of elements,

m

, in the first dimension (number of rows of each two-dimensional frame) of the input data,

A

Constraint:

m \geq 2

5: $n$ – Integer Input

On entry: the number of elements,

n

, in the second dimension (number of columns of each two-dimensional frame) of the input data,

A

Constraint:

n \geq 2

6: $fr$ – Integer Input

On entry: the number of elements,

fr

, in the third dimension (number of frames) of the input data,

A

Constraint:

fr \geq 2

7: $nwlmax$ – Integer Output

On exit: the maximum number of levels of resolution,

l_{\max}

, that can be computed if a multi-level discrete wavelet transform is applied (

wtrans ='M'

). It is such that

2^{l_{\max}} \leq \min (m, n, fr) < 2^{l_{\max} + 1}

, for

l_{\max}

an integer.

wtrans ='S'

, nwlmax is not set.

8: $nf$ – Integer Output

On exit: the filter length,

n_{f}

, for the supplied mother wavelet. This is used to determine the number of coefficients to be generated by the chosen transform.

9: $nwct$ – Integer Output

On exit: the total number of wavelet coefficients,

n_{ct}

, that will be generated. When

wtrans ='S'

the number of rows required (i.e., the first dimension of each two-dimensional frame) in each of the output coefficient arrays can be calculated as

n_{cm} = n_{ct} / (8 \times n_{cn} \times n_{cfr})

. When

wtrans ='M'

the length of the array used to store all of the coefficient matrices must be at least

n_{ct}

10: $nwcn$ – Integer Output

On exit: for a single-level transform (

wtrans ='S'

), the number of coefficients that would be generated in the second dimension,

n_{cn}

, for each coefficient type. For a multi-level transform (

wtrans ='M'

) this is set to

1

11: $nwcfr$ – Integer Output

On exit: for a single-level transform (

wtrans ='S'

), the number of coefficients that would be generated in the third dimension,

n_{cfr}

, for each coefficient type. For a multi-level transform (

wtrans ='M'

) this is set to

1

12: $icomm (260)$ – Integer array Communication Array

On exit: contains details of the wavelet transform and the problem dimension which is to be communicated to the three-dimensional discrete transform routines in this chapter.

13: $ifail$ – Integer Input/Output

On entry: ifail must be set to

0

−1

1

to set behaviour on detection of an error; these values have no effect when no error is detected.

A value of

0

causes the printing of an error message and program execution will be halted; otherwise program execution continues. A value of

−1

means that an error message is printed while a value of

1

means that it is not.

If halting is not appropriate, the value

−1

1

is recommended. If message printing is undesirable, then the value

1

is recommended. Otherwise, the value

0

is recommended. When the value $- 1$ or $1$ is used it is essential to test the value of ifail on exit.

On exit:

ifail = 0

unless the routine detects an error or a warning has been flagged (see Section 6).

6 Error Indicators and Warnings

If on entry

ifail = 0

−1

, explanatory error messages are output on the current error message unit (as defined by x04aaf).

Errors or warnings detected by the routine:

$ifail = 1$: On entry, $wavnam = ⟨ value ⟩$ was an illegal value.

$ifail = 2$: On entry, $wtrans = ⟨ value ⟩$ was an illegal value.

$ifail = 3$: On entry, $mode = ⟨ value ⟩$ was an illegal value.

$ifail = 4$: On entry, $fr = ⟨ value ⟩$ .
Constraint: $fr \geq 2$ .

On entry, $m = ⟨ value ⟩$ .
Constraint: $m \geq 2$ .

On entry, $n = ⟨ value ⟩$ .
Constraint: $n \geq 2$ .

$ifail = - 99$: An unexpected error has been triggered by this routine. Please contact NAG.
See Section 7 in the Introduction to the NAG Library FL Interface for further information.

$ifail = - 399$: Your licence key may have expired or may not have been installed correctly.
See Section 8 in the Introduction to the NAG Library FL Interface for further information.

$ifail = - 999$: Dynamic memory allocation failed.
See Section 9 in the Introduction to the NAG Library FL Interface for further information.

7 Accuracy

Not applicable.

8 Parallelism and Performance

Background information to multithreading can be found in the Multithreading documentation.

c09acf is not threaded in any implementation.

9 Further Comments

None.

10 Example

This example computes the three-dimensional multi-level resolution for

8 \times 8 \times 8

input data by a discrete wavelet transform using the Daubechies wavelet with four vanishing moments (see

wavnam ='DB4'

in c09acf) and zero end extension. The number of levels of transformation actually performed is one less than the maximum possible. This number of levels, the length of the wavelet filter, the total number of coefficients and the number of coefficients in each dimension for each level are printed along with the approximation coefficients before a reconstruction is performed. This example also demonstrates in general how to access any set of coefficients at any level following a multi-level transform.

c09ac: FL CL CPP AD PY MB

NAG FL Interfacec09acf (dim3_​init)

▸▿ Contents

1 Purpose

2 Specification

3 Description

4 References

5 Arguments

6 Error Indicators and Warnings

7 Accuracy

8 Parallelism and Performance

9 Further Comments

10 Example

10.1 Program Text

10.2 Program Data

10.3 Program Results

NAG FL Interface
c09acf (dim3_init)