Fortran-arrays, Unterprogramme (sub-arrays)

Werde ich durch ein Fortran-code, und ein bisschen hat mich ein wenig verwirrt.

Es ist ein Unterprogramm, sagen

SUBROUTINE SSUB(X,...)
REAL*8 X(0:N1,1:N2,0:N3-1),...
...
RETURN 
END

Die aufgerufen wird, in einem anderen Unterprogramm von:

CALL SSUB(W(0,1,0,1),...)

wo W ist eine "Arbeits-array'. Es scheint, dass Sie einen bestimmten Wert von W wird an den X, aber X ist so dimensioniert wie ein array. Was ist Los?

InformationsquelleAutor ccook | 2010-03-22
  1. 8

    Dies ist nicht ungewöhnlich idiom für das Unterprogramm arbeiten auf einem (rechteckigen in N-Dimensionen) Teilmenge des ursprünglichen Arrays.

    Alle Parameter in Fortran (zumindest vor Fortran 90) werden per Referenz übergeben, so dass die eigentliche array-argument ist aufgelöst, von einem Speicherort im Arbeitsspeicher. Wählen Sie einen Speicherort innen der Speicherplatz für das ganze array und der subroutine manipuliert nur ein Teil des Arrays.

    Größte Problem: müssen Sie sich bewusst sein, wie das array im Speicher angeordnet und wie Fortran ' s array-Indizierung Schema funktioniert. Fortran verwendet Spalte große Auswahl bestellen, das ist das Gegenteil Konvention von c. Betrachten Sie ein array mit 5x5 in der Größe (und index beide Richtungen von 0 zu machen, der Vergleich mit c einfacher). In beiden Sprachen 0,0 ist das erste element im Speicher. In c das nächste element im Speicher ist [0][1] aber in Fortran ist es (1,0). Dies wirkt sich auf die Indizes fallen Sie bei der Wahl eines Subraum: wenn das original-array mit A(i,j,k,l), und das Unterprogramm arbeitet an drei-dimensionaler Unterraum (wie in deinem Beispiel), in c funktioniert es auf Aprime[i=constant][j][k][l], aber in Fortran arbeiten auf Aprime(i,j,k,l=constant).

    Die andere Gefahr ist die wrap-around. Die Abmessungen der (sub)array im Unterprogramm müssen denen entsprechen, die in der aufrufenden routine, oder seltsame, merkwürdige Dinge passieren (denken Sie). Also, wenn Einem erklärt wird, von der Größe (0:4,0:5,0:6,0:7), und wir rufen Sie mit element A(0,1,0,1) die empfangsroutine ist frei, start der index jeder dimension, wo auch immer es sein mag, aber müssen die Größen (4,5,6) oder anderes; aber das bedeutet, dass das Letzte element in der j-Richtung eigentlich umschlingt! Die Sache zu tun, über das nicht das Letzte element. Machen Sie sicher, dass das passiert, ist die Programmierer-job, und ist ein Schmerz in den Hintern. Kümmern. Viel Sorgfalt.

    • Danke, das Beispiel und den link besonders nützlich sind. Macht es Sinn für die übergabe sagen 0,0 oder 0,5, aber in diesem Fall, wo es 0,1,0,1, wie geht das zweite 1 ins Spiel kommen? Würde die nächsten Elemente werden 1,1,0,1 und dann nach dem ersten dimension 0,2,0,1?
    • Ihr Kommentar gekreuzt mit meinem edit. Der Letzte index in der aufrufenden routine ist effektiv Konstante in der aufgerufenen routine.
    • Ah, tut mir Leid, dass. Die Bearbeitungen, die sicherlich helfen, ty wieder. Ich nehme an, ich bin nur Probleme mit dem zweiten 1. Bringen, dass zu 2d -, ist, zu sagen, der start aus dem index von (0,1) die erste Spalte, die zweite Zeile, und der Reihe nach "Karte" in das array?
    • Okay, ich bin ein Narr. Die zweite dimension, beginnt, von ein und (0,1,0,1) ist (0,0,0,const). Danke!!!
    • Tatsächlich, auch wenn es nicht gilt, diese Zeit, diese Frage nicht kommen.
    • Meine computational methods professor wäre so stolz auf mich. Er war der Meister von diesem Unsinn. Zum Glück ist er nur uns beweisen, dass wir wickeln konnte aus wenig Köpfe um ihn herum einmal.
    • Nochmals vielen Dank 🙂 ich Liebe die Tangenten, die einen code senden uns auf.

    InformationsquelleAutor dmckee
  2. 3

    Dieser nennt sich "Sequenz " Verein". In diesem Fall, was zu sein scheint ein scaler, ein element eines Arrays (tatsächliche argument im Aufrufer) mit einem zugeordneten array (implizit beim ersten element), die dummy-argument in die Unterroutine . Danach werden die Elemente des arrays verbunden sind, die durch die Lagerung um, bekannt als "Sequenz". Dies wurde in Fortran 77 und früher aus verschiedenen Gründen, die hier offenbar für einen Arbeitsbereich array-vielleicht die Programmierer tun Ihr eigenes memory-management. Dies ist aufbewahrt in Fortran >=90 für die rückwärts-Kompatibilität, aber IMO gehört nicht in neuem code.

    • Was ist der moderne Ansatz?
    InformationsquelleAutor M. S. B.
  3. 2

    in fortran werden Variablen übergeben-Adresse.
    So W(0,1,0,1) ist Wert und Adresse. so im Grunde you passieren subarray ab W(0,1,0,1).

    • Vielen Dank für das feedback. Ich vermutete dabei, aber könnten Sie erklären, wie das subarray geprägt ist? Das ist, wo es beginnt, und wie es Karten für die neue Dimensionierung. Es tut Lesen/anzeigen der Reihe nach in das neue sizing?
    InformationsquelleAutor Anycorn
Schreibe einen Kommentar