     float *p = (float*)malloc(nn*sizeof(float));
     float *u = (float*)malloc(nn*sizeof(float));
 
+    // Pré-condicionamento.
     for (i = 0; i < nn; i++) {
         ri = nodes[i].calc ? rsum[i] - dsum[i]*V[i] : 0.0;
         p[i] = ri;

    P[i] = ri;
}

__global__ void kernel_iter_element(int ne) { }

// Função externa que processa o problema na GPU.
//    ne: número de elementos.

                      elementri *elements, node *nodes, float *V, bool verbose,
                      float *bench) {
    int i, k;

    // Pre-processamento. Calcula as matrizes de contribuição dos elementos.
    float J1, J2, J3, J4, dJ;
