Skip to contents

Fast version of Matrix::.bdiag() – for the case of many (k x k) matrices: Copyright (C) 2016 Martin Maechler, ETH Zurich

Usage

bdiag_m(lmat)

Arguments

lmat

list(<mat1>, <mat2>, ....., <mat_N>) where each mat_j is a k x k 'matrix'

Value

a sparse (Nk x Nk) matrix of class Matrix::dgCMatrix.