Yinjian Zhao

@RemiLehe Hi Remi, thanks for your detailed reply! I have a few more questions. 1. Within the above defined function `phi_fix()`, when I use domain decomposition, each MPI rank will...

Well, I made a mistake I guess, this may be better. ``` #elif defined(WARPX_DIM_XZ) || defined(WARPX_DIM_RZ) int rank,nprocs,n,nc,ns; MPI_Comm_size(MPI_COMM_WORLD,&nprocs); MPI_Comm_rank(MPI_COMM_WORLD,&rank); n = 256/nprocs; nc = 0; Real phi0 = 0._rt,...

> Hi Yin - For you question 1, when you do `phi = fields.PhiFPWrapper(0, True)[...]`, the phi that is returned is a full sized array covering the full domain and...

I will close this issue, because I have had a working scipt already. Thanks for all the help!

So far, this version of the PR still has some issues, the results do not match that well. Smilei's plot: Ours: