pbdDMAT (version 0.4-2)

Hilbert: Generate Hilbert Matrices

Description

Methods for constructing Hilbert matrices: H[i,j] = 1/(i+j-1)

Usage

Hilbert(n, type = "matrix", ..., bldim = .pbd_env$BLDIM,
  ICTXT = .pbd_env$ICTXT)

Arguments

n

number of rows and columns.

type

"matrix" or "ddmatrix".

...

Additional arguments.

bldim

blocking dimension.

ICTXT

BLACS context number.

Value

Returns a matrix or a distributed matrix.

Details

This constructs the square Hilbert matrix of order n. The return is either a matrix or a distributed matrix depending on the argument type=.

Examples

# NOT RUN {
# Save code in a file "demo.r" and run with 2 processors by
# > mpiexec -np 2 Rscript demo.r

library(pbdDMAT, quiet = TRUE)
init.grid()

dx <- Hilbert(100, type="ddmatrix")

print(dx)

finalize()
# }
# NOT RUN {
# }