         kUpdateParticleVelocity<<<CONFIG(host_settings.num_mpm_markers)>>>(
             grid_vel.data_d, old_vel_node_mpm.data_d, pos.data_d, vel.data_d, marker_Fe.data_d, marker_Fp.data_d);
     }
-    printf("kIncrementVelocity: %f\n", time_measured);
+    printf("kUpdateParticleVelocity: %f\n", time_measured);
     time_measured = 0;
 
     vel.copyDeviceToHost();

        PolarR[p] = RE;
        PolarR[p][3] = JP;

        Mat33 PED = real(2.) * current_mu * (FE_hat - RE) +
                    current_lambda * JE_hat * (JE_hat - real(1.)) * AdjointTranspose(FE_hat) * real(1.0 / JE_hat);

        Mat33 vPEDFepT = marker_volume[p] * MultTranspose(PED, FE);

        LOOP_TWO_RING_GPU(                                                  //
            real3 d_weight = dN(xi - current_node_location, inv_bin_edge);  //
            real3 force = device_settings.dt * (vPEDFepT * d_weight) / (JE * JP);

            real mass = node_mass[current_node];  //
            if (mass > 0) {

        real3 V_flip = vel_marker[p];
        real3 V_pic = real3(0.0);

        const real bin_edge = device_settings.bin_edge;
        const real inv_bin_edge = 1.f / bin_edge;
        int cx, cy, cz;

        LOOP_TWO_RING_GPUSP(

            real weight = N((xi.x - current_node_locationx) * inv_bin_edge) *
                          N((xi.y - current_node_locationy) * inv_bin_edge) *
                          N((xi.z - current_node_locationz) * inv_bin_edge);

            real vnx = grid_vel[current_node * 3 + 0];  //
            real vny = grid_vel[current_node * 3 + 1];  //
            real vnz = grid_vel[current_node * 3 + 2];

            V_pic.x += vnx * weight;                                              //
            V_pic.y += vny * weight;                                              //
            V_pic.z += vnz * weight;                                              //
            V_flip.x += (vnx - old_vel_node_mpm[current_node * 3 + 0]) * weight;  //
            V_flip.y += (vny - old_vel_node_mpm[current_node * 3 + 1]) * weight;  //
            V_flip.z += (vnz - old_vel_node_mpm[current_node * 3 + 2]) * weight;  //
            )
        real3 new_vel = (1.0 - alpha) * V_pic + alpha * V_flip;

