Как векторизовать петлю ARM Neon

У меня есть цикл, как указано ниже в классе C ++, который выполняет разложение по Холески, где A, B, C, D и E являются целочисленными массивами.

int i = 10;
for (int j = A [i] + 1 ; j < C [i] ; j++)
{
E[D [j]] -= B[j] * val ;
}

На каждой итерации один элемент массива обновляется и накапливается в другой индекс E. Как я могу эффективно векторизовать этот цикл для оптимизации с использованием встроенных функций ARM Neon?