fixed dimension of nodevals output

WIAS-PDELib · Nov 4, 2024 · 463acd2 · 463acd2
1 parent 711c691
commit 463acd2
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/interpolations.jl b/src/interpolations.jl
@@ -819,7 +819,7 @@ function nodevalues(source::FEVectorBlock{T, Tv, Ti, FEType, APT}, operator::Typ
 	if abs
 		nvals = 1
 	else
-		xdim = size(source.FES.dofgrid[Coordinates], 2)
+		xdim = size(source.FES.dofgrid[Coordinates], 1)
 		ncomponents = get_ncomponents(eltype(source.FES))
 		nvals = Length4Operator(operator, xdim, ncomponents)
 	end