C-dll Speicheradressen

[Fritzew]

nur auf den ersten Blick:

Das ist falsch, ändere

Delphi-Quellcode:

Tarray = array[0..1] of double;
TMatrix = array[0..1,0..1] of Pdouble;

zu

Delphi-Quellcode:

Tarray = array[0..1] of double;
TMatrix = array[0..1,0..1] of double; // ist auch Double

Hintergrund die funktion erwartet einen Pointer auf ein Array of Double nicht einen Pointer auf ein Array of Pointer

[Johannes G.]

Hallo Fritz,

danke für deine schnelle Antwort!

Leider funktioniert das ganze noch nicht. Wenn ich aus dem Pdouble ein double mache, also die Matrix zurück gebe, kommt die Fehlermeldung 'Access Violation at ....: read of the address ...'

Eine andere Funktion der gleichen DLL ist auch so aufgebaut, nur mit dem Unterschied keine Matrix(*df[]) sondern eine normale Variable (*y). Ich habe es bei dieser Funktion auch so gemacht, dass ich auf diese Adresse schreibe mit y^:=, siehe unten. Hat alles wunderbar funktioniert!

Code:

function f(k:integer; x: Tarray; y: Pdouble): integer; cdecl;
 begin
    if k=0 then
    begin
      y^ := x[0]*x[0]/4 + x[1]*x[1]/5 + x[0] - 0.5;
    end else
    begin
      y^ := x[0]*x[0]/4 + x[1]*x[1]/2 + x[1] - 0.5;
    end;
    result := 0;
 end;

C-Funktionsaufruf:

Code:

int f( int k, double x[], double *y)

Mit der Beschreibung der Funktion:

Code:

Funktion, welche das n × n–System beschreibt.
In
int k;
double x[], y;
...
f( k, x, &y);
...
wird der Wert der k-ten Komponentenfunktion des Systems an der Stelle x
(Vektor der Länge n) berechnet und das Ergebnis (also eine Zahl) nach y geschrieben.
(Im Gegensatz zum Verfahren newt berechnet hier die Funktion f
nur eine Komponente des nicht linearen Systems!)

Ich stehe im Moment ziemlich auf der Leitung wie ich das mit dem Pointer auf das Array (bzw. Matrix) realisieren soll.
Vielleicht hat noch jemand eine Idee?

Gruß Johannes

[Fritzew]

Code:

"int jaco( double x[], double *df[])

das ist auf Delphi Seite eigentlich

Delphi-Quellcode:

type
pDoubleArray = ^TdoubleArray;
tDoubleArray = array[0..0] of Double;

function Jaco(x : tDoubleArray;  df : pDoubleArray) : integer; cdecl;

[Johannes G.]

Danke, jetzt habe ich das mit dem Pointer auf ein Array verstanden

Ich habe mein Bsp. einwenig abgeändert, da der Pointer df auf die Matrix df[][] und nicht auf das Array x[] zeigen soll!

Code:

pMatrix = ^TMatrix;
  tMatrix = array[0..1,0..1] of Double;

Funktion jac:

Code:

function Jac(var x: Tarray; var df: pMatrix) : integer; cdecl;
begin

  df[0][0] := x[0] / 2 + 1;
  df[0][1] := 2 * x[1] / 5;
  df[1][0] := x[0] / 2;
  df[1][1] := x[1] + 1;

  result :=  0;
end;

Der ptr df zeigt jetzt auf den ersten Eintrag der Matrix df[0][0], die Adressen der anderen 3 Einträge sind jeweils um 8Byte dahinter (aneinander gefolgt). Ich denke das passt soweit!

Irgendwas stimmt aber trotzdem noch nicht. Wenn ich durch die Funktion debugge, öffnet sich am Ende der Reiter "CPU" mit folgender Zeile:

7775067D C605A592787700 mov byte ptr [$777892a5],$00 was hat das zu bedeuten?

Auch die Funktion newt der DLL gibt mir einen Fehler zurück: Matrix singulär. Heißt also, dass die Einträge der df-Matrix doch nicht am richtigen Platz landen

[Fritzew]

Zitat von Johannes G.:

Danke, jetzt habe ich das mit dem Pointer auf ein Array verstanden

Ich habe mein Bsp. einwenig abgeändert, da der Pointer df auf die Matrix df[][] und nicht auf das Array x[] zeigen soll!

Code:

pMatrix = ^TMatrix;
  tMatrix = array[0..1,0..1] of Double;

Funktion jac:

Code:

function Jac(var x: Tarray; var df: pMatrix) : integer; cdecl;
begin

  df[0][0] := x[0] / 2 + 1;
  df[0][1] := 2 * x[1] / 5;
  df[1][0] := x[0] / 2;
  df[1][1] := x[1] + 1;

  result :=  0;
end;

Der ptr df zeigt jetzt auf den ersten Eintrag der Matrix df[0][0], die Adressen der anderen 3 Einträge sind jeweils um 8Byte dahinter (aneinander gefolgt). Ich denke das passt soweit!

Irgendwas stimmt aber trotzdem noch nicht. Wenn ich durch die Funktion debugge, öffnet sich am Ende der Reiter "CPU" mit folgender Zeile:

7775067D C605A592787700 mov byte ptr [$777892a5],$00 was hat das zu bedeuten?

Auch die Funktion newt der DLL gibt mir einen Fehler zurück: Matrix singulär. Heißt also, dass die Einträge der df-Matrix doch nicht am richtigen Platz landen

Deine Jac stimmt nicht der erste Parameter ist kein var..........

[Johannes G.]

Das Stimmt, habe ich ausgebessert. Allerdings leider das gleiche Ergebnis...

[Zacherl]

Zitat von Johannes G.:

Das Stimmt, habe ich ausgebessert. Allerdings leider das gleiche Ergebnis...

Es sollte meiner Meinung nach entweder nur pMatrix oder var TMatrix sein. Edit: Ah nvm, habe die falsche Signatur angeschaut. Diese ist richtig? int jaco( double x[], double *df[])

Arrays in C sind automatisch immer Zeiger. Also statt dobule bla[] könnte man in C auch einfach double* bla schreiben. Durch die spezielle Zeigerarithmetik in C darf der Array Operator auf jeden Zeiger angewendet werden und greift dabei auf den Speicher bei Index * sizeof(ElementType) zu. Ähnlich funktioniert ja auch Delphis Inc , wenn du es auf eine "typisierte Zeigervariable" z.b. PDobule anwendest. Der Zeiger würde dann automatisch um SizeOf(Double) erhöht werden.

Wenn du deinen eigenen Matrix Typ erstmal weglässt, funktioniert es dann folgendermaßen:

Delphi-Quellcode:

type
  PDoubleArray = ^TDoubleArray;
  TDoubleArray = array[0..0] of Double;

function Jaco(X: PDoubleArray; var DF: PDoubleArray): Integer; cdecl;