Skip to content

Commit 18a2304

Browse files
authored
Add files via upload
1 parent a29f20e commit 18a2304

1 file changed

Lines changed: 14 additions & 33 deletions

File tree

kernel/arm/zscal.c

Lines changed: 14 additions & 33 deletions
Original file line numberDiff line numberDiff line change
@@ -48,46 +48,27 @@ int CNAME(BLASLONG n, BLASLONG dummy0, BLASLONG dummy1, FLOAT da_r,FLOAT da_i, F
4848

4949
if ((n <= 0) || (inc_x <= 0))
5050
return(0);
51-
printf("arm zscal kernel\n");
51+
5252
inc_x2 = 2 * inc_x;
5353
if (dummy2 == 0) {
54-
printf("arm zscal kernel, dummy2 is zero\n");
55-
for ( i=0; i<n; i++ )
54+
for (i = 0; i < n; i++)
5655
{
57-
if ( da_r == 0.0 )
58-
{
59-
if ( da_i == 0.0 )
60-
{
61-
temp = 0.0;
62-
x[ip+1] = 0.0 ;
63-
}
64-
else
65-
{
66-
temp = - da_i * x[ip+1] ;
67-
x[ip+1] = da_i * x[ip] ;
68-
}
69-
}
70-
else
71-
{
72-
if ( da_i == 0.0 )
73-
{
74-
temp = da_r * x[ip] ;
75-
x[ip+1] = da_r * x[ip+1];
76-
}
77-
else
78-
{
79-
temp = da_r * x[ip] - da_i * x[ip+1] ;
80-
x[ip+1] = da_r * x[ip+1] + da_i * x[ip] ;
81-
}
82-
}
83-
x[ip] = temp;
56+
if (da_r == 0.0 && da_i == 0.0)
57+
{
58+
x[ip] = 0.0;
59+
x[ip+1] = 0.0;
60+
}
61+
else
62+
{
63+
temp = da_r * x[ip] - da_i * x[ip+1];
64+
x[ip+1] = da_r * x[ip+1] + da_i * x[ip] ;
65+
x[ip] = temp;
66+
}
8467

85-
ip += inc_x2;
68+
ip += inc_x2;
8669
}
87-
8870
return(0);
8971
}
90-
printf("arm zscal kernel in nan-propagating mode\n");
9172
for (i = 0; i < n; i++)
9273
{
9374
temp = da_r * x[ip] - da_i * x[ip+1];

0 commit comments

Comments
 (0)