Expressing Modal C0 shapefuns as linear combination

Author: ericneiva

src/ReferenceFEs/ModalC0RefFEs.jl

@@ -69,7 +69,7 @@ function ModalC0RefFE(
 end
 
 function get_orders(reffe::GenericRefFE{ModalC0,D}) where{D}
-  get_orders(get_shapefuns(reffe))
+  get_orders(reffe.prebasis)
 end
 
 function compute_reffe_data(::Type{T},

src/ReferenceFEs/ReferenceFEInterfaces.jl

@@ -488,7 +488,7 @@ struct GenericRefFE{T,D} <: ReferenceFE{D}
       conformity,
       metadata,
       face_dofs,
-      shapefuns)
+      linear_combination(Eye{Int}(ndofs),shapefuns)) # Trick to be able to eval dofs af shapefuns
   end
 end