Merge pull request #259 from xtalax/hotfx

hotfix

Author: xtalax

docs/src/tutorials/PIDE.md

@@ -55,7 +55,7 @@ plot(sol[x], transpose(solu))
 To have an integral over the whole domain, be sure to wrap the integral in an auxiliary variable.
 Due to a limitation, the whole domain integral needs to have the same arguments as the integrand, but is constant in x. To use it in an equation one dimension lower, use a boundary value like integral(t, 0.0)
 
-```@example integrals2
+``` julia
 using MethodOfLines, ModelingToolkit, DomainSets, OrdinaryDiffEq, Plots
 
 @parameters t, x

src/discretization/generate_bc_eqs.jl

@@ -209,7 +209,7 @@ end
 
 @inline function generate_corner_eqs!(bceqs, s, interiormap, N, u)
     interior = interiormap.I[interiormap.pde[u]]
-    N == 0 || ndims(u, s) == 0 && return
+    ndims(u, s) == 0 && return
     sd(i, j) = selectdim(interior, j, i)
     domain = setdiff(s.Igrid[u], interior)
     II1 = unitindices(N)

src/system_parsing/interior_map.jl

@@ -59,7 +59,7 @@ end
 function generate_interior(lower, upper, u, s, disc::MOLFiniteDifference{G,D}) where {G, D<:ScalarizedDiscretization}
     args = remove(arguments(u), s.time)
 
-    ret = collect(s.Igrid[u][[(1+lower[x2i(s, u, x)]:length(s.grid[x])-upper[x2i(s, u, x)]) for x in args]...])
+    ret = s.Igrid[u][[(1+lower[x2i(s, u, x)]:length(s.grid[x])-upper[x2i(s, u, x)]) for x in args]...]
     return ret
 end