A note on the consistent initialization for nonlinear index-2 differential-algebraic equations in Hessenberg form
2004 (English)Report (Other academic)
In a previous paper, we developed a new algorithm for the consistent initialization of general index-2 differential-algebraic equations arising within the method of lines for solving partial differential equations. In the case of Hessenberg systems, the structural information allows for a lot of simplifications thus allowing for much larger systems to be solved. The crucial point consists of providing sparse projections by the use of sparse approximations to the inverse of the mass matrix. We obtain almost linear computational complexity with respect to the number of degrees of freedom.
Place, publisher, year, edition, pages
Stockholm: KTH and Parallel and Scientific Computing Institute , 2004. , 10 p.
, PSCI Report, 04:02
differential-algebraic equations, consistent initial values, consistent initialization, method of lines, MATLAB
IdentifiersURN: urn:nbn:se:kth:diva-48499OAI: oai:DiVA.org:kth-48499DiVA: diva2:457763
QC 201111222011-11-222011-11-192011-11-22Bibliographically approved