//-----------------------------------------------------------------------------
UINT64 Crc64( UINT64 start_crc , UINT8 *addr ,UINT32 size)
{
UINT8 * limit;
UINT32 u32tmp ;
PAGED_CODE ();
start_crc = ~start_crc;
if (size > 4)
{
while ((uintptr_t)(addr) & 3)
{
start_crc = crc64_table[0][*addr++ ^ A1(start_crc)] ^ S8(start_crc);
--size;
}
limit = addr + (size & ~(size_t)(3));
size &= (size_t)(3);
while (addr < limit)
{
u32tmp = (UINT32)start_crc ^ (*(UINT32*)(addr));
addr += 4;
start_crc = crc64_table[3][A(u32tmp)] ^ crc64_table[2][B(u32tmp)] ^ S32(start_crc)^ crc64_table[1][C(u32tmp)] ^ crc64_table[0][D(u32tmp)];
}
}
while (size-- != 0)
start_crc = crc64_table[0][*addr++ ^ A1(start_crc)] ^ S8(start_crc);
return ~start_crc;
}
///
/// Generates source code for the king moves
///
void CFiestyGen::genKingAttacks()
{
std::cout << "const YBitBoard CGen::mbbKingAttacks[CSqix::kNumSquares] = {";
for ( U8 sq = 0; sq < CSqix::kNumSquares; sq++ )
{
CSqix sqix( sq );
S8 r = S8( sqix.getRank().get() );
S8 f = S8( sqix.getFile().get() );
CBitBoard bb( 0ULL );
bb.setIfValid( r + 1, f - 1 );
bb.setIfValid( r + 1, f );
bb.setIfValid( r + 1, f + 1 );
bb.setIfValid( r, f - 1 );
bb.setIfValid( r, f + 1 );
bb.setIfValid( r - 1, f - 1 );
bb.setIfValid( r - 1, f );
bb.setIfValid( r - 1, f + 1 );
#ifdef BBTRACE
printBitBoardDiagram( sqix.asAbbr(), bb );
#endif
std::cout << "\n /* " << sqix.asAbbr() << " */ ";
std::cout << bb.asAbbr() << "ULL";
if( sq != CSqix::kNumSquares - 1 )
std::cout << ",";
}
std::cout << " };\n" << std::flush;
}
// If you make any changes, do some bench marking! Seemingly unrelated
// changes can very easily ruin the performance (and very probably is
// very compiler dependent).
extern uint32_t
lzma_crc32(const uint8_t *buf, size_t size, uint32_t crc)
{
crc = ~crc;
#ifdef WORDS_BIGENDIAN
crc = bswap32(crc);
#endif
if (size > 8) {
// Fix the alignment, if needed. The if statement above
// ensures that this won't read past the end of buf[].
while ((uintptr_t)(buf) & 7) {
crc = lzma_crc32_table[0][*buf++ ^ A(crc)] ^ S8(crc);
--size;
}
// Calculate the position where to stop.
const uint8_t *const limit = buf + (size & ~(size_t)(7));
// Calculate how many bytes must be calculated separately
// before returning the result.
size &= (size_t)(7);
// Calculate the CRC32 using the slice-by-eight algorithm.
while (buf < limit) {
crc ^= *(const uint32_t *)(buf);
buf += 4;
crc = lzma_crc32_table[7][A(crc)]
^ lzma_crc32_table[6][B(crc)]
^ lzma_crc32_table[5][C(crc)]
^ lzma_crc32_table[4][D(crc)];
const uint32_t tmp = *(const uint32_t *)(buf);
buf += 4;
// At least with some compilers, it is critical for
// performance, that the crc variable is XORed
// between the two table-lookup pairs.
crc = lzma_crc32_table[3][A(tmp)]
^ lzma_crc32_table[2][B(tmp)]
^ crc
^ lzma_crc32_table[1][C(tmp)]
^ lzma_crc32_table[0][D(tmp)];
}
}
while (size-- != 0)
crc = lzma_crc32_table[0][*buf++ ^ A(crc)] ^ S8(crc);
#ifdef WORDS_BIGENDIAN
crc = bswap32(crc);
#endif
return ~crc;
}
extern uint64_t
lzma_crc64_8bchk(const uint8_t *buf, size_t size, uint64_t crc, uint64_t *cnt)
{
crc = ~crc;
uint64_t cnt1 = 0;
#ifdef WORDS_BIGENDIAN
crc = bswap64(crc);
#endif
if (size > 4) {
while ((uintptr_t)(buf) & 3) {
crc = lzma_crc64_table[0][*buf++ ^ A1(crc)] ^ S8(crc);
--size;
}
const uint8_t *const limit = buf + (size & ~(size_t)(3));
size &= (size_t)(3);
while (buf < limit) {
uint32_t val = *(const uint32_t *)(buf);
#ifdef WORDS_BIGENDIAN
const uint32_t tmp = (crc >> 32) ^ val;
#else
const uint32_t tmp = crc ^ val;
#endif
buf += 4;
val >>= 7; cnt1 += (val & 1);
val >>= 8; cnt1 += (val & 1);
crc = lzma_crc64_table[3][A(tmp)]
^ lzma_crc64_table[2][B(tmp)]
^ S32(crc)
^ lzma_crc64_table[1][C(tmp)]
^ lzma_crc64_table[0][D(tmp)];
}
}
while (size-- != 0)
crc = lzma_crc64_table[0][*buf++ ^ A1(crc)] ^ S8(crc);
#ifdef WORDS_BIGENDIAN
crc = bswap64(crc);
#endif
*cnt = cnt1;
return ~crc;
}
uint64_t crc64_update(uint8_t* buf, size_t size, uint64_t crc)
{
crc = ~crc;
#ifdef WORDS_BIGENDIAN
crc = bswap64(crc);
#endif
if (size > 4) {
while ((uintptr_t)(buf) & 3) {
crc = crc64_table[0][*buf++ ^ A1(crc)] ^ S8(crc);
--size;
}
const uint8_t *const limit = buf + (size & ~(size_t)(3));
size &= (size_t)(3);
while (buf < limit) {
#ifdef WORDS_BIGENDIAN
const uint32_t tmp = (crc >> 32)
^ *(const uint32_t *)(buf);
#else
const uint32_t tmp = crc ^ *(const uint32_t *)(buf);
#endif
buf += 4;
crc = crc64_table[3][A(tmp)]
^ crc64_table[2][B(tmp)]
^ S32(crc)
^ crc64_table[1][C(tmp)]
^ crc64_table[0][D(tmp)];
}
}
while (size-- != 0)
crc = crc64_table[0][*buf++ ^ A1(crc)] ^ S8(crc);
#ifdef WORDS_BIGENDIAN
crc = bswap64(crc);
#endif
return ~crc;
return 0;
}
///
/// Generates source code for the knight moves
///
void CFiestyGen::genKnightAttacks()
{
std::cout << "const YBitBoard CGen::mbbKnightAttacks[CSqix::kNumSquares] = {";
for ( U8 sq = 0; sq < CSqix::kNumSquares; sq++ )
{
CSqix sqix( sq );
S8 r = S8( sqix.getRank().get() );
S8 f = S8( sqix.getFile().get() );
std::cout << "\n /* " << sqix.asAbbr() << " */ ";
CBitBoard bb( 0ULL );
bb.setIfValid( r + 2, f + 1 );
bb.setIfValid( r + 2, f - 1 );
bb.setIfValid( r - 2, f + 1 );
bb.setIfValid( r - 2, f - 1 );
bb.setIfValid( r + 1, f + 2 );
bb.setIfValid( r + 1, f - 2 );
bb.setIfValid( r - 1, f + 2 );
bb.setIfValid( r - 1, f - 2 );
std::cout << bb.asAbbr() << "ULL";
if( sq != CSqix::kNumSquares - 1 )
std::cout << ",";
}
std::cout << " };\n" << std::flush;
}
AbstractModel::SufficientStatisticsVector
MultivariateModel
::GetSufficientStatistics(const Realizations &R, const Observations& Obs)
{
/// Computation of the suffisient statistics of the model
/// Basically, it is a vector (named S1 to S9) of vectors
/// S1 <- y_ij * eta_ij & S2 <- eta_ij * eta_ij
VectorType S1(m_NbTotalOfObservations), S2(m_NbTotalOfObservations);
auto itS1 = S1.begin(), itS2 = S2.begin();
for(size_t i = 0; i < m_NumberOfSubjects; ++i)
{
for(size_t j = 0; j < Obs.GetNumberOfTimePoints(i); ++j)
{
VectorType PC = ComputeParallelCurve(i, j);
*itS1 = dot_product(PC, Obs.GetSubjectCognitiveScore(i, j));
*itS2 = PC.squared_magnitude();
++itS1, ++itS2;
}
}
/// S3 <- Ksi_i & S4 <- Ksi_i * Ksi_i
VectorType S3 = R.at("Ksi");
VectorType S4 = R.at("Ksi") % R.at("Ksi");
/// S5 <- Tau_i & S6 <- Tau_i * Tau_i
VectorType S5 = R.at("Tau");
VectorType S6 = R.at("Tau") % R.at("Tau");
/// S7 <- G
VectorType S7(1, R.at("G", 0));
/// S8 <- beta_k
VectorType S8((m_ManifoldDimension-1) * m_NbIndependentSources);
ScalarType * itS8 = S8.memptr();
for(size_t i = 0; i < S8.size(); ++i)
itS8[i] = R.at("Beta#" + std::to_string(i), 0);
/// S8 <- delta_k
VectorType S9(m_ManifoldDimension - 1);
ScalarType * itS9 = S9.memptr();
for(size_t i = 1; i < S9.size(); ++i)
itS9[i] = R.at("Delta#" + std::to_string(i), 0);
return {S1, S2, S3, S4, S5, S6, S7, S8, S9};
}
static void restart(struct net_device *dev)
{
struct fs_enet_private *fep = netdev_priv(dev);
const struct fs_platform_info *fpi = fep->fpi;
fcc_t __iomem *fccp = fep->fcc.fccp;
fcc_c_t __iomem *fcccp = fep->fcc.fcccp;
fcc_enet_t __iomem *ep = fep->fcc.ep;
dma_addr_t rx_bd_base_phys, tx_bd_base_phys;
u16 paddrh, paddrm, paddrl;
const unsigned char *mac;
int i;
C32(fccp, fcc_gfmr, FCC_GFMR_ENR | FCC_GFMR_ENT);
/* clear everything (slow & steady does it) */
for (i = 0; i < sizeof(*ep); i++)
out_8((u8 __iomem *)ep + i, 0);
/* get physical address */
rx_bd_base_phys = fep->ring_mem_addr;
tx_bd_base_phys = rx_bd_base_phys + sizeof(cbd_t) * fpi->rx_ring;
/* point to bds */
W32(ep, fen_genfcc.fcc_rbase, rx_bd_base_phys);
W32(ep, fen_genfcc.fcc_tbase, tx_bd_base_phys);
/* Set maximum bytes per receive buffer.
* It must be a multiple of 32.
*/
W16(ep, fen_genfcc.fcc_mrblr, PKT_MAXBLR_SIZE);
W32(ep, fen_genfcc.fcc_rstate, (CPMFCR_GBL | CPMFCR_EB) << 24);
W32(ep, fen_genfcc.fcc_tstate, (CPMFCR_GBL | CPMFCR_EB) << 24);
/* Allocate space in the reserved FCC area of DPRAM for the
* internal buffers. No one uses this space (yet), so we
* can do this. Later, we will add resource management for
* this area.
*/
W16(ep, fen_genfcc.fcc_riptr, fpi->dpram_offset);
W16(ep, fen_genfcc.fcc_tiptr, fpi->dpram_offset + 32);
W16(ep, fen_padptr, fpi->dpram_offset + 64);
/* fill with special symbol... */
memset_io(fep->fcc.mem + fpi->dpram_offset + 64, 0x88, 32);
W32(ep, fen_genfcc.fcc_rbptr, 0);
W32(ep, fen_genfcc.fcc_tbptr, 0);
W32(ep, fen_genfcc.fcc_rcrc, 0);
W32(ep, fen_genfcc.fcc_tcrc, 0);
W16(ep, fen_genfcc.fcc_res1, 0);
W32(ep, fen_genfcc.fcc_res2, 0);
/* no CAM */
W32(ep, fen_camptr, 0);
/* Set CRC preset and mask */
W32(ep, fen_cmask, 0xdebb20e3);
W32(ep, fen_cpres, 0xffffffff);
W32(ep, fen_crcec, 0); /* CRC Error counter */
W32(ep, fen_alec, 0); /* alignment error counter */
W32(ep, fen_disfc, 0); /* discard frame counter */
W16(ep, fen_retlim, 15); /* Retry limit threshold */
W16(ep, fen_pper, 0); /* Normal persistence */
/* set group address */
W32(ep, fen_gaddrh, fep->fcc.gaddrh);
W32(ep, fen_gaddrl, fep->fcc.gaddrh);
/* Clear hash filter tables */
W32(ep, fen_iaddrh, 0);
W32(ep, fen_iaddrl, 0);
/* Clear the Out-of-sequence TxBD */
W16(ep, fen_tfcstat, 0);
W16(ep, fen_tfclen, 0);
W32(ep, fen_tfcptr, 0);
W16(ep, fen_mflr, PKT_MAXBUF_SIZE); /* maximum frame length register */
W16(ep, fen_minflr, PKT_MINBUF_SIZE); /* minimum frame length register */
/* set address */
mac = dev->dev_addr;
paddrh = ((u16)mac[5] << 8) | mac[4];
paddrm = ((u16)mac[3] << 8) | mac[2];
paddrl = ((u16)mac[1] << 8) | mac[0];
W16(ep, fen_paddrh, paddrh);
W16(ep, fen_paddrm, paddrm);
W16(ep, fen_paddrl, paddrl);
W16(ep, fen_taddrh, 0);
W16(ep, fen_taddrm, 0);
W16(ep, fen_taddrl, 0);
W16(ep, fen_maxd1, 1520); /* maximum DMA1 length */
W16(ep, fen_maxd2, 1520); /* maximum DMA2 length */
/* Clear stat counters, in case we ever enable RMON */
W32(ep, fen_octc, 0);
W32(ep, fen_colc, 0);
W32(ep, fen_broc, 0);
W32(ep, fen_mulc, 0);
W32(ep, fen_uspc, 0);
W32(ep, fen_frgc, 0);
W32(ep, fen_ospc, 0);
W32(ep, fen_jbrc, 0);
W32(ep, fen_p64c, 0);
W32(ep, fen_p65c, 0);
W32(ep, fen_p128c, 0);
W32(ep, fen_p256c, 0);
W32(ep, fen_p512c, 0);
W32(ep, fen_p1024c, 0);
W16(ep, fen_rfthr, 0); /* Suggested by manual */
W16(ep, fen_rfcnt, 0);
W16(ep, fen_cftype, 0);
fs_init_bds(dev);
/* adjust to speed (for RMII mode) */
if (fpi->use_rmii) {
if (fep->phydev->speed == 100)
C8(fcccp, fcc_gfemr, 0x20);
else
S8(fcccp, fcc_gfemr, 0x20);
}
fcc_cr_cmd(fep, CPM_CR_INIT_TRX);
/* clear events */
W16(fccp, fcc_fcce, 0xffff);
/* Enable interrupts we wish to service */
W16(fccp, fcc_fccm, FCC_ENET_TXE | FCC_ENET_RXF | FCC_ENET_TXB);
/* Set GFMR to enable Ethernet operating mode */
W32(fccp, fcc_gfmr, FCC_GFMR_TCI | FCC_GFMR_MODE_ENET);
/* set sync/delimiters */
W16(fccp, fcc_fdsr, 0xd555);
W32(fccp, fcc_fpsmr, FCC_PSMR_ENCRC);
if (fpi->use_rmii)
S32(fccp, fcc_fpsmr, FCC_PSMR_RMII);
/* adjust to duplex mode */
if (fep->phydev->duplex)
S32(fccp, fcc_fpsmr, FCC_PSMR_FDE | FCC_PSMR_LPB);
else
C32(fccp, fcc_fpsmr, FCC_PSMR_FDE | FCC_PSMR_LPB);
/* Restore multicast and promiscuous settings */
set_multicast_list(dev);
S32(fccp, fcc_gfmr, FCC_GFMR_ENR | FCC_GFMR_ENT);
}
void test(int M, int N)
{
/* Scattering iterators. */
int c2, c4, c6;
/* Original iterators. */
int i, j, k;
if ((M <= -1) && (N >= 1)) {
S1(0) ;
}
if ((M >= 0) && (N >= 1)) {
S1(0) ;
}
if ((M >= 1) && (N >= 2)) {
for (c4=0;c4<=M-1;c4++) {
S2(0,c4) ;
}
S3(0) ;
for (c4=0;c4<=M-1;c4++) {
S4(0,c4) ;
}
S10(0) ;
S1(1) ;
S5(0) ;
}
if ((M <= 0) && (N >= 2)) {
S3(0) ;
S10(0) ;
S1(1) ;
S5(0) ;
}
if ((M >= 1) && (N == 1)) {
for (c4=0;c4<=M-1;c4++) {
S2(0,c4) ;
}
S3(0) ;
for (c4=0;c4<=M-1;c4++) {
S4(0,c4) ;
}
S10(0) ;
S5(0) ;
}
if ((M <= 0) && (N == 1)) {
S3(0) ;
S10(0) ;
S5(0) ;
}
for (c2=2;c2<=min(N-1,M);c2++) {
for (c4=c2-1;c4<=N-1;c4++) {
i = c2-2 ;
S6(c2-2,c4) ;
for (c6=c2-2;c6<=M-1;c6++) {
i = c2-2 ;
S7(c2-2,c4,c6) ;
}
i = c2-2 ;
S8(c2-2,c4) ;
for (c6=c2-2;c6<=M-1;c6++) {
i = c2-2 ;
S9(c2-2,c4,c6) ;
}
}
for (c4=c2-1;c4<=M-1;c4++) {
i = c2-1 ;
S2(c2-1,c4) ;
}
i = c2-1 ;
S3(c2-1) ;
for (c4=c2-1;c4<=M-1;c4++) {
i = c2-1 ;
S4(c2-1,c4) ;
}
i = c2-1 ;
S10(c2-1) ;
S1(c2) ;
i = c2-1 ;
S5(c2-1) ;
}
if ((M >= 1) && (M <= N-2)) {
c2 = M+1 ;
for (c4=M;c4<=N-1;c4++) {
i = M-1 ;
S6(M-1,c4) ;
c6 = M-1 ;
i = M-1 ;
k = M-1 ;
S7(M-1,c4,M-1) ;
i = M-1 ;
S8(M-1,c4) ;
c6 = M-1 ;
i = M-1 ;
k = M-1 ;
S9(M-1,c4,M-1) ;
}
S3(M) ;
S10(M) ;
i = M+1 ;
S1(M+1) ;
S5(M) ;
}
if ((M >= N) && (N >= 2)) {
c4 = N-1 ;
i = N-2 ;
j = N-1 ;
S6(N-2,N-1) ;
for (c6=N-2;c6<=M-1;c6++) {
i = N-2 ;
j = N-1 ;
S7(N-2,N-1,c6) ;
}
i = N-2 ;
j = N-1 ;
S8(N-2,N-1) ;
for (c6=N-2;c6<=M-1;c6++) {
i = N-2 ;
j = N-1 ;
S9(N-2,N-1,c6) ;
}
for (c4=N-1;c4<=M-1;c4++) {
i = N-1 ;
S2(N-1,c4) ;
}
i = N-1 ;
S3(N-1) ;
for (c4=N-1;c4<=M-1;c4++) {
i = N-1 ;
S4(N-1,c4) ;
}
i = N-1 ;
S10(N-1) ;
i = N-1 ;
S5(N-1) ;
}
if ((M == N-1) && (N >= 2)) {
c4 = N-1 ;
i = N-2 ;
j = N-1 ;
S6(N-2,N-1) ;
c6 = N-2 ;
i = N-2 ;
j = N-1 ;
k = N-2 ;
S7(N-2,N-1,N-2) ;
i = N-2 ;
j = N-1 ;
S8(N-2,N-1) ;
c6 = N-2 ;
i = N-2 ;
j = N-1 ;
k = N-2 ;
S9(N-2,N-1,N-2) ;
i = N-1 ;
S3(N-1) ;
i = N-1 ;
S10(N-1) ;
i = N-1 ;
S5(N-1) ;
}
for (c2=max(M+2,2);c2<=N-1;c2++) {
for (c4=c2-1;c4<=N-1;c4++) {
i = c2-2 ;
S6(c2-2,c4) ;
i = c2-2 ;
S8(c2-2,c4) ;
}
i = c2-1 ;
S3(c2-1) ;
i = c2-1 ;
S10(c2-1) ;
S1(c2) ;
i = c2-1 ;
S5(c2-1) ;
}
if ((M <= N-2) && (N >= 2)) {
c4 = N-1 ;
i = N-2 ;
j = N-1 ;
S6(N-2,N-1) ;
i = N-2 ;
j = N-1 ;
S8(N-2,N-1) ;
i = N-1 ;
S3(N-1) ;
i = N-1 ;
S10(N-1) ;
i = N-1 ;
S5(N-1) ;
}
}
///
/// Generates source code for the bishop rays.
///
void CFiestyGen::genBishopRays()
{
std::cout << "const SBishopRays CGen::mbbBishopRays[CSqix::kNumSquares] = {";
for ( U8 sq = 0; sq < CSqix::kNumSquares; sq++ )
{
CSqix sqix( sq );
U8 bishopRank = U8( sqix.getRank().get() );
U8 bishopFile = U8( sqix.getFile().get() );
CBitBoard bbNorthEast( 0ULL );
CBitBoard bbSouthEast( 0ULL );
CBitBoard bbSouthWest( 0ULL );
CBitBoard bbNorthWest( 0ULL );
//
// Loop setting the north-east bits
//
S8 r = bishopRank + 1;
S8 f = bishopFile + 1;
while ( r <= S8( ERank::kRank8 ) && f <= S8( EFile::kFileH ) )
{
bbNorthEast.setSquare( CSqix( ERank( r ), EFile( f ) ).get() );
r++;
f++;
}
//
// Loop setting the south-east bits
//
r = bishopRank - 1;
f = bishopFile + 1;
while ( r >= S8( ERank::kRank1 ) && f <= S8( EFile::kFileH ) )
{
bbSouthEast.setSquare( CSqix( ERank( r ), EFile( f ) ).get() );
r--;
f++;
}
//
// Loop setting the south-west bits
//
r = bishopRank - 1;
f = bishopFile - 1;
while ( r >= S8( ERank::kRank1 ) && f >= S8( EFile::kFileA ) )
{
bbSouthWest.setSquare( CSqix( ERank( r ), EFile( f ) ).get() );
r--;
f--;
}
//
// Loop setting the north-west bits
//
r = bishopRank + 1;
f = bishopFile - 1;
while ( r <= S8( ERank::kRank8 ) && f >= S8( EFile::kFileA ) )
{
bbNorthWest.setSquare( CSqix( ERank( r ), EFile( f ) ).get() );
r++;
f--;
}
#ifdef BBTRACE
printBitBoardDiagram( sqix.asAbbr() + " NorthEast", bbNorthEast );
printBitBoardDiagram( sqix.asAbbr() + " SouthEast", bbSouthEast );
printBitBoardDiagram( sqix.asAbbr() + " SouthWest", bbSouthWest );
printBitBoardDiagram( sqix.asAbbr() + " NorthWest", bbNorthWest );
#endif
std::cout << "\n /* " << sqix.asAbbr() << " */ { "
<< bbNorthEast.asAbbr() << "ULL" << ", "
<< bbSouthEast.asAbbr() << "ULL" << ", "
<< bbSouthWest.asAbbr() << "ULL" << ", "
<< bbNorthWest.asAbbr() << "ULL" << " }";
if ( sq != CSqix::kNumSquares - 1 )
std::cout << ", ";
}
std::cout << "};\n" << std::flush;
}
void test(int M, int N, int O, int P, int Q, int R)
{
/* Scattering iterators. */
int p1, p3, p5;
/* Original iterators. */
int i, j, k;
if (M == 1) {
S1() ;
S2() ;
S3() ;
S4() ;
S5() ;
S6() ;
S7() ;
S8() ;
S9() ;
S10() ;
S11() ;
S12() ;
S13() ;
S14() ;
S15() ;
S16() ;
S17() ;
S18() ;
S19() ;
S20() ;
S21() ;
S22() ;
S23() ;
S24() ;
S25() ;
S26() ;
S27() ;
}
if (M == 1) {
for (p1=1;p1<=N;p1++) {
for (p3=1;p3<=N;p3++) {
S28(p1,p3) ;
S29(p1,p3) ;
S30(p1,p3) ;
}
S31(p1) ;
}
}
if (M == 1) {
S32() ;
S33() ;
S34() ;
}
if ((M == 1) && (O <= 1)) {
S35() ;
}
if (M == 1) {
S36() ;
S37() ;
}
if ((M == 1) && (N >= 1) && (Q >= 1) && (R >= 1)) {
for (p1=2;p1<=P;p1++) {
S38(p1) ;
S39(p1) ;
for (p3=1;p3<=Q;p3++) {
for (p5=1;p5<=R;p5++) {
S40(p1,p3,p5) ;
S41(p1,p3,p5) ;
S42(p1,p3,p5) ;
S43(p1,p3,p5) ;
}
}
for (p3=1;p3<=Q;p3++) {
S44(p1,p3) ;
S45(p1,p3) ;
S46(p1,p3) ;
S47(p1,p3) ;
}
for (p3=1;p3<=R;p3++) {
S48(p1,p3) ;
S49(p1,p3) ;
S50(p1,p3) ;
S51(p1,p3) ;
}
S52(p1) ;
S53(p1) ;
S54(p1) ;
S55(p1) ;
S56(p1) ;
S57(p1) ;
S58(p1) ;
for (p3=1;p3<=Q;p3++) {
for (p5=1;p5<=R;p5++) {
S59(p1,p3,p5) ;
S60(p1,p3,p5) ;
S61(p1,p3,p5) ;
}
}
for (p3=1;p3<=Q;p3++) {
S62(p1,p3) ;
S63(p1,p3) ;
S64(p1,p3) ;
}
for (p3=1;p3<=R;p3++) {
S65(p1,p3) ;
S66(p1,p3) ;
S67(p1,p3) ;
}
S68(p1) ;
S69(p1) ;
S70(p1) ;
S71(p1) ;
S72(p1) ;
S73(p1) ;
S74(p1) ;
S75(p1) ;
S76(p1) ;
S77(p1) ;
S78(p1) ;
S79(p1) ;
S80(p1) ;
S81(p1) ;
S82(p1) ;
S83(p1) ;
S84(p1) ;
S85(p1) ;
S86(p1) ;
S87(p1) ;
S88(p1) ;
S89(p1) ;
S90(p1) ;
S91(p1) ;
S92(p1) ;
S93(p1) ;
S94(p1) ;
for (p3=1;p3<=N;p3++) {
for (p5=1;p5<=N;p5++) {
S95(p1,p3,p5) ;
S96(p1,p3,p5) ;
S97(p1,p3,p5) ;
}
S98(p1,p3) ;
}
S99(p1) ;
S100(p1) ;
S101(p1) ;
for (p3=1;p3<=Q;p3++) {
for (p5=1;p5<=R;p5++) {
S102(p1,p3,p5) ;
S103(p1,p3,p5) ;
S104(p1,p3,p5) ;
S105(p1,p3,p5) ;
S106(p1,p3,p5) ;
S107(p1,p3,p5) ;
}
}
for (p3=1;p3<=Q;p3++) {
S108(p1,p3) ;
S109(p1,p3) ;
S110(p1,p3) ;
S111(p1,p3) ;
S112(p1,p3) ;
S113(p1,p3) ;
}
for (p3=1;p3<=R;p3++) {
S114(p1,p3) ;
S115(p1,p3) ;
S116(p1,p3) ;
S117(p1,p3) ;
S118(p1,p3) ;
S119(p1,p3) ;
}
S120(p1) ;
S121(p1) ;
S122(p1) ;
S123(p1) ;
S124(p1) ;
S125(p1) ;
}
}
if ((M == 1) && (N <= 0) && (Q >= 1) && (R >= 1)) {
for (p1=2;p1<=P;p1++) {
S38(p1) ;
S39(p1) ;
for (p3=1;p3<=Q;p3++) {
for (p5=1;p5<=R;p5++) {
S40(p1,p3,p5) ;
S41(p1,p3,p5) ;
S42(p1,p3,p5) ;
S43(p1,p3,p5) ;
}
}
for (p3=1;p3<=Q;p3++) {
S44(p1,p3) ;
S45(p1,p3) ;
S46(p1,p3) ;
S47(p1,p3) ;
}
for (p3=1;p3<=R;p3++) {
S48(p1,p3) ;
S49(p1,p3) ;
S50(p1,p3) ;
S51(p1,p3) ;
}
S52(p1) ;
S53(p1) ;
S54(p1) ;
S55(p1) ;
S56(p1) ;
S57(p1) ;
S58(p1) ;
for (p3=1;p3<=Q;p3++) {
for (p5=1;p5<=R;p5++) {
S59(p1,p3,p5) ;
S60(p1,p3,p5) ;
S61(p1,p3,p5) ;
}
}
for (p3=1;p3<=Q;p3++) {
S62(p1,p3) ;
S63(p1,p3) ;
S64(p1,p3) ;
}
for (p3=1;p3<=R;p3++) {
S65(p1,p3) ;
S66(p1,p3) ;
S67(p1,p3) ;
}
S68(p1) ;
S69(p1) ;
S70(p1) ;
S71(p1) ;
S72(p1) ;
S73(p1) ;
S74(p1) ;
S75(p1) ;
S76(p1) ;
S77(p1) ;
S78(p1) ;
S79(p1) ;
S80(p1) ;
S81(p1) ;
S82(p1) ;
S83(p1) ;
S84(p1) ;
S85(p1) ;
S86(p1) ;
S87(p1) ;
S88(p1) ;
S89(p1) ;
S90(p1) ;
S91(p1) ;
S92(p1) ;
S93(p1) ;
S94(p1) ;
S99(p1) ;
S100(p1) ;
S101(p1) ;
for (p3=1;p3<=Q;p3++) {
for (p5=1;p5<=R;p5++) {
S102(p1,p3,p5) ;
S103(p1,p3,p5) ;
S104(p1,p3,p5) ;
S105(p1,p3,p5) ;
S106(p1,p3,p5) ;
S107(p1,p3,p5) ;
}
}
for (p3=1;p3<=Q;p3++) {
S108(p1,p3) ;
S109(p1,p3) ;
S110(p1,p3) ;
S111(p1,p3) ;
S112(p1,p3) ;
S113(p1,p3) ;
}
for (p3=1;p3<=R;p3++) {
S114(p1,p3) ;
S115(p1,p3) ;
S116(p1,p3) ;
S117(p1,p3) ;
S118(p1,p3) ;
S119(p1,p3) ;
}
S120(p1) ;
S121(p1) ;
S122(p1) ;
S123(p1) ;
S124(p1) ;
S125(p1) ;
}
}
if ((M == 1) && (N >= 1) && (Q <= 0) && (R >= 1)) {
for (p1=2;p1<=P;p1++) {
S38(p1) ;
S39(p1) ;
for (p3=1;p3<=R;p3++) {
S48(p1,p3) ;
S49(p1,p3) ;
S50(p1,p3) ;
S51(p1,p3) ;
}
S52(p1) ;
S53(p1) ;
S54(p1) ;
S55(p1) ;
S56(p1) ;
S57(p1) ;
S58(p1) ;
for (p3=1;p3<=R;p3++) {
S65(p1,p3) ;
S66(p1,p3) ;
S67(p1,p3) ;
}
S68(p1) ;
S69(p1) ;
S70(p1) ;
S71(p1) ;
S72(p1) ;
S73(p1) ;
S74(p1) ;
S75(p1) ;
S76(p1) ;
S77(p1) ;
S78(p1) ;
S79(p1) ;
S80(p1) ;
S81(p1) ;
S82(p1) ;
S83(p1) ;
S84(p1) ;
S85(p1) ;
S86(p1) ;
S87(p1) ;
S88(p1) ;
S89(p1) ;
S90(p1) ;
S91(p1) ;
S92(p1) ;
S93(p1) ;
S94(p1) ;
for (p3=1;p3<=N;p3++) {
for (p5=1;p5<=N;p5++) {
S95(p1,p3,p5) ;
S96(p1,p3,p5) ;
S97(p1,p3,p5) ;
}
S98(p1,p3) ;
}
S99(p1) ;
S100(p1) ;
S101(p1) ;
for (p3=1;p3<=R;p3++) {
S114(p1,p3) ;
S115(p1,p3) ;
S116(p1,p3) ;
S117(p1,p3) ;
S118(p1,p3) ;
S119(p1,p3) ;
}
S120(p1) ;
S121(p1) ;
S122(p1) ;
S123(p1) ;
S124(p1) ;
S125(p1) ;
}
}
if ((M == 1) && (N <= 0) && (Q <= 0) && (R >= 1)) {
for (p1=2;p1<=P;p1++) {
S38(p1) ;
S39(p1) ;
for (p3=1;p3<=R;p3++) {
S48(p1,p3) ;
S49(p1,p3) ;
S50(p1,p3) ;
S51(p1,p3) ;
}
S52(p1) ;
S53(p1) ;
S54(p1) ;
S55(p1) ;
S56(p1) ;
S57(p1) ;
S58(p1) ;
for (p3=1;p3<=R;p3++) {
S65(p1,p3) ;
S66(p1,p3) ;
S67(p1,p3) ;
}
S68(p1) ;
S69(p1) ;
S70(p1) ;
S71(p1) ;
S72(p1) ;
S73(p1) ;
S74(p1) ;
S75(p1) ;
S76(p1) ;
S77(p1) ;
S78(p1) ;
S79(p1) ;
S80(p1) ;
S81(p1) ;
S82(p1) ;
S83(p1) ;
S84(p1) ;
S85(p1) ;
S86(p1) ;
S87(p1) ;
S88(p1) ;
S89(p1) ;
S90(p1) ;
S91(p1) ;
S92(p1) ;
S93(p1) ;
S94(p1) ;
S99(p1) ;
S100(p1) ;
S101(p1) ;
for (p3=1;p3<=R;p3++) {
S114(p1,p3) ;
S115(p1,p3) ;
S116(p1,p3) ;
S117(p1,p3) ;
S118(p1,p3) ;
S119(p1,p3) ;
}
S120(p1) ;
S121(p1) ;
S122(p1) ;
S123(p1) ;
S124(p1) ;
S125(p1) ;
}
}
if ((M == 1) && (N >= 1) && (Q <= 0) && (R <= 0)) {
for (p1=2;p1<=P;p1++) {
S38(p1) ;
S39(p1) ;
S52(p1) ;
S53(p1) ;
S54(p1) ;
S55(p1) ;
S56(p1) ;
S57(p1) ;
S58(p1) ;
S68(p1) ;
S69(p1) ;
S70(p1) ;
S71(p1) ;
S72(p1) ;
S73(p1) ;
S74(p1) ;
S75(p1) ;
S76(p1) ;
S77(p1) ;
S78(p1) ;
S79(p1) ;
S80(p1) ;
S81(p1) ;
S82(p1) ;
S83(p1) ;
S84(p1) ;
S85(p1) ;
S86(p1) ;
S87(p1) ;
S88(p1) ;
S89(p1) ;
S90(p1) ;
S91(p1) ;
S92(p1) ;
S93(p1) ;
S94(p1) ;
for (p3=1;p3<=N;p3++) {
for (p5=1;p5<=N;p5++) {
S95(p1,p3,p5) ;
S96(p1,p3,p5) ;
S97(p1,p3,p5) ;
}
S98(p1,p3) ;
}
S99(p1) ;
S100(p1) ;
S101(p1) ;
S120(p1) ;
S121(p1) ;
S122(p1) ;
S123(p1) ;
S124(p1) ;
S125(p1) ;
}
}
if ((M == 1) && (N <= 0) && (Q <= 0) && (R <= 0)) {
for (p1=2;p1<=P;p1++) {
S38(p1) ;
S39(p1) ;
S52(p1) ;
S53(p1) ;
S54(p1) ;
S55(p1) ;
S56(p1) ;
S57(p1) ;
S58(p1) ;
S68(p1) ;
S69(p1) ;
S70(p1) ;
S71(p1) ;
S72(p1) ;
S73(p1) ;
S74(p1) ;
S75(p1) ;
S76(p1) ;
S77(p1) ;
S78(p1) ;
S79(p1) ;
S80(p1) ;
S81(p1) ;
S82(p1) ;
S83(p1) ;
S84(p1) ;
S85(p1) ;
S86(p1) ;
S87(p1) ;
S88(p1) ;
S89(p1) ;
S90(p1) ;
S91(p1) ;
S92(p1) ;
S93(p1) ;
S94(p1) ;
S99(p1) ;
S100(p1) ;
S101(p1) ;
S120(p1) ;
S121(p1) ;
S122(p1) ;
S123(p1) ;
S124(p1) ;
S125(p1) ;
}
}
if ((M == 1) && (N >= 1) && (Q >= 1) && (R <= 0)) {
for (p1=2;p1<=P;p1++) {
S38(p1) ;
S39(p1) ;
for (p3=1;p3<=Q;p3++) {
S44(p1,p3) ;
S45(p1,p3) ;
S46(p1,p3) ;
S47(p1,p3) ;
}
S52(p1) ;
S53(p1) ;
S54(p1) ;
S55(p1) ;
S56(p1) ;
S57(p1) ;
S58(p1) ;
for (p3=1;p3<=Q;p3++) {
S62(p1,p3) ;
S63(p1,p3) ;
S64(p1,p3) ;
}
S68(p1) ;
S69(p1) ;
S70(p1) ;
S71(p1) ;
S72(p1) ;
S73(p1) ;
S74(p1) ;
S75(p1) ;
S76(p1) ;
S77(p1) ;
S78(p1) ;
S79(p1) ;
S80(p1) ;
S81(p1) ;
S82(p1) ;
S83(p1) ;
S84(p1) ;
S85(p1) ;
S86(p1) ;
S87(p1) ;
S88(p1) ;
S89(p1) ;
S90(p1) ;
S91(p1) ;
S92(p1) ;
S93(p1) ;
S94(p1) ;
for (p3=1;p3<=N;p3++) {
for (p5=1;p5<=N;p5++) {
S95(p1,p3,p5) ;
S96(p1,p3,p5) ;
S97(p1,p3,p5) ;
}
S98(p1,p3) ;
}
S99(p1) ;
S100(p1) ;
S101(p1) ;
for (p3=1;p3<=Q;p3++) {
S108(p1,p3) ;
S109(p1,p3) ;
S110(p1,p3) ;
S111(p1,p3) ;
S112(p1,p3) ;
S113(p1,p3) ;
}
S120(p1) ;
S121(p1) ;
S122(p1) ;
S123(p1) ;
S124(p1) ;
S125(p1) ;
}
}
if ((M == 1) && (N <= 0) && (Q >= 1) && (R <= 0)) {
for (p1=2;p1<=P;p1++) {
S38(p1) ;
S39(p1) ;
for (p3=1;p3<=Q;p3++) {
S44(p1,p3) ;
S45(p1,p3) ;
S46(p1,p3) ;
S47(p1,p3) ;
}
S52(p1) ;
S53(p1) ;
S54(p1) ;
S55(p1) ;
S56(p1) ;
S57(p1) ;
S58(p1) ;
for (p3=1;p3<=Q;p3++) {
S62(p1,p3) ;
S63(p1,p3) ;
S64(p1,p3) ;
}
S68(p1) ;
S69(p1) ;
S70(p1) ;
S71(p1) ;
S72(p1) ;
S73(p1) ;
S74(p1) ;
S75(p1) ;
S76(p1) ;
S77(p1) ;
S78(p1) ;
S79(p1) ;
S80(p1) ;
S81(p1) ;
S82(p1) ;
S83(p1) ;
S84(p1) ;
S85(p1) ;
S86(p1) ;
S87(p1) ;
S88(p1) ;
S89(p1) ;
S90(p1) ;
S91(p1) ;
S92(p1) ;
S93(p1) ;
S94(p1) ;
S99(p1) ;
S100(p1) ;
S101(p1) ;
for (p3=1;p3<=Q;p3++) {
S108(p1,p3) ;
S109(p1,p3) ;
S110(p1,p3) ;
S111(p1,p3) ;
S112(p1,p3) ;
S113(p1,p3) ;
}
S120(p1) ;
S121(p1) ;
S122(p1) ;
S123(p1) ;
S124(p1) ;
S125(p1) ;
}
}
}
SolidesP ConstruireNormale (FacesP Face)
// Construit un solide symbolisant la normale à la face.
// Le solide est placé sur le barycentre de la face et peut donc
// être hors de la face.
// Sa longueur est égale à 5% de la dimension du solide à laquelle la
// face appartient s'il existe. Sinon, c'est 10 % de la dimension de la face.
{
SolidesP Sol = Face.Solide () ;
double HauteurTot ;
Englobants3D EnglobantFace = Face.Englobant () ;
if (Sol.EstInitialise ()) {
Englobants3D E = Sol.Englobant () ;
Vecteurs3D V (E.BasGauche (),E.HautDroit ()) ;
HauteurTot = 0.05*V.Norme () ; ;
}
else {
Vecteurs3D V (EnglobantFace.BasGauche (),EnglobantFace.HautDroit ()) ;
HauteurTot = 0.1*V.Norme () ; ;
}
double LargeurPied = HauteurTot / 24,
LargeurTete = HauteurTot / 4,
HauteurPied = (7./8.)*HauteurTot ;
// Construction d'un repère centré sur le "barycentre" de la face, dont
// le (xOy) correspond au plan de la face et dont le z est la normale à
// la face.
Points3D Origine = 0.5*(EnglobantFace.BasGauche ()+EnglobantFace.HautDroit ()) ;
Vecteurs3D K = Face.VecteurNormal ().VecteurNorme (),
U (Face.ContourExterieur ().IemeAreteOrientee (0).Origine ().Point3D (),
Face.ContourExterieur ().IemeAreteOrientee (0).Extremite ().Point3D ()),
I = (K & U).VecteurNorme (),
J = (K & I).VecteurNorme () ;
SommetsP S1 (Origine-LargeurPied/2*I - LargeurPied/2*J) ;
SommetsP S2 (Origine-LargeurPied/2*I + LargeurPied/2*J) ;
SommetsP S3 (Origine+LargeurPied/2*I + LargeurPied/2*J) ;
SommetsP S4 (Origine+LargeurPied/2*I - LargeurPied/2*J) ;
SommetsP S5 (Origine-LargeurPied/2*I - LargeurPied/2*J + HauteurPied*K) ;
SommetsP S6 (Origine-LargeurPied/2*I + LargeurPied/2*J + HauteurPied*K) ;
SommetsP S7 (Origine+LargeurPied/2*I + LargeurPied/2*J + HauteurPied*K) ;
SommetsP S8 (Origine+LargeurPied/2*I - LargeurPied/2*J + HauteurPied*K) ;
SommetsP S9 (Origine-LargeurTete/2*I - LargeurTete/2*J + HauteurPied*K) ;
SommetsP S10 (Origine-LargeurTete/2*I + LargeurTete/2*J + HauteurPied*K) ;
SommetsP S11 (Origine+LargeurTete/2*I + LargeurTete/2*J + HauteurPied*K) ;
SommetsP S12 (Origine+LargeurTete/2*I - LargeurTete/2*J + HauteurPied*K) ;
SommetsP S13 (Origine+HauteurTot*K) ;
AretesP A1 (S1,S2) ;
AretesP A2 (S2,S3) ;
AretesP A3 (S3,S4) ;
AretesP A4 (S4,S1) ;
AretesP A5 (S5,S6) ;
AretesP A6 (S6,S7) ;
AretesP A7 (S7,S8) ;
AretesP A8 (S8,S5) ;
AretesP A9 (S9,S10) ;
AretesP A10 (S10,S11) ;
AretesP A11 (S11,S12) ;
AretesP A12 (S12,S9) ;
AretesP A13 (S1,S5) ;
AretesP A14 (S2,S6) ;
AretesP A15 (S3,S7) ;
AretesP A16 (S4,S8) ;
AretesP A17 (S9,S13) ;
AretesP A18 (S10,S13) ;
AretesP A19 (S11,S13) ;
AretesP A20 (S12,S13) ;
Listes <AretesP> L1 ;
L1.InsertionEnQueue (A1) ;
L1.InsertionEnQueue (A2) ;
L1.InsertionEnQueue (A3) ;
L1.InsertionEnQueue (A4) ;
Listes <AretesP> L2 ;
L2.InsertionEnQueue (A1) ;
L2.InsertionEnQueue (A13) ;
L2.InsertionEnQueue (A5) ;
L2.InsertionEnQueue (A14) ;
Listes <AretesP> L3 ;
L3.InsertionEnQueue (A2) ;
L3.InsertionEnQueue (A14) ;
L3.InsertionEnQueue (A6) ;
L3.InsertionEnQueue (A15) ;
Listes <AretesP> L4 ;
L4.InsertionEnQueue (A3) ;
L4.InsertionEnQueue (A15) ;
L4.InsertionEnQueue (A7) ;
L4.InsertionEnQueue (A16) ;
Listes <AretesP> L5 ;
L5.InsertionEnQueue (A4) ;
L5.InsertionEnQueue (A16) ;
L5.InsertionEnQueue (A8) ;
L5.InsertionEnQueue (A13) ;
Listes <AretesP> L6 ;
L6.InsertionEnQueue (A9) ;
L6.InsertionEnQueue (A10) ;
L6.InsertionEnQueue (A11) ;
L6.InsertionEnQueue (A12) ;
Listes <AretesP> L7 ;
L7.InsertionEnQueue (A8) ;
L7.InsertionEnQueue (A7) ;
L7.InsertionEnQueue (A6) ;
L7.InsertionEnQueue (A5) ;
Listes <AretesP> L8 ;
L8.InsertionEnQueue (A9) ;
L8.InsertionEnQueue (A17) ;
L8.InsertionEnQueue (A18) ;
Listes <AretesP> L9 ;
L9.InsertionEnQueue (A10) ;
L9.InsertionEnQueue (A18) ;
L9.InsertionEnQueue (A19) ;
Listes <AretesP> L10 ;
L10.InsertionEnQueue (A11) ;
L10.InsertionEnQueue (A19) ;
L10.InsertionEnQueue (A20) ;
Listes <AretesP> L11 ;
L11.InsertionEnQueue (A12) ;
L11.InsertionEnQueue (A20) ;
L11.InsertionEnQueue (A17) ;
Listes <FacesP> ListeFaces ;
AttributsFaces *PAF = new AttributsFaces (PMateriauRouge) ;
FacesP F1 (ContoursP (L1),PAF) ;
ListeFaces.InsertionEnQueue (F1) ;
ListeFaces.InsertionEnQueue (FacesP (ContoursP (L2),PAF)) ;
ListeFaces.InsertionEnQueue (FacesP (ContoursP (L3),PAF)) ;
ListeFaces.InsertionEnQueue (FacesP (ContoursP (L4),PAF)) ;
ListeFaces.InsertionEnQueue (FacesP (ContoursP (L5),PAF)) ;
ContoursP C (L6) ;
FacesP F (C,PAF) ;
F.AjouterContourInterieur (ContoursP (L7)) ;
ListeFaces.InsertionEnQueue (F) ;
ListeFaces.InsertionEnQueue (FacesP (ContoursP (L8),PAF)) ;
ListeFaces.InsertionEnQueue (FacesP (ContoursP (L9),PAF)) ;
ListeFaces.InsertionEnQueue (FacesP (ContoursP (L10),PAF)) ;
ListeFaces.InsertionEnQueue (FacesP (ContoursP (L11),PAF)) ;
return SolidesP (ListeFaces) ;
}
void test(int m)
{
/* Scattering iterators. */
int glT1, rp1, local1, local2;
/* Original iterators. */
int coordT1, coordP1, other1, other2, other3, other4, other5, other6;
if (m >= 2) {
S1(0,1,1,1) ;
S2(0,1,1,1,1,1,2,1) ;
S3(0,1,1,2,1,1,1,2) ;
S4(0,1,2,2,1,1,2,2) ;
S8(0,1) ;
}
if (m == 1) {
S1(0,1,1,1) ;
S8(0,1) ;
}
if (m >= 3) {
S5(0,1,1,1,1,1,2,1) ;
S1(1,1,2,1) ;
S2(1,1,2,1,2,1,3,1) ;
S3(1,1,2,2,2,1,2,2) ;
S4(1,1,3,2,2,1,3,2) ;
S6(0,1,1,2,1,1,1,2) ;
S7(0,1,2,2,1,1,2,2) ;
S1(1,2,1,2) ;
S2(1,2,2,2,1,2,2,2) ;
S3(1,2,2,3,1,2,1,3) ;
S4(1,2,3,3,1,2,2,3) ;
for (coordP1=1;coordP1<=2;coordP1++) {
S8(1,coordP1) ;
}
}
for (glT1=2;glT1<=m-2;glT1++) {
coordT1 = glT1-1 ;
other5 = glT1+1 ;
S5(glT1-1,1,glT1,1,glT1,1,glT1+1,1) ;
other1 = glT1+1 ;
S1(glT1,1,glT1+1,1) ;
local1 = glT1+1 ;
other1 = glT1+1 ;
other3 = glT1+1 ;
other5 = glT1+2 ;
S2(glT1,1,glT1+1,1,glT1+1,1,glT1+2,1) ;
other3 = glT1+1 ;
other5 = glT1+1 ;
S3(glT1,1,glT1+1,2,glT1+1,1,glT1+1,2) ;
other1 = glT1+2 ;
other3 = glT1+1 ;
other5 = glT1+2 ;
S4(glT1,1,glT1+2,2,glT1+1,1,glT1+2,2) ;
for (rp1=2;rp1<=glT1;rp1++) {
local1 = glT1-rp1+1 ;
coordT1 = glT1-1 ;
other3 = glT1-rp1+1 ;
other5 = glT1-rp1+2 ;
S5(glT1-1,rp1,glT1,rp1,glT1-rp1+1,rp1,glT1-rp1+2,rp1) ;
local1 = glT1-rp1+2 ;
local2 = rp1-1 ;
coordT1 = glT1-1 ;
coordP1 = rp1-1 ;
other3 = glT1-rp1+2 ;
other4 = rp1-1 ;
other5 = glT1-rp1+2 ;
S6(glT1-1,rp1-1,glT1,rp1,glT1-rp1+2,rp1-1,glT1-rp1+2,rp1) ;
other1 = glT1+1 ;
other3 = glT1-rp1+2 ;
other4 = rp1-1 ;
other5 = glT1-rp1+3 ;
S7(glT1-1,rp1-1,glT1+1,rp1,glT1-rp1+2,rp1-1,glT1-rp1+3,rp1) ;
other1 = glT1-rp1+2 ;
S1(glT1,rp1,glT1-rp1+2,rp1) ;
local1 = glT1-rp1+2 ;
other1 = glT1+1 ;
other3 = glT1-rp1+2 ;
other5 = glT1-rp1+3 ;
S2(glT1,rp1,glT1+1,rp1,glT1-rp1+2,rp1,glT1-rp1+3,rp1) ;
other2 = rp1+1 ;
other3 = glT1-rp1+2 ;
other5 = glT1-rp1+2 ;
other6 = rp1+1 ;
S3(glT1,rp1,glT1+1,rp1+1,glT1-rp1+2,rp1,glT1-rp1+2,rp1+1) ;
other1 = glT1+2 ;
other2 = rp1+1 ;
other3 = glT1-rp1+2 ;
other5 = glT1-rp1+3 ;
other6 = rp1+1 ;
S4(glT1,rp1,glT1+2,rp1+1,glT1-rp1+2,rp1,glT1-rp1+3,rp1+1) ;
}
rp1 = glT1+1 ;
coordT1 = glT1-1 ;
other2 = glT1+1 ;
other6 = glT1+1 ;
S6(glT1-1,glT1,glT1,glT1+1,1,glT1,1,glT1+1) ;
other1 = glT1+1 ;
other2 = glT1+1 ;
other6 = glT1+1 ;
S7(glT1-1,glT1,glT1+1,glT1+1,1,glT1,2,glT1+1) ;
coordP1 = glT1+1 ;
other2 = glT1+1 ;
S1(glT1,glT1+1,1,glT1+1) ;
local2 = glT1+1 ;
coordP1 = glT1+1 ;
other1 = glT1+1 ;
other2 = glT1+1 ;
other4 = glT1+1 ;
other6 = glT1+1 ;
S2(glT1,glT1+1,glT1+1,glT1+1,1,glT1+1,2,glT1+1) ;
other2 = glT1+2 ;
other4 = glT1+1 ;
other6 = glT1+2 ;
S3(glT1,glT1+1,glT1+1,glT1+2,1,glT1+1,1,glT1+2) ;
other1 = glT1+2 ;
other2 = glT1+2 ;
other4 = glT1+1 ;
other6 = glT1+2 ;
S4(glT1,glT1+1,glT1+2,glT1+2,1,glT1+1,2,glT1+2) ;
for (coordP1=1;coordP1<=glT1+1;coordP1++) {
S8(glT1,coordP1) ;
}
}
if (m >= 3) {
glT1 = m-1 ;
local1 = m-1 ;
coordT1 = m-2 ;
other1 = m-1 ;
other3 = m-1 ;
S5(m-2,1,m-1,1,m-1,1,m,1) ;
coordT1 = m-1 ;
S1(m-1,1,m,1) ;
coordT1 = m-1 ;
S3(m-1,1,m,2,m,1,m,2) ;
for (rp1=2;rp1<=m-1;rp1++) {
local1 = -rp1+m ;
coordT1 = m-2 ;
other1 = m-1 ;
other3 = -rp1+m ;
other5 = -rp1+m+1 ;
S5(m-2,rp1,m-1,rp1,-rp1+m,rp1,-rp1+m+1,rp1) ;
local1 = -rp1+m+1 ;
local2 = rp1-1 ;
coordT1 = m-2 ;
coordP1 = rp1-1 ;
other1 = m-1 ;
other3 = -rp1+m+1 ;
other4 = rp1-1 ;
other5 = -rp1+m+1 ;
S6(m-2,rp1-1,m-1,rp1,-rp1+m+1,rp1-1,-rp1+m+1,rp1) ;
other3 = -rp1+m+1 ;
other4 = rp1-1 ;
other5 = -rp1+m+2 ;
S7(m-2,rp1-1,m,rp1,-rp1+m+1,rp1-1,-rp1+m+2,rp1) ;
coordT1 = m-1 ;
other1 = -rp1+m+1 ;
S1(m-1,rp1,-rp1+m+1,rp1) ;
local1 = -rp1+m+1 ;
coordT1 = m-1 ;
other3 = -rp1+m+1 ;
other5 = -rp1+m+2 ;
S2(m-1,rp1,m,rp1,-rp1+m+1,rp1,-rp1+m+2,rp1) ;
other2 = rp1+1 ;
other3 = -rp1+m+1 ;
other5 = -rp1+m+1 ;
other6 = rp1+1 ;
S3(m-1,rp1,m,rp1+1,-rp1+m+1,rp1,-rp1+m+1,rp1+1) ;
other1 = m+1 ;
other2 = rp1+1 ;
other3 = -rp1+m+1 ;
other5 = -rp1+m+2 ;
other6 = rp1+1 ;
S4(m-1,rp1,m+1,rp1+1,-rp1+m+1,rp1,-rp1+m+2,rp1+1) ;
}
local2 = m-1 ;
coordT1 = m-2 ;
coordP1 = m-1 ;
other1 = m-1 ;
other4 = m-1 ;
S6(m-2,m-1,m-1,m,1,m-1,1,m) ;
other4 = m-1 ;
S7(m-2,m-1,m,m,1,m-1,2,m) ;
coordT1 = m-1 ;
S1(m-1,m,1,m) ;
coordT1 = m-1 ;
S2(m-1,m,m,m,1,m,2,m) ;
coordT1 = m-1 ;
for (coordP1=1;coordP1<=m;coordP1++) {
S8(m-1,coordP1) ;
}
}
for (glT1=m;glT1<=2*m-4;glT1++) {
rp1 = glT1-m+2 ;
local1 = m-1 ;
local2 = glT1-m+2 ;
coordT1 = glT1-1 ;
coordP1 = glT1-m+2 ;
other2 = glT1-m+2 ;
other3 = m-1 ;
other4 = glT1-m+2 ;
other6 = glT1-m+2 ;
S5(glT1-1,glT1-m+2,glT1,glT1-m+2,m-1,glT1-m+2,m,glT1-m+2) ;
local2 = glT1-m+1 ;
coordT1 = glT1-1 ;
coordP1 = glT1-m+1 ;
other2 = glT1-m+2 ;
other4 = glT1-m+1 ;
other6 = glT1-m+2 ;
S6(glT1-1,glT1-m+1,glT1,glT1-m+2,m,glT1-m+1,m,glT1-m+2) ;
coordP1 = glT1-m+2 ;
other2 = glT1-m+2 ;
S1(glT1,glT1-m+2,m,glT1-m+2) ;
local2 = glT1-m+2 ;
coordP1 = glT1-m+2 ;
other1 = glT1+1 ;
other2 = glT1-m+3 ;
other4 = glT1-m+2 ;
other6 = glT1-m+3 ;
S3(glT1,glT1-m+2,glT1+1,glT1-m+3,m,glT1-m+2,m,glT1-m+3) ;
for (rp1=glT1-m+3;rp1<=m-1;rp1++) {
local1 = glT1-rp1+1 ;
coordT1 = glT1-1 ;
other3 = glT1-rp1+1 ;
other5 = glT1-rp1+2 ;
S5(glT1-1,rp1,glT1,rp1,glT1-rp1+1,rp1,glT1-rp1+2,rp1) ;
local1 = glT1-rp1+2 ;
local2 = rp1-1 ;
coordT1 = glT1-1 ;
coordP1 = rp1-1 ;
other3 = glT1-rp1+2 ;
other4 = rp1-1 ;
other5 = glT1-rp1+2 ;
S6(glT1-1,rp1-1,glT1,rp1,glT1-rp1+2,rp1-1,glT1-rp1+2,rp1) ;
other1 = glT1+1 ;
other3 = glT1-rp1+2 ;
other4 = rp1-1 ;
other5 = glT1-rp1+3 ;
S7(glT1-1,rp1-1,glT1+1,rp1,glT1-rp1+2,rp1-1,glT1-rp1+3,rp1) ;
other1 = glT1-rp1+2 ;
S1(glT1,rp1,glT1-rp1+2,rp1) ;
local1 = glT1-rp1+2 ;
other1 = glT1+1 ;
other3 = glT1-rp1+2 ;
other5 = glT1-rp1+3 ;
S2(glT1,rp1,glT1+1,rp1,glT1-rp1+2,rp1,glT1-rp1+3,rp1) ;
other2 = rp1+1 ;
other3 = glT1-rp1+2 ;
other5 = glT1-rp1+2 ;
other6 = rp1+1 ;
S3(glT1,rp1,glT1+1,rp1+1,glT1-rp1+2,rp1,glT1-rp1+2,rp1+1) ;
other1 = glT1+2 ;
other2 = rp1+1 ;
other3 = glT1-rp1+2 ;
other5 = glT1-rp1+3 ;
other6 = rp1+1 ;
S4(glT1,rp1,glT1+2,rp1+1,glT1-rp1+2,rp1,glT1-rp1+3,rp1+1) ;
}
local1 = glT1-m+1 ;
coordT1 = glT1-1 ;
other3 = glT1-m+1 ;
other5 = glT1-m+2 ;
S5(glT1-1,m,glT1,m,glT1-m+1,m,glT1-m+2,m) ;
local1 = glT1-m+2 ;
local2 = m-1 ;
coordT1 = glT1-1 ;
coordP1 = m-1 ;
other3 = glT1-m+2 ;
other4 = m-1 ;
other5 = glT1-m+2 ;
S6(glT1-1,m-1,glT1,m,glT1-m+2,m-1,glT1-m+2,m) ;
other1 = glT1+1 ;
other3 = glT1-m+2 ;
other4 = m-1 ;
other5 = glT1-m+3 ;
S7(glT1-1,m-1,glT1+1,m,glT1-m+2,m-1,glT1-m+3,m) ;
other1 = glT1-m+2 ;
S1(glT1,m,glT1-m+2,m) ;
local1 = glT1-m+2 ;
other1 = glT1+1 ;
other3 = glT1-m+2 ;
other5 = glT1-m+3 ;
S2(glT1,m,glT1+1,m,glT1-m+2,m,glT1-m+3,m) ;
for (coordP1=glT1-m+2;coordP1<=m;coordP1++) {
S8(glT1,coordP1) ;
}
}
if (m >= 3) {
glT1 = 2*m-3 ;
rp1 = m-1 ;
local1 = m-1 ;
local2 = m-1 ;
coordT1 = 2*m-4 ;
coordP1 = m-1 ;
other1 = 2*m-3 ;
other2 = m-1 ;
other3 = m-1 ;
other4 = m-1 ;
other6 = m-1 ;
S5(2*m-4,m-1,2*m-3,m-1,m-1,m-1,m,m-1) ;
local2 = m-2 ;
coordT1 = 2*m-4 ;
coordP1 = m-2 ;
other1 = 2*m-3 ;
other2 = m-1 ;
other4 = m-2 ;
other6 = m-1 ;
S6(2*m-4,m-2,2*m-3,m-1,m,m-2,m,m-1) ;
coordT1 = 2*m-3 ;
coordP1 = m-1 ;
other2 = m-1 ;
S1(2*m-3,m-1,m,m-1) ;
local2 = m-1 ;
coordT1 = 2*m-3 ;
coordP1 = m-1 ;
other1 = 2*m-2 ;
other4 = m-1 ;
S3(2*m-3,m-1,2*m-2,m,m,m-1,m,m) ;
local1 = m-2 ;
coordT1 = 2*m-4 ;
other1 = 2*m-3 ;
other3 = m-2 ;
other5 = m-1 ;
S5(2*m-4,m,2*m-3,m,m-2,m,m-1,m) ;
local1 = m-1 ;
local2 = m-1 ;
coordT1 = 2*m-4 ;
coordP1 = m-1 ;
other1 = 2*m-3 ;
other3 = m-1 ;
other4 = m-1 ;
other5 = m-1 ;
S6(2*m-4,m-1,2*m-3,m,m-1,m-1,m-1,m) ;
other1 = 2*m-2 ;
other3 = m-1 ;
other4 = m-1 ;
S7(2*m-4,m-1,2*m-2,m,m-1,m-1,m,m) ;
coordT1 = 2*m-3 ;
other1 = m-1 ;
S1(2*m-3,m,m-1,m) ;
local1 = m-1 ;
coordT1 = 2*m-3 ;
other1 = 2*m-2 ;
other3 = m-1 ;
S2(2*m-3,m,2*m-2,m,m-1,m,m,m) ;
coordT1 = 2*m-3 ;
for (coordP1=m-1;coordP1<=m;coordP1++) {
S8(2*m-3,coordP1) ;
}
}
if (m == 2) {
S5(0,1,1,1,1,1,2,1) ;
S1(1,1,2,1) ;
S3(1,1,2,2,2,1,2,2) ;
S6(0,1,1,2,1,1,1,2) ;
S7(0,1,2,2,1,1,2,2) ;
S1(1,2,1,2) ;
S2(1,2,2,2,1,2,2,2) ;
for (coordP1=1;coordP1<=2;coordP1++) {
S8(1,coordP1) ;
}
}
if (m >= 2) {
glT1 = 2*m-2 ;
local1 = m-1 ;
coordT1 = 2*m-3 ;
other1 = 2*m-2 ;
other3 = m-1 ;
S5(2*m-3,m,2*m-2,m,m-1,m,m,m) ;
local2 = m-1 ;
coordT1 = 2*m-3 ;
coordP1 = m-1 ;
other1 = 2*m-2 ;
other4 = m-1 ;
S6(2*m-3,m-1,2*m-2,m,m,m-1,m,m) ;
coordT1 = 2*m-2 ;
S1(2*m-2,m,m,m) ;
coordT1 = 2*m-2 ;
S8(2*m-2,m) ;
}
}
void test(int M, int N, int O, int P, int Q, int R, int S, int T, int U)
{
/* Scattering iterators. */
int c2, c4, c6;
/* Original iterators. */
int i, j, k;
if ((M >= 2) && (N >= 4)) {
for (c2=1;c2<=O-1;c2++) {
for (c6=1;c6<=M;c6++) {
S1(c2,1,c6) ;
}
for (c6=1;c6<=M-1;c6++) {
S6(c2,1,c6) ;
S7(c2,1,c6) ;
}
for (c6=1;c6<=M;c6++) {
S3(c2,1,c6) ;
}
for (c6=1;c6<=M-1;c6++) {
S1(c2,2,c6) ;
}
S1(c2,2,M) ;
for (c6=1;c6<=M-1;c6++) {
S6(c2,2,c6) ;
S7(c2,2,c6) ;
}
for (c6=1;c6<=M-1;c6++) {
S11(c2,1,c6) ;
}
for (c4=3;c4<=2*N-5;c4++) {
for (c6=1;c6<=M-1;c6++) {
if ((c4+1)%2 == 0) {
j = (c4-1)/2 ;
S10(c2,(c4-1)/2,c6) ;
}
}
for (c6=1;c6<=M;c6++) {
if ((c4+1)%2 == 0) {
j = (c4+1)/2 ;
S3(c2,(c4+1)/2,c6) ;
}
}
for (c6=1;c6<=M-1;c6++) {
if (c4%2 == 0) {
j = (c4+2)/2 ;
S6(c2,(c4+2)/2,c6) ;
S7(c2,(c4+2)/2,c6) ;
}
if ((c4+1)%2 == 0) {
j = (c4+3)/2 ;
S1(c2,(c4+3)/2,c6) ;
}
}
if ((c4+1)%2 == 0) {
j = (c4+3)/2 ;
S1(c2,(c4+3)/2,M) ;
}
for (c6=1;c6<=M-1;c6++) {
if (c4%2 == 0) {
S11(c2,c4/2,c6) ;
}
}
}
c4 = 2*N-4 ;
for (c6=1;c6<=M-1;c6++) {
j = N-1 ;
S6(c2,N-1,c6) ;
S7(c2,N-1,c6) ;
}
for (c6=1;c6<=M-1;c6++) {
j = N-2 ;
S11(c2,N-2,c6) ;
}
c4 = 2*N-3 ;
for (c6=1;c6<=M-1;c6++) {
j = N-2 ;
S10(c2,N-2,c6) ;
}
for (c6=1;c6<=M;c6++) {
j = N-1 ;
S3(c2,N-1,c6) ;
}
c4 = 2*N-2 ;
for (c6=1;c6<=M-1;c6++) {
j = N-1 ;
S11(c2,N-1,c6) ;
}
c4 = 2*N-1 ;
for (c6=1;c6<=M-1;c6++) {
j = N-1 ;
S10(c2,N-1,c6) ;
}
}
}
if ((M >= 2) && (N == 3)) {
for (c2=1;c2<=O-1;c2++) {
for (c6=1;c6<=M;c6++) {
S1(c2,1,c6) ;
}
for (c6=1;c6<=M-1;c6++) {
S6(c2,1,c6) ;
S7(c2,1,c6) ;
}
for (c6=1;c6<=M;c6++) {
S3(c2,1,c6) ;
}
for (c6=1;c6<=M-1;c6++) {
S1(c2,2,c6) ;
}
S1(c2,2,M) ;
for (c6=1;c6<=M-1;c6++) {
S6(c2,2,c6) ;
S7(c2,2,c6) ;
}
for (c6=1;c6<=M-1;c6++) {
S11(c2,1,c6) ;
}
for (c6=1;c6<=M-1;c6++) {
S10(c2,1,c6) ;
}
for (c6=1;c6<=M;c6++) {
S3(c2,2,c6) ;
}
for (c6=1;c6<=M-1;c6++) {
S11(c2,2,c6) ;
}
for (c6=1;c6<=M-1;c6++) {
S10(c2,2,c6) ;
}
}
}
if ((M >= 2) && (N == 2)) {
for (c2=1;c2<=O-1;c2++) {
for (c6=1;c6<=M;c6++) {
S1(c2,1,c6) ;
}
for (c6=1;c6<=M-1;c6++) {
S6(c2,1,c6) ;
S7(c2,1,c6) ;
}
for (c6=1;c6<=M;c6++) {
S3(c2,1,c6) ;
}
for (c6=1;c6<=M-1;c6++) {
S11(c2,1,c6) ;
}
for (c6=1;c6<=M-1;c6++) {
S10(c2,1,c6) ;
}
}
}
if ((M == 1) && (N >= 3)) {
for (c2=1;c2<=O-1;c2++) {
for (c4=-1;c4<=0;c4++) {
if ((c4+1)%2 == 0) {
j = (c4+3)/2 ;
S1(c2,(c4+3)/2,1) ;
}
}
for (c4=1;c4<=2*N-5;c4++) {
if ((c4+1)%2 == 0) {
j = (c4+1)/2 ;
S3(c2,(c4+1)/2,1) ;
}
if ((c4+1)%2 == 0) {
j = (c4+3)/2 ;
S1(c2,(c4+3)/2,1) ;
}
}
for (c4=2*N-4;c4<=2*N-3;c4++) {
if ((c4+1)%2 == 0) {
j = (c4+1)/2 ;
S3(c2,(c4+1)/2,1) ;
}
}
}
}
if ((M == 1) && (N == 2)) {
for (c2=1;c2<=O-1;c2++) {
S1(c2,1,1) ;
S3(c2,1,1) ;
}
}
if ((M >= 2) && (N >= 3)) {
for (c2=1;c2<=O-1;c2++) {
for (c6=1;c6<=M;c6++) {
S2(c2,1,c6) ;
}
for (c6=1;c6<=M-1;c6++) {
S8(c2,1,c6) ;
}
for (c4=3;c4<=2*N-2;c4++) {
for (c6=1;c6<=M;c6++) {
if (c4%2 == 0) {
S2(c2,c4/2,c6) ;
}
}
for (c6=1;c6<=M-1;c6++) {
if (c4%2 == 0) {
S8(c2,c4/2,c6) ;
}
}
for (c6=1;c6<=M-1;c6++) {
if ((c4+1)%2 == 0) {
j = (c4-1)/2 ;
S9(c2,(c4-1)/2,c6) ;
}
}
}
c4 = 2*N-1 ;
for (c6=1;c6<=M-1;c6++) {
j = N-1 ;
S9(c2,N-1,c6) ;
}
}
}
if ((M >= 2) && (N == 2)) {
for (c2=1;c2<=O-1;c2++) {
for (c6=1;c6<=M;c6++) {
S2(c2,1,c6) ;
}
for (c6=1;c6<=M-1;c6++) {
S8(c2,1,c6) ;
}
for (c6=1;c6<=M-1;c6++) {
S9(c2,1,c6) ;
}
}
}
if ((M == 1) && (N >= 2)) {
for (c2=1;c2<=O-1;c2++) {
for (c4=2;c4<=2*N-2;c4++) {
if (c4%2 == 0) {
S2(c2,c4/2,1) ;
}
}
}
}
if ((M >= 2) && (N >= 2)) {
for (c2=1;c2<=O-1;c2++) {
for (c4=1;c4<=N-1;c4++) {
for (c6=1;c6<=M-1;c6++) {
S4(c2,c4,c6) ;
}
}
}
}
if ((M >= 2) && (N >= 2)) {
for (c2=1;c2<=O-1;c2++) {
for (c4=1;c4<=N-1;c4++) {
for (c6=1;c6<=M-1;c6++) {
S5(c2,c4,c6) ;
}
}
}
}
if ((M >= P+1) && (N >= Q+1)) {
for (c2=R;c2<=O-1;c2++) {
for (c4=Q;c4<=N-1;c4++) {
for (c6=P;c6<=M-1;c6++) {
S12(c2,c4,c6) ;
}
}
}
}
if ((M >= 2) && (N >= Q+1)) {
for (c2=R;c2<=O-1;c2++) {
for (c4=Q;c4<=N-1;c4++) {
for (c6=1;c6<=M-1;c6++) {
S13(c2,c4,c6) ;
}
}
}
}
if ((M >= P+1) && (N >= 2)) {
for (c2=R;c2<=O-1;c2++) {
for (c4=1;c4<=N-1;c4++) {
for (c6=P;c6<=M-1;c6++) {
S14(c2,c4,c6) ;
}
}
}
}
if ((M >= 2) && (N >= 2)) {
for (c2=R;c2<=O-1;c2++) {
for (c4=1;c4<=N-1;c4++) {
for (c6=1;c6<=M-1;c6++) {
S15(c2,c4,c6) ;
}
}
}
}
}
