On 8/31/22 11:09, Lucas Myers wrote:
Thanks for the help! I printed out the number of constraints after each
of your suggestions, and after setting the hanging node constraints the
number of constraints is 0 (as it should be, I'm working with a 1D mesh
that has been globally refined some number
Dear all,
I want to differentiate my pde solution U with respect to some design
variables E, that is, computing the Jacobian dU/dE.
I already implemented the Jacobian analytically, however, I would like to
double check it using AD.
" (In theory an entire program could be made differentiable
Thanks for the help! I printed out the number of constraints after each of
your suggestions, and after setting the hanging node constraints the number
of constraints is 0 (as it should be, I'm working with a 1D mesh that has
been globally refined some number of times), but after interpolating th
On 8/31/22 06:59, Daniel Brauss wrote:
#1 typename SUNDIALS::KINSOL>::AdditionalData
additional_data;
#2 additional_data.strategy = newton;
You probably need to write this as
additional_data.strategy = KINSOL>::AdditionalData::newton;
The SolutionStrategy enum is declare
Dear All,
I am writing for help with SUNDIALS KINSOL in deal.II.
In deal.II’s step-77 tutorial, it looks like KINSOL’s linesearch algorithm
is used and linesearch appears to be the default SolutionStrategy. I am
trying to use a different solution strategy (like newton, picard, or
fixed_point).
