  We construct finite time blow-up solutions to the Landau-Lifshitz-Gilbert
equation (LLG) from ${\mathbb R}^2$ into $S^2$ \begin{equation*} \begin{cases}
u_t= a(\Delta u+|\nabla u|^2u) -b u\wedge \Delta u &\ \mbox{ in }\ {\mathbb
R}^2\times(0,T), u(\cdot,0) = u_0\in S^2 &\ \mbox{ in }\ {\mathbb R}^2,
\end{cases} \end{equation*} where $a^2+b^2=1,~a > 0,~ b\in {\mathbb R}$. Given
any prescribed $N$ points in $\mathbb{R}^2$ and small $T>0$, we prove that
there exists regular initial data such that the solution blows up precisely at
these points at finite time $t=T$, taking around each point the profile of
sharply scaled degree 1 harmonic map with the type II blow-up speed
\begin{equation*} \| \nabla u\|_{L^\infty } \sim \frac{|\ln(T-t)|^2}{ T-t } \
\mbox{ as } \ t\to T. \end{equation*} The proof is based on the {\em parabolic
inner-outer gluing method}, developed in \cite{17HMF} for Harmonic Map Flow
(HMF). However, a direct consequence of the presence of dispersion is the {\em
lack of maximum principle} for suitable quantities, which makes the analysis
more delicate even at the linearized level. To overcome this difficulty, we
make use of two key technical ingredients: first, for the inner problem we
employ the tool of {\em distorted Fourier transform}, as developed by Krieger,
Miao, Schlag and Tataru \cite{Krieger09Duke,KMS20WM}. Second, the linear theory
for the outer problem is achieved by means of the sub-Gaussian estimate for the
fundamental solution of parabolic system in non-divergence form with
coefficients of Dini mean oscillation in space ($\mathsf{DMO_x}$), which was
proved by Dong, Kim and Lee \cite{dong22-non-divergence} recently.
