Troubleshoot tri6 scheme in Norma

[ikalash]

When running NO-SAM with a tet10 mesh, which invokes the tri6 interpolation schemes for the Schwarz BCs, there are issues/errors. With a 3-point quadrature rule, negative Jacobians are encountered for the cuboid problem. With a 4-point quadrature rule (the default and preferred scheme), there seem to be singular matrices:

ERROR: LoadError: LinearAlgebra.SingularException(25)
Stacktrace:
  [1] checknonsingular
    @ ~/Install/julia-1.11.1/share/julia/stdlib/v1.11/LinearAlgebra/src/factorization.jl:69 [inlined]
  [2] _check_lu_success
    @ ~/Install/julia-1.11.1/share/julia/stdlib/v1.11/LinearAlgebra/src/lu.jl:84 [inlined]
  [3] #lu!#182
    @ ~/Install/julia-1.11.1/share/julia/stdlib/v1.11/LinearAlgebra/src/lu.jl:92 [inlined]
  [4] lu!
    @ ~/Install/julia-1.11.1/share/julia/stdlib/v1.11/LinearAlgebra/src/lu.jl:90 [inlined]
  [5] lu!
    @ ~/Install/julia-1.11.1/share/julia/stdlib/v1.11/LinearAlgebra/src/lu.jl:89 [inlined]
  [6] _lu
    @ ~/Install/julia-1.11.1/share/julia/stdlib/v1.11/LinearAlgebra/src/lu.jl:347 [inlined]
  [7] lu(::Matrix{Float64}; kwargs::@Kwargs{})
    @ LinearAlgebra ~/Install/julia-1.11.1/share/julia/stdlib/v1.11/LinearAlgebra/src/lu.jl:341
  [8] lu
    @ ~/Install/julia-1.11.1/share/julia/stdlib/v1.11/LinearAlgebra/src/lu.jl:341 [inlined]
  [9] inv(A::Matrix{Float64})
    @ LinearAlgebra ~/Install/julia-1.11.1/share/julia/stdlib/v1.11/LinearAlgebra/src/dense.jl:993
 [10] \(A::Matrix{Float64}, J::LinearAlgebra.UniformScaling{Bool})
    @ LinearAlgebra ~/Install/julia-1.11.1/share/julia/stdlib/v1.11/LinearAlgebra/src/uniformscaling.jl:277
 [11] compute_neumann_projector(dst_model::Main.Norma.SolidMechanics, dst_bc::Main.Norma.SolidMechanicsNonOverlapSchwarzBoundaryCondition)
    @ Main.Norma ~/Norma.jl-tri6/src/ics_bcs.jl:674
 [12] initialize_bc_projectors(sim::Main.Norma.MultiDomainSimulation)
    @ Main.Norma ~/Norma.jl-tri6/src/simulation.jl:812
 [13] initialize(sim::Main.Norma.MultiDomainSimulation)
    @ Main.Norma ~/Norma.jl-tri6/src/simulation.jl:363
 [14] evolve(sim::Main.Norma.MultiDomainSimulation)
    @ Main.Norma ~/Norma.jl-tri6/src/simulation.jl:222
 [15] run(sim::Main.Norma.MultiDomainSimulation)
    @ Main.Norma ~/Norma.jl-tri6/src/Norma.jl:27
 [16] run(input_file::String)
    @ Main.Norma ~/Norma.jl-tri6/src/Norma.jl:14
 [17] top-level scope
    @ ~/Norma.jl-tri6/src/Norma.jl:39
in expression starting at /home/ikalash/Norma.jl-tri6/src/Norma.jl:7

Need to troubleshoot.