On 8/25/21 9:37 PM, Anton Ermakov wrote:
>>> You can do that, but I believe that it's also possible to work with the
original system where one of the two matrices is only semidefinite.
I have redone it with the original system, but it seems that PETSc is
complaining about missing diagonal ent
Dear all
I am using TimeStepper with classical RK4 which works fine. I am facing some
issue with using embedded RK schemes. My usage is like this.
typedef LinearAlgebra::distributed::Vector PVector;
const TimeStepping::runge_kutta_method method = TimeStepping::DOPRI;
If I do
TimeStepping::Emb
I modified step-52 in following way
Change this line
TimeStepping::EmbeddedExplicitRungeKutta>
embedded_explicit_runge_kutta(method,
coarsen_param,
refine_param,
min_delta,
Dear all
In TimeStepping methods
We have functions taking these arguments
https://www.dealii.org/current/doxygen/deal.II/classTimeStepping_1_1ExplicitRungeKutta.html#ad43855dc35abf17e9c707444c0901790
https://www.dealii.org/current/doxygen/deal.II/classTimeStepping_1_1EmbeddedExplicitRungeKutta.
On 8/26/21 10:22 AM, Praveen C wrote:
ExplicitRungeKutta and EmbeddedExplicitRungeKutta are derived from RungeKutta
but there is no virtual function with these arguments in RungeKutta.
Is it possible to declare such a virtual function in the base class ? Seems
useful to have this.
Do all of
On 8/26/21 10:14 AM, Praveen C wrote:
Embedded explicit methods:
An error occurred in line <356> of file
in function
void dealii::Vector::sadd(const Number, const Number, const
Vector &) [number = double]
The violated condition