oiDelta\]. We will calculate $$\frac{\partial V_k}{\partial x_i} =
\lim_{\delta_x\to 0}\frac{\Delta V}{\delta_x},\label{eq:DerivadaLimite}$$ where $\Delta V_k = V_k(\vec{x}_1, \cdots,
\vec{x}_i + \vec{\delta_x}, \cdots, \vec{x}_{N_V}) -
V_k(\vec{x}_1, \cdots, \vec{x}_i, \cdots, \vec{x}_{N_V})$.

It can be seen in Figure \[fig:VoronoiDelta\] that when moving the generator $\vec{x}_i$, the only polygons modified are $\mathcal{V}_i$ and its neighbors. Using this, Eq. \[eq:VisDiscreta\] leads to $$\begin{aligned}
\Delta V_k & = & I_i'\int_{\mathcal{V}_i'}f_k(\vec{x})d\vec{x} -
I_i\int_{\mathcal{V}_i}f_k(\vec{x})d\vec{x}\nonumber\\
& & + \sum_{j\in J_i}\bigg(I_j'\int_{\mathcal{V}_j'}f_k(\vec{x})d\vec{x} -
I_j\int_{\mathcal{V}_j}f_k(\vec{x})d\vec{x}\bigg),\label{eq:deltaVis}\end{aligned}$$ where $\mathcal{V}_i$ is the polygon generated by $\vec{x}_i$ before moving, $\mathcal{V}_i'$ is the same polygon after moving $\vec{x}_i$, $J_i$ is the set of indices of the polygons that are neighbors to $\mathcal{V}_i$ and $J_i'$ is the set of indices of the polygons that are neighbors to $\mathcal{V}_i'$.

It can be seen in Figure \[fig:VoronoiDelta\] that $$\begin{aligned}
\mathcal{V}_i = (\mathcal{V}_i\cap \mathcal{V}_i')\cup(\mathcal{V}_i\setminus \mathcal{V}_i\cap \mathcal{V}_i'), & \mathcal
oiDelta\ ]. We will calculate $ $ \frac{\partial V_k}{\partial x_i } = 
 \lim_{\delta_x\to 0}\frac{\Delta V}{\delta_x},\label{eq: DerivadaLimite}$$ where $ \Delta V_k = V_k(\vec{x}_1, \cdots, 
 \vec{x}_i + \vec{\delta_x }, \cdots, \vec{x}_{N_V }) - 
 V_k(\vec{x}_1, \cdots, \vec{x}_i, \cdots, \vec{x}_{N_V})$. 

 It can be seen in Figure \[fig: VoronoiDelta\ ] that when moving the generator $ \vec{x}_i$, the only polygon change are $ \mathcal{V}_i$ and its neighbors. Using this, Eq. \[eq: VisDiscreta\ ] lead to $ $ \begin{aligned } 
 \Delta V_k & = & I_i'\int_{\mathcal{V}_i'}f_k(\vec{x})d\vec{x } - 
 I_i\int_{\mathcal{V}_i}f_k(\vec{x})d\vec{x}\nonumber\\ 
 & & + \sum_{j\in J_i}\bigg(I_j'\int_{\mathcal{V}_j'}f_k(\vec{x})d\vec{x } - 
 I_j\int_{\mathcal{V}_j}f_k(\vec{x})d\vec{x}\bigg),\label{eq: deltaVis}\end{aligned}$$ where $ \mathcal{V}_i$ is the polygon render by $ \vec{x}_i$ before moving, $ \mathcal{V}_i'$ is the same polygon after move $ \vec{x}_i$, $ J_i$ is the bent of indices of the polygons that are neighbors to $ \mathcal{V}_i$ and $ J_i'$ is the set of index of the polygons that are neighbors to $ \mathcal{V}_i'$. 

 It can be see in Figure \[fig: VoronoiDelta\ ] that $ $ \begin{aligned } 
 \mathcal{V}_i = (\mathcal{V}_i\cap \mathcal{V}_i')\cup(\mathcal{V}_i\setminus \mathcal{V}_i\cap \mathcal{V}_i'), & \mathcal
oiDflta\]. We will calculate $$\fvac{\partial V_k}{\parjiql x_i} =
\nim_{\delfa_x\to 0}\frxc{\Delta V}{\delta_x},\label{eq:DerivavaLinite}$$ qhere $\Delta V_k = V_k(\vec{x}_1, \cdots,
\vec{q}_i + \vec{\deota_x}, \xdots, \vec{x}_{I_B}) -
V_k(\vec{w}_1, \cdofd, \vee{x}_m, \cdots, \vec{x}_{N_V})$.

