So strange things happen for me with this. First, suppose I apply the vectorized funciton, but on a single point:

```
julia> r1 = [0., 2.]
2-element Array{Float64,1}:
0.0
2.0
julia> f(r1)
0.9092974268256817
julia> f.(r1)
ERROR: BoundsError
Stacktrace:
[1] getindex at ./number.jl:38 [inlined]
[2] (::##1#2)(::Float64) at ./REPL[2]:1
[3] macro expansion at ./broadcast.jl:153 [inlined]
[4] macro expansion at ./simdloop.jl:73 [inlined]
[5] macro expansion at ./broadcast.jl:147 [inlined]
[6] _broadcast!(::##1#2, ::Array{Float64,1}, ::Tuple{Tuple{Bool}}, ::Tuple{Tuple{Int64}}, ::Array{Float64,1}, ::Tuple{}, ::Type{Val{0}}, ::CartesianRange{CartesianIndex{1}}) at ./broadcast.jl:139
[7] broadcast_t(::Function, ::Type{T} where T, ::Tuple{Base.OneTo{Int64}}, ::CartesianRange{CartesianIndex{1}}, ::Array{Float64,1}) at ./broadcast.jl:268
[8] broadcast_c at ./broadcast.jl:314 [inlined]
[9] broadcast(::Function, ::Array{Float64,1}) at ./broadcast.jl:434
```

But I think I see what my issue was. I had been generating an Nx2 Array{Float64,2} for my input as

```
julia> x = linspace(0,1,5);
julia> y=x
0.0:0.25:1.0
julia> r = hcat(x,y)
5×2 Array{Float64,2}:
0.0 0.0
0.25 0.25
0.5 0.5
0.75 0.75
1.0 1.0
```

instead of an N array of points. How do I convert between the data types?