On 3/13/19 9:18 AM, Apurva Tiwari wrote:
>
> So far I'd been using an RK3 integrator that I had written myself. To
> rectify
> the problem I faced (mentioned in the mail earlier), I thought of, was using
> the built-in explicit RK4 integrator TimeStepping::ExplicitRungeKutta
> <https://dealii.org/developer/doxygen/deal.II/classTimeStepping_1_1ExplicitRungeKutta.html>.
> Now, I have a function called "assemble_rhs" that uses MeshWorker to evaluate
> the rhs of the explicit time dependent equation. Its signature is:
> Step12<dim>::assemble_rhs (RHSIntegrator<dim>& rhs_integrator)
> (I'm modifying part of the step-12 tutorial program). But the evolve_one_step
> function is of format:
> TimeStepping::ExplicitRungeKutta
> <https://dealii.org/developer/doxygen/deal.II/classTimeStepping_1_1ExplicitRungeKutta.html><
>
> VectorType >::evolve_one_time_step ( const std::function<
> VectorType(const
> double, const VectorType &)> & f,
> double t,
> double delta_t,
> VectorType & y
> )
>
>
> To this, I tried to pass assemble_rhs as the first arg and ofcourse it failed
> to compile. I've not been able to figure out how do I use my rhs function to
> get to evaluate the next time-step. I'm probably making quite a fundamental
> mistake, can you please guide me?
Yes, the function you have is not the function that works in the context of
the time integrator.
You will want to take a look at the tutorial program that shows how to use the
time integrator to see what kind of function can be passed as first argument.
Best
W.
--
------------------------------------------------------------------------
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because you are subscribed to the Google Groups
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to dealii+unsubscr...@googlegroups.com.
For more options, visit https://groups.google.com/d/optout.
