Samenvatting van bijeenkomst

Bespreking van het algorithme om pseudospectra uit te rekenen.

Aanwezigen:

• Petra Kok
• Tijmen Collignon,
• Joost Rommes,
• Ildiko Flesch

Wij bekeken, hoe het algorithme van de berekenen van pseudospectra in elkaar zit. Voor het berekenen kunnen we dus gebruiken algorithme in figuur 2 uit raport 'Large-scale computation of speudospectra using ARPACK and eigs'.

In dit algorithme zitten al een aantal "optimaliseringstrukjes" om een snellere methode te krijgen (zie algorithme). Er wordt Lanczos gebruikt om de systeem iteratief op te lossen.

Er zijn twee belangrijkere vragen:

• hoe ver moeten we nog optimalizeren? Er zijn nog meerdere ideeen uitgelegd in de copies, die we kunnen gebruiken. Deze optimalisaties zijn makkelijk te implementeren.
• hoe gaan we aanpakken? Welke taal zullen we gebruiken? Wij kunnen C++ of Matlab gebruiken, maar alle twee hebben voor- en nadelen. Wij kunnen C++ waarschinlijk (niet zeker dus) sneller krijgen en parallel maken, en het is wel overdraagbaar. Maar de programma moet dan wel helemaal vanaf het begin gemaakt worden. Aan de andere kant er is al een Matlab programma, die aangepasd en uitgebreid zou kunnen worden. Een heel groot voordeel is van Matlab, dat heel veel functies zitten al erin, en moeten we niet hun schrijven, of op Internet zoeken (bij voorbeeld eigenwaarden berekenen, verschillende eigenschappen van de matrix bekijken, etc.).

Tot dinsdag zal nog iedereen doen:

• meer informatie zoeken over uitbreidingen,
• een goed uitgaan algorithme zoeken, d.w.z. dat het duidelijk moet worden, waar wij beginnen.