[cfe-dev] [test-suite] making polybench/symm succeed with "-Ofast" and "-ffp-contract=on"

[Hal Finkel via cfe-dev]

----- Original Message -----
> From: "Sebastian Pop" <sebpop.llvm at gmail.com>
> To: "Hal Finkel" <hfinkel at anl.gov>
> Cc: "Sebastian Paul Pop" <s.pop at samsung.com>, "llvm-dev" <llvm-dev at lists.llvm.org>, "Matthias Braun"
> <matze at braunis.de>, "Clang Dev" <cfe-dev at lists.llvm.org>, "nd" <nd at arm.com>, "Abe Skolnik" <a.skolnik at samsung.com>,
> "Renato Golin" <renato.golin at linaro.org>
> Sent: Monday, October 10, 2016 9:10:01 AM
> Subject: [test-suite] making polybench/symm succeed with "-Ofast" and "-ffp-contract=on"
> 
> Hi,
> 
> I would need some help to fix polybench/symm:
> 
> void kernel_symm(int ni, int nj,
> DATA_TYPE alpha,
> DATA_TYPE beta,
> DATA_TYPE POLYBENCH_2D(C,NI,NJ,ni,nj),
> DATA_TYPE POLYBENCH_2D(A,NJ,NJ,nj,nj),
> DATA_TYPE POLYBENCH_2D(B,NI,NJ,ni,nj))
> {
>   int i, j, k;
>   DATA_TYPE acc;
> 
>   /*  C := alpha*A*B + beta*C, A is symetric */
>   for (i = 0; i < _PB_NI; i++)
>     for (j = 0; j < _PB_NJ; j++)
>       {
>         acc = 0;
>         for (k = 0; k < j - 1; k++)
>           {
>              C[k][j] += alpha * A[k][i] * B[i][j];
>              acc += B[k][j] * A[k][i];
>           }
>         C[i][j] = beta * C[i][j] + alpha * A[i][i] * B[i][j] + alpha
>         * acc;
>       }
> }
> 
> Compiling this kernel with __attribute__((optnone)) and outputing the
> contents of the C[][] array does not match the reference output.

Why is this? What compiler are you using? Are we not using IEEE FP @ -O0 (e.g. using x87 floating point)? IEEE FP, without FMA, should be completely deterministic. Sounds like a bug.

> Furthermore, compiling this kernel at -Ofast and comparing against
> -O0
> only passes for FP_ABSTOLERANCE=10.
> All the 10 other polybench tests that I have transformed to check FP
> are passing at FP_ABSTOLERANCE=1e-5 (and most likely they could pass
> at an even more reduced tolerance.)
> 
> The symm benchmark seems to accumulate all the errors as it is a big
> reduction from the first elements of the C[][] array into the last
> elements.
> I'm not sure we can rely on this benchmark to check FP correctness.
> 
> One option is to completely specify which optimization flags have
> been
> used to compute the reference output and only use that to compile
> this
> benchmark.
> 
> Please share your ideas on how to deal with this particular test.

If the test is not numerically stable, we can:

 1. Only test the non-FP-contracted output
 2. Run the FP-contracted test only for a very small size (so that we'll stay within some reasonable tolerance of the reference output)
 3. Change the matrix to something that will make the test numerically stable (it does not look like the matrix itself matters to the performance; where do the values come from?).

 -Hal

> 
> Thanks,
> Sebastian
> 

-- 
Hal Finkel
Lead, Compiler Technology and Programming Languages
Leadership Computing Facility
Argonne National Laboratory