Viewing file:  pr104152.c (1.05 KB) -rw-r--r--
Select action/file-type:
 (+) |  (+) |  (+) | Code (+) | Session (+) |  (+) | SDB (+) |  (+) |  (+) |  (+) |  (+) |  (+) |
/* { dg-do compile } */
/* { dg-additional-options "-Ofast" } */
/* { dg-additional-options "-march=armv8.2-a+sve" { target aarch64-*-* } } */
#define M_PI 3.14f
#define NK 24
#define DIM 3
typedef float rvec_ [3];
static rvec_ v0[NK] = {
{ 1, 0, 0 }, { 0, 1, 0 }, { 0, 0, 1 }, { 1, 1, 0 }, { 1, -1, 0 },
{ 1, 0, 1 }, { 1, 0, -1 }, { 0, 1, 1 }, { 0, 1, -1 }, { 1, 1, 1 },
{ 1, 1, -1 }, { 1, -1, 1 }, { -1, 1, 1 }, { 2, 0, 0 }, { 0, 2, 0 },
{ 0, 0, 2 }, { 3, 0, 0 }, { 0, 3, 0 }, { 0, 0, 3 }, { 4, 0, 0 },
{ 0, 4, 0 }, { 0, 0, 4 } };
static inline float iprod__(const rvec_ a, const rvec_ b)
{
return (a[0] * b[0] + a[1] * b[1] + a[2] * b[2]);
}
int badaboum(rvec_ cm_mol, float **tc)
{
float sx;
int k, d;
rvec_ kfac[3];
for (k = 0; k < DIM; k++)
for (d = 0; d < DIM; d++)
kfac[k][d] = M_PI * v0[k][d] / v0[d][d];
for (k = 0; k < DIM; k++)
{
sx = __builtin_sinf(iprod__(kfac[k], cm_mol));
tc[k][0] += sx * iprod__(v0[k], cm_mol);
}
return 0;
}
|
