Export of super parametric elements #167
Comments
julia> grid = generate_grid(QuadraticTriangle, (5, 5));
julia> dh = DofHandler(grid);
julia> push!(dh, :u, 1, Lagrange{2,RefTetrahedron,1}());
julia> close!(dh);
julia> u = collect(1:ndofs(dh));
julia> vtk = vtk_grid("test", grid);
julia> vtk_point_data(vtk, dh, u);
ERROR: BoundsError: attempt to access 3-element Array{Int64,1} at index [4]
Stacktrace: |
|
Always just export it as linear? |
|
Yes, or if both fields are quadratic I guess we can export that |
fredrikekre
changed the title
7 tasks
fredrikekre
added a commit
that referenced
this issue
May 10, 2023
This patch deprecates `reshape_to_nodes` in favor of the new `evaluate_at_grid_nodes` function which more or less does the same thing. `evaluate_at_grid_nodes` uses the interpolation of the field to evaluate the approximated function in the cell coordinates (by using `reference_coordinates` of the cell as quadrature points). Fixes #167.
fredrikekre
added a commit
that referenced
this issue
May 10, 2023
This patch deprecates `reshape_to_nodes` in favor of the new `evaluate_at_grid_nodes` function which more or less does the same thing. `evaluate_at_grid_nodes` uses the interpolation of the field to evaluate the approximated function in the cell coordinates (by using `reference_coordinates` of the cell as quadrature points). This is much simpler, and more correct, compared to figuring out which Dofs live on which node etc. In particular, this now support non-standard interpolations which may not even have dofs in the nodes. Fixes #167.
fredrikekre
added a commit
that referenced
this issue
May 11, 2023
This patch deprecates `reshape_to_nodes` in favor of the new `evaluate_at_grid_nodes` function which more or less does the same thing. `evaluate_at_grid_nodes` uses the interpolation of the field to evaluate the approximated function in the cell coordinates (by using `reference_coordinates` of the cell as quadrature points). This is much simpler, and more correct, compared to figuring out which Dofs live on which node etc. In particular, this now support non-standard interpolations which may not even have dofs in the nodes. Fixes #167.
fredrikekre
added a commit
that referenced
this issue
May 11, 2023
This patch deprecates `reshape_to_nodes` in favor of the new `evaluate_at_grid_nodes` function which more or less does the same thing. `evaluate_at_grid_nodes` uses the interpolation of the field to evaluate the approximated function in the cell coordinates (by using `reference_coordinates` of the cell as quadrature points). This is much simpler, and more correct, compared to figuring out which Dofs live on which node etc. In particular, this now support non-standard interpolations which may not even have dofs in the nodes. Fixes #167.
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
I don't think it works
.
The text was updated successfully, but these errors were encountered: