003
07.07.2008, 16:35 Uhr
Arno
|
Aus Langeweile mache ich es nicht, sondern weil ich es lernen soll Mit einer lib(rary) möchte ich nur ungern arbeiten - CBLAS sagt mir z.B. auch nichts.
Was den Rest deiner Antwort angeht bin ich mir nicht sicher ob ich richtig verstanden wurde, deshalb nochmal ein paar zusätzliche Worte:
F ist wie gesagt eine n x m Matrix, wobei die Spaltenanzahl um eins größer ist als die Zeilenzahl, also m = n+1 u ist eine einspaltige Matrix mit m Zeilen.
Wie ich die Matrix F mit der Matrix u multipliziere ist mir schon klar, es geht mir um den Aufbau/die Programmierung der Matrix F. Dafür benötige ich die Werte b1 bis bn die ich vorher einlese (siehe Programmcode oben). In der Matrix F tauchen die Werte dann entweder als 1/b, -1/b oder gar nicht auf, daher komme ich auch auf die Nullen bzw. die Divisionen.
Ich schreibe den Aufbau der Matrix F nochmal in Worten:
1. Zeile, 1. Spalte: 1/b1 1. Zeile, 2. Spalte: -1/b1 Alle anderen Spalten der 1. Zeile sind 0
2. Zeile, 1. Spalte: 0 2. Zeile, 2. Spalte: 1/b2 2. Zeile, 3. Spalte: -1/b2 Alle weiteren Spalten der 2. Zeile sind 0
3. Zeile, 1. Spalte: 0 3. Zeile, 2. Spalte: 0 3. Zeile, 3. Spalte: 1/b3 3. Zeile, 4. Spalte: -1/b3 Alle weiteren Spalten der 3. Zeile sind 0
So setzt sich das Schema bis zur letzten (n-ten) Zeile fort. Die Belegung der n-ten Zeile lautet dann: n-te Zeile, 1. Spalte bis n-1. Spalte: 0 n-te Zeile, n. Spalte: 1/bn n-te Zeile, n+1. Spalte: -1/bn
Ich hoffe es ist jetzt klarer um was es mir geht |