I have the following code in Julia

```
julia> x = [1 2];
julia> y = reshape([3,4], (2, 1));
julia> x .+ y
2×2 Array{Int64,2}:
4 5
5 6
```

In julia, the `.`

is broadcasting and is equivalent to

```
julia> x_expanded = vcat(x, x)
2×2 Array{Int64,2}:
1 2
1 2
julia> y_expanded = hcat(y, y)
2×2 Array{Int64,2}:
3 3
4 4
julia> x_expanded + y_expanded
2×2 Array{Int64,2}:
4 5
5 6
```

How would I do the same in R? After some benchmarking in Julia the broadcasted summation looks like the fastest option for my code. I am wondering if there is some sort of highly optimized equivalent in R.

```
r$> x = matrix(c(1, 2), nrow = 1, ncol = 2)
r$> y = matrix(c(3, 4), nrow = 2, ncol = 1)
r$> x + y
Error in x + y : non-conformable arrays
```

Can someone point me in the right direction?