inst/doc/examples/blockmatrix2.R

In blockmatrix: blockmatrix: Tools to solve algebraic systems with partitioned matrices

# file blockmatrix2.R
#
# This file contains a script building a 2 x 2 block-matrix from 4 3 x 3 real matrices. 
# The lines coded in this script are lumbed in the function 'blockmatrix'
#
# author: Emanuele Cordano on 22-02-2012

#This program is free software: you can redistribute it and/or modify
#it under the terms of the GNU General Public License as published by
#the Free Software Foundation, either version 3 of the License, or
#(at your option) any later version.
#
#This program is distributed in the hope that it will be useful,
#but WITHOUT ANY WARRANTY; without even the implied warranty of
#MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
#GNU General Public License for more details.
#
#You should have received a copy of the GNU General Public License
#along with this program.  If not, see <http://www.gnu.org/licenses/>.

###############################################################################
rm(list=ls())
library(blockmatrix) 

A <- array(rnorm(9,mean=1),c(3,3))
B <- 0 #array(rnorm(9,mean=2),c(3,3))
C <- 0
D <- array(rnorm(9,mean=4),c(3,3))

list <- list(A=A,B=B,C=C,D=D)
names <- names(list)
zeros <- names(list[lapply(FUN=length,list)==1]==0)
 
names[names %in% zeros] <- "0"
M <- list[names!="0"]
X <- array(names,c(2,2))
M$value <- X

class(M) <- "blockmatrix"