Magic Cube

///////////////////////////////////////////////////////////////////////////////////////////////////

// Copyright 2020 Leszek Koltunski                                                               //

//                                                                                               //

// This file is part of Magic Cube.                                                              //

//                                                                                               //

// Magic Cube is free software: you can redistribute it and/or modify                            //

// it under the terms of the GNU General Public License as published by                          //

// the Free Software Foundation, either version 2 of the License, or                             //

// (at your option) any later version.                                                           //

//                                                                                               //

// Magic Cube is distributed in the hope that it will be useful,                                 //

// but WITHOUT ANY WARRANTY; without even the implied warranty of                                //

// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the                                 //

// GNU General Public License for more details.                                                  //

//                                                                                               //

// You should have received a copy of the GNU General Public License                             //

// along with Magic Cube.  If not, see <http://www.gnu.org/licenses/>.                           //

///////////////////////////////////////////////////////////////////////////////////////////////////

package org.distorted.objects;

import android.content.res.Resources;

import org.distorted.helpers.FactoryCubit;

import org.distorted.helpers.ObjectShape;

import org.distorted.helpers.ObjectSticker;

import org.distorted.helpers.QuatHelper;

import org.distorted.library.effect.MatrixEffectQuaternion;

import org.distorted.library.main.DistortedEffects;

import org.distorted.library.main.DistortedTexture;

import org.distorted.library.mesh.MeshBase;

import org.distorted.library.mesh.MeshSquare;

import org.distorted.library.type.Static3D;

import org.distorted.library.type.Static4D;

import org.distorted.main.R;

///////////////////////////////////////////////////////////////////////////////////////////////////

public class TwistyKilominx extends TwistyMinx

{

  private static final int[] mCenterFaceMap = new int[]

      {

        0,0,0,0,1,

        1,0,1,1,0,

        2,0,1,1,0,

        2,2,1,0,2,

        2,1,0,0,1,

        1,2,0,1,0,

        0,1,0,1,1,

        0,1,0,2,0,

        2,1,2,2,2,

        1,0,2,1,2,

        2,1,0,1,2,

        2,2,2,2,2

      };

  private static final float CENTER_CORR = 0.87f;

  private static final float[][] STICKERS = new float[][]

      {

        { -0.36616942f, -0.36327124f, 0.5f, -0.36327124f, 0.23233888f, 0.4605048f, -0.36616942f, 0.26603764f },

        { -0.36327127f, -0.5f, 0.36327127f, -0.26393202f, 0.36327127f, 0.5f, -0.36327127f, 0.26393202f },

        { -0.3249197f, -0.39442718f, 0.3249197f, -0.39442718f, 0.3249197f, 0.5f, -0.3249197f, 0.2888544f }

      };

  static

    {

    final float C = 1.14f; // make the 'center' sticker artificially larger, so that we paint

                           // over the area in the center of the face.

    STICKERS[0][0] *= C;

    STICKERS[0][1] *= C;

    STICKERS[0][2] *= C;

    STICKERS[0][3] *= C;

    STICKERS[0][4] *= C;

    STICKERS[0][5] *= C;

    STICKERS[0][6] *= C;

    STICKERS[0][7] *= C;

    STICKERS[0][2] *= CENTER_CORR;

    STICKERS[0][3] *= CENTER_CORR;

    }

  private static final ObjectSticker[] mStickers;

  static

    {

    mStickers = new ObjectSticker[STICKERS.length];

    float R = 0.10f;

    final float[][] radii = { {R,R,R,R},{R,R,R,R},{R,R,R,R} };

    final float[] strokes = { 0.20f, 0.11f, 0.10f };

    for(int s=0; s<STICKERS.length; s++)

      {

      mStickers[s] = new ObjectSticker(STICKERS[s],null,radii[s],strokes[s]);

      }

    }

///////////////////////////////////////////////////////////////////////////////////////////////////

  TwistyKilominx(int size, Static4D quat, DistortedTexture texture, MeshSquare mesh,

                 DistortedEffects effects, int[][] moves, Resources res, int scrWidth)

    {

    super(size, size, quat, texture, mesh, effects, moves, ObjectList.KILO, res, scrWidth);

    }

///////////////////////////////////////////////////////////////////////////////////////////////////

  private int numCubitsPerCorner(int numLayers)

    {

    return 3*((numLayers-3)/2)*((numLayers-5)/2) + (numLayers<5 ? 0:1);

    }

///////////////////////////////////////////////////////////////////////////////////////////////////

