//VTK::System::Dec
// SPDX-FileCopyrightText: Copyright (c) Ken Martin, Will Schroeder, Bill Lorensen
// SPDX-License-Identifier: BSD-3-Clause
/// Model-to-device coordinate transform
uniform mat4 MCDCMatrix;
/// Model-to-view coordinate transform
uniform mat4 MCVCMatrix;
/// matrix used to transform normals from model to view coordinate
uniform mat3 normalMatrix;
/// View coordinate normal for this vertex.
smooth out vec3 vertexNormalVCVS;

uniform samplerBuffer vertices;

// The parametric coordinates of each shape-attribute integration point.
uniform samplerBuffer cell_parametrics;

// Tuples of (cell-id, side-id) to render.
uniform highp isamplerBuffer sides;
// Look up the offset at which local cell-side connectivity starts in \a side_local.
uniform highp isamplerBuffer side_offsets;
// The subset of DOFs in shape_conn that correspond to each side in cells of this shape.
uniform highp isamplerBuffer side_local;
// The indices of all DOFs for all cells of this shape.
uniform highp isamplerBuffer shape_conn;
// The (x,y,z) points of each DOF in the entire mesh.
uniform samplerBuffer shape_vals;
// We may need this if the color and shape have different interpolation orders
// **and** the color attribute is continuous (i.e., shares coefficients at cell boundaries).
uniform highp isamplerBuffer color_conn;
// The coefficient of each basis function at each integration point in each cell.
uniform samplerBuffer color_vals;
// For vector- or tensor-valued field attributes, which should determine the color?
// If -1 or -2, the L1-norm or L2-norm are used.
uniform int field_component;

flat out int cellIdVS;
flat out int sideIdVS;

{commonDefs}
{cellEval}
{cellUtil}

/// Parametric coordinates of this vertex.
smooth out vec3 pcoordVS;
/// View coordinate normal for this vertex.
smooth out vec3 normalVCVS;
/// View coordinate position for this vertex.
smooth out vec4 vertexPositionVCVS;

/// Store coefficients for one finite element's shape and color attributes.
///
/// These are used by downstream shader(s) to interpolate values at
/// intermediate vertices (i.e., tessellation control shader) and fragments
/// for pixel-accurate renderings. Fetching these in the vertex shader
/// saves the work of fetching them for each fragment.
flat out float shapeValuesVS[{ShapeCoeffPerCell}];
flat out float colorValuesVS[{ColorCoeffPerCell}];

