1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
|
unit pseudohadamard;
{$mode objfpc}{$H+}
interface
uses
Classes, SysUtils, Gmp, Math;
procedure pseudoHadamardTransformation(var daten: array of mpz_t; modul: mpz_t; invers: boolean); overload;
procedure pseudoHadamardTransformation(var daten: array of mpz_t; Ordnung, Verschiebung: longint; modul: mpz_t; var buff: mpz_t; invers: boolean); overload;
implementation
procedure pseudoHadamardTransformation(var daten: array of mpz_t; modul: mpz_t; invers: boolean);
var
buff: mpz_t;
begin
if length(daten)<>round(power(2,round(ln(length(daten))/ln(2)))) then
raise exception.create('Die Länge der Daten für eine Pseudohadamardtransformation muss eine Zweierpotenz sein (und nicht '+intToStr(length(daten))+')!');
mpz_init(buff);
pseudoHadamardTransformation(daten,length(daten),0,modul,buff,invers);
mpz_clear(buff);
end;
procedure pseudoHadamardTransformation(var daten: array of mpz_t; Ordnung, Verschiebung: longint; modul: mpz_t; var buff: mpz_t; invers: boolean);
var
i: longint;
begin
Ordnung:=Ordnung div 2;
if Ordnung>1 then
for i:=0 to 1 do
pseudoHadamardTransformation(daten,Ordnung,Verschiebung + i*Ordnung,modul,buff,invers);
if invers then begin
for i:=0 to Ordnung-1 do begin
mpz_set(buff,daten[Verschiebung+i]);
mpz_sub(daten[Verschiebung+i],daten[Verschiebung+i],daten[Verschiebung+i+Ordnung]); // a = a' - b
mpz_mod(daten[Verschiebung+i],daten[Verschiebung+i],modul);
mpz_mul_ui(daten[Verschiebung+i+Ordnung],daten[Verschiebung+i+Ordnung],2); // b = 2*b'
mpz_sub(daten[Verschiebung+i+Ordnung],daten[Verschiebung+i+Ordnung],buff); // - a'
mpz_mod(daten[Verschiebung+i+Ordnung],daten[Verschiebung+i+Ordnung],modul);
end;
end
else
for i:=0 to Ordnung-1 do begin
mpz_set(buff,daten[Verschiebung+i]);
mpz_mul_ui(daten[Verschiebung+i],daten[Verschiebung+i],2); // a' = 2 * a
mpz_add(daten[Verschiebung+i],daten[Verschiebung+i],daten[Verschiebung+i+Ordnung]); // + b
mpz_mod(daten[Verschiebung+i],daten[Verschiebung+i],modul);
mpz_add(daten[Verschiebung+i+Ordnung],daten[Verschiebung+i+Ordnung],buff); // b' = a + b
mpz_mod(daten[Verschiebung+i+Ordnung],daten[Verschiebung+i+Ordnung],modul);
end;
end;
end.
|