  private int numCubitsPerEdge(int numLayers)

    {

    return numLayers<5 ? 0 : 2*(numLayers-4);

    }

///////////////////////////////////////////////////////////////////////////////////////////////////

  int getNumStickerTypes(int numLayers)

    {

    return numLayers<5 ? 1 : numLayers/2 + 1;

    }

///////////////////////////////////////////////////////////////////////////////////////////////////

  float getScreenRatio()

    {

    return 1.00f;

    }

///////////////////////////////////////////////////////////////////////////////////////////////////

  float[][] getCuts(int numLayers)

    {

    float[][] cuts = new float[6][numLayers-1];

    float D = numLayers*MovementMinx.DIST3D;

    float E = 2*SIN54;

    float X = 2*D*E/(1+2*E);  // height of the 'upper' part of a dodecahedron, i.e. put it on a table,

                              // its height is then D*2*DIST3D, it has one 'lower' part of height X, one

                              // 'middle' part of height Y and one upper part of height X again.

                              // It's edge length = numLayers/3.0f.

    int num = (numLayers-1)/2;

    float G = X*0.5f/num;     // height of one Layer

    for(int i=0; i<num; i++)

      {

      float cut = -D + (i+0.5f)*G;

      int j = 2*num-1-i;

      cuts[0][i] = +cut;

      cuts[0][j] = -cut;

      cuts[1][i] = +cut;

      cuts[1][j] = -cut;

      cuts[2][i] = +cut;

      cuts[2][j] = -cut;

      cuts[3][i] = +cut;

      cuts[3][j] = -cut;

      cuts[4][i] = +cut;

      cuts[4][j] = -cut;

      cuts[5][i] = +cut;

      cuts[5][j] = -cut;

      }

    return cuts;

    }

///////////////////////////////////////////////////////////////////////////////////////////////////

// Fill out mCurrCorner{X,Y,Z} by applying appropriate Quat to mBasicCorner{X,Y,Z}

// Appropriate one: QUATS[QUAT_INDICES[corner]].

  private void computeBasicCornerVectors(int corner)

    {

    Static4D quat = QUATS[QUAT_CORNER_INDICES[corner]];

    mCurrCornerV[0] = QuatHelper.rotateVectorByQuat(mBasicCornerV[0],quat);

    mCurrCornerV[1] = QuatHelper.rotateVectorByQuat(mBasicCornerV[1],quat);

    mCurrCornerV[2] = QuatHelper.rotateVectorByQuat(mBasicCornerV[2],quat);

    }

///////////////////////////////////////////////////////////////////////////////////////////////////

  private float[] computeCorner(int numCubitsPerCorner, int numLayers, int corner, int part)

    {

    float D = numLayers/3.0f;

    float[] corn = CORNERS[corner];

    if( part==0 )

      {

      return new float[] { corn[0]*D, corn[1]*D, corn[2]*D };

      }

    else

      {

      float E = D/(0.5f*(numLayers-1));   // ?? maybe 0.5*

      int N = (numCubitsPerCorner-1)/3;

      int block = (part-1) % N;

      int index = (part-1) / N;

      Static4D pri = mCurrCornerV[index];

      Static4D sec = mCurrCornerV[(index+2)%3];

      int layers= (numLayers-5)/2;

      int multP = (block % layers) + 1;

      int multS = (block / layers);

      return new float[] {

                          corn[0]*D + (pri.get0()*multP + sec.get0()*multS)*E,

                          corn[1]*D + (pri.get1()*multP + sec.get1()*multS)*E,

                          corn[2]*D + (pri.get2()*multP + sec.get2()*multS)*E

                         };

      }

    }

///////////////////////////////////////////////////////////////////////////////////////////////////

  private float[] computeCenter(int numLayers, int center, int part)

    {

    int corner = mCenterMap[center][part];

    float[] cent = mCenterCoords[center];

    float[] corn = CORNERS[corner];

    float D = numLayers/3.0f;

    float F = 1.0f - (2.0f*numLayers-6.0f)/(numLayers-1)*COS54*COS54;

    return new float[]

      {

        D * ( cent[0] + (corn[0]-cent[0])*F),

        D * ( cent[1] + (corn[1]-cent[1])*F),

        D * ( cent[2] + (corn[2]-cent[2])*F)

      };

    }

///////////////////////////////////////////////////////////////////////////////////////////////////

