Ich habe noch den Code angeschaut, welcher mit der Originalarbeit veröffentlicht worden ist.
Dort steht "delphisch" geschrieben:
Delphi-Quellcode:
y2 := int_state.coef[i + 1][0] * x_4 +
int_state.coef[i + 1][1] * x_3 +
int_state.coef[i + 1][2] * x_2 +
int_state.coef[i + 1][3] * x_1 + 1;
In deinem
DEC Code steht nach deiner Korrektur von gestern genau das, was du im Originalcode auch findest. Wenn ich nicht was voll falsch interpretiere/lese, dann war der SCOP Code in
DEC bisher nicht OK.
Wenn du die Originalarbeit liest, siehst du, dass bei GP8 acht Polynome vierten Grades vom Typ Y = a*X4 + b*X'3 + c*X^2 + d*X + 1 an gewissen Stellen ausgewertet werden. Der bisherige Code machte keinen Sinn. Im korrigierten Code (wie oben), wird Polynom p[i+1] mit den Koeffizienten A=int_state.coef[i+1][0], B=int_state.coef[i+1][1], C=int_state.coef[i+1][2], D=int_state.coef[i+1][3] an der Stelle x_1 ausgewertet. Und das ist gut so.
Tests könntest du mit dem C Code aus der Originalarbeit erstellen (?).