New user self-registration is disabled due to spam. Please email eigen-core-team @ if you need an account.
Before reporting a bug, please make sure that your Eigen version is up-to-date!
Bug 1555 - Eigen 3.3.4 compilation fails on xlc
Summary: Eigen 3.3.4 compilation fails on xlc
Alias: None
Product: Eigen
Classification: Unclassified
Component: Core - vectorization (show other bugs)
Version: 3.3 (current stable)
Hardware: All Other
: Normal Compilation Problem
Assignee: Nobody
Depends on: 1494
Blocks: 3.4
  Show dependency treegraph
Reported: 2018-06-14 07:09 UTC by Kate Clark
Modified: 2018-10-08 20:59 UTC (History)
10 users (show)


Description Kate Clark 2018-06-14 07:09:55 UTC
Compiling Eigen 3.3.4 (the latest release) with xlc++ 16.0.3 gives this error message:

In file included from .../eigen/src/Eigen/Eigen/Core:376:
.../eigen/src/Eigen/Eigen/src/Core/arch/AltiVec/PacketMath.h:1025:13: error:
cannot initialize a variable of type 'Packet2bl' (aka 'Packet2ul') with an
rvalue of type '__vector __bool unsigned long long' (vector of 2
'unsigned long long' values)
...mask = vec_cmpeq(reinterpret_cast(select), reinterpret_cast(p2l_ONE));
1 error generated.

However, compiling an earlier version of Eigen (3.2.10) has no such error.  Can this error be fixed for compilation with Eigen?

Now that the POWER9 Summit/Sierra supercomputer has gone online, the use of XLC for HPC codes is expected to increase dramatically.

Thank you
Comment 1 Konstantinos Margaritis 2018-06-14 07:22:19 UTC

I just got my Power9 Talos II board and will set it up in the next days. I will investigate the xlc issue, though this one looks similar to the CLANG issue, in the code line ~780 there is the code:

typedef Packet2ul                    Packet2bl;
typedef __vector __bool long         Packet2bl;

Could you try adding a check for XLC in the first case?
Comment 2 Kate Clark 2018-06-14 17:41:45 UTC
(Passing on comments from one of my collaborators who looked at this)

After inspection of PacketMath.h from eigen-3.3.4, the suggested code snippet is already there. Macro EIGEN_COMP_CLANG is defined for XLC compiler, along with EIGEN_COMP_IBM macro.  (So I extended with an additional check for EIGEN_COMP_IBM, for consistency).

The problem was with line 1025, namely I added a reinterpret cast to type Packet2bl:
Packet2bl mask = reinterpret_cast<Packet2bl>( vec_cmpeq(reinterpret_cast<Packet2d>(select), reinterpret_cast<Packet2d>(p2l_ONE));
With this modification I was able to compile the code. I'm still curious why the problem did not arise with clang.
Comment 3 Nicole Trudeau 2018-06-20 21:41:28 UTC
One of my colleagues has been able to synthesize a short test, which he expects is representative of the Eigen diagnostic:


void foo(__vector double v) {
__vector bool long long r = vec_cmpeq(v, v);

Clang and g++ accept this because they allow lax vector conversions while the XL compiler does not. Clang emits a diagnostic as well if -fno-lax-vector-conversions is specified. The type of 'r' differs from the result type of the function. The diagnostic is valid and he would suggest the Eigen code be updated to improve portability.

Nicole Trudeau, Customer Support for IBM Compilers
Comment 4 Gael Guennebaud 2018-06-21 08:30:06 UTC
I applied the reinterpret_cast fix: (default) (3.3)
Comment 5 Samuel Matzek 2018-06-22 22:03:57 UTC
Line 1057 in the default commit is missing a lost paren at the end of the line:
Comment 6 Rasmus Munk Larsen 2018-06-22 22:23:41 UTC
I submitted a PR to fix the typo. Can somebody please merge it?
Comment 7 Rasmus Munk Larsen 2018-06-22 22:24:25 UTC
Comment 8 Kate Clark 2018-06-25 21:41:05 UTC
Even with this PR in place, there is still in issue with XLC using the latest version of Eigen:

/gpfs/gpfs_gl4_16mb/b7p319/b7p319ab/install/eigen-devel-2/include/eigen3/Eigen/src/Core/arch/AltiVec/PacketMath.h:945:71: error: 1540-0557 The asm operand does not match the specified constraint. __asm__ ("xvcmpgtdp %x0,%x2,%x1\n\txxsel %x0,%x1,%x2,%x0" : "=&wa" (ret) : "wa" (a), "wa" (b));

Note that the original patch does work on the latest release of Eigen, however, with above problem only appearing on the latest repository version of Eigen.
Comment 9 Gael Guennebaud 2018-10-08 20:59:37 UTC
This issue is related to bug 1494.

Note You need to log in before you can comment on or make changes to this bug.