--- a/src/Grids/tensor_grid.jl	Tue Sep 17 09:12:52 2024 +0200
+++ b/src/Grids/tensor_grid.jl	Tue Sep 17 11:19:08 2024 +0200
@@ -84,17 +84,28 @@
     return TensorGrid(new_grids...)
 end
 
+
+function boundary_indices(g::TensorGrid{T,1} where T, id::TensorGridBoundary)
+    return boundary_indices(g.grids[grid_id(id)], boundary_id(id))
+end
 function boundary_indices(g::TensorGrid, id::TensorGridBoundary)
-    per_grid_ind = map(g.grids) do g
-        ntuple(i->:, ndims(g))
-    end
+    all_indices = map(eachindex, g.grids)
 
     local_b_ind = boundary_indices(g.grids[grid_id(id)], boundary_id(id))
-    b_ind = Base.setindex(per_grid_ind, local_b_ind, grid_id(id))
+
+    b_ind = Base.setindex(all_indices, local_b_ind, grid_id(id))
+
+    return view(_combine_indices(all_indices...), LazyTensors.concatenate_tuples(bla.(b_ind)...)...)
+end
 
-    return LazyTensors.concatenate_tuples(b_ind...)
+# function _combine_indices(Is::Vararg{Union{Int, <:AbstractRange}})
+function _combine_indices(Is...)
+    return CartesianIndices(LazyTensors.concatenate_tuples(bla.(Is)...))
 end
 
+bla(a) = (a,)
+bla(a::CartesianIndices) = a.indices
+
 function combined_coordinate_vector_type(coordinate_types...)
     combined_coord_length = mapreduce(_ncomponents, +, coordinate_types)
     combined_coord_type = mapreduce(eltype, promote_type, coordinate_types)