 {
     // A material constant
     int Eta = 13136;    // 13136 liver
-    cudaMemcpyToSymbol("Eta_gpu", &Eta, sizeof(int));
+    cudaMemcpyToSymbol("Eta_hex_gpu", &Eta, sizeof(int));
 
     // The structure tensor (a defines the preferred material direction)
     float a[3] = {0, 0.707f, 0.707f};


}

__global__ void CudaHexahedronTLEDForceField3f_calcForce_kernel2(float Lambda, float Mu, int nbElem, float4 * Di1, float4 * Di2, float4 * Dv1, float4 * Dv2, float4* F0_gpu, float4* F1_gpu, float4* F2_gpu, float4* F3_gpu, float4* F4_gpu, float4* F5_gpu, float4* F6_gpu, float4* F7_gpu)
{
    int index0 = umul24(blockIdx.x,BSIZE); //blockDim.x;
    int index1 = threadIdx.x;


/// GPU pointers
// List of nodes for each element
int4* NodesPerElement_hex_gpu = 0;
// Shape function derivatives arrays
float4* DhC0_hex_gpu = 0;
float4* DhC1_hex_gpu = 0;
float4* DhC2_hex_gpu = 0;
// Hourglass control
float4* HG_hex_gpu = 0;
// Force coordinates for each node
int2* FCrds_hex_gpu = 0;
// Jacobian determinant array
float* DetJ_hex_gpu;
// Element nodal force contribution
float4* F0_hex_gpu = 0;
float4* F1_hex_gpu = 0;
float4* F2_hex_gpu = 0;
float4* F3_hex_gpu = 0;
float4* F4_hex_gpu = 0;
float4* F5_hex_gpu = 0;
float4* F6_hex_gpu = 0;
float4* F7_hex_gpu = 0;
// Displacements
float4* Disp = 0;

// Viscoelasticity
float4 * Di1_hex_gpu = 0;
float4 * Di2_hex_gpu = 0;
float4 * Dv1_hex_gpu = 0;
float4 * Dv2_hex_gpu = 0;

/// CPU pointers
float* test = 0;
