include("nodecalc.jl")

order = parse(ARGS[1])
sbp = TetSBP{Float64}(degree=order, reorder=false, internal=true)
#ref_verts = [0. 1 0; 0 0 1]

xi, coords = nodecalc(sbp, true)

nnodes = sbp.numnodes
println("size(coords) = ", size(coords))
for i=1:nnodes
  println("node $i xi coordinates = ", xi[:, i])
end
print("\n")
for i=1:nnodes
  println("node $i xy coordinates = ", coords[:, i])
end

println("coords = ", coords)

