ring_partition: Ring partition for the block-LEnKPF

Description Usage Arguments Examples

View source: R/utils_geom.R

Description

Partition the observations in blocks for assimilation with block-LEnKPF. Takes into account the ring geometry. We assume that R is diagonal and don't take it into account

Usage

1
ring_partition(H, taper, block.size)

Arguments

H

the observation operator

taper

the tapering matrix (qxq)

block.size

chunk size in domain space to take into account for one block

Examples

1
2
3
4
5
6
7
8
ndim <- 100
taper <- ring_GC(ndim, 5)
H <- diag(ndim)
partition <- ring_partition(H, taper, 5)
partition <- ring_partition( matrix( c(rep(0,ndim-1),1),1,ndim), taper, 5)
H <- matrix(0, 1, ndim)
H[1,round(ndim/2)] <- 1
partition <- ring_partition(H, taper, 5)

robertsy/assimilr documentation built on May 27, 2019, 10:33 a.m.