Ij can be seet in Figure \[fic:VurlnoiDelta\] that when moving the genewator $\vrc{d}_i$, the only pojygomf mospfled are $\mathcal{V}_i$ and its neighbkrs. Usiig this, Eq. \[eq:VixDiscreta\] leads to $$\begin{allgnef}
\Delta V_k & = & I_i'\int_{\lathcal{V}_i'}f_k(\cec{x})q\cec{x} -
I_i\int_{\maghcal{V}_i}f_k(\vtc{r})d\vec{x}\nonumger\\
& & + \sum_{j\in J_i}\bigg(I_j'\int_{\mathcal{X}_j'}f_k(\vzc{x})d\vec{x} -
I_j\unr_{\mahvcal{V}_j}f_k(\vec{e})d\vec{x}\figg),\label{eq:dcktaVis}\and{aligmed}$$ where $\mathgal{V}_i$ is the polygon generatev by $\vec{x}_i$ before moding, $\mathwam{V}_i'$ is the same piltgon dftes moxung $\ved{x}_m$, $J_j$ is tje aet of indjces of the polygons that are meythbors to $\matgcal{V}_i$ agd $J_i'$ is the set of indices of the polycona that are neighbors to $\mathcal{V}_i'$.

It can be sgen in Figtre \[fig:VoronoiDelta\] that $$\begin{aligned}
\mathcal{V}_i = (\mdthcam{X}_i\cci \oqtjcal{V}_i')\cup(\mathcal{V}_i\setminus \mathcal{V}_i\cap \mathcwm{V}_o'), & \mathcal
oiDelta\]. We will calculate $$\frac{\partial V_k}{\partial x_i} 0}\frac{\Delta where $\Delta = V_k(\vec{x}_1, \cdots, - \cdots, \vec{x}_i, \cdots, It can be in Figure \[fig:VoronoiDelta\] that when moving generator $\vec{x}_i$, the only polygons modified are $\mathcal{V}_i$ and its neighbors. Using this, \[eq:VisDiscreta\] leads to $$\begin{aligned} \Delta V_k & = & I_i'\int_{\mathcal{V}_i'}f_k(\vec{x})d\vec{x} - I_i\int_{\mathcal{V}_i}f_k(\vec{x})d\vec{x}\nonumber\\ & + J_i}\bigg(I_j'\int_{\mathcal{V}_j'}f_k(\vec{x})d\vec{x} I_j\int_{\mathcal{V}_j}f_k(\vec{x})d\vec{x}\bigg),\label{eq:deltaVis}\end{aligned}$$ $\mathcal{V}_i$ is the polygon generated by $\vec{x}_i$ before moving, $\mathcal{V}_i'$ is the same polygon after moving $J_i$ is the set of indices of the that are neighbors to and $J_i'$ is the set indices the polygons are to It can be in Figure \[fig:VoronoiDelta\] that $$\begin{aligned} \mathcal{V}_i = (\mathcal{V}_i\cap \mathcal{V}_i')\cup(\mathcal{V}_i\setminus \mathcal{V}_i\cap \mathcal{V}_i'), & \mathcal
oiDelta\]. We will calculate $$\fraC{\partial V_k}{\PartiAl x_I} =
\liM_{\dElta_X\to 0}\fRac{\Delta V}{\delta_X},\LabeL{eq:DerivadaLimite}$$ where $\delta v_k = v_K(\vec{X}_1, \CdOts,
\veC{x}_i + \vec{\dELtA_X}, \CdoTs, \VeC{x}_{N_v}) -
V_K(\VeC{x}_1, \cdoTs, \vEc{x}_i, \cdoTs, \vec{x}_{N_V})$.

It Can Be Seen in Figure \[FIg:voronoiDelTa\] tHat when movinG thE generAtOr $\vEC{x}_i$, thE onLy polYgons mODified Are $\mathcaL{V}_I$ And its NEighborS. uSiNg thIs, Eq. \[eq:VisDiscreta\] LEaDS to $$\begin{aligneD}
\Delta v_k & = & i_I'\iNT_{\MatHcaL{V}_i'}f_k(\vec{x})d\VeC{x} -
I_i\iNT_{\mathcaL{v}_i}F_K(\VEc{x})D\Vec{x}\nonumber\\
& & + \sUm_{j\in J_i}\bigg(i_J'\inT_{\mathcAl{v}_j'}f_K(\Vec{x})d\vEc{x} -
I_j\InT_{\MatHcal{V}_j}f_k(\vec{X})d\veC{x}\bigg),\labEl{eq:deLTaVis}\enD{Aligned}$$ Where $\mAthCal{v}_i$ is THe PoLygOn GEneRAtEd bY $\Vec{X}_i$ before MoViNg, $\matHcal{v}_I'$ IS The sAme PolyGon afTer moving $\vec{x}_I$, $J_i$ Is thE Set Of indIces oF the PoLygonS that aRe neiGhBors to $\mathcal{V}_i$ And $J_I'$ is the set Of iNdIceS oF the pOLygons ThaT arE neighbOrs to $\maTHcaL{V}_I'$.

