Frama-C

Informations légales et droit de diffusion

L'ensemble du site relève de la législation française et internationale sur le droit d'auteur et de la propriété intellectuelle. Tous les droits de reproduction sont réservés, y compris pour les documents téléchargeables et les représentations iconographiques et photographiques.

La reproduction de tout ou partie de ce site sur un support électronique quel qu'il soit est interdite sauf autorisation préalable du directeur de publication.

Directeur de la publication : Claude Marché
Entreprise éditrice : INRIA, centre de recherche de Saclay - Île-de-France
4, rue Jacques Monod, 91893 Orsay cedex - France
Hébergeur : OVH, 2 rue Kellermann - 59100 Roubaix - France

Aucune information personnelle n'est collectée à votre insu, aucune information personnelle n'est cédée à des tiers, aucune information personnelle n'est utilisée à des fins non prévues.

Legal information and diffusion rights

All materials and graphical elements contained on this site are protected by French and international copyright law and intellectual property rights. Reproduction in part or in whole of materials and graphical elements of this site is prohibited without the express permission of the Director of publication.

Déclaration CNIL

Conformément à la loi n° 78-17 du 6 janvier 1978, relative à l'Informatique, aux Fichiers et aux Libertés, vous disposez d'un droit d'accès et de rectification des informations nominatives vous concernant. Vous pouvez l'exercer en vous adressant à :

Florent.Kirchner

cea.fr
F. Kirchner
CEA Saclay
F91191 Gif sur Yvette cedex

static int x = 1; static long mc1[8][8], mc2[8][8]; /* if (init) { for (;;) } (long m1[8][8], long m2[8][8]) { long i, j, k, tmp1[8][8], tmp2[8][8]; /* Loops indexes and temporary matrices. */ double ftmp1, ftmp2; static int init = 1; static long mc1[8][8], mc2[8][8]; /* hard-coded cosines matrices. */if (init) { for (i = 0; i < 8; i++) for (j = 0; j < 8; j++) { ftmp1 = ((j == 0) ? 0; static long mc1[8][8], mc2[8][8]; /* hard-coded cosines matrices. */if (init) { for (i = 0; i < 8; i++) for (j = 0; j < 8; j++) { ftmp1 = ((j == 0) ? 0.5 / sqrt (2.0) : 0.5) * cos ((2.0 * i + 1.0) * j * TH); ftmp2 = ftmp1; /* The well known formula. The max absolute value for ftmp1 and ftmp2 is 0.5. */ ftmp1 *= (1 << NBC1); if (ftmp1 < 0) ftmp1 -= 0.5; else ftmp1 += 0.5; mc1[i][j] = ftmp1; ftmp2 *= (1 << NBC2); /* Multiply the cosine coefficient by 2^NBC2. The max absolute value for * ftmp2 is 2^(NBC2-1). */ if (ftmp2 < 0) ftmp2 –= 0.5; else ftmp2 += 0.5; /* For symetrical rounding. */ mc2[i][j] = ftmp2; } init = 0; } /* Then the first pass. */ for(i = 0; i < 8; i++) for(j = 0; j < 8; j++) { for(k = 0, tmp1[i][j] = 0; k < 8; k++)tmp1[i][j] += mc1[i][k] * m1[k][j]; /* The [i,j] coefficient of the matrix product MC1*M1. */ tmp1[i][j] >>= (NBC1 + 10 – NBI); tmp1[i][j] += 1; /* For rounding purpose. */ tmp1[i][j] >>= 1; /* Final rounding. tmp1[i][j] is now represented on NBI bits. */if (tmp1[i][j] < -(1 <<

Informations légales et droit de diffusion

Legal information and diffusion rights

Contacts

Contacts

Déclaration CNIL