ich habe mal wieder eine Frage in Bezug auf die Übersetzung von Fortran Code in C++.

In Fortran sind Konstrukte wie

integer grdpoi(ngrdp,nbl,ngrid)
...
call foo(grdpoi(1,1,n))
...

üblich. In foo kommt dann nur noch ein 2D Array an.

subroutine foo(grdpoi)
integer grdpoi(ngrdp,nbl)
...

Wenn man wegen der Übersetzung die Paramterreihenfolge wie in Fortran belassen will. Wie kann man aus einer n-Dimensionalen Matrix Klasse eine n-1 Dimensionale Matrix übergeben?

Gruß Andreas Georg

aber ich verstehe valarray nicht so richtig. So wie ich es verstehe kann ich mit slices auf meine Elemente zugreifen. Aber wenn ich das valarray an foo übergebe, übergebe ich doch den ganzen Vektor. Wie reduziert sich die Dimension von z.B. 3D in der aufrufenden Funktion zu 2D in foo?

In Fortran verstehe ich was passiert, dort läuft ja der linke Index am schnellsten, so das grdpoi(1,1,n) mir die Adresse für das n-te 2D Array übergibt. Und in foo muß ich die Dimension für das 2D Array halt wissen (woher auch immer).

Ich behelf mir bisher damit

void foo(double* GRDPOI)
{
  auto grdpoi = [] (int i, int j) -> double { return GRDPOI[(i-1)+(j-1)*ngrdp];};

...

}

So ein Vorspann am Anfang der Funktion funktioniert zwar, ist aber genauso unansehnlich wie der FORTRAN Code und sollte vieleicht eleganter zu realisieren sein.

Gruß Andreas Georg