indexing_mod Module

private elemental function fuse_symm_idx_int32(p, q) result(pq)

Fuse an index

create a composite index out of two indices, assuming they are unordered i.e. their ordering does not matter.

The ordering looks like this:

1 2 4 7 3 5 8 6 9 10

Arguments

Return Value integer(kind=int32)

1d-array index assuming the array is symmetric w.r. p<->q

private elemental function fuse_symm_idx_int64(p, q) result(pq)

Fuse an index

create a composite index out of two indices, assuming they are unordered i.e. their ordering does not matter.

The ordering looks like this:

1 2 4 7 3 5 8 6 9 10

Arguments

Return Value integer(kind=int64)

1d-array index assuming the array is symmetric w.r. p<->q

private elemental subroutine inv_fuse_symm_idx_int32(pq, p, q)

Invert the fuse_symm_idx routine

Returns index pair x <= y for the fused index xy. This assumes that the underlying matrix is symmetric, i.e. (x, y) and (y, x) map to the same xy.

Arguments

private elemental subroutine inv_fuse_symm_idx_int64(pq, p, q)

Invert the fuse_symm_idx routine

Returns index pair x <= y for the fused index xy. This assumes that the underlying matrix is symmetric, i.e. (x, y) and (y, x) map to the same xy.

Arguments

private elemental function fuse_idx_int32(p, q, n_dim, fortran) result(pq)

Fuse p, q into one contiguous index

If fortran is true, then this assumes columnwise Fortran order, i.e. the index looks like this:

1 5 2 6 3 7 4 8

Otherwise it is C-order, i.e. row major (but still one-indexed).

Arguments

Return Value integer(kind=int32)

contiguous 1d-array index

private elemental function fuse_idx_int64(p, q, n_dim, fortran) result(pq)

Fuse p, q into one contiguous index

If fortran is true, then this assumes columnwise Fortran order, i.e. the index looks like this:

1 5 2 6 3 7 4 8

Otherwise it is C-order, i.e. row major (but still one-indexed).

Arguments

Return Value integer(kind=int64)

contiguous 1d-array index

private elemental subroutine inv_fuse_idx_int32(pq, n_dim, p, q, fortran)

Fuse p, q into one contiguous index

This assumes columnwise Fortran order, i.e. the index looks like this:

1 5 2 6 3 7 4 8

Arguments

private elemental subroutine inv_fuse_idx_int64(pq, n_dim, p, q, fortran)

Fuse p, q into one contiguous index

This assumes columnwise Fortran order, i.e. the index looks like this:

1 5 2 6 3 7 4 8

Arguments

private elemental function new_fuse_symm_idx_int32(p, q, n_dim) result(pq)

Fuse an index

create a composite index out of two indices, assuming they are unordered i.e. their ordering does not matter.

The ordering looks like this:

1 2 3 4 5 6 7 8 9 10

Arguments

Return Value integer(kind=int32)

1d-array index assuming the array is symmetric w.r. p<->q

private elemental function new_fuse_symm_idx_int64(p, q, n_dim) result(pq)

Fuse an index

create a composite index out of two indices, assuming they are unordered i.e. their ordering does not matter.

The ordering looks like this:

1 2 3 4 5 6 7 8 9 10

Arguments

Return Value integer(kind=int64)

1d-array index assuming the array is symmetric w.r. p<->q

private elemental function triangle_number_int32(n) result(res)

Arguments

Return Value integer(kind=int32)

private elemental function triangle_number_int64(n) result(res)

Arguments

Return Value integer(kind=int64)

private elemental function inv_triangle_number_int32(n) result(res)

Arguments

Return Value integer(kind=int32)

private elemental function inv_triangle_number_int64(n) result(res)

Arguments

Return Value integer(kind=int64)