Re: [deal.II] Re: DG explicit time integration for linear advection equation with MeshWorker (suggestions)

2019-09-17 Thread Praveen C
> 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

[deal.II] Re: DG explicit time integration for linear advection equation with MeshWorker (suggestions)

2019-09-17 Thread vachan potluri
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

Re: [deal.II] Instantiation of Vector< Sacado::Fad::DFad >

2019-09-17 Thread Jean-Paul Pelteret
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

Re: [deal.II] DG explicit time integration for linear advection equation with MeshWorker (suggestions)

2019-09-17 Thread Praveen C
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

Re: [deal.II] Instantiation of Vector< Sacado::Fad::DFad >

2019-09-17 Thread Doug
> > 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

Re: [deal.II] Problem warning in eclipse: invalid argument for the distance function

2019-09-17 Thread Daniel Arndt
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

Re: [deal.II] Cmake-GUI error

2019-09-17 Thread Wolfgang Bangerth
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: > > >

Re: [deal.II] Instantiation of Vector< Sacado::Fad::DFad >

2019-09-17 Thread Wolfgang Bangerth
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,

Re: [deal.II] Instantiation of Vector< Sacado::Fad::DFad >

2019-09-17 Thread Doug
> > 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? > >

[deal.II] Instantiation MappingFEField with hp::DoFHandler

2019-09-17 Thread Doug
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

[deal.II] Re: DG explicit time integration for linear advection equation with MeshWorker (suggestions)

2019-09-17 Thread Doug
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

[deal.II] Problem warning in eclipse: invalid argument for the distance function

2019-09-17 Thread Brian Zhang
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"

[deal.II] DG explicit time integration for linear advection equation with MeshWorker (suggestions)

2019-09-17 Thread vachan potluri
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,