Viewing file:  lxc1-sxc1-1.c (1.16 KB) -rw-r--r--
Select action/file-type:
 (+) |  (+) |  (+) | Code (+) | Session (+) |  (+) | SDB (+) |  (+) |  (+) |  (+) |  (+) |  (+) |
/* { dg-options "(HAS_LXC1) -mhard-float -mno-lxc1-sxc1" } */
/* { dg-final { scan-assembler-not "\tldxc1\t" } } */
/* { dg-final { scan-assembler-not "\tsdxc1\t" } } */
#ifndef __mips_no_lxc1_sxc1
#error missing definition of __mips_no_lxc1_sxc1
#endif
double ldexp(double x, int exp);
typedef struct
{
double** rows;
} d_mat_struct;
typedef d_mat_struct d_mat_t[1];
#define d_mat_entry(mat,i,j) (*((mat)->rows[i] + (j)))
double __attribute__((noinline))
ldxc1_test (int kappa, int zeros, double ctt, int* expo, d_mat_t r, double* s)
{
int kappa2 = kappa;
double tmp = 0.0;
do
{
kappa--;
if (kappa > zeros + 1)
{
tmp = d_mat_entry(r, kappa - 1, kappa - 1) * ctt;
tmp = ldexp(tmp, (expo[kappa - 1] - expo[kappa2]));
}
}
while ((kappa >= zeros + 2) && (s[kappa - 1] <= tmp));
return tmp;
}
#define SIZE 20
int main(void)
{
int kappa = SIZE - 1;
int zeros = 1;
double ctt = 2;
int expo[SIZE] = {0};
double s[SIZE] = {0};
double rows_data[SIZE][SIZE] = {0};
double* rows[SIZE];
for (int i = 0; i < SIZE; i++)
rows[i] = rows_data[i];
d_mat_t r = { rows };
ldxc1_test(kappa, zeros, ctt, expo, r, s);
return 0;
}
|
