Danke, hat jetzt geklappt mit dem Vektor x[] als Pointer.
Allerdings habe ich immer noch das Problem, dass mein Übergabepointer
df nicht stimmt.
Dazu habe ich nach längerem stöbern noch folgende nützliche Information gefunden:
Code:
int jaco(); *
* jaco hat die Form: *
* *
* int jaco (int n, REAL x[], REAL *mem[]) *
* { *
* REAL **
df; *
*
df = mem; *
* for (i = 0; i < n; i++) *
* for (j = 0; j < n; j++) *
*
df[i][j] = ...; *
* return (0); *
* } *
* *
* Dabei ist ... durch die partielle Ableitung der *
* i-ten Funktionskomponente nach der j-ten x-Komponen- *
* te zu ersetzten. mem ist hierbei die Lokation fuer *
* den Speicherbereich, der nach Ausfuehrung die *
* Jacobimatrix beinhaltet[/B].
und folgender Hinweis zu Übergabe von Matrizen:
Zitat:
Da die Übergabe von zweidimensionalen Feldern (Matrizen) an Funktionen in AnsiC für Anwendungen
ziemlich ungeeignet ist, werden in der Bibliothek nicht die Matrizen selbst an die
Funktionen übergeben, sondern mittels eines Feldes von Zeigern die Anfangsadressen der einzelnen
Zeilen der Matrix.
Man muss also etwa neben einer m × n–Matrix selbst zusätzlich ein Feld der Länge m von
Zeigern definieren, wobei der i–te Zeiger auf die i–te Zeile der Matrix zeigt.
Es wird also das Feld von Zeigern an die Funktion übergeben und es wird davon ausgegangen,
dass die Zeiger der Reihe nach auf die Zeilen der Matrix zeigen."
Ich habe mal folgende Varianten ausprobiert, bin mir nicht sicher welche Variante die richtige ist:
1. Übergabe der Matrix
df als Array von Pointern mit
df:= [(AdresseZeile1,AdresseZeile2)]
und
2. Übergabe als Pointer auf das Array of Pointer.
Leider haben beide Varianten nicht funktioniert. Die
DLL-Funktion gab mir bei Variante 1 als Rückgabewert 3 (Matrix sei singulär, was nicht sein kann) an.
Bei Variante 2 springt die
DLL einen Iterationsschritt weiter, doch als nächster Startwert (Vektor x[]) kommt so ein Blödsinn raus, dass ich in Gleitkomma überlauf habe (irgendwas ^300)
Hier noch meine Funktion:
Delphi-Quellcode:
PTArray = ^Tarray; //Pointer auf Tarray
Tarray = array[0..1] of double;
PTArray2 = ^Tarray2;
Tarray2 = array[0..1] of double;
PTarrayofPointer = ^TarrayofPointer; //Pointer auf array of Pointer
TarrayofPointer = Array[0..1] of Pdouble; //Array of Pointer;
// pMatrix = ^TMatrix;
TMatrix = array[0..1,0..1] of PDouble;
Delphi-Quellcode:
function Jac( x: PTarray;
df: TarrayofPointer) : integer;
cdecl;
var
JacMat: TMatrix;
i: integer;
Pinhalt:double;
begin
//Speicher für Jacobimatrix anfordern
new(JacMat[0][0]);
new(JacMat[0][1]);
new(JacMat[1][0]);
new(JacMat[1][1]);
//Jacobimatrix beschreiben
JacMat[0][0]^ := x[0] / 2 + 1;
JacMat[0][1]^ := 2 * x[1] / 5;
JacMat[1][0]^ := x[0] / 2;
JacMat[1][1]^ := x[1] + 1;
df[0] := @jacMat[0];
//Adresse 1. Zeile der Jacobimatrix in ein Array of Pointer
df[1] := @jacMat[1];
//Adresse 2. Zeile
result := 0;
end;
Wahrscheinlich mache ich bei der Adresszuweisung hier irgendwas falsch...
