Distorted Android

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

// Copyright 2020 Leszek Koltunski  leszek@koltunski.pl                                          //

//                                                                                               //

// This file is part of Distorted.                                                               //

//                                                                                               //

// This library is free software; you can redistribute it and/or                                 //

// modify it under the terms of the GNU Lesser General Public                                    //

// License as published by the Free Software Foundation; either                                  //

// version 2.1 of the License, or (at your option) any later version.                            //

//                                                                                               //

// This library 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                             //

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

//                                                                                               //

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

// License along with this library; if not, write to the Free Software                           //

// Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301  USA                //

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

package org.distorted.library.mesh;

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

import org.distorted.library.main.DistortedLibrary;

/**

 * Create a polygon of any shape and varying elevations from the edges towards the center.

 * <p>

 * Specify a list of vertices. Any two adjacent vertices + the center (0,0,0) form a triangle. The

 * polygon is going to be split into such triangles, and each triangle is split into adjustable number

 * of 'bands' form the outer edge towards the center. Edges of each band can can at any elevation.

 */

public class MeshPolygon extends MeshBase

  {

  private static final float NOT_DONE_YET = -1000;

  private static final int SHAPE_DD = 0;

  private static final int SHAPE_DU = 1;

  private static final int SHAPE_UD = 2;

  private static final int SHAPE_UU = 3;

  private static final int SHAPE_DUD= 4;

  private float[] mPolygonVertices;

  private int mNumPolygonVertices;

  private float[] mPolygonBands;

  private int mNumPolygonBands;

  private boolean[] mEdgeUp;

  private int remainingVert;

  private int numVertices;

  private int extraIndex, extraVertices;

  private float[] mBandQuot;

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

// polygonVertices>=3 , polygonBands>=2

  private void computeNumberOfVertices()

     {

     if( mNumPolygonBands==2 && extraIndex>0 )

       {

       numVertices = 1 + 2*mNumPolygonVertices*(1+extraIndex+2*extraVertices);

       }

     else

       {

       numVertices = (mNumPolygonVertices*mNumPolygonBands+2)*(mNumPolygonBands-1) - 1;

       numVertices+= 2*mNumPolygonVertices*(2*extraIndex*extraVertices);

       }

     remainingVert = numVertices;

     }

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

  private void computeCache()

    {

    mBandQuot = new float[mNumPolygonBands];

    int next, prev;

    for(int band=0; band<mNumPolygonBands; band++)

      {

      next = (band==mNumPolygonBands-1 ? band : band+1);

      prev = (band==                 0 ? band : band-1);

      mBandQuot[band] = (mPolygonBands[2*prev+1]-mPolygonBands[2*next+1]) / (mPolygonBands[2*next]-mPolygonBands[2*prev]);

      }

    }

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

  private float getSpecialQuot(int index)

    {

    int num = 1 + extraIndex + 2*extraVertices;

    int change1 = extraVertices+1;

    int change2 = num-change1;

    float quot = 1.0f/(extraIndex+1);

    if( index<change1 )      return index*quot/change1;

    else if( index>change2 ) return 1-quot + (index-change2)*quot/change1;

    else                     return (index-change1+1)*quot;

    }

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

  private float getQuot(int index, int band, boolean isExtra)

    {

    int num = mNumPolygonBands-1-band;

    if( num>0 )

      {

      if( isExtra )

        {

        int extra = extraIndex-band+extraVertices;

        if( index < extra )

          {

          float quot = ((float)extraIndex-band)/(extra*num);

          return index*quot;

          }

        else if( index > num+2*extraVertices-extra )

          {

          float quot = ((float)extraIndex-band)/(extra*num);

          return (1.0f-((float)extraIndex-band)/num) + (index-num-2*extraVertices+extra)*quot;

          }

        else

          {

          return ((float)(index-extraVertices))/num;

          }

        }

      return (float)index/num;

      }

    return 1.0f;

    }

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

  private float computeZEdge(float quot)

    {

    if( quot>=1.0f ) return 0.0f;

    for(int band=1; band<mNumPolygonBands; band++)

      {

      float curr = mPolygonBands[2*band];

      if( curr<=quot )

        {

        float prev = mPolygonBands[2*band-2];

        float prevH= mPolygonBands[2*band-1];

        float currH= mPolygonBands[2*band+1];

        float A = (prev-quot)/(prev-curr);

        return A*currH + (1-A)*prevH;

        }

      }

    return 0.0f;

    }

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

  private int addVertex(int vertex, int polyBand, int polyVertex, int polyEndVer, float quot,

                        int edgeShape, float[] attribs1, float[] attribs2)

    {

    remainingVert--;

    float Xfirst= mPolygonVertices[2*polyVertex  ];

    float Yfirst= mPolygonVertices[2*polyVertex+1];

    float Xlast = mPolygonVertices[2*polyEndVer  ];

    float Ylast = mPolygonVertices[2*polyEndVer+1];

    float xEdge = Xfirst + quot*(Xlast-Xfirst);

    float yEdge = Yfirst + quot*(Ylast-Yfirst);

    float zEdge;

    switch(edgeShape)

      {

      case SHAPE_DD : zEdge = 0.0f;

                      break;

      case SHAPE_UU : zEdge = mPolygonBands[2*mNumPolygonBands-1];

                      break;

      case SHAPE_DU : zEdge = quot>=0.5f ? mPolygonBands[2*mNumPolygonBands-1] : computeZEdge(1-2*quot);

                      break;

      case SHAPE_UD : zEdge = quot<=0.5f ? mPolygonBands[2*mNumPolygonBands-1] : computeZEdge(2*quot-1);

                      break;

      default       : zEdge = quot<=0.5f ? computeZEdge(1-2*quot) : computeZEdge(2*quot-1);

                      break;

      }

    float q =  mPolygonBands[2*polyBand];

    float x = q*xEdge;

    float y = q*yEdge;

    float z = q*zEdge + mPolygonBands[2*polyBand+1];

    attribs1[VERT1_ATTRIBS*vertex + POS_ATTRIB  ] = x;

    attribs1[VERT1_ATTRIBS*vertex + POS_ATTRIB+1] = y;

    attribs1[VERT1_ATTRIBS*vertex + POS_ATTRIB+2] = z;

    int index = VERT1_ATTRIBS*vertex + NOR_ATTRIB;

    if( quot==0.0f || quot==1.0f )

      {

      float vx = mBandQuot[polyBand]*xEdge;

      float vy = mBandQuot[polyBand]*yEdge;

      float vz = xEdge*xEdge + yEdge*yEdge;

      float len = (float)Math.sqrt(vx*vx + vy*vy + vz*vz);

      attribs1[index  ] = vx/len;

      attribs1[index+1] = vy/len;

      attribs1[index+2] = vz/len;

      }

    else

      {

      attribs1[index  ] = NOT_DONE_YET;

      }

    attribs2[VERT2_ATTRIBS*vertex + TEX_ATTRIB  ] = x+0.5f;

    attribs2[VERT2_ATTRIBS*vertex + TEX_ATTRIB+1] = y+0.5f;

    return vertex+1;

    }

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

  private void addVertexNormal(int vertex, int end1, int end2, float[] attribs1)

    {

    // android.util.Log.e("D", "vertex="+vertex+" end1="+end1+" end2="+end2);

    int index = VERT1_ATTRIBS*vertex + NOR_ATTRIB;

    if( attribs1[index] == NOT_DONE_YET)

      {

      int iv = VERT1_ATTRIBS*vertex + POS_ATTRIB;

      int i1 = VERT1_ATTRIBS*end1   + POS_ATTRIB;

      int i2 = VERT1_ATTRIBS*end2   + POS_ATTRIB;

      float vx = attribs1[iv  ];

      float vy = attribs1[iv+1];

      float vz = attribs1[iv+2];

      float x1 = attribs1[i1  ];

      float y1 = attribs1[i1+1];

      float z1 = attribs1[i1+2];

      float x2 = attribs1[i2  ];

      float y2 = attribs1[i2+1];

      float z2 = attribs1[i2+2];

      float dx1 = vx-x1;

      float dy1 = vy-y1;

      float dz1 = vz-z1;

      float dx2 = vx-x2;

      float dy2 = vy-y2;

      float dz2 = vz-z2;

      float cx = dy1*dz2 - dz1*dy2;

      float cy = dz1*dx2 - dx1*dz2;

      float cz = dx1*dy2 - dy1*dx2;

      float len = (float)Math.sqrt(cx*cx + cy*cy + cz*cz);

      attribs1[index  ] = cx/len;

      attribs1[index+1] = cy/len;

      attribs1[index+2] = cz/len;

      }

    }

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

  private int createBandStrip(int vertex, int polyBand, int polyVertex, int edgeShape, float[] attribs1, float[] attribs2)

    {

    int initVertex = vertex;

    if( polyVertex==0 )

      {

      vertex = addVertex(vertex,polyBand,0,1,0,edgeShape,attribs1,attribs2);

      if( polyBand>0 ) vertex = addVertex(vertex,polyBand,0,1,0,edgeShape,attribs1,attribs2);

      }

    boolean specialCase = mNumPolygonBands==2 && polyBand==0 && extraIndex>0;

    boolean isExtra = polyBand<extraIndex;

    int numPairs = specialCase ? extraIndex+1 : mNumPolygonBands-1-polyBand;

    if( isExtra ) numPairs += 2*extraVertices;

    int polyEndVer = polyVertex==mNumPolygonVertices-1 ? 0 : polyVertex+1;

    float quot1, quot2;

    for(int index=0; index<numPairs; index++)

      {

      if( specialCase )

        {

        quot1 = 1.0f;

        quot2 = getSpecialQuot(index+1);

        }

      else

        {

        quot1 = getQuot(index  ,polyBand+1, isExtra);

        quot2 = getQuot(index+1,polyBand  , isExtra);

        }

      vertex = addVertex(vertex,polyBand+1,polyVertex,polyEndVer,quot1,edgeShape,attribs1,attribs2);

      vertex = addVertex(vertex,polyBand  ,polyVertex,polyEndVer,quot2,edgeShape,attribs1,attribs2);

      }

    if( polyVertex==0 )

      {

      if( polyBand>0 ) addVertexNormal(initVertex,initVertex+3,initVertex+2,attribs1);

      else             addVertexNormal(initVertex,initVertex+2,initVertex+1,attribs1);

      }

    else

      {

      addVertexNormal(initVertex,initVertex-1,initVertex+1,attribs1);

      }

    boolean lower = (polyVertex>0 || polyBand>0);

    for(int index=initVertex+1; index<vertex-1; index++)

      {

      addVertexNormal(index,(lower ? index+2 : index-1),index+1,attribs1);

      lower = !lower;

      }

    addVertexNormal(vertex-1,vertex-2,vertex-3,attribs1);

    return vertex;

    }

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

  private int computeEdgeShape(int curr)

    {

    if( mEdgeUp==null || !mEdgeUp[curr] ) return SHAPE_DD;

    int prev = (curr==0 ? mNumPolygonVertices-1 : curr-1);

    int next = (curr==mNumPolygonVertices-1 ? 0 : curr+1);

    boolean l = mEdgeUp[prev];

    boolean r = mEdgeUp[next];

    return l ? (r ? SHAPE_UU : SHAPE_UD) : (r ? SHAPE_DU : SHAPE_DUD);

    }

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

  private void buildGrid(float[] attribs1, float[] attribs2)

    {

    int vertex=0;

    int[] edgeShape = new int[mNumPolygonVertices];

    for(int polyVertex=0; polyVertex<mNumPolygonVertices; polyVertex++)

      edgeShape[polyVertex] = computeEdgeShape(polyVertex);

    for(int polyBand=0; polyBand<mNumPolygonBands-1; polyBand++)

      for(int polyVertex=0; polyVertex<mNumPolygonVertices; polyVertex++)

        {

       // android.util.Log.e("D", "creating strip polyBand="+polyBand+" polyVertex="+polyVertex+" : "+vertex);

        vertex = createBandStrip(vertex,polyBand,polyVertex,edgeShape[polyVertex],attribs1,attribs2);

       // android.util.Log.e("D", "  "+vertex);

        }

    }

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

// PUBLIC API

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

/**

 * Create a polygon of any shape and varying elevations from the edges towards the center.

 * Optionally make it more dense at the vertices.

 *

 * @param verticesXY 2N floats - packed description of polygon vertices. N pairs (x,y).

 *                   Vertices HAVE TO be specified in a COUNTERCLOCKWISE order (starting from any).

 * @param bands      2K floats; K pairs of two floats each describing a single band.

 *                   From (1.0,Z[0]) (outer edge, its Z elevation) to (0.0,Z[K]) (the center,

 *                   its elevation). The polygon is split into such concentric bands.

 *                   Must be band[2*i] > band[2*(i+1)] !

 * @param edgeUp     N booleans - one for each edge; the first connects vertices (0,1) and (2,3);

 *                   then just like 'verticesXY', i.e. counterclockwise.

 *                   If this is null, all edges are by default 'down'.

 *                   'Down' means that edge's Z is equal to 0; 'up' means that its Z, at least in its

 *                   middle, is equal to the highest elevation in the middle of the mesh.

 *                   If the 'previous' edge is also up, then the Z is up horizontally from its middle

 *                   to the left vertex; else it goes down along the same function it goes along the bands.

 *                   Same with the right half of the edge.

 * @param exIndex    This and the next parameter describe how to make the mesh denser at the

 *                   polyVertices. If e.g. exIndex=3 and exVertices=2, then 3 triangles of the

 *                   outermost band (and 2 triangles of the next band, and 1 triangle of the third

 *                   band) get denser - the 3 triangles become 3+2 = 5.

 * @param exVertices See above.

 * @param centerX    the X coordinate of the 'center' of the Polygon, i.e. point of the mesh

 *                   all bands go to.

 * @param centerY    Y coordinate of the center.

 */

  public MeshPolygon(float[] verticesXY, float[] bands, boolean[] edgeUp, int exIndex, int exVertices, float centerX, float centerY)

    {

    super();

    mPolygonVertices   = verticesXY;

    mPolygonBands      = bands;

    mNumPolygonVertices= mPolygonVertices.length /2;

    mNumPolygonBands   = mPolygonBands.length /2;

    mEdgeUp            = edgeUp;

    extraIndex         = exIndex;

    extraVertices      = exVertices;

    if( centerX!=0.0f || centerY!=0.0f )

      {

      for(int v=0; v<mNumPolygonVertices; v++)

        {

        mPolygonVertices[2*v  ] -= centerX;

        mPolygonVertices[2*v+1] -= centerY;

        }

      }

    computeNumberOfVertices();

    computeCache();

    float[] attribs1= new float[VERT1_ATTRIBS*numVertices];

    float[] attribs2= new float[VERT2_ATTRIBS*numVertices];

    buildGrid(attribs1,attribs2);

    if( remainingVert!=0 )

      DistortedLibrary.logMessage("MeshPolygon: remainingVert " +remainingVert );

    if( centerX!=0.0f || centerY!=0.0f )

      {

      for(int v=0; v<numVertices; v++)

        {

        attribs1[VERT1_ATTRIBS*v + POS_ATTRIB  ] += centerX;

        attribs1[VERT1_ATTRIBS*v + POS_ATTRIB+1] += centerY;

        attribs2[VERT2_ATTRIBS*v + TEX_ATTRIB  ] += centerX;

        attribs2[VERT2_ATTRIBS*v + TEX_ATTRIB+1] += centerY;

        }

      }

    setAttribs(attribs1,attribs2);

    }

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

  public MeshPolygon(float[] verticesXY, float[] bands, int exIndex, int exVertices)

    {

    this(verticesXY,bands,null,exIndex,exVertices,0.0f,0.0f);

    }

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

/**

 * Create a polygon of any shape and varying elevations from the edges towards the center.

 * Equivalent of the previous with exIndex=0 or exVertices=0.

 */

  public MeshPolygon(float[] verticesXY, float[] bands)

    {

    this(verticesXY,bands,null,0,0,0.0f,0.0f);

    }

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

/**

 * Copy constructor.

 */

  public MeshPolygon(MeshPolygon mesh, boolean deep)

    {

    super(mesh,deep);

    }

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

/**

 * Copy the Mesh.

 *

 * @param deep If to be a deep or shallow copy of mVertAttribs1, i.e. the array holding vertices,

 *             normals and inflates (the rest, in particular the mVertAttribs2 containing texture

 *             coordinates and effect associations, is always deep copied)

 */

  public MeshPolygon copy(boolean deep)

    {

    return new MeshPolygon(this,deep);

    }

 }