void generate_conic_constraint(CoinPackedMatrix *& A, CoinPackedVector *& b,
CoinPackedVector *& d, double & h) {
fill_A(A);
b = new CoinPackedVector();
int b_indices[] = {0,1};
double b_values[] = {1.0, 2.0};
b->setVector(2, b_indices, b_values);
d = new CoinPackedVector();
int d_indices[] = {0,1,3};
double d_values[] = {-1.0, -1.0, 1.0};
d->setVector(3, d_indices, d_values);
h = 2.0;
}
int Choose_Intra_Mode (int *pcoeff, int *store_coeff, int xpos, int ypos, int newgob)
{
int E0[6], E1[8][6], E2[8][6], A[8][6], B[8][6];
int i, j;
int SAD0 = 0, SAD1 = 0, SAD2 = 0;
int Intra_Mode;
for (i = 0; i < 6; i++)
{
if (xpos == 0 && ypos == 0)
{
/* top left corner */
((i == 2 || i == 3)) ? fill_A (A, store_coeff, xpos, ypos, i - 2, i) : fill_null (A, i);
((i == 1 || i == 3)) ? fill_B (B, store_coeff, xpos, ypos, i - 1, i) : fill_null (B, i);
} else
{
if (xpos == 0)
{
/* left edge of the picture */
(i == 2 || i == 3) ? fill_A (A, store_coeff, xpos, ypos, i - 2, i) :
((i == 0 || i == 1) && !(newgob)) ? fill_A (A, store_coeff, xpos, ypos - 1, i + 2, i) :
((i == 4 || i == 5) && !(newgob)) ? fill_A (A, store_coeff, xpos, ypos - 1, i, i) :
fill_null (A, i);
(i == 1 || i == 3) ? fill_B (B, store_coeff, xpos, ypos, i - 1, i) : fill_null (B, i);
} else
{
if (ypos == 0)
{
/* top border of picture */
(i == 2 || i == 3) ? fill_A (A, store_coeff, xpos, ypos, i - 2, i) : fill_null (A, i);
((i == 4 || i == 5) && coded_map[ypos + 1][xpos] == 2) ?
fill_B (B, store_coeff, xpos - 1, ypos, i, i) :
(i == 1 || i == 3) ? fill_B (B, store_coeff, xpos, ypos, i - 1, i) :
((i == 0 || i == 2) && coded_map[ypos + 1][xpos] == 2) ?
fill_B (B, store_coeff, xpos - 1, ypos, i + 1, i) : fill_null (B, i);
} else
{
/* anywhere else in the picture, do not
* cross GOB boundary */
(i == 2 || i == 3) ? fill_A (A, store_coeff, xpos, ypos, i - 2, i) :
((i == 0 || i == 1) && !(newgob) && coded_map[ypos][xpos + 1] == 2) ?
fill_A (A, store_coeff, xpos, ypos - 1, i + 2, i) :
((i == 4 || i == 5) && !(newgob) && coded_map[ypos][xpos + 1] == 2) ?
fill_A (A, store_coeff, xpos, ypos - 1, i, i) :
fill_null (A, i);
((i == 4 || i == 5) && coded_map[ypos + 1][xpos] == 2) ?
fill_B (B, store_coeff, xpos - 1, ypos, i, i) :
(i == 1 || i == 3) ? fill_B (B, store_coeff, xpos, ypos, i - 1, i) :
((i == 0 || i == 2) && coded_map[ypos + 1][xpos] == 2) ?
fill_B (B, store_coeff, xpos - 1, ypos, i + 1, i) : fill_null (B, i);
}
}
}
/* DC prediction */
E0[i] = store_coeff[(xpos + ypos * pels / MB_SIZE) * 384 + i * 64]
- (A[0][i] + B[0][i] + 1) / 2;
for (j = 0; j < 8; j++)
{
/* DC/AC prediction from above (block A) */
E1[j][i] = store_coeff[(xpos + ypos * pels / MB_SIZE) * 384 + j + i * 64]
- A[j][i];
/* DC/AC prediction from left (block B) */
E2[j][i] = store_coeff[(xpos + ypos * pels / MB_SIZE) * 384 + j * 8 + i * 64]
- B[j][i];
}
/* Choose Intra Mode based on SAD of the prediction error */
SAD0 += E0[i];
SAD1 += E1[0][i];
SAD2 += E2[0][i];
for (j = 1; j < 8; j++)
{
SAD0 += 32 * abs (store_coeff[(xpos + ypos * pels / MB_SIZE) * 384 + j + i * 64])
+ 32 * abs (store_coeff[(xpos + ypos * pels / MB_SIZE) * 384 + j * 8 + i * 64]);
SAD1 += 32 * abs (E1[j][i]) + 32 * abs (store_coeff[(xpos + ypos * pels / MB_SIZE) * 384 + j * 8 + i * 64]);
SAD2 += 32 * abs (store_coeff[(xpos + ypos * pels / MB_SIZE) * 384 + j + i * 64]) + 32 * abs (E2[j][i]);
}
}
if (SAD0 <= SAD1)
if (SAD0 <= SAD2)
Intra_Mode = INTRA_MODE_DC;
else
Intra_Mode = INTRA_MODE_HORI_AC;
else if (SAD1 <= SAD2)
Intra_Mode = INTRA_MODE_VERT_AC;
else
Intra_Mode = INTRA_MODE_HORI_AC;
#ifdef PRINTBLOCKS
for (i = 0; i < 6; i++)
{
printf ("\nMB: %2i Block: %1i\n", xpos + ypos * pels / MB_SIZE, i);
printf ("qcoeff: \n");
for (j = 0; j < 8; j++)
{
printf ("%4i%4i%4i%4i%4i%4i%4i%4i\n", pcoeff[j * 8 + i * 64], pcoeff[j * 8 + 1 + i * 64], pcoeff[j * 8 + 2 + i * 64], pcoeff[j * 8 + 3 + i * 64], pcoeff[j * 8 + 4 + i * 64], pcoeff[j * 8 + 5 + i * 64], pcoeff[j * 8 + 6 + i * 64], pcoeff[j * 8 + 7 + i * 64]);
}
printf ("\nA[%1i]: ", i);
for (j = 0; j < 8; j++)
{
printf ("%4i", A[j][i]);
}
printf ("\nB[%1i]: ", i);
for (j = 0; j < 8; j++)
{
printf ("%4i", B[j][i]);
}
printf ("\nE0[%1i]: %3i", i, E0[i]);
printf ("\nE1[%1i]: ", i);
for (j = 0; j < 8; j++)
{
printf ("%4i", E1[j][i]);
}
printf ("\nE2[%1i]: ", i);
for (j = 0; j < 8; j++)
{
printf ("%4i", E2[j][i]);
}
}
printf ("\nSAD0: %4i, SAD1: %4i, SAD2: %4i\n", SAD0, SAD1, SAD2);
#endif
return Intra_Mode;
}
void Intra_AC_DC_Decode (int *rcoeff, int *store_rcoeff, int INTRA_AC_DC, int xpos, int ypos, int newgob, int i)
{
int A[8][6], B[8][6];
int j, k, tempDC;
if (xpos == 0 && ypos == 0)
{
/* top left corner */
(i == 2 || i == 3) ? fill_A (A, store_rcoeff, xpos, ypos, i - 2, i) : fill_null (A, i);
(i == 1 || i == 3) ? fill_B (B, store_rcoeff, xpos, ypos, i - 1, i) : fill_null (B, i);
} else
{ /* left border of picture */
if (xpos == 0)
{ /* left edge of the picture */
(i == 2 || i == 3) ? fill_A (A, store_rcoeff, xpos, ypos, i - 2, i) :
((i == 0 || i == 1) && !(newgob)) ? fill_A (A, store_rcoeff, xpos, ypos - 1, i + 2, i) :
((i == 4 || i == 5) && !(newgob)) ? fill_A (A, store_rcoeff, xpos, ypos - 1, i, i) : fill_null (A, i);
(i == 1 || i == 3) ? fill_B (B, store_rcoeff, xpos, ypos, i - 1, i) : fill_null (B, i);
} else
{
if (ypos == 0)
{ /* top border of picture */
(i == 2 || i == 3) ? fill_A (A, store_rcoeff, xpos, ypos, i - 2, i) : fill_null (A, i);
(i == 4 || i == 5) ? fill_B (B, store_rcoeff, xpos - 1, ypos, i, i) :
(i == 1 || i == 3) ? fill_B (B, store_rcoeff, xpos, ypos, i - 1, i) :
fill_B (B, store_rcoeff, xpos - 1, ypos, i + 1, i);
} else
{
/* anywhere else in the picture, do not
* cross GOB boundary */
(i == 2 || i == 3) ? fill_A (A, store_rcoeff, xpos, ypos, i - 2, i) :
((i == 0 || i == 1) && !(newgob)) ? fill_A (A, store_rcoeff, xpos, ypos - 1, i + 2, i) :
((i == 4 || i == 5) && !(newgob)) ? fill_A (A, store_rcoeff, xpos, ypos - 1, i, i) : fill_null (A, i);
(i == 4 || i == 5) ? fill_B (B, store_rcoeff, xpos - 1, ypos, i, i) :
(i == 1 || i == 3) ? fill_B (B, store_rcoeff, xpos, ypos, i - 1, i) :
fill_B (B, store_rcoeff, xpos - 1, ypos, i + 1, i);
}
}
}
switch (INTRA_AC_DC)
{
case INTRA_MODE_DC:
/* It is OK to use 1024 to check the DC direction since the DC *
* reconstructed coefficient will always be odd */
tempDC = rcoeff[i * 64]
+ ((A[0][i] == 1024 && B[0][i] == 1024) ? 1024 :
(A[0][i] == 1024) ? B[0][i] :
(B[0][i] == 1024) ? A[0][i] : (A[0][i] + B[0][i]) / 2);
for (j = 0; j < 8; j++)
for (k = 0; k < 8; k++)
rcoeff[i * 64 + j * 8 + k] = clipAC (rcoeff[i * 64 + j * 8 + k]);
rcoeff[i * 64] = oddifyclipDC (tempDC);
break;
case INTRA_MODE_VERT_AC:
tempDC = rcoeff[i * 64] + A[0][i];
for (j = 1; j < 8; j++)
{
rcoeff[i * 64 + j] = clipAC (rcoeff[i * 64 + j] + A[j][i]);
}
for (j = 1; j < 8; j++)
for (k = 0; k < 8; k++)
rcoeff[i * 64 + j * 8 + k] = clipAC (rcoeff[i * 64 + j * 8 + k]);
rcoeff[i * 64] = oddifyclipDC (tempDC);
break;
case INTRA_MODE_HORI_AC:
tempDC = rcoeff[i * 64] + B[0][i];
for (j = 1; j < 8; j++)
rcoeff[i * 64 + j * 8] = clipAC (rcoeff[i * 64 + j * 8] + B[j][i]);
for (j = 0; j < 8; j++)
for (k = 1; k < 8; k++)
rcoeff[i * 64 + j * 8 + k] = clipAC (rcoeff[i * 64 + j * 8 + k]);
rcoeff[i * 64] = oddifyclipDC (tempDC);
break;
default:
printf ("Error in Prediction in Advanced Intra Coding\n");
exit (-1);
break;
}
/* store the reconstructed DCT coefficients */
memcpy ((void *) (store_rcoeff + i * 64 + (xpos + ypos * pels / MB_SIZE) * 384), (void *) (rcoeff + i * 64), sizeof (int) * 64);
#ifdef PRINTBLOCKS
printf ("MB: %2i Block: %1i\n", xpos + ypos * pels / MB_SIZE, i);
printf ("rcoeff: \n");
for (j = 0; j < 8; j++)
{
printf ("%4i%4i%4i%4i%4i%4i%4i%4i\n", rcoeff[j * 8 + i * 64], rcoeff[j * 8 + 1 + i * 64], rcoeff[j * 8 + 2 + i * 64], rcoeff[j * 8 + 3 + i * 64], rcoeff[j * 8 + 4 + i * 64], rcoeff[j * 8 + 5 + i * 64], rcoeff[j * 8 + 6 + i * 64], rcoeff[j * 8 + 7 + i * 64]);
}
#endif
return;
}
void Intra_AC_DC_Encode (int *pcoeff, int *store_rcoeff, int Intra_Mode, int xpos, int ypos, int newgob, int i)
{
int A[8][6], B[8][6];
int j;
if (xpos == 0 && ypos == 0)
{
/* top left corner */
(i == 2 || i == 3) ? fill_A (A, store_rcoeff, xpos, ypos, i - 2, i) : fill_null (A, i);
(i == 1 || i == 3) ? fill_B (B, store_rcoeff, xpos, ypos, i - 1, i) : fill_null (B, i);
} else
{
/* left border of picture */
if (xpos == 0)
{ /* left edge of the picture */
(i == 2 || i == 3) ? fill_A (A, store_rcoeff, xpos, ypos, i - 2, i) :
((i == 0 || i == 1) && !(newgob)) ? fill_A (A, store_rcoeff, xpos, ypos - 1, i + 2, i) :
((i == 4 || i == 5) && !(newgob)) ? fill_A (A, store_rcoeff, xpos, ypos - 1, i, i) : fill_null (A, i);
(i == 1 || i == 3) ? fill_B (B, store_rcoeff, xpos, ypos, i - 1, i) : fill_null (B, i);
} else
{
if (ypos == 0)
{
/* top border of picture */
(i == 2 || i == 3) ? fill_A (A, store_rcoeff, xpos, ypos, i - 2, i) : fill_null (A, i);
(i == 4 || i == 5) ? fill_B (B, store_rcoeff, xpos - 1, ypos, i, i) :
(i == 1 || i == 3) ? fill_B (B, store_rcoeff, xpos, ypos, i - 1, i) :
fill_B (B, store_rcoeff, xpos - 1, ypos, i + 1, i);
} else
{
/* anywhere else in the picture, do not
* cross GOB boundary */
(i == 2 || i == 3) ? fill_A (A, store_rcoeff, xpos, ypos, i - 2, i) :
((i == 0 || i == 1) && !(newgob)) ? fill_A (A, store_rcoeff, xpos, ypos - 1, i + 2, i) :
((i == 4 || i == 5) && !(newgob)) ? fill_A (A, store_rcoeff, xpos, ypos - 1, i, i) : fill_null (A, i);
(i == 4 || i == 5) ? fill_B (B, store_rcoeff, xpos - 1, ypos, i, i) :
(i == 1 || i == 3) ? fill_B (B, store_rcoeff, xpos, ypos, i - 1, i) :
fill_B (B, store_rcoeff, xpos - 1, ypos, i + 1, i);
}
}
}
#ifdef PRINTBLOCKS
printf ("Intra_Mode: %i", Intra_Mode);
printf ("\nA[%1i]: ", i);
for (j = 0; j < 8; j++)
{
printf ("%4i", A[j][i]);
}
printf ("\nB[%1i]: ", i);
for (j = 0; j < 8; j++)
{
printf ("%4i", B[j][i]);
}
#endif
/* replace the qcoeff with the predicted values pcoeff */
switch (Intra_Mode)
{
case INTRA_MODE_DC:
/* It is OK to use 1024 to check the DC direction since the DC *
* reconstructed coefficient will always be odd */
pcoeff[i * 64] -=
(A[0][i] == 1024 && B[0][i] == 1024) ? 1024 :
(A[0][i] == 1024) ? B[0][i] :
(B[0][i] == 1024) ? A[0][i] : ((A[0][i] + B[0][i]) / 2);
break;
case INTRA_MODE_VERT_AC:
for (j = 0; j < 8; j++)
{
pcoeff[j + i * 64] -= A[j][i];
}
break;
case INTRA_MODE_HORI_AC:
for (j = 0; j < 8; j++)
{
pcoeff[j * 8 + i * 64] -= B[j][i];
}
break;
default:
printf ("Error in Prediction in Advanced Intra Coding\n");
exit (-1);
break;
}
#ifdef PRINTBLOCKS
printf ("MB: %2i Block: %1i\n", xpos + ypos * pels / MB_SIZE, i);
printf ("pcoeff: \n");
for (j = 0; j < 8; j++)
{
printf ("%4i%4i%4i%4i%4i%4i%4i%4i\n", pcoeff[j * 8 + i * 64], pcoeff[j * 8 + 1 + i * 64], pcoeff[j * 8 + 2 + i * 64], pcoeff[j * 8 + 3 + i * 64], pcoeff[j * 8 + 4 + i * 64], pcoeff[j * 8 + 5 + i * 64], pcoeff[j * 8 + 6 + i * 64], pcoeff[j * 8 + 7 + i * 64]);
}
#endif
return;
}