~alexander-erlich/schroedingersolver/main

Viewing all changes in revision 5.

  • Committer: Alexander Erlich
  • Date: 2009-08-30 22:06:27 UTC
  • Revision ID: alexander.erlich@gmail.com-20090830220627-zw2ekzc6lm2mscj7
In diesem commit funktioniert prinzipiell sowohl die lineare, als auch die polynomiale Interpolation. Die Version ist sehr chaotisch, ein großer Teil ist auskommentiert. Als nächstes sollte der Code gesäubert und modularisiert werden. 

Die lineare Interpolation kann noch keine Fälle behandeln, wo mehrere Punkte denselben x-Wert haben bzw. x-Werte, die sich um weniger als die Diskretisierungs-Länge delta unterscheiden. Diesem Problem könnte man begegnen, indem man im bubble-Sort gleiche x-Werte detektiert und sie 'unterschiedlich macht', indem man eine (oder mehrere) Diskretisuerungs-Längen draufaddiert oder abzieht, so dass das Intervall [xMin, xMax] nicht überschritten wird. 

Die Polynomialinterpolation nutzt dgesv und funktioniert schon ganz gut. Es ist aber noch kein interface geschrieben worden, in kein module ausgelagert worden etc., dies 'Code-Kosmetik' steht noch an. Der Interpolationstyp kann auch noch nicht vom user gewählt werden (noch ist keine entsprechende Abfrage eingebaut), aber das sollte kein großes Problem sein. 

Beim Installieren der Lapack-Bibliothek unter Ubuntu sollte man aufpassen, dass man nicht nur liblapack aus den repositories installiert, sondern auch liblapack-dev. Ohne das -dev Paket funktioniert nämlich das Linken nicht.

Für das Einbinden von Lapack war neben der Übung zum Gauss-Eliminationsprogramm noch diese Seite sehr nützlich: http://de.wikibooks.org/wiki/Fortran:_LAPACK

expand all expand all

Show diffs side-by-side

added added

removed removed

Lines of Context: