> I now have one more question. In assemble_face_term function of step-33, the
> normal numerical flux is calculated. This function is called for every
> cell-neighbor pair from the assemble_system function. This means, if I am not
> wrong, at every interior face quadrature point, the
Doug and Praveen,
Thanks for your answers. I had a look at step-33. As far as I understand,
although the looping through cells is not through MeshWorker, the assembly
is still global! So, for a non-cartesian mesh, I think you are suggesting
using such a loop over all cells to calculate local
H Doug and Wolfgang,
At the time when we introduced the AD framework I took a stab at quickly adding
support for to our Vector class. To summarise, I hit the same issue and
couldn’t find an elegant solution past to get past it (so hopefully this is the
only outstanding problem to add support
If you are working only on Cartesian grids, these matrices will be same for all
cells, except maybe for scalar factors. In this case, you can first compute
them on a reference cell. For assembly, you loop over cells and perform a small
matrix-vector product on each cell. Similar things can be
>
> so 'a' is of type Sacado::Fad::DFad. It then needs to call
>
>numbers::is_finite (Sacado::Fad::DFad), which doesn't exist.
> It probably just tries to go through all of the overloads of
> numbers::is_finite() and wants to see whether it can convert the
> argument Sacado::Fad::DFad to
Brian,
The warnings you are seeing come from a static analyzer in Eclipse and not
from the compiler. It seems like Eclipse does not have full access to the
source code
or doesn't understand it correctly. E.g., passing a non-const variable to a
function that doesn't modify its parameter, is never
On 9/17/19 4:13 PM, GaryR wrote:
> I keep getting the following error when running cmake-gui.
> |
>
>
> CMakeErrorat
> /home/gary/Deal.II/share/deal.II/macros/macro_deal_ii_setup_target.cmake:64(INCLUDE):
> INCLUDE could notfind load file:
>
>
>
On 9/17/19 6:58 PM, Doug wrote:
>
> /home/ddong/Libraries/dealii/include/deal.II/base/numbers.h:583:3:
> note: no known conversion for argument 1 from ‘const
> Sacado::Fad::DFad’ to ‘const std::complex&’
> In file included from
> /home/ddong/Libraries/dealii/include/deal.II/base/cuda.h:21:0,
>
> What I meant is: Can you show the compiler error message that illustrates
> where the assertion is located, what the template arguments are, how it
> came
> that we called that function with these template arguments, etc?
>
>
Hello again,
I am trying to use MappingFEField to represent my high-order mesh points
(which will then be deformed through optimization). Similarly to this
thread High order mesh from Gmsh
Hello Vachan,
What you are describing is very similar to Hesthaven's framework of build
matrix operators and apply them to your degrees of freedom. The main
advantage of this operator form is to have a common mass (if linear
elements), differentiation, flux, and boundary operators for all your
Hello, everyone,
I am pretty new to deal.ii.
I have studied the tutorial step-1, where it is recommended to learn a
debugger ASAP.
So I put it into the Eclipse IDE.
However, I found there is an error as shown in the bellow picture. It is
related to that the argument of the "distance"
Hello all,
I am a beginner in dealii. I want to solve a linear, transient advection
equation explicitly in two dimensions using DG. The resulting discrete
equation will have a mass matrix as the system matrix and a sum of terms
which depend on previous solution (multiplied by mass,
13 matches
Mail list logo