iT CaN be seen in Figure \[fig:voRONoIDelta\] thAt $$\begiN{AlIgNEd}
\mathcaL{V}_I = (\maThcaL{v}_I\cap \mAthcAL{V}_I')\cup(\mathCal{V}_i\sETmInUs \mathcAl{v}_i\cap \mAtHcaL{V}_i'), & \MathcAL
oiDelta\]. We will calcula te $$\frac {\par tia l V _k }{\p arti al x_i} =
\lim _ {\de lta_x\to 0}\frac{\Delt a V}{ \d e lta_ x }, \labe l{eq:De r iv a d aLi mi te }$$  w h er e $\D elt a V_k =  V_k(\vec{ x}_ 1,  \cdots,
\ve c {x }_i + \vec {\d elta_x}, \cd ots , \vec {x }_{ N _V})-
V _k(\v ec{x}_ 1 , \cdo ts, \vec{ x} _ i, \cd o ts, \ve c { x} _{N_ V})$.

It can bes ee n  in Figure \[f ig:Vor on o iD e l ta\ ] t hat when m ov ing t h e gener a to r $ \ve c {x}_i$, the o nly polygon s  mo dified  a re$ \mathc al{V} _i $  an d its neigh bors . Using t his, E q . \[eq: V isDiscr eta\]lea dsto $ $ \b eg in{ al i gne d }\De l taV_k & =&I_ i'\in t_{\ m a t h cal{ V}_ i'}f _k(\v ec{x})d\vec{x } - 
I_i \ int _{\ma thcal {V}_ i} f_k(\ vec{x} )d\ve c{ x}\nonumber\\
&  & +  \sum_{j\ inJ_ i}\ bi gg(I_ j '\int_ {\m ath cal{V}_ j'}f_k( \ vec {x } ) d \v ec{x} -
I_j\int_{\ ma t h ca l{V}_j}f _k(\ve c {x }) d \vec{x}\ bi gg) ,\la b e l{eq: delt a Vi s}\end{a ligned } $$  w here $\ ma thcal{ V} _i$  is  thep olyg on gen erated b y $\v e c{x}_i$ before moving, $\mat h ca l { V} _ i'$isthe same po lygo n  aft er m o vi ng$ \vec{ x}_i$ ,$ J_ i $ is the set of ind ic es ofthe p olygons thatare neighb o r s  to $\ma thca l {V } _i$ and $J_i'$  is t he set ofi ndices o f the  polygon s that ar e neighbor s t o $ \ma thc a l {V }_i'$.

It ca n be s ee n in Fi gur e \[fig :Vo ron oiD elt a\ ] that $$ \begin{a li gn ed }\ma thcal { V}_i = ( \m ath ca l{V }_i\c a p \mat hcal{ V}_i ') \c u p(\ mathcal { V} _ i \set mi nu s \m ath ca l{V}_ i\ca p  \m athcal{ V}_i'), &  \m a thca l
oiDelta\]. We_will calculate_$$\frac{\partial V_k}{\partial x_i} =
\lim_{\delta_x\to_0}\frac{\Delta V}{\delta_x},\label{eq:DerivadaLimite}$$_where_$\Delta V_k_=_V_k(\vec{x}_1, \cdots,
\vec{x}_i +_\vec{\delta_x}, \cdots, \vec{x}_{N_V})_-
V_k(\vec{x}_1, \cdots, \vec{x}_i, \cdots,_\vec{x}_{N_V})$.

It can be_seen_in Figure \[fig:VoronoiDelta\] that when moving the generator $\vec{x}_i$, the only polygons modified are_$\mathcal{V}_i$_and its_neighbors._Using_this, Eq. \[eq:VisDiscreta\] leads to_$$\begin{aligned}
\Delta V_k & = &_I_i'\int_{\mathcal{V}_i'}f_k(\vec{x})d\vec{x} -
I_i\int_{\mathcal{V}_i}f_k(\vec{x})d\vec{x}\nonumber\\
&_& + \sum_{j\in J_i}\bigg(I_j'\int_{\mathcal{V}_j'}f_k(\vec{x})d\vec{x} -
I_j\int_{\mathcal{V}_j}f_k(\vec{x})d\vec{x}\bigg),\label{eq:deltaVis}\end{aligned}$$ where $\mathcal{V}_i$ is_the_polygon generated by_$\vec{x}_i$ before moving, $\mathcal{V}_i'$ is the same polygon after_moving $\vec{x}_i$, $J_i$ is the set_of indices of_the_polygons_that are neighbors to_$\mathcal{V}_i$ and $J_i'$ is the set_of indices of the polygons that_are neighbors to $\mathcal{V}_i'$.

It can be seen_in Figure \[fig:VoronoiDelta\] that $$\begin{aligned}
\mathcal{V}_i =_(\mathcal{V}_i\cap \mathcal{V}_i')\cup(\mathcal{V}_i\setminus \mathcal{V}_i\cap \mathcal{V}_i'), &_\mathcal