void main()
{{
  // Identify the cell ID and side ID (if applicable) we are to render.
  ivec2 cellAndSide;
  if (DrawingCellsNotSides)
  {{
     cellAndSide = ivec2(gl_InstanceID, 0);
     // These two VS outputs are to be used for picking:
     cellIdVS = cellAndSide.s;
     sideIdVS = -1;
     // Note sideIdVS != cellAndSide.t; this is so that we can properly
     // compute the cell connectivity (using the zero offset in cellAndSide.t)
     // while also computing the normal (using sideIdVS).
  }}
  else
  {{
     cellAndSide = texelFetchBuffer(sides, gl_InstanceID).st;
     // These two VS outputs are to be used for picking:
     cellIdVS = cellAndSide.s;
     sideIdVS = cellAndSide.t;
  }}

  // Fetch the offset into the (ragged) side_local table:
  int sideRaggedOffset = texelFetchBuffer(side_offsets, {ShapeIndex}).s;

  // Fetch point coordinates and per-integration-point field values for
  // the entire cell, passing them in bulk to the fragment shader using
  // a "flat" keyword.
  //
  // Doing this lookup in the vertex shader saves lots of work in the
  // fragment shader (usu. invoked many times per vertex).
  //
  // NB: Currently, a cell-grid's shape attribute *must* be continuous
  // (i.e., share degrees of freedom at cell boundaries). This makes
  // fetching shapeValuesVS much simpler than fetching colorValuesVS.
  for (int ii = 0; ii < {ShapeNumBasisFun}; ++ii)
  {{
    int vertexId = texelFetchBuffer(shape_conn, cellIdVS * {ShapeNumBasisFun} + ii).s;
    for (int jj = 0; jj < {ShapeMultiplicity}; ++jj)
    {{
      shapeValuesVS[ii * {ShapeMultiplicity} + jj] = texelFetchBuffer(shape_vals, vertexId * {ShapeMultiplicity} + jj).x;
    }}
  }}

  if ({HaveColors})
  {{
    if ({ColorContinuous})
    {{
      // Continuous (shared) field values
      for (int ii = 0; ii < {ColorNumBasisFun}; ++ii)
      {{
        int dofId = texelFetchBuffer(color_conn, cellIdVS * {ColorNumBasisFun} + ii).s;
        for (int jj = 0; jj < {ColorMultiplicity}; ++jj)
        {{
          colorValuesVS[ii * {ColorMultiplicity} + jj] = texelFetchBuffer(color_vals, dofId * {ColorMultiplicity} + jj).x;
        }}
      }}
    }}
    else
    {{
      // Discontinuous field values
      float cv;
      for (int ii = 0; ii < {ColorNumBasisFun}; ++ii)
      {{
        for (int jj = 0; jj < {ColorMultiplicity}; ++jj)
        {{
          cv = texelFetchBuffer(color_vals, (cellIdVS * {ColorNumBasisFun} + ii) * {ColorMultiplicity} + jj).x;
          colorValuesVS[ii * {ColorMultiplicity} + jj] = cv;
        }}
      }}
    }}
  }}

  // Now compute the location of the current vertex inside the current side or cell.
  // NB: {{SideOffset}} = {SideOffset} is the offset to apply to side IDs so we can
  // look up their connectivity offset from sideRaggedOffset (the start in side_local
  // for connectivity of sides of our type).
  int sideVertexIndex =
    texelFetchBuffer(side_local,
      sideRaggedOffset + (cellAndSide.t - {SideOffset}) * NumPtsPerSide + gl_VertexID).s;
  int vertexId = texelFetchBuffer(shape_conn, cellIdVS * NumPtsPerCell + sideVertexIndex).s;
  // Parametric coordinate for this vertex.
  pcoordVS = texelFetchBuffer(cell_parametrics, sideVertexIndex).xyz;
  // position for this vertex as defined in vtk data model.
  vec4 vertexMC = vec4(
      texelFetchBuffer(shape_vals, vertexId * 3).x,
      texelFetchBuffer(shape_vals, vertexId * 3 + 1).x,
      texelFetchBuffer(shape_vals, vertexId * 3 + 2).x,
      1.0f);
  // default eye direction in model coordinates.
  vec3 eyeNormalMC = vec3(0.0f, 0.0f, 1.0f);
  vec3 vertexNormalMC = normalToSideAt(sideIdVS, shapeValuesVS, pcoordVS, -eyeNormalMC);

  // Transform the vertex by the model-to-device coordinate matrix.
  // This matrix must be the result of the following multiplication:
  // MCDCMatrix = ModelToWorld X WorldToView X ViewToDisplay
  gl_Position = MCDCMatrix * vertexMC;

  // Transform vertex posittion to view coordinate.
  // Some operations in fragment shader want it that way.
  vertexPositionVCVS = MCVCMatrix * vertexMC;

  // Normal vectors are transformed in a different way than vertices.
  // Instead of pre-multiplying with MCDCMatrix, a different matrix is used.
  // This `normalMatrix` is computed on the CPU. It must be the result of the following:
  // normalMatrix = inverse(ModelToWorld X WorldToView)
  // Read more about normal matrix at http://www.songho.ca/opengl/gl_normaltransform.html
  if ((DrawingCellsNotSides && NumPtsPerCell <= 2) || (!DrawingCellsNotSides && NumPtsPerSide <= 2))
  {{
    // for lines or vertices, the normal will always face the camera.
    vertexNormalVCVS = vec3(vertexNormalMC.xy, 1.0f);
  }}
  else
  {{
    vertexNormalVCVS = normalMatrix * vertexNormalMC;
  }}
}}
