source edits

Author: eddelbuettel

_posts/2014-03-21-simple-array-class.md

@@ -28,7 +28,10 @@ partial model matrices where the entries in the vector `X%*%theta` (with
 `nd*...*n2*n1` entries) are the same as the entries in the array with
 dimension `c(n1, n2, ..., nd)` returned by `RH(Xd, ... , RH(X2, RH(X1,
 Theta))...)`.  
-`Theta` is an array with dimensions `c(c1, c2, ..., cd)` containing `theta` and `RH(X, A)` -- the "rotated H-transform" -- is an operation generalizing transposed pre-multiplication `t(X %*% A)` of a matrix `A` by a matrix `X` to the case of higher dimensional array-valued `A`. 
+`Theta` is an array with dimensions `c(c1, c2, ..., cd)` containing `theta`
+and `RH(X, A)` -- the "rotated H-transform" -- is an operation generalizing
+transposed pre-multiplication `t(X %*% A)` of a matrix `A` by a matrix `X` to
+the case of higher dimensional array-valued `A`.  
 
 The code below implements a simple array class for numeric arrays and the
 rotated H-transform in `RcppArmadillo` and compares the performance to both
@@ -280,7 +283,7 @@ benchmark(
 
 <pre class="output">
                                          test elapsed relative
-1 array(X %*% theta_vec, dim = c(n1, n2, n3))  27.159  198.241
+1 array(X %*% theta_vec, dim = c(n1, n2, n3))  27.252  198.920
 3   Reduce(RH, list(X3, X2, X1), Theta, TRUE)   0.137    1.000
 2 Reduce(RH_r, list(X3, X2, X1), Theta, TRUE)   0.325    2.372
 </pre>

src/2014-03-21-simple-array-class.Rmd

@@ -26,7 +26,10 @@ partial model matrices where the entries in the vector `X%*%theta` (with
 `nd*...*n2*n1` entries) are the same as the entries in the array with
 dimension `c(n1, n2, ..., nd)` returned by `RH(Xd, ... , RH(X2, RH(X1,
 Theta))...)`.  
-`Theta` is an array with dimensions `c(c1, c2, ..., cd)` containing `theta` and `RH(X, A)` -- the "rotated H-transform" -- is an operation generalizing transposed pre-multiplication `t(X %*% A)` of a matrix `A` by a matrix `X` to the case of higher dimensional array-valued `A`. 
+`Theta` is an array with dimensions `c(c1, c2, ..., cd)` containing `theta`
+and `RH(X, A)` -- the "rotated H-transform" -- is an operation generalizing
+transposed pre-multiplication `t(X %*% A)` of a matrix `A` by a matrix `X` to
+the case of higher dimensional array-valued `A`.  
 
 The code below implements a simple array class for numeric arrays and the
 rotated H-transform in `RcppArmadillo` and compares the performance to both