Moment fitting with simplices

[pmartorell]

Dear @ericneiva ,

the moment fitting machinery is not working for simplex background meshes.

Since the JacobiPolynomialBasis are a tensor product, the following line tries to invert a rectangular matrix.

GridapEmbedded.jl/src/MomentFittedQuadratures/MomentFittedQuadratures.jl

Line 322 in 4cf152d

change = transpose(inv(change))

I guess that we should use pinv instead. What do you think?

Here a MWE:

using GridapEmbedded
using GridapEmbedded.LevelSetCutters
using GridapEmbedded.Interfaces
using Gridap
using Gridap.ReferenceFEs
geo = sphere(0.5,x0=Point(0.0,0.0))
domain = (-1,1,-1,1)
partition = (5,5)
model = CartesianDiscreteModel(domain,partition) |> simplexify
cutgeo = cut(model,geo)
act = Triangulation(cutgeo,ACTIVE)
Measure(act,Quadrature(momentfitted,cutgeo,2))