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.library.main.DistortedTexture;

import org.distorted.library.mesh.MeshSquare;

import org.distorted.library.type.Static4D;

import org.distorted.main.R;

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

  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);

    }

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

  void initializeCenterFaceMap()

    {

    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

      };

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

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

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

  float getScreenRatio()

    {

    return 1.00f;

    }

                              // 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++)

      {

      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)

    {

    if( mQuats==null ) initializeQuats();

    if( mCurrCornerV==null || mBasicCornerV==null ) initializeCornerV();

    Static4D quat = mQuats[mQuatCornerIndices[corner]];

    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)

    {

    if( mCurrCornerV==null || mBasicCornerV==null ) initializeCornerV();

    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)

    {

    if( mCorners     ==null ) initializeCorners();

    if( mCenterMap   ==null ) initializeCenterMap();

    float[] cent = mCenterCoords[center];

    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;

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

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

    {

    if( mCorners==null ) initializeCorners();

    if( mEdgeMap==null ) initializeEdgeMap();

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

    if( part==0 )

      {

      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]);

      }

    else

      {

      int mult = (part+1)/2;

      int dir  = (part+1)%2;

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

      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;

      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);

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

  float[][] getCubitPositions(int numLayers)

    {

    if( numLayers<5 ) return mCorners;

    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 part=0; part<numCubitsPerCenter; part++, index++)

        {

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

        }

      }

    return CENTERS;

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

    if( mCenterMap==null ) initializeCenterMap();

      {

      int corner = cubit/numCubitsPerCorner;

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

      {

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

    if( numCubitsPerCorner==0 )

      {

    else

      {

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

      int numCubitsPerCenter = 5;

      int face = cubit/numCubitsPerCenter;

      int index= cubit%numCubitsPerCenter;

      int corner=mCenterMap[face][index];

    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[][]

            {   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 },

            {4,5,1,0},

            {7,4,0,3},

            {0,1,2,3},

        {

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

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

        };

      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 };

      }

    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 },

        };

            {4,5,1,0},

            {7,4,0,3},

            {7,6,2,3},

            {6,5,1,2},

            {0,1,2,3},

            {4,5,6,7}

        };

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

        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 numBands1 = numLayers<=5 ? 3 : 2;

        {

         {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);

      }

    int numCubitsPerEdge   = numCubitsPerEdge(numLayers);

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

    int variants = sizes.length;

    int highestSize = sizes[variants-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 )

      return type+1;

    int variants = sizes.length;

    int highestSize = sizes[variants-1];

    }

    int part  = cubit % numCubitsPerCorner;

    int corner= cubit / numCubitsPerCorner;

    if( part==0 )

      {

      int N = (numCubitsPerCorner-1)/3;

      int block = (part-1) % N;

      int index = (part-1) / N;