Вариация функционала по вектору

[L-sky]

Добрый день.
При рассмотрении уравнения Ландау-Лифшица столкнулся с выражением: $\frac{\delta E}{\delta\vec{M}}$, в котором происходит вариация энергии (некоторый известный интегральный функционал) по вектору намагниченности (3-мерному), причем $E$ не содержит явно $\vec{M}$, а зависит от угловых координат $(\theta, \varphi)$. $\vec{M}$ также зависит от $(\theta, \varphi)$.
То есть, имеем выражение вида:

$$\frac{\delta E(\theta,\varphi)}{\delta \vec{M}(\theta,\varphi)}$$

В примерах, которые удалось найти, выполнение вариации рассматриваются только по скалярным функциям. Можете подсказать, как действовать в приведенном выше случае?