Fix IntegrationMap on empty subds with physical quads

gridap · Jul 8, 2024 · 78a70fc · 78a70fc
1 parent 9079dd1
commit 78a70fc
Showing 1 changed file with 18 additions and 0 deletions.
diff --git a/src/Fields/FieldsInterfaces.jl b/src/Fields/FieldsInterfaces.jl
@@ -491,6 +491,15 @@ end
 
 struct IntegrationMap <: Map end
 
+function return_value(k::IntegrationMap,ax::AbstractVector,w)
+  if size(ax,1) == length(w)
+    evaluate(k,ax,w)
+  else
+    T = typeof( testitem(ax)*testitem(w) + testitem(ax)*testitem(w) )
+    zero(T)
+  end
+end
+
 function evaluate!(cache,k::IntegrationMap,ax::AbstractVector,w)
   T = typeof( testitem(ax)*testitem(w) + testitem(ax)*testitem(w) )
   z = zero(T)
@@ -513,6 +522,15 @@ function evaluate!(cache,k::IntegrationMap,aq::AbstractVector,w,jq::AbstractVect
   z
 end
 
+function return_value(k::IntegrationMap,ax::AbstractArray,w)
+  if size(ax,1) == length(w)
+    evaluate(k,ax,w)
+  else
+    c = return_cache(k,ax,w)
+    c.array
+  end
+end
+
 function return_cache(k::IntegrationMap,ax::AbstractArray,w)
   T = typeof( testitem(ax)*testitem(w) + testitem(ax)*testitem(w) )
   r = zeros(T,size(ax)[2:end])