[Top][All Lists]

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: DAE daspk solver problem

From: SunilShah
Subject: Re: DAE daspk solver problem
Date: Fri, 28 Oct 2016 10:17:52 -0700 (PDT)

As I understand it, Radau5.f demands a structure of DAEs with Linearly
implicit form:

M Y' = F(Y,T)

DASPK admits allows a more general form F(Y,Y',T)=0.  There can be problems
that can not be transformed to the linearly implicit form.  

For large scale problems, (> 1e5 variables), retaining sparsity is an
important issue, so DASPK structure can offer some advantages.  Even when
the equations can be symbolically transformed, symbolic manipulation may not
scale well (perhaps the problem with scaling in Modelica),  or may destroy
sparsity and therefore have inefficient discrete propagation.

Lastly, the fortran DASPK / DASKR can take advantage of implicit sparse
solvers.  Do octave implementations allow this?

View this message in context:
Sent from the Octave - General mailing list archive at

reply via email to

[Prev in Thread] Current Thread [Next in Thread]