int DoTest(int N, TYPE *alpha0, int incX, int incY)
{
int iret;
const int npad=Mmax(4*Mabs(incY), 16);
const TYPE padval=(-2271.0);
TYPE *Yg, *Yt, *X, *x, *y;
#ifdef TREAL
TYPE alpha = *alpha0;
#else
TYPE *alpha = alpha0;
#endif
Yg = getvec(npad, padval, N, incY);
Yt = dupvec(npad, N, Yg, incY);
X = getvec(0, padval, N, incX); /* no padding for read-only X */
x = X;
y = Yg + (npad SHIFT);
if (incX < 1) x -= ((N-1)SHIFT) * incX;
if (incY < 1) y -= ((N-1)SHIFT) * incY;
good_axpy(N, alpha, x, incX, y, incY);
y = Yt + (npad SHIFT);
if (incY < 1) y -= ((N-1)SHIFT) * incY;
TEST_AXPY(N, alpha, x, incX, y, incY);
iret = CheckY(npad, padval, N, Yg, incY, Yt, incY);
free(X);
free(Yg);
free(Yt);
return(iret);
}
int DoTest(int N, TYPE *alpha0, int incY)
{
int iret;
const int npad=Mmax(4*Mabs(incY), 16);
const TYPE padval=(-2271.0);
TYPE *Yg, *Yt, *y;
#ifdef TREAL
TYPE alpha = *alpha0;
#else
TYPE *alpha = alpha0;
#endif
Yg = getvec(npad, padval, N, incY);
Yt = dupvec(npad, N, Yg, incY);
y = Yg + (npad SHIFT);
if (incY < 1) y -= ((N-1)SHIFT) * incY;
good_set(N, alpha, y, incY);
y = Yt + (npad SHIFT);
if (incY < 1) y -= ((N-1)SHIFT) * incY;
TEST_SET(N, alpha, Yt+(npad SHIFT), incY);
iret = CheckY(npad, padval, N, Yg, incY, Yt, incY);
free(Yg);
free(Yt);
return(iret);
}
void IsIntersection(int x, int y)
{
int win_x = CheckX(x);
int win_y = CheckY(y);
if(win_x == YES && win_y == YES)
{
WinMessage();
display_type = PAUSE;
NewGame();
}
}
