     int haveHelper = (numObjectsH->numHelperMarkers > 0) ? 1 : 0;
 
     int4 updatePortion = mI4(fsiSystem.fsiGeneralData->referenceArray[haveHelper].x,
-                             fsiSystem.fsiGeneralData->referenceArray[haveHelper + haveGhost].y, 0, 0);
+        fsiSystem.fsiGeneralData->referenceArray[haveHelper + haveGhost].y, 0, 0);
 
     cout << "time step in UpdateFluid_Implicit " << paramsH->dT << endl;
     bool *isErrorH, *isErrorD;

// Author: Milad Rakhsha, Arman Pazouki, Wei Hu
// =============================================================================
//
// Class for performing time integration in fluid system.//
// =============================================================================

#include "chrono_fsi/physics/ChFluidDynamics.cuh"

    uint numBlocks, numThreads;
    computeGridSize(updatePortion.y - updatePortion.x, 256, numBlocks, numThreads);
    UpdateActivityD<<<numBlocks, numThreads>>>(
        mR4CAST(sphMarkersD2->posRadD), mR3CAST(sphMarkersD1->velMasD), 
        mR3CAST(fsiBodiesD->posRigid_fsiBodies_D),
        U1CAST(fsiSystem.fsiGeneralData->activityIdentifierD), 
        U1CAST(fsiSystem.fsiGeneralData->extendedActivityIdD),
        updatePortion, numObjectsH->numRigidBodies, Time, isErrorD);
    cudaDeviceSynchronize();
    cudaCheckError();
    //------------------------

    cudaMemcpy(isErrorH, isErrorD, sizeof(bool), cudaMemcpyDeviceToHost);
    if (*isErrorH == true)
        throw std::runtime_error("Error! program crashed in UpdateActivityD!\n");
    cudaFree(isErrorD);
    free(isErrorH);
}
