permutedims for Number block arrays (#164)

* permutedims for Number block arrays

* more layout_getindex overloads

* Add mortar example to README

* test on 1.6

Author: dlfivefifty

.github/workflows/ci.yml

@@ -11,6 +11,7 @@ jobs:
       matrix:
         version:
           - '1.5'
+          - '^1.6.0-0'
           - 'nightly'
         os:
           - ubuntu-latest

Project.toml

@@ -1,6 +1,6 @@
 name = "BlockArrays"
 uuid = "8e7c35d0-a365-5155-bbbb-fb81a777f24e"
-version = "0.15"
+version = "0.15.1"
 
 [deps]
 ArrayLayouts = "4c555306-a7a7-4459-81d9-ec55ddd5c99a"
@@ -9,7 +9,7 @@ LinearAlgebra = "37e2e46d-f89d-539d-b4ee-838fcccc9c8e"
 
 [compat]
 ArrayLayouts = "0.5, 0.6"
-FillArrays = "0.11"
+FillArrays = "0.11.6"
 julia = "1.5"
 
 [extras]

README.md

@@ -8,6 +8,27 @@ A block array is a partition of an array into blocks or subarrays, see [wikipedi
 
 Secondly, it also implements two different type of block arrays that follow the `AbstractBlockArray` interface. The type `BlockArray` stores each block contiguously while the type `PseudoBlockArray` stores the full matrix contiguously. This means that `BlockArray` supports fast non copying extraction and insertion of blocks while `PseudoBlockArray` supports fast access to the full matrix to use in in for example a linear solver.
 
+A simple way to produce `BlockArray`s is via `mortar`, which combines an array of arrays into a `BlockArray`:
+```julia
+julia> mortar([randn(3), randn(4)])
+2-blocked 7-element BlockArray{Float64,1}:
+ -0.19808699390960527 
+  0.04711385377738941 
+ -0.6308529482215658  
+ ─────────────────────
+ -0.021279626465135287
+ -1.0991149020591062  
+  1.0817971931026398  
+ -0.012442892450142308
+
+julia> mortar(reshape([randn(2,2), randn(1,2), randn(2,3), randn(1,3)],2,2))
+2×2-blocked 3×5 BlockArray{Float64,2}:
+ -1.17797    0.359738   │   0.87676    -2.06495    1.74256 
+  1.54787    1.64133    │  -0.0416484  -2.00241   -0.522441
+ ───────────────────────┼──────────────────────────────────
+  0.430093  -0.0263753  │  -1.31275     0.278447  -0.139579
+```
+
 ## Documentation
 
 - [**STABLE**][docs-stable-url] — **most recently tagged version of the documentation.**

src/abstractblockarray.jl

@@ -178,3 +178,16 @@ end
 @inline Base.getindex(A::AbstractArray{T,N}, kr::Block{1}, jrs...) where {T,N} = ArrayLayouts.layout_getindex(A, kr, jrs...)
 @inline Base.getindex(A::AbstractArray{T,N}, block::Block{N}) where {T,N} = ArrayLayouts.layout_getindex(A, block)
 @inline Base.getindex(A::AbstractMatrix, kr::AbstractVector, jr::Block) = ArrayLayouts.layout_getindex(A, kr, jr)
+@inline Base.getindex(A::AbstractMatrix, kr::BlockRange{1}, jr::BlockRange{1}) = ArrayLayouts.layout_getindex(A, kr, jr)
+@inline Base.getindex(A::LayoutMatrix, kr::BlockRange{1}, jr::BlockRange{1}) = ArrayLayouts.layout_getindex(A, kr, jr)
+
+###
+# permutedims
+#
+# just use transpose for now
+###
+
+Base.permutedims(A::AbstractBlockVector{<:Number}) = transpose(A)
+Base.permutedims(A::AbstractBlockMatrix{<:Number}) = transpose(A)
+
+