  private int computeEdgeType(int cubit, int numCubitsPerCorner, int numCubitsPerEdge)

    {

    int part = (cubit - NUM_CORNERS*numCubitsPerCorner) % numCubitsPerEdge;

    return part - 2*(part/4);

    }

///////////////////////////////////////////////////////////////////////////////////////////////////

  private float[] computeEdge(int numLayers, int edge, int part)

    {

    float D = numLayers/3.0f;

    float[] c1 = CORNERS[ mEdgeMap[edge][0] ];

    float[] c2 = CORNERS[ mEdgeMap[edge][1] ];

    int leftRight = 2*(part%2) -1;

    part /= 2;

    if( part==0 )

      {

      float T = 0.5f + leftRight/(numLayers-1.0f);

      float x = D * (T*c1[0]+(1.0f-T)*c2[0]);

      float y = D * (T*c1[1]+(1.0f-T)*c2[1]);

      float z = D * (T*c1[2]+(1.0f-T)*c2[2]);

      return new float[] { x, y, z };

      }

    else

      {

      int mult = (part+1)/2;

      int dir  = (part+1)%2;

      float[] center = mCenterCoords[ mEdgeMap[edge][dir+2] ];

      float x = 0.5f * D * (c1[0]+c2[0]);

      float y = 0.5f * D * (c1[1]+c2[1]);

      float z = 0.5f * D * (c1[2]+c2[2]);

      float vX = D*center[0] - x;

      float vY = D*center[1] - y;

      float vZ = D*center[2] - z;

      float T = 0.5f + leftRight*(mult*SIN18 + 1.0f)/(numLayers-1);

      x = D * (T*c1[0]+(1.0f-T)*c2[0]);

      y = D * (T*c1[1]+(1.0f-T)*c2[1]);

      z = D * (T*c1[2]+(1.0f-T)*c2[2]);

      float H = mult*D*COS18/(numLayers-1);

      H /= (float)Math.sqrt(vX*vX+vY*vY+vZ*vZ);

      return new float[] { x + H*vX, y + H*vY, z + H*vZ };

      }

    }

///////////////////////////////////////////////////////////////////////////////////////////////////

  float[][] getCubitPositions(int numLayers)

    {

    if( numLayers<5 ) return CORNERS;

    int numCubitsPerCorner = numCubitsPerCorner(numLayers);

    int numCubitsPerEdge   = numCubitsPerEdge(numLayers);

    int numCubitsPerCenter = 5;

    int numCubits = NUM_CORNERS*numCubitsPerCorner + NUM_EDGES*numCubitsPerEdge + NUM_CENTERS*numCubitsPerCenter;

    int index=0;

    final float[][] CENTERS = new float[numCubits][];

    for(int corner=0; corner<NUM_CORNERS; corner++)

      {

      computeBasicCornerVectors(corner);

      for(int part=0; part<numCubitsPerCorner; part++, index++)

        {

        CENTERS[index] = computeCorner(numCubitsPerCorner,numLayers,corner,part);

        }

      }

    for(int edge=0; edge<NUM_EDGES; edge++)

      {

      for(int part=0; part<numCubitsPerEdge; part++, index++)

        {

        CENTERS[index] = computeEdge(numLayers, edge, part );

        }

      }

    for(int center=0; center<NUM_CENTERS; center++)

      {

      for(int part=0; part<numCubitsPerCenter; part++, index++)

        {

        CENTERS[index] = computeCenter(numLayers,center, part);

        }

      }

    return CENTERS;

    }

///////////////////////////////////////////////////////////////////////////////////////////////////

  private int getQuat(int cubit, int numCubitsPerCorner, int numCubitsPerEdge)

    {

    if( cubit < NUM_CORNERS*numCubitsPerCorner )

      {

      int corner = cubit/numCubitsPerCorner;

      return QUAT_CORNER_INDICES[corner];

      }

    if( cubit < NUM_CORNERS*numCubitsPerCorner + NUM_EDGES*numCubitsPerEdge )

      {

      int edge = (cubit-NUM_CORNERS*numCubitsPerCorner)/numCubitsPerEdge;

      return QUAT_EDGE_INDICES[edge];

      }

    if( numCubitsPerCorner==0 )

      {

      return QUAT_CORNER_INDICES[cubit];

      }

    else

      {

      cubit -= (NUM_CORNERS*numCubitsPerCorner + NUM_EDGES*numCubitsPerEdge);

      int numCubitsPerCenter = 5;

      int face = cubit/numCubitsPerCenter;

      int index= cubit%numCubitsPerCenter;

      int corner=mCenterMap[face][index];

      return QUAT_CORNER_INDICES[corner];

      }

    }

///////////////////////////////////////////////////////////////////////////////////////////////////

