search

OptimalDesignLab / SummationByParts.jl

Other
6 stars 4 forks source link

Interpolation Operator Accuracy #15

Open JaredCrean2 opened 8 years ago

JaredCrean2 commented 8 years ago

Using this code:

      sbp = TetSBP{Float64}(degree=order, reorder=false, internal=false)
      ref_verts = sbp.vtx
      println("ref_verts = \n", ref_verts)
      interp_op = SummationByParts.buildinterpolation(sbp, ref_verts.')
      println("interp_op = \n", interp_op)

produces an interpolation operator that is not the identity matrix for a SBP Gamma element of degree 3 or greater.