Re: [deal.II] Automatic Differentiation of the Finite Element Jacobian

Dear Doug, > Looking at the current implementation in fe_values, mapping_q_generic, and > local integrators, I see that it looks relatively easy to template the value > type since it is already partially templated. > > Is there a reason why those functions were not templates? No, its likely

[deal.II] Automatic Differentiation of the Finite Element Jacobian

Hello, We had chosen to use deal.ii since most of it looked templated in terms of accepting the float type. However, I noticed that the FEValuesBase returns a 'double', instead of a generic 'value_type' for the Jacobian term. Our goal would be to differentiate through the entire residual (and

Re: [deal.II] Mistake in Step 16?

On 2/5/19 3:29 AM, ky...@math.uh.edu wrote: > > Hi All, > > I was working through example step 16 and I think I found a mistake. > > In| > LaplaceProblem::setup_system() > | > > |We have the following| > > || > | > | > std::set         dirichlet_boundary_ids; > Functions::ZeroFunction >

[deal.II] Mistake in Step 16?

Hi All, I was working through example step 16 and I think I found a mistake. In LaplaceProblem::setup_system() We have the following std::set dirichlet_boundary_ids; Functions::ZeroFunction