  ObjectShape getObjectShape(int cubit, int numLayers)

    {

    int variant = getCubitVariant(cubit,numLayers);

    int numVariants = getNumCubitVariants(numLayers);

    if( variant==0 )

      {

      float width = (numLayers/3.0f)/(numLayers-1);

      float A = (2*SQ3/3)*SIN54;

      float B = 0.4f;

      double X = width*COS18*SIN_HALFD;

      double Y = width*SIN18;

      double Z = width*COS18*COS_HALFD;

      double[][] vertices = new double[][]

        {

            { 0.0, 0.0      , 0.0 },

            {   X,   Y      ,  -Z },

            { 0.0, 2*Y      ,-2*Z },

            {  -X,   Y      ,  -Z },

            { 0.0, 0.0-width, 0.0 },

            {   X,   Y-width,  -Z },

            { 0.0, 2*Y-width,-2*Z },

            {  -X,   Y-width,  -Z },

        };

      int[][] vertIndexes = new int[][]

        {

            {4,5,1,0},

            {7,4,0,3},

            {0,1,2,3},

            {4,5,6,7},

            {6,5,1,2},

            {7,6,2,3}

        };

      float[][] bands     = new float[][]

        {

         {0.04f,34,0.3f,0.2f, 3, 1, 0},

         {0.00f, 0,0.0f,0.0f, 2, 1, 0}

        };

      int[] bandIndices   = new int[] { 0,0,0,1,1,1};

      float[][] corners   = new float[][] { {0.04f,0.10f} };

      int[] cornerIndices = new int[] { 0,-1,-1,-1,-1,-1,-1,-1 };

      float[][] centers   = new float[][] { {0.0f, -(float)Math.sqrt(1-A*A)*B,-A*B} };

      int[] centerIndices = new int[] { 0,-1,-1,-1,-1,-1,-1,-1 };

      return new ObjectShape(vertices,vertIndexes,bands,bandIndices,corners,cornerIndices,centers,centerIndices,getNumCubitFaces(), null);

      }

    if( variant<numVariants-1 )

      {

      int numCubitsPerCorner = numCubitsPerCorner(numLayers);

      int numCubitsPerEdge   = numCubitsPerEdge(numLayers);

      int type = computeEdgeType(cubit,numCubitsPerCorner,numCubitsPerEdge);

      float tmpVal= (numLayers/3.0f)/(numLayers-1);

      float height= tmpVal*COS18;

      float width = tmpVal + (type/2)*tmpVal*SIN18;

      boolean left = (type%2)==0;

      double X = height*SIN_HALFD;

      double Y = height*SIN18/COS18;

      double Z = height*COS_HALFD;

      double[][] vertices = new double[][]

        {

            { 0.0, 0.0   , 0.0 },

            {   X,   Y   ,  -Z },

            { 0.0, 2*Y   ,-2*Z },

            {  -X,   Y   ,  -Z },

            { 0.0, -width, 0.0 },

            {   X, -width,  -Z },

            { 0.0, -width,-2*Z },

            {  -X, -width,  -Z },

        };

      int[][] vertIndexes = new int[][]

        {

            {4,5,1,0},

            {7,4,0,3},

            {7,6,2,3},

            {6,5,1,2},

            {0,1,2,3},

            {4,5,6,7}

        };

      if( !left )

        {

        int tmp, len = vertices.length;

        for(int i=0; i<len; i++) vertices[i][1] = -vertices[i][1];

        len = vertIndexes.length;

        for(int i=0; i<len; i++)

          {

          tmp = vertIndexes[i][0];

          vertIndexes[i][0] = vertIndexes[i][3];

          vertIndexes[i][3] = tmp;

          tmp = vertIndexes[i][1];

          vertIndexes[i][1] = vertIndexes[i][2];

          vertIndexes[i][2] = tmp;

          }

        }

      int numBands0 = numLayers<=5 ? 4 : 3;

      int numBands1 = numLayers<=5 ? 3 : 2;

      float[][] bands     = new float[][]

        {

         {0.04f,34,0.2f,0.2f,numBands0,1,1},

         {0.00f, 0,0.0f,0.0f,numBands1,0,0}

        };

      int[] bandIndices   = new int[] { 0,0,1,1,1,1};

      float[][] corners   = new float[][] { {0.04f,0.10f} };

      int[] cornerIndices = new int[] { 0,-1,-1,-1, 0,-1,-1,-1 };

      float[][] centers   = new float[][] { {0.0f, -width/2, (float)(-2*Z)} };

      int[] centerIndices = new int[] { 0,-1,-1,-1, 0,-1,-1,-1 };

      return new ObjectShape(vertices,vertIndexes,bands,bandIndices,corners,cornerIndices,centers,centerIndices,getNumCubitFaces(), null);

      }

    else

      {

      float width = (1+0.5f*(numLayers-3)*SIN18)*(numLayers/3.0f)/(numLayers-1);

      double X = width*COS18*SIN_HALFD;

      double Y = width*SIN18;

      double Z = width*COS18*COS_HALFD;

      double H = width*(SIN54/COS54);

      double H3= H/COS_HALFD;

      double X3= H*SIN_HALFD;

      double Z3= H*COS_HALFD;

      double C = 1/(COS54*Math.sqrt(2-2*SIN18));

      int N = numLayers==3 ? 4 : 3;

      int E = numLayers==3 ? 1 : 0;

      double[][] vertices = new double[][]

        {

            { 0.0, 0.0  ,   0.0 },

            {   X,   Y  ,    -Z },

            { 0.0,C*2*Y ,-2*C*Z },

            {  -X,   Y  ,    -Z },

            { 0.0,-width,   0.0 },

            {  X3,-width,   -Z3 },

            { 0.0,-width,   -H3 },

            { -X3,-width,   -Z3 }

        };

      int[][] vertIndexes = new int[][]

        {

            {4,5,1,0},

            {7,4,0,3},

            {0,1,2,3},

            {7,6,2,3},

            {6,5,1,2},

            {4,5,6,7}

        };

      float[][] bands = new float[][]

        {

         {0.04f,17,0.3f,0.2f,N,1,E},

         {0.00f,17,0.3f,0.2f,N,1,E}

        };

      float A = (2*SQ3/3)*SIN54;

      float B = 0.4f;

      int[] bandIndices   = new int[] { 0,0,0,1,1,1};

      float[][] corners   = new float[][] { {0.03f,0.10f} };

      int[] cornerIndices = new int[] { 0, 0,-1, 0, 0,-1,-1,-1 };

      float[][] centers   = new float[][] { {0.0f, -(float)Math.sqrt(1-A*A)*B,-A*B} };

      int[] centerIndices = new int[] { 0, 0,-1, 0, 0,-1,-1,-1 };

      return new ObjectShape(vertices,vertIndexes,bands,bandIndices,corners,cornerIndices,centers,centerIndices,getNumCubitFaces(), null);

      }

    }

///////////////////////////////////////////////////////////////////////////////////////////////////

