solve.HDF5Matrix

solve.HDF5Matrix

DECOMPOSITIONS

1 Description

Computes the matrix inverse of a symmetric positive-definite HDF5Matrix using Cholesky decomposition + back-substitution. Equivalent to base::solve(A) for SPD matrices.

2 Usage

solve.HDF5Matrix(...)

3 Arguments

Parameter	Description
`a`	An (square, symmetric positive-definite).
`b`	Not supported for ; must be missing.
`full_matrix`	Logical. Return full symmetric inverse. Default .
`overwrite`	Logical. Overwrite existing result. Default .
`threads`	Integer. OpenMP threads (-1 = auto).
`block_size`	Integer or NULL. Elements per block. NULL = auto.
`compression`	Integer (0-9) or NULL. gzip compression level for the result dataset. NULL uses the global option set by (default 6). Use to disable compression (faster for benchmarks).

4 Value

containing the matrix inverse.

5 Examples

\donttest{
tmp <- tempfile(fileext = ".h5")

X  <- hdf5_create_matrix(tmp, "data/X", data = matrix(rnorm(10000), 100, 100))

X   <- hdf5_matrix(tmp, "data/X")
AtA <- crossprod(X)              # HDF5Matrix, square SPD
inv <- solve(AtA)               # inverse of AtA

hdf5_close_all()
unlink(tmp)
}

---
title: "solve.HDF5Matrix"
subtitle: "solve.HDF5Matrix"
---

<span class="category-badge decompositions">DECOMPOSITIONS</span>

## Description

Computes the matrix inverse of a symmetric positive-definite HDF5Matrix
using Cholesky decomposition + back-substitution.  Equivalent to
`base::solve(A)` for SPD matrices.

## Usage

```r
solve.HDF5Matrix(...)
```
## Arguments

::: {.param-table}
| Parameter | Description |
|-----------|-------------|
| `a` | An \code{HDF5Matrix} (square, symmetric positive-definite). |
| `b` | Not supported for \code{HDF5Matrix}; must be missing. |
| `full_matrix` | Logical. Return full symmetric inverse.  Default \code{TRUE}. |
| `overwrite` | Logical. Overwrite existing result.  Default \code{FALSE}. |
| `threads` | Integer. OpenMP threads (-1 = auto). |
| `block_size` | Integer or NULL. Elements per block.  NULL = auto. |
| `compression` | Integer (0-9) or NULL. gzip compression level for the   result dataset.  NULL uses the global option set by   \code{\link{hdf5matrix_options}} (default 6).  Use \code{0} to disable   compression (faster for benchmarks). |
:::

## Value

::: {.return-value}
\code{HDF5Matrix} containing the matrix inverse.

:::
## Examples
```{r}
#| eval: false
#| warning: false

\donttest{
tmp <- tempfile(fileext = ".h5")

X  <- hdf5_create_matrix(tmp, "data/X", data = matrix(rnorm(10000), 100, 100))

X   <- hdf5_matrix(tmp, "data/X")
AtA <- crossprod(X)              # HDF5Matrix, square SPD
inv <- solve(AtA)               # inverse of AtA

hdf5_close_all()
unlink(tmp)
}

```