Comments (2)
You could average the solution over each cell. There may be a cleaner way, but the internal routine Trixi.wrap_array
will reshape the solution vector into an array of size [num_fields, N+1, N+1, num_elements]
in 2D.
The best way to average the solution is to use the quadrature nodes available in dg.basis
or solver.basis
to compute integral averages of the solution over each cell, but you could take the arithmetic average over each cell u[i, :, :, cell_index]
as well.
from trixi.jl.
If you use trixi2vtk
to convert the solution files it, by default, reinterpolates to the solution from the nnodes
Legendre-Gauss-Lobatto points (which is polydeg+1
) onto 2*nnodes
that are uniform. This would explain your situation.
-
If you want to change the number of uniform nodes onto which you reinterpolate then you can pass the additional argument
nvisnodes
to yourtrixi2vtk
call, e.g.,trixi2vtk("out/solution_files", nvisnodes=4)
. -
To deactivate this reinterpolation completely and make the solution on the original solver nodes available you can set additional argument
reinterpolate
to be false in yourtrixi2vtk
call, e.g.,trixi2vtk("out/solution_files", reinterpolate=false)
.
from trixi.jl.
Related Issues (20)
- Add Positivity Preserving Feature for GaussSBP nodes HOT 2
- `AnalysisSurfaceIntegral` not ready for AMR HOT 1
- Filter out dependabot commits in release changelog HOT 2
- Parallel p4est simulation crashes when ranks have zero elements HOT 3
- Missing tests for boundary conditions HOT 1
- Finalize formatting
- Release v0.9 HOT 1
- RFC: alternative computation of source terms based on indices HOT 2
- The GLL method cannot perform exact numerical integration for DGSEM ? HOT 2
- A few questions about the formulas in the flux differencing tutorial HOT 2
- Spatially-varying number of fields HOT 2
- hamilton-jacobi / level-set equations HOT 1
- Unused parameter in `AnalysisSurfaceIntegral` constructor
- `trixi_include` does not expose variables inside a Pluto notebook HOT 4
- Viscous terms using FDSBP discretizations
- Variable namings within `MPIInterfaceContainer2D`, `MPIL2MortarContainer2D` HOT 1
- saving solution while using T8code meshes HOT 1
- Remove `normal_direction_ll` in nonconservative terms HOT 8
- Subcell limiting: Improving interface for `get_boundary_outer_state()`
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from trixi.jl.