@@ -124,28 +124,18 @@ end
124
124
125
125
# ###########################################################################
126
126
"""
127
- $(SIGNATURES )
127
+ nodeflux(system, F,U; data )
128
128
129
129
Reconstruction of an edge function as vector function on the nodes of the
130
130
triangulation. The result can be seen as a piecewiesw linear vector
131
131
function in the FEM space spanned by the discretization nodes exhibiting
132
132
the flux density.
133
133
134
134
The reconstruction is based on the "magic formula"
135
- R. Eymard, T. Gallouet , R. Herbin, IMA Journal of Numerical Analysis (2006)
136
- 26, 326−353, Lemma 2.4 .
135
+ R. Eymard, T. Gallouët , R. Herbin, IMA Journal of Numerical Analysis (2006)
136
+ 26, 326−353, Lemma 2.4 (also: [hal.science/hal-00004840](https://hal.science/hal-00004840) ) .
137
137
138
- The return value is a `dim x nspec x nnodes` vector. The flux of species i
139
- can e.g. plotted via GridVisualize.vectorplot.
140
-
141
- Example:
142
- ```julia
143
- ispec=3
144
- vis=GridVisualizer(Plotter=Plotter)
145
- scalarplot!(vis,grid,solution[ispec,:],clear=true,colormap=:summer)
146
- vectorplot!(vis,grid,nf[:,ispec,:],clear=false)
147
- reveal(vis)
148
- ```
138
+ The return value is a `dim x nspec x nnodes` vector.
149
139
150
140
CAVEAT: there is a possible unsolved problem with the values at domain
151
141
corners in the code. Please see any potential boundary artifacts as a manifestation
@@ -206,10 +196,24 @@ function nodeflux(
206
196
end
207
197
208
198
"""
209
- $(SIGNATURES )
199
+ nodeflux(system, U; data )
210
200
211
201
Reconstruction of the edge flux as vector function on the nodes of the
212
202
triangulation. See [`nodeflux(system,F,U;data)`](@ref).
203
+
204
+ The flux of species i
205
+ can e.g. plotted via GridVisualize.vectorplot.
206
+
207
+ Example:
208
+ ```julia
209
+ ispec=3
210
+ vis=GridVisualizer(Plotter=Plotter)
211
+ scalarplot!(vis,grid,solution[ispec,:],clear=true,colormap=:summer)
212
+ vectorplot!(vis,grid,nf[:,ispec,:],clear=false)
213
+ reveal(vis)
214
+ ```
215
+
216
+
213
217
"""
214
218
function nodeflux (
215
219
system:: AbstractSystem{Tv, Tc, Ti, Tm} , U:: AbstractArray{Tu, 2} ;
0 commit comments