lasq6

slasq6

Functions

void slasq6(
    const INT           i0,
    const INT           n0,
          f32* restrict Z,
    const INT           pp,
          f32*          dmin,
          f32*          dmin1,
          f32*          dmin2,
          f32*          dn,
          f32*          dnm1,
          f32*          dnm2
);

void slasq6(const INT i0, const INT n0, f32 *restrict Z, const INT pp, f32 *dmin, f32 *dmin1, f32 *dmin2, f32 *dn, f32 *dnm1, f32 *dnm2)

SLASQ6 computes one dqd (shift equal to zero) transform in ping-pong form, with protection against underflow and overflow.

Parameters

in

i0

First index (0-based).

in

n0

Last index (0-based).

inout

Z

Double precision array, dimension (4*N). Z holds the qd array. EMIN is stored in Z(4*N0) to avoid an extra argument.

in

pp

PP=0 for ping, PP=1 for pong.

out

dmin

Minimum value of d.

out

dmin1

Minimum value of d, excluding D(N0).

out

dmin2

Minimum value of d, excluding D(N0) and D(N0-1).

out

dn

d(N0), the last value of d.

out

dnm1

d(N0-1).

out

dnm2

d(N0-2).

dlasq6

Functions

void dlasq6(
    const INT           i0,
    const INT           n0,
          f64* restrict Z,
    const INT           pp,
          f64*          dmin,
          f64*          dmin1,
          f64*          dmin2,
          f64*          dn,
          f64*          dnm1,
          f64*          dnm2
);

void dlasq6(const INT i0, const INT n0, f64 *restrict Z, const INT pp, f64 *dmin, f64 *dmin1, f64 *dmin2, f64 *dn, f64 *dnm1, f64 *dnm2)

DLASQ6 computes one dqd (shift equal to zero) transform in ping-pong form, with protection against underflow and overflow.

Parameters

in

i0

First index (0-based).

in

n0

Last index (0-based).

inout

Z

Double precision array, dimension (4*N). Z holds the qd array. EMIN is stored in Z(4*N0) to avoid an extra argument.

in

pp

PP=0 for ping, PP=1 for pong.

out

dmin

Minimum value of d.

out

dmin1

Minimum value of d, excluding D(N0).

out

dmin2

Minimum value of d, excluding D(N0) and D(N0-1).

out

dn

d(N0), the last value of d.

out

dnm1

d(N0-1).

out

dnm2

d(N0-2).