fastfloat_maketestvectors_random.m

来自「ADI BF 16位定点DSP的快速浮点仿真的汇编代码」· M 代码 · 共 57 行

% ff16_makeTestVectors_random.m

VEC_SIZE = 64;

a = (rand(VEC_SIZE,1)-0.5).*2.^(128*(rand(VEC_SIZE,1)-0.5));

% a = [0, -0, 2, 6.5, -6.5, 2^-126, 2^-127, 2^-149; \
0, -0, 2, 6.5, -6.5, 2^-126, 2^-127, 2^-149; \
0, -0, 2, 6.5, -6.5, 2^-126, 2^-127, 2^-149; \
0, -0, 2, 6.5, -6.5, 2^-126, 2^-127, 2^-149; \
0, -0, 2, 6.5, -6.5, 2^-126, 2^-127, 2^-149; \
0, -0, 2, 6.5, -6.5, 2^-126, 2^-127, 2^-149; \
0, -0, 2, 6.5, -6.5, 2^-126, 2^-127, 2^-149; \
0, -0, 2, 6.5, -6.5, 2^-126, 2^-127, 2^-149]

b = (rand(VEC_SIZE,1)-0.5).*2.^(128*(rand(VEC_SIZE,1)-0.5));
% b = a'

% a = reshape(a, 64, 1);
% b = reshape(b, 64, 1);

fprintf("min(a) = %10g, min(abs(a)) = %10g, max(a) = %10g\n", min(a), min(abs(a)), max(a));
fid = fopen("vec_a.dat", "w");
fprintf(fid, "%30g,\n", a(1:VEC_SIZE-1));
fprintf(fid, "%30g\n", a(VEC_SIZE));
fclose(fid);

fprintf("min(b) = %10g, min(abs(b)) = %10g, max(b) = %10g\n", min(b), min(abs(b)), max(b));
fid = fopen("vec_b.dat", "w");
fprintf(fid, "%30g,\n", b(1:VEC_SIZE-1));
fprintf(fid, "%30g\n", b(VEC_SIZE));
fclose(fid);

% add

ans_add = a+b;
fid = fopen("add.ans", "w");
fprintf(fid, "%30g,\n", ans_add(1:VEC_SIZE-1));
fprintf(fid, "%30g\n", ans_add(VEC_SIZE));
fclose(fid);

% subtract

ans_sub = a-b;
fid = fopen("sub.ans", "w");
fprintf(fid, "%30g,\n", ans_sub(1:VEC_SIZE-1));
fprintf(fid, "%30g\n", ans_sub(VEC_SIZE));
fclose(fid);

% multiply

ans_mult = a.*b;
fid = fopen("mult.ans", "w");
fprintf(fid, "%30g,\n", ans_mult(1:VEC_SIZE-1));
fprintf(fid, "%30g\n", ans_mult(VEC_SIZE));
fclose(fid);