void	e02zac (Integer px, Integer py, const double lamda[], const double mu[], Integer m, const double x[], const double y[], Integer point[], NagError *fail)

The function may be called by the names: e02zac, nag_fit_dim2_spline_sort or nag_2d_panel_sort.

3 Description

A set of

m

data points with rectangular Cartesian coordinates

x_{r}, y_{r}

are sorted into panels defined by lines parallel to the

y

and

x

axes. The intercepts of these lines on the

x

and

y

axes are given in

lamda [i - 1]

, for

i = 5, 6, \dots, px - 4

and

mu [j - 1]

, for

j = 5, 6, \dots, py - 4

, respectively. The function orders the data so that all points in a panel occur before data in succeeding panels, where the panels are numbered from bottom to top and then left to right, with the usual arrangement of axes, as shown in the diagram. Within a panel the points maintain their original order.

Figure 1

A data point lying exactly on one or more panel sides is taken to be in the highest-numbered panel adjacent to the point. The function does not physically rearrange the data, but provides the array point which contains a linked list for each panel, pointing to the data in that panel. The total number of panels is

(px - 7) \times (py - 7)

4 References

None.

5 Arguments

1: $px$ – Integer Input
2: $py$ – Integer Input: On entry: px and py must specify eight more than the number of intercepts on the $x$ axis and $y$ axis, respectively.

Constraint: $px \geq 8$ and $py \geq 8$ .
3: $lamda [px]$ – const double Input: On entry: $lamda [4]$ to $lamda [px - 5]$ must contain, in nondecreasing order, the intercepts on the $x$ axis of the sides of the panels parallel to the $y$ axis.
4: $mu [py]$ – const double Input: On entry: $mu [4]$ to $mu [py - 5]$ must contain, in nondecreasing order, the intercepts on the $y$ axis of the sides of the panels parallel to the $x$ axis.
5: $m$ – Integer Input: On entry: the number $m$ of data points.
6: $x [m]$ – const double Input
7: $y [m]$ – const double Input: On entry: the coordinates of the $r$ th data point $(x_{r}, y_{r})$ , for $r = 1, 2, \dots, m$ .
8: $point [\dim]$ – Integer Output: Note: the dimension, dim, of the array point must be at least $(m + (px - 7) \times (py - 7))$ .

On exit: for $i = 1, 2, \dots, (m + (px - 7) \times (py - 7))$ , $point [m + i - 1] = I1$ is the index of the first point in panel $i$ , $point [I1 - 1] = I2$ is the index of the second point in panel $i$ and so on.
$point [In - 1] = 0$ indicates that $x [In - 1], y [In - 1]$ was the last point in the panel.
9: $fail$ – NagError * Input/Output: The NAG error argument (see Section 7 in the Introduction to the NAG Library CL Interface).

6 Error Indicators and Warnings

NE_ALLOC_FAIL: Dynamic memory allocation failed.
See Section 3.1.2 in the Introduction to the NAG Library CL Interface for further information.
NE_BAD_PARAM: On entry, argument $⟨ value ⟩$ had an illegal value.
NE_INT: On entry, $m = ⟨ value ⟩$ .
Constraint: $m > 0$ .

On entry, $px = ⟨ value ⟩$ .
Constraint: $px \geq 8$ .

On entry, $py = ⟨ value ⟩$ .
Constraint: $py \geq 8$ .
NE_INT_2: On entry, $px = ⟨ value ⟩$ and $py = ⟨ value ⟩$ .
Constraint: $px \geq 8$ and $py \geq 8$ .
NE_INTERNAL_ERROR: An internal error has occurred in this function. Check the function call and any array sizes. If the call is correct then please contact NAG for assistance.
See Section 7.5 in the Introduction to the NAG Library CL Interface for further information.
NE_NO_LICENCE: Your licence key may have expired or may not have been installed correctly.
See Section 8 in the Introduction to the NAG Library CL Interface for further information.
NE_NOT_NON_DECREASING: On entry, $I = ⟨ value ⟩$ , $lamda [I - 1] = ⟨ value ⟩$ and $lamda [I - 2] = ⟨ value ⟩$ .
Constraint: $lamda [I - 1] \geq lamda [I - 2]$ .

On entry, $I = ⟨ value ⟩$ , $mu [I - 1] = ⟨ value ⟩$ and $mu [I - 2] = ⟨ value ⟩$ .
Constraint: $mu [I - 1] \geq mu [I - 2]$ .

7 Accuracy

Not applicable.

8 Parallelism and Performance

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

e02zac is not threaded in any implementation.

9 Further Comments

The time taken is approximately proportional to

m \times \log ((m + (px - 7) \times (py - 7)))

This function was written to sort two-dimensional data in the manner required by function e02dac. The first

9

arguments of e02zac are the same as the arguments in e02dac which have the same name.

10 Example

This example reads in data points and the intercepts of the panel sides on the

x

and

y

axes; it calls e02zac to set up the index array point; and finally it prints the data points in panel order.

e02za: FL CL CPP AD PY MB

NAG CL Interfacee02zac (dim2_​spline_​sort)

▸▿ 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 CL Interface
e02zac (dim2_spline_sort)