 #define    epsilon0     8.8542e-12
 
 __global__ void spUpdateFieldYeeKernel(Real dt, Real3 dt_inv,
-                                       dim3 N, dim3 I, int rank,
+                                       dim3 N, dim3 I,
                                        Real const *Rho,
                                        Real const *Jx,
                                        Real const *Jy,

                B[0], B[1], B[2]
    );

//    if (spMPIRank() == 0)
//    {
//        Real tmp[64];
//        spParallelMemcpy(tmp, B[1], 64 * sizeof(Real));
//
//        for (int i = 0; i < 8; ++i)
//        {
//            for (int j = 0; j < 8; ++j)
//            {
//                printf("%d \t", (int) tmp[i * 8 + j]);
//            }
//            printf("\n");
//
//        }
//        printf("---------------------------\n");
//    }
    spFieldSync(fE);
    spFieldSync(fB);
//    if (spMPIRank() == 0)
//    {
//        Real tmp[64];
//        spParallelMemcpy(tmp, B[1], 64 * sizeof(Real));
//
//        for (int i = 0; i < 8; ++i)
//        {
//            for (int j = 0; j < 8; ++j)
//            {
//                printf("%d \t", (int) tmp[i * 8 + j]);
//            }
//            printf("\n");
//        }
//        printf("---------------------------\n");
//    }
    return SP_SUCCESS;
}

//
//#include <assert.h>
//#include <math.h>
//
//extern "C" {
//#include </usr/local/cuda/include/cuda_runtime.h>
//#include </usr/local/cuda/include/device_launch_parameters.h>
//#include </usr/local/cuda/include/device_functions_decls.h>
//
//#include "sp_lite_def.h"
//
//#include "spParallel.h"
//
//#include "spParticle.h"
//#include "BorisYee.h"
//
//}
//
//__global__ void
//spBorisInitializeParticleKernel(boris_particle *d, spParticleFiber **bucket, size_type PIC)
//{

    spFieldSubArray(fB, (void **) B);

    LOAD_KERNEL(spUpdateFieldYeeKernel, sizeType2Dim3(dims), 1,
                dt, real2Real3(dt_inv), sizeType2Dim3(dims), sizeType2Dim3(strides),
                (const Real *) rho,
                (const Real *) J[0],
                (const Real *) J[1],
