LU-Zerlegung in C++

Sourcen:

• Vector.h, Vector.cc	dynamische Vektoren, Indizes starten ab 1
  Matrix.h, Matrix.cc	dynamische Matrizen, Fortran-Indizierung
  Random.h, Random.cc	Wrapper für Routine dran, versteckt seed
  Lutest.h, Lutest.cc	enthält die Lapack-Calls
  main.cc	einfaches Startprogramm
  veclib.h, lapack.h	extern-``C''-Wrapper um MLIB-Header

Bemerkungen:

• Matrizen und Vektoren werden als Handle (Referenzen) weitergereicht
  Kopien (incl. Daten) explizit mit copy
• Matrizen werden in Unterroutinen erzeugt und als Handle hochgereicht
  Vorteil: keine Kopien, bequem zu nutzen
  Nachteil: problematisch für Operator-Funktionen

Performance:

• erstaunlicherweise schneller als Fortran 90:

  
  

  

  
  

• und besser parallel: