|Subject:||Re: DAE daspk solver problem|
|Date:||Fri, 28 Oct 2016 14:00:01 -0400|
As I understand it, Radau5.f demands a structure of DAEs with Linearly
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: http://octave.1599824.n4.
Sent from the Octave - General mailing list archive at Nabble.com.
Help-octave mailing list
|[Prev in Thread]||Current Thread||[Next in Thread]|