struct GLvector {
 GLfloat fX, fY, fZ; };
...
GLvector asEdgeVertex[12];
GLvector asEdgeNorm[12];
...
vGetNormal(asEdgeNorm[iEdge], asEdgeVertex[iEdge].fX, asEdgeVertex[iEdge].fY, asEdgeVertex[iEdge].fZ);

  //Make a local copy of the values at the cube's corners
  for(iVertex = 0; iVertex < 8; iVertex++)
  {
    afCubeValue[iVertex] = fSample(fX + a2fVertexOffset[iVertex][0]*fScale,
               fY + a2fVertexOffset[iVertex][1]*fScale,
               fZ + a2fVertexOffset[iVertex][2]*fScale);
  }

 
// 3D (4*4*4) массив представленный как одномерный
static const GLfloat MyValues[64]={
20,20,20,20,20,20,20,20,20,20,20,20,20,20,20,20,
20,20,20,20,20, 0, 0,20,20, 0, 0,20,20,20,20,20,
20,20,20,20,20, 0, 0,20,20, 0, 0,20,20,20,20,20,
20,20,20,20,20,20,20,20,20,20,20,20,20,20,20,20
};
...
void __fastcall TForm1::vMarchCube1(GLint iX, GLint iY, GLint iZ, GLfloat fIsoLevel)
{
  ...
  // Делаем локальную копию данных кубика из большого массива
  afCubeValue[0] = MyValues [iZ+4*iX+iY*4*4];
  afCubeValue[1] = MyValues [iZ+4*(iX+1)+iY*4*4];
  afCubeValue[2] = MyValues [iZ+4*(iX+1)+(iY+1)*4*4];
  afCubeValue[3] = MyValues [iZ+4*iX+(iY+1)*4*4];
  afCubeValue[4] = MyValues [(iZ+1)+4*iX+iY*4*4];
  afCubeValue[5] = MyValues [(iZ+1)+4*(iX+1)+iY*4*4];
  afCubeValue[6] = MyValues [(iZ+1)+4*(iX+1)+(iY+1)*4*4];
  afCubeValue[7] = MyValues [(iZ+1)+4*iX+(iY+1)*4*4];
  ...
}

// линейная интерполяция
GLfloat fBlend( GLfloat a, GLfloat b, GLfloat t )
{
  return a * ( 1.0 - t ) + b * t;
}

// важно: fX, fY, fZ должны соответствовать размерам имеющегося массива
GLfloat fSample5( GLfloat fX, GLfloat fY, GLfloat fZ )
{
  GLfloat dX, dY, dZ;
  GLint iX, iY, iZ;

  GLfloat fCV[8];

  iX = (int)fX;
  iY = (int)fY;
  iZ = (int)fZ;

  // это коеффициенты для интерполяции
  dX = fX - iX;
  dY = fY - iY;
  dZ = fZ - iZ;

  // Здесь мы заполняем куб, точку в котором мы вычисляем
  fCV[0] = MyValues [iZ+4*iX+iY*4*4];
  fCV[1] = MyValues [iZ+4*(iX+1)+iY*4*4];
  fCV[2] = MyValues [iZ+4*(iX+1)+(iY+1)*4*4];
  fCV[3] = MyValues [iZ+4*iX+(iY+1)*4*4];
  fCV[4] = MyValues [(iZ+1)+4*iX+iY*4*4];
  fCV[5] = MyValues [(iZ+1)+4*(iX+1)+iY*4*4];
  fCV[6] = MyValues [(iZ+1)+4*(iX+1)+(iY+1)*4*4];
  fCV[7] = MyValues [(iZ+1)+4*iX+(iY+1)*4*4];

  return fBlend( fBlend( fBlend( fCV[0], fCV[1], dX ),
                         fBlend( fCV[3], fCV[2], dX ), dY ),
                 fBlend( fBlend( fCV[4], fCV[5], dX ),
                         fBlend( fCV[7], fCV[6], dX ), dY ), dZ );

}