stable release
updated documentation
replace all previous functions with fatdiag
add fatdig <- function
add support for replacement other than 0 (zero)
break out rectangles to separate functions: rectangle_matrix() and minus_rectangle_matrix()
block_matrix(): for extracting block diagonals
minus_block_matrix(): for dropping block diagonals
add support for non-square matrices