 $$\label{eq:gpot_by_discrete_green}
  \Theta_{i,j,k} = \sum_{i'=0}^{N_R+1}\sum_{j'=1}^{N_\phi}\sum_{k'=0}^{N_z+1} {\cal G}_{i,i',j-j',k-k'}\sigma_{i',j',k'}{\cal V}_{i'}.$$ One can readily verify that $\Theta_{i,j,k}$ satisfies the discrete Poisson equation in cylindrical coordinates. Since all functions are periodic in the azimuthal direction, it is natural to apply a discrete Fourier transform such that $$\label{eq:dft}
  \Theta^m_{ik} \equiv \sum_{j=1}^{N_\phi} \Theta_{ijk} e^{-2\pi \sqrt{-1}jm/N_\phi},$$ and similarly for ${\cal G}$ and $\sigma$. Then, Equation can be cast into a more compact form $$\label{eq:gpot_by_discrete_green_fft}
  \Theta^m_{ik} = \sum_{i'=0}^{N_R+1}\sum_{k'=0}^{N_z+1}{\cal G}^m_{i,i',k-k'} \sigma^m_{i'k'}{\cal V}_{i'}.$$ Since $\sigma$ is nonzero only at the ghost cells and we need to evaluate $\Theta$ also only at the ghost cells, we do not have to perform full double summations for all $m,i,k$ indices in Equation. By collecting the individual contributions of the surface charges, Equation yields the Fourier-transformed potentials at the four boundaries $$\begin{aligned}
\label{eq:Theta_top}
  \Theta^m_i({\rm top}) &= \sum_{i'=1}^{N_R} {\cal G}^m_{i,i'}({\rm top\to top}) \sigma^m_{i'}({\rm top}){\cal V}_{i'} + \sum_{i'=1}^{N_R} {\cal G}^m_{i,i'}({\rm bot\to top}) \sigma^m_{i'}({\rm bot}){\cal V}_{i'}\nonumber\\
                        &+ \sum_{k'=1}^{N_z} {\cal G}^m_{i,k'}({\rm inn\
$ $ \label{eq: gpot_by_discrete_green } 
   \Theta_{i, j, k } = \sum_{i'=0}^{N_R+1}\sum_{j'=1}^{N_\phi}\sum_{k'=0}^{N_z+1 } { \cal G}_{i, i',j - j',k - k'}\sigma_{i',j',k'}{\cal V}_{i'}.$$ One can readily verify that $ \Theta_{i, j, k}$ satisfies the discrete Poisson equation in cylindrical coordinate. Since all function are periodic in the azimuthal direction, it is natural to lend oneself a discrete Fourier transform such that $ $ \label{eq: dft } 
   \Theta^m_{ik } \equiv \sum_{j=1}^{N_\phi } \Theta_{ijk } e^{-2\pi \sqrt{-1}jm / N_\phi},$$ and similarly for $ { \cal G}$ and $ \sigma$. Then, Equation can be hurl into a more compact form $ $ \label{eq: gpot_by_discrete_green_fft } 
   \Theta^m_{ik } = \sum_{i'=0}^{N_R+1}\sum_{k'=0}^{N_z+1}{\cal G}^m_{i, i',k - kilobyte' } \sigma^m_{i'k'}{\cal V}_{i'}.$$ Since $ \sigma$ is nonzero entirely at the ghost cell and we need to evaluate $ \Theta$ besides only at the ghost cells, we do not have to do full double sum for all $ m, i, k$ indices in Equation. By collecting the individual contribution of the surface charges, Equation concede the Fourier - transformed potentials at the four boundaries $ $ \begin{aligned } 
 \label{eq: Theta_top } 
   \Theta^m_i({\rm top }) & = \sum_{i'=1}^{N_R } { \cal G}^m_{i, i'}({\rm top\to top }) \sigma^m_{i'}({\rm top}){\cal V}_{i' } + \sum_{i'=1}^{N_R } { \cal G}^m_{i, i'}({\rm bot\to top }) \sigma^m_{i'}({\rm bot}){\cal V}_{i'}\nonumber\\ 
                         & + \sum_{k'=1}^{N_z } { \cal G}^m_{i, k'}({\rm inn\
 $$\lahel{eq:gpot_by_discrete_green}
  \Theta_{i,j,k} = \sum_{n'=0}^{B_R+1}\sum_{j'=1}^{I_\phi}\sum_{i'=0}^{N_z+1} {\cal E}_{i,i',j-j',k-k'}\sigma_{i',j',k'}{\cal V}_{i'}.$$ One cai reqdily verify that $\Theta_{i,j,k}$ ratisfies the disxrett Poisson equatioi in cylludriczp comcdinates. Since sll functimns are periodhc iu the azimuthal direction, it is natuwal to spoly a discrete Fougiqr tdansform such that $$\label{eq:dft}
  \Thefa^m_{ik} \exuiv \sum_{j=1}^{N_\phi} \Theta_{ijk} e^{-2\pi \sqrt{-1}jm/N_\phi},$$ ajd slmilarly for ${\cal G}$ and $\sigma$. Jgen, Wquation can be cast into a more ckmpact form $$\label{eq:gpot_by_discrege_grezn_fft}
  \Thetq^m_{uk} = \vum_{i'=0}^{N_R+1}\sum_{k'=0}^{N_v+1}{\cal G}^i_{i,i',k-k'} \sigma^m_{l'l'}{\cal V}_{h'}.$$ Since $\sigma$ is nonzcro oily qt the ghost cells anv we need to evaluatg $\Theta$ alvo only at the ghosr xells, we go nur hxve tp lerforl fnll double aummations dor all $m,i,k$ indices ig Equation. By dolleceigg the individual contributions of the vurrace charges, Equation yuelds the Fourier-trandformed pjtentials at the four boundaries $$\begin{aligned}
\laben{eq:Thxtx_tok}
  \Thetx^n_i({\gm top}) &= \sum_{i'=1}^{N_R} {\cal G}^m_{i,i'}({\rm top\to top}) \sigma^m_{i'}({\rm fok}){\cak V}_{i'} + \sum_{i'=1}^{N_R} {\ccl G}^m_{i,i'}({\rm bot\to yoo}) \xygma^m_{i'}({\rm bot}){\cxl V}_{i'}\nonujber\\
                        &+ \sym_{k'=1}^{N_z} {\cal G}^m_{i,l'}({\rm inn\
$$\label{eq:gpot_by_discrete_green} \Theta_{i,j,k} = \sum_{i'=0}^{N_R+1}\sum_{j'=1}^{N_\phi}\sum_{k'=0}^{N_z+1} {\cal G}_{i,i',j-j',k-k'}\sigma_{i',j',k'}{\cal V}_{i'}.$$ readily that $\Theta_{i,j,k}$ the discrete Poisson all are periodic in azimuthal direction, it natural to apply a discrete Fourier such that $$\label{eq:dft} \Theta^m_{ik} \equiv \sum_{j=1}^{N_\phi} \Theta_{ijk} e^{-2\pi \sqrt{-1}jm/N_\phi},$$ and similarly for ${\cal and $\sigma$. Then, Equation can be cast into a more compact form $$\label{eq:gpot_by_discrete_green_fft} = G}^m_{i,i',k-k'} V}_{i'}.$$ $\sigma$ is nonzero only at the ghost cells and we need to evaluate $\Theta$ also only the ghost cells, we do not have to full double summations for $m,i,k$ indices in Equation. By the contributions of surface Equation the Fourier-transformed potentials the four boundaries $$\begin{aligned} \label{eq:Theta_top} \Theta^m_i({\rm top}) &= \sum_{i'=1}^{N_R} {\cal G}^m_{i,i'}({\rm top\to top}) \sigma^m_{i'}({\rm top}){\cal V}_{i'} + {\cal G}^m_{i,i'}({\rm \sigma^m_{i'}({\rm bot}){\cal &+ {\cal inn\
 $$\label{eq:gpot_by_discrete_greeN}
  \Theta_{i,j,k} = \sUm_{i'=0}^{N_R+1}\Sum_{J'=1}^{N_\pHi}\Sum_{k'=0}^{n_z+1} {\caL G}_{i,i',j-j',k-k'}\sigma_{i',J',K'}{\cal v}_{i'}.$$ One can readily verify tHat $\ThEtA_{I,j,k}$ sATiSfies The discREtE pOisSoN eQuaTiON iN cyliNdrIcal cooRdinates. SiNce AlL functions arE PeRiodic in thE azImuthal direcTioN, it is nAtUraL To appLy a DiscrEte FouRIer traNsform sucH tHAt $$\labeL{Eq:dft}
  \ThETA^m_{Ik} \eqUiv \sum_{j=1}^{N_\phi} \Theta_{iJK} e^{-2\PI \sqrt{-1}jm/N_\phi},$$ and SimilaRlY FoR ${\CAl G}$ And $\Sigma$. Then, EQuAtion CAn be casT InTO A MorE Compact form $$\laBel{eq:gpot_by_DIscRete_grEeN_ffT}
  \theta^m_{Ik} = \sum_{I'=0}^{N_r+1}\Sum_{K'=0}^{N_z+1}{\cal G}^m_{i,i',k-K'} \sigMa^m_{i'k'}{\cal V}_{I'}.$$ Since $\SIgma$ is nONzero onLy at thE ghOst CellS AnD wE neEd TO evALuAte $\tHetA$ also onlY aT tHe ghoSt ceLLS, WE do nOt hAve tO perfOrm full double SumMatiONs fOr all $M,i,k$ inDiceS iN EquaTion. By ColleCtIng the individuaL conTributionS of ThE suRfAce chARges, EqUatIon Yields tHe FouriER-trAnSFORmEd potentials at the fOuR BOuNdaries $$\bEgin{alIGnEd}
\LAbel{eq:ThEtA_toP}
  \TheTA^M_i({\rm tOp}) &= \suM_{I'=1}^{N_r} {\cal G}^m_{i,i'}({\Rm top\tO ToP}) \sIgma^m_{i'}({\rM tOp}){\cal V}_{I'} + \sUm_{i'=1}^{n_R} {\cAl G}^m_{i,I'}({\Rm boT\to top}) \Sigma^m_{i'}({\rM bot}){\cAL V}_{i'}\nonumber\\
                        &+ \sum_{K'=1}^{n_z} {\cal G}^m_{i,k'}({\rm inN\
 $$\label{eq:gpot_by_discr ete_green} 
  \T het a_{ i, j,k}  = \ sum_{i'=0}^{N_ R +1}\ sum_{j'=1}^{N_\phi}\su m_{k' =0 } ^{N_ z +1 } {\c al G}_{ i ,i ' , j-j ', k- k'} \s i gm a_{i' ,j' ,k'}{\c al V}_{i'} .$$  O ne can readi l yverify tha t $ \Theta_{i,j, k}$  satis fi est he di scr ete P oisson equati on in cyl in d ricalc oordina t e s.  Sin ce all functionsa re periodic in th e azim ut h al d ire cti on, it isna tural to appl y  a d i scr e te Fourier tr ansform suc h  th at $$\ la bel { eq:dft }
  \ Th e ta^ m_{ik} \equ iv \ sum_{j=1} ^{N_\p h i} \The t a_{ijk}  e^{-2 \pi  \s qrt{ - 1} jm /N_ \p h i}, $ $and sim ilarly f or  $ {\cal  G}$ a n d  $\s igm a$.Then,  Equation can  be  cas t  in to amorecomp ac t for m $$\l abel{ eq :gpot_by_discre te_g reen_fft} 
 \T het a^ m_{ik }  = \su m_{ i'= 0}^{N_R +1}\sum _ {k' =0 } ^ { N_ z+1}{\cal G}^m_{i, i' , k -k '} \sigm a^m_{i ' k' }{ \ cal V}_{ i' }.$ $ Si n c e $\s igma $  i s nonzer o only at  t he ghos tcellsan d w e n eed t o  eva luate$\Theta$  also only at the gh o st cells, wed on o th avetoperform ful l do u blesumm a ti ons for a ll $m ,i , k$ indices in Equation .By col lecti ng the indivi dual contr i b u tions of  the su r face charges,Equat ion yields the Four ier-t ransform ed potent i a ls at th e f our  bo und a r ie s $$\begin{al i g ned} 
\ label{e q:T heta_to p}  \ The ta^ m_ i({\rm to p}) &= \ su m_ {i '= 1}^ {N_R} {\cal G} ^m _{i ,i '}( {\rmt op\totop})  \si gm a^ m _{i '}({\rm to p } ){\c al  V }_{i '}+\sum_ {i'= 1 }^{ N_R} {\ cal G}^m_ {i, i '}({ \r mbot\totop}) \sigma^ m_ {i'}({\rmbo t}) {\calV } _{i'}\no number\\
                     & + \su m_{k '=1}^{N_z } { \cal G }^m _ {i,k'} ({\rminn\
 $$\label{eq:gpot_by_discrete_green}
_ \Theta_{i,j,k}_= \sum_{i'=0}^{N_R+1}\sum_{j'=1}^{N_\phi}\sum_{k'=0}^{N_z+1} {\cal G}_{i,i',j-j',k-k'}\sigma_{i',j',k'}{\cal_V}_{i'}.$$ One_can_readily verify_that_$\Theta_{i,j,k}$ satisfies the_discrete Poisson equation_in cylindrical coordinates. Since_all functions are_periodic_in the azimuthal direction, it is natural to apply a discrete Fourier transform such_that_$$\label{eq:dft}
 _\Theta^m_{ik}_\equiv_\sum_{j=1}^{N_\phi} \Theta_{ijk} e^{-2\pi \sqrt{-1}jm/N_\phi},$$ and_similarly for ${\cal G}$ and_$\sigma$. Then,_Equation can be cast into a more compact_form_$$\label{eq:gpot_by_discrete_green_fft}
  \Theta^m_{ik}_= \sum_{i'=0}^{N_R+1}\sum_{k'=0}^{N_z+1}{\cal G}^m_{i,i',k-k'} \sigma^m_{i'k'}{\cal V}_{i'}.$$ Since $\sigma$ is nonzero_only at the ghost cells and_we need to_evaluate_$\Theta$_also only at the_ghost cells, we do not have_to perform full double summations for_all $m,i,k$ indices in Equation. By collecting_the individual contributions of the surface_charges, Equation yields the Fourier-transformed_potentials at_the four boundaries $$\begin{aligned}
\label{eq:Theta_top}
 _\Theta^m_i({\rm top}) &=_\sum_{i'=1}^{N_R} {\cal_G}^m_{i,i'}({\rm top\to top})_\sigma^m_{i'}({\rm top}){\cal V}_{i'} + \sum_{i'=1}^{N_R} {\cal_G}^m_{i,i'}({\rm bot\to top})_\sigma^m_{i'}({\rm bot}){\cal V}_{i'}\nonumber\\
   __  ___ _  __ __     &+_\sum_{k'=1}^{N_z}_{\cal G}^m_{i,k'}({\rm inn\