  private Static4D getQuat(int cubit, int numLayers)

    {

    int numCubitsPerCorner = numCubitsPerCorner(numLayers);

    int numCubitsPerEdge   = numCubitsPerEdge(numLayers);

    return QUATS[getQuat(cubit,numCubitsPerCorner,numCubitsPerEdge)];

    }

///////////////////////////////////////////////////////////////////////////////////////////////////

  private int getNumCubitVariants(int numLayers)

    {

    int[] sizes = ObjectList.KILO.getSizes();

    int variants = sizes.length;

    int highestSize = sizes[variants-1];

    return highestSize-1;

    }

///////////////////////////////////////////////////////////////////////////////////////////////////

  int getCubitVariant(int cubit, int numLayers)

    {

    int numCubitsPerCorner = numCubitsPerCorner(numLayers);

    if( cubit<NUM_CORNERS*numCubitsPerCorner ) return 0;

    int numCubitsPerEdge   = numCubitsPerEdge(numLayers);

    if( cubit<NUM_CORNERS*numCubitsPerCorner + NUM_EDGES*numCubitsPerEdge )

      {

      int type = computeEdgeType(cubit,numCubitsPerCorner,numCubitsPerEdge);

      return type+1;

      }

    int[] sizes = ObjectList.KILO.getSizes();

    int variants = sizes.length;

    int highestSize = sizes[variants-1];

    return highestSize-2;

    }

///////////////////////////////////////////////////////////////////////////////////////////////////

  MeshBase createCubitMesh(int cubit, int numLayers)

