Структура данных для ОЦК и ГЦК решеток

// ...
typedef vector<int> vec_1d;
typedef vector<vec_1d> vec_2d;
typedef vector<vec_2d> vec_3d;
/* ...
Вводим 3 состояния:
0 - ячейка доступна, но пуста
1 - ячейка доступна и заполнена атомом
2 - ячейка недоступна
... */
vec_3d frames(z, vec_2d(x, vec_1d(y, 2))); // делаем все ячейки недоступными
 
/* Открываем доступ ячейкам в соответствии с выбранной структурой решетки,
параллельно заполняя подложку (первые 3 слоя) атомами */
for( int i = 0; i < z; i++ )
{
    for( int j = 0; j < x; j++ )
    {
        for( int k = 0; k < y; k++ )
        {
            switch (structType)
            {
            case 0: // ОЦК-структура
                if ( ( i % 2 == 0 && j % 2 == 0 && k % 2 == 0 ) ||
                     ( i % 2 != 0 && j % 2 != 0 && k % 2 != 0 ) )
                    if (i < 3)
                        frames[i][j][k] = 1;
                    else
                        frames[i][j][k] = 0;
                break;
            case 1: // ГЦК-структура
                if ( ( i % 2 == 0 && j % 2 == 0 && k % 2 == 0 ) ||
                     ( i % 2 == 0 && j % 2 != 0 && k % 2 != 0 ) ||
                     ( i % 2 != 0 && j % 2 == 0 && k % 2 != 0 ) ||
                     ( i % 2 != 0 && j % 2 != 0 && k % 2 == 0 ) )
                    if (i < 3)
                            frames[i][j][k] = 1;
                    else
                        frames[i][j][k] = 0;
                break;
            default: // простая структура
                if (i < 3)
                    frames[i][j][k] = 1;
                else
                    frames[i][j][k] = 0;
                break;
            }
        }
    }
}
// ...