Zum Problem:
Ich habe eine Matrix (das ist wie eine 2D Tabelle) die entweder nur mit double-werten gefüllt ist ODER pro element mit 2 Double werten (wenn sie komplexwertig ist).
Bsp:
1.0 0.0 4.0
4.4 0.0 0.0
0.0 0.0 1.5

oder
1.0 +9.5i 0.0 +0.0i ....

Desweiteren wird aus der Ausgangsmatrix immer nur ein Teil betrachtet der sich auch ändert. D.h einmal sind es die spalten 1-4 dann mal wieder 7-9 und so weiter. Ich speichere also nur einen Teil der Matrix.
Ich will da über Templates gehen wegen den daten.

Ich dachte ich mache folgendes:

class Matrix
{    
      //hier einige public/private sachen die zur ganzen Matrix gehören

       //das sind die gestauchten Spalten der ausgansmatrix
       Lines lines;
};

class lines
{
   //hier einige public/privae sachen UND
   // ein 2Darray aus lauter werten
   rcv **arr; //2D array referenziert lauter arr_buf's
   rcv *arr_buf; //array das die templates trägt ,

};

//template struct
template< type T > struct rcv < T >
{
   int r; //Zeile
   int c; //Spalte
   T val; //Wert der ja auch complex sein kann.
};

Für weitere Vorschläge wäre ich sehr dankbar

also ganz ehrlich gesagt weiß ich nicht genau was du willst.
eine matrix koenntest du auch mit zwei vektoren anlegen

template<class T>
class Matrix
{
   public:
   typedef std::vector<T> line;

   line& operator[](std::size_t pos) { return matrix[pos]; }

   private:
      std::vector<line> matrix;
};

Matrix<double> matrix1;

ist das sowas was du suchst ?

Meep Meep

ich versuchs mal an nem beispiel:

matrix :
0 0 0
3 0 0
0 1 0
6 0 0

2D array A =
-> 3 6 //1D array
-> 1 //1D array

//trägt nur die Zeilenindizes wo die werte stehen, das erste array gehört zur
// ersten spalte, das zweite zur zweiten
2D array ptrs =
-> 2 4 //1D array
-> 3 //1D array

bisschen freaky....

ich glaub das ich dich nun verstanden habe.
fuer die wahl der richtigen datenstruktur bzw. kontainer, muesstest du jedoch auch mal offen legen welchen kriterien die matrix unterlegen ist.
brauchst du da auch direkten zugriff auf die daten ? wird haeufig geloescht oder hinzugefuegt. wie gross wird die matrix ueberhaupt sein ? kann man das von vornherein schon abschaetzen ?
was willst du ueberhaupt machen ? welchen grund hat dein vorgehen ?

Meep Meep

brauchst du da auch direkten zugriff auf die daten?

Falls Du meinst ob ich direkt auf die Werte zugreifen will und nicht nur über die Indizes: Dann ja

wird haeufig geloescht oder hinzugefuegt.

eher selten

wie gross wird die matrix ueberhaupt sein ?

Die matrix ist dünnbesetzt , d.h. die meisten einträge sind null- deswegen werden nur die nicht -Null einträge gespeichert. das führt dann zu array-grössen von z.B 100 - 50000. (Bei 10000x10000 wäre es ja insgesamt 10^8 !) Also kann man sagen ja: die arrays können groß werden.

kann man das von vornherein schon abschaetzen ?

Ja - der input ist eben eine solche matrix. Ich muss zuerst mal übers file laufen um mir die infos zu holen. Die dimension steht schon vorher fest. Das
file beinhaltet also eine input matrix.

was willst du ueberhaupt machen ? welchen grund hat dein vorgehen ?

Puh - da müsste ich jetzt gewaltig ausholen: Prinzipiell hole ich mir eine Matrix und speichere mir nur die nicht-null einträge - das steht auch in dem file so drin - da sind nur nicht-null einträge vorhanden. Diese zeilen/spalten/value werte puffere ich mir also um damit arbeiten zu können.
Aus meinem Werten (das sind wohlgemerkt nur die nicht-null einträge) extrahiere ich in einer for-schleife immer nur einen Teil der Matrix. Also z.B in iteration 1: spalten 1-4, dann in iteration 2: spalten 5-9 etc...
Mit diesen Submatrizen hantiere ich jetzt rum - berechne ne QR-Zerlegung und pflantsche das Ergebnis in eine neue Matrix. Also mathematisches Zeug mache ich.
Wichtig ist: es geht um "optimalen" code. D.h die Laufzeit ist erstes Kriterium.
Und noch etwas: Der input kann sowohl komplex sein als auch nur real. Das was ich am anfang gepostet habe in meinem ersten post.

Ich weiß nicht ob ihr/Du mir da helfen könnt - aber über weitere Hilfe bin ich sehr dankbar.

so mal gucken ob ich das nun richtig verstanden habe.
du hast eine datei, in welcher eine matrix bestimmter groesse (sprich bekannter groesse) gespeichert ist. jedoch nur die werte welche != 0 sind.
dafuer haben sie jeweils zeile und spalte als angabe dabei.
nun willst du die datei in den speicher laden um damit zu arbeiten. die komplette matrix wuerde jedoch zuviel speicher in anspruch nehmen, weshalb du auch im speicher nur die werte != 0 halten willst.
damit du jedoch zuegig damit arbeiten kannst, brauchst du eine struktur, welche nur die werte != 0 enthaelt, aber nach aussen hin wie eine komplette matrix aussieht, damit du mit normalen indizies arbeiten kannst. ist das soweit korrekt ?
sind die daten der matrix schon in einer bestimmten art vorsortiert ?