    {

    int variant = getCubitVariant(cubit,numLayers);

    if( mMeshes==null )

      {

      FactoryCubit factory = FactoryCubit.getInstance();

      factory.clear();

      mMeshes = new MeshBase[getNumCubitVariants(numLayers)];

      }

    if( mMeshes[variant]==null )

      {

      ObjectShape shape = getObjectShape(cubit,numLayers);

      FactoryCubit factory = FactoryCubit.getInstance();

      factory.createNewFaceTransform(shape);

      mMeshes[variant] = factory.createRoundedSolid(shape);

      }

    MeshBase mesh = mMeshes[variant].copy(true);

    MatrixEffectQuaternion quat = new MatrixEffectQuaternion( getQuat(cubit,numLayers), new Static3D(0,0,0) );

    mesh.apply(quat,0xffffffff,0);

    return mesh;

    }

///////////////////////////////////////////////////////////////////////////////////////////////////

  int getCornerColor(int cubit, int cubitface, int numLayers, int numCubitsPerCorner)

    {

    if( cubitface<0 || cubitface>2 ) return NUM_TEXTURES;

    int part  = cubit % numCubitsPerCorner;

    int corner= cubit / numCubitsPerCorner;

    if( part==0 )

      {

      return mCornerFaceMap[corner][cubitface];

      }

    else

      {

      int N = (numCubitsPerCorner-1)/3;

      int block = (part-1) % N;

      int index = (part-1) / N;

      if( block< (numLayers-3)/2 )

        {

        switch(index)

          {

          case 0: return cubitface==1 ? NUM_TEXTURES : mCornerFaceMap[corner][cubitface];

          case 1: return cubitface==0 ? NUM_TEXTURES : mCornerFaceMap[corner][cubitface];

          case 2: return cubitface==2 ? NUM_TEXTURES : mCornerFaceMap[corner][cubitface];

          }

        }

      else

        {

        switch(index)

          {

          case 0: return cubitface==0 ? mCornerFaceMap[corner][cubitface] : NUM_TEXTURES;

          case 1: return cubitface==2 ? mCornerFaceMap[corner][cubitface] : NUM_TEXTURES;

          case 2: return cubitface==1 ? mCornerFaceMap[corner][cubitface] : NUM_TEXTURES;

          }

        }

      }

    return NUM_TEXTURES;

    }

///////////////////////////////////////////////////////////////////////////////////////////////////

  int getEdgeColor(int edge, int cubitface, int numCubitsPerEdge)

    {

    if( cubitface<0 || cubitface>1 ) return NUM_TEXTURES;

    int part    = edge % numCubitsPerEdge;

    int variant = edge / numCubitsPerEdge;

    part /=2;

    return (part==0 || cubitface==((part+1)%2)) ? mEdgeMap[variant][cubitface+2] + ((part+3)/2)*NUM_FACES : NUM_TEXTURES;

    }

///////////////////////////////////////////////////////////////////////////////////////////////////

  int getCenterColor(int center, int cubitface, int numLayers)

    {

     if( numLayers==3 )

      {

      return cubitface>=0 && cubitface<3 ? mCornerFaceMap[center][cubitface] : NUM_TEXTURES;

      }

    return cubitface==mCenterFaceMap[center] ? center/5 + NUM_FACES*(numLayers-1)/2 : NUM_TEXTURES;

    }

///////////////////////////////////////////////////////////////////////////////////////////////////

  int getFaceColor(int cubit, int cubitface, int numLayers)

    {

    int numCubitsPerCorner = numCubitsPerCorner(numLayers);

    int numCubitsPerEdge   = numCubitsPerEdge(numLayers);

    if( cubit < NUM_CORNERS*numCubitsPerCorner )

      {

      return getCornerColor(cubit,cubitface,numLayers,numCubitsPerCorner);

      }

    else if( cubit<NUM_CORNERS*numCubitsPerCorner + NUM_EDGES*numCubitsPerEdge )

      {

      int edge = cubit - NUM_CORNERS*numCubitsPerCorner;

      return getEdgeColor(edge,cubitface,numCubitsPerEdge);

      }

    else

      {

      int center = cubit-NUM_CORNERS*numCubitsPerCorner-NUM_EDGES*numCubitsPerEdge;

      return getCenterColor( center, cubitface, numLayers);

      }

    }

///////////////////////////////////////////////////////////////////////////////////////////////////

  int getColor(int face)

    {

    return FACE_COLORS[face];

    }

///////////////////////////////////////////////////////////////////////////////////////////////////