
From:  Bill Greene 
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
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: http://octave.1599824.n4.nabble.com/DAEdaspksolver problemtp4680377p4680409.html
Sent from the Octave  General mailing list archive at Nabble.com.
_______________________________________________
Helpoctave mailing list
address@hidden
https://lists.gnu.org/mailman/listinfo/helpoctave
[Prev in Thread]  Current Thread  [Next in Thread] 