Meep Meep

Puh - da müsste ich jetzt gewaltig ausholen: Prinzipiell hole ich mir eine Matrix und speichere mir nur die nicht-null einträge - das steht auch in dem file so drin - da sind nur nicht-null einträge vorhanden. Diese zeilen/spalten/value werte puffere ich mir also um damit arbeiten zu können.
Aus meinem Werten (das sind wohlgemerkt nur die nicht-null einträge) extrahiere ich in einer for-schleife immer nur einen Teil der Matrix. Also z.B in iteration 1: spalten 1-4, dann in iteration 2: spalten 5-9 etc...
Mit diesen Submatrizen hantiere ich jetzt rum - berechne ne QR-Zerlegung und pflantsche das Ergebnis in eine neue Matrix. Also mathematisches Zeug mache ich.
Wichtig ist: es geht um "optimalen" code. D.h die Laufzeit ist erstes Kriterium.
Und noch etwas: Der input kann sowohl komplex sein als auch nur real. Das was ich am anfang gepostet habe in meinem ersten post.

Ich weiß nicht ob ihr/Du mir da helfen könnt - aber über weitere Hilfe bin ich sehr dankbar.

Wenn es nur darum geht, ein dünn besetztes LGS schnell zu Lösen, verwende einfach

UMFPACK (http://www.cise.ufl.edu/research/sparse/umfpack/)

oder Sparselib.

Wenn es jedoch um die eigentliche implementierung geht, google mach nach "Compressed Row / Compressed Column Matrix"

Gruss, Gast

nun willst du die datei in den speicher laden um damit zu arbeiten. die komplette matrix wuerde jedoch zuviel speicher in anspruch nehmen, weshalb du auch im speicher nur die werte != 0 halten willst.
damit du jedoch zuegig damit arbeiten kannst, brauchst du eine struktur, welche nur die werte != 0 enthaelt, aber nach aussen hin wie eine komplette matrix aussieht, damit du mit normalen indizies arbeiten kannst. ist das soweit korrekt ?
sind die daten der matrix schon in einer bestimmten art vorsortiert ?

jein: in der datei sind auch nur die nicht-nuller drin - aber die komplette dimension der matrix schon. Daraus weiß ich schon mal wie groß sie ist. Alle anderen Indizes sind ja Null die nicht in der datei über row column value angegeben sind.
Nach außen hin muss es gar nicht so aussehen. Also intern muss ich halt aufpassen: Es gibt da z.B ne stelle wo ich eine kleine Submatrix aufbauen muss z.B aus den spalten 2-4 und da muss ich dann halt Nuller in das array an den richtigen positionen einfügen und die nicht-nuller halt anhand meiner gespeicherten Indizes. Also wenn du das mit nach außen hin meinst - dann ja->nach außen muss es als komplette Matrix gelten.
Die Daten sind leider in keiner bestimmten Weise vorsortiert. Die können willkürlich im file rumliegen: Prinzipiell ist das file so aufebaut:

Im folgenden ist die erste Zahl der zeilenindex, die zweite zahl der spaltenindex und die dritte zahl der wert . Also r, c, v

dim_x dim_y
2 3 5.99
2 4 9,7
4 7 8
... ... ...

Ich nutze natürlich die Lapack-routinen die schneller gar nicht implementiert existieren - aber ich muss gewisse dinge mit meinen daten machen für die es noch nichts gibt Deswegen das ganze tratra

wuerde es reichen wenn du alle deiner werte in ein array rein tust.
dann sortieren und dann mittels std::lower_bound oder so danach suchen wuerdest ?
oder wie schnell waere fuer dich schnell genug ?

Meep Meep

Hmm...naja es könnte evtl. alles in ein array rein. das sollte dann aber ein 2D array sein - pro spalte ein array im array.
ich überlege gerade weiter ob ich irgendwo in teufels küche käme so....

Gibt es einen bestimmten grund dafuer ?

http://www.intel.com/software/products/mkl/docs/webhelp/appendices/mkl_appA_SMSF.html

die section CSC Format!!!

so stehen da ja 3 arrays - jetzt weiß ich nicht ob es sinn machen würde nicht 3 einzelne arrays zu nehmen sondern z.B ein 2D array. Performance-technisch bin ich da unsicher.

Meine Idee wäre jetzt mal ein 2D array zu bauen, wo jedes element ein array ist das insgesamt eine Spalte der Matrix trägt. D.h ein struct mit r,c,v bauen und das dann ins array. Dieses array an position 1 des "Hauptarrays" wenn es die erste spalte der matrix ist.
Zusätzlich könnte man noch ein 2tes array mitschleifen welches nur die Länge der einzelnen arrays wiedergibt.

Ich muss halt noch das mit dem complex-wertigen eingaben mittels templates verbauen...

Wahrscheinlich hängt das alles von meiner Anwendung ab und ihr könnt mir da wohl nicht wirklich weiter helfen - nichts desto trotz-danke euch!