Трехмерный массив в одномерный

[КРИЖ]

Надо было сделать трехмерную модель Изинга, с примесями, но так как там приходиться бегать по 3 вложеным циклам, а тут ещё и вдобавок бегать по ячейкам трехмерого массива где хранятся значения спина, то захотелось сделать массив одномерный. И появилась такая вот проблема:в трехмерном массиве двигаюсь по нему таким вот образом

for(a=0;a<L*L*L;++a)//пробегаем по решётки и ищем случайный спин
{
	j = Mersenne.IRandomX(0,L-1);
	i = Mersenne.IRandomX(0,L-1);
	k = Mersenne.IRandomX(0,L-1);

	r = Mersenne.Random();

	if(sp[i+L*j+L*k])
	{
		if(r<0.5)
		{
			sp_new=sp[i+L*j+L*k];
		}
		else
		{
			sp_new=-sp[i+L*j+L*k];
			dE=(double)2*sp[i+L*j+L*k]*Sosedu(sp,i,j,k);
							
			//тепловая баня
			if(dE<=0.0)//if((dE<=0.0)||(dE<=-0.0))
			{
				//sp[i+L*j+L*k]=sp_new;
				r = Mersenne.Random();
				if(w(dE,T)>r)
				{
					sp[i+L*j+L*k]=sp_new;
				}
			}
			else
			{
				r = Mersenne.Random();
				if(w(dE,T)>r)
				{
					sp[i+L*j+L*k]=sp_new;
				}
			}
		}
	}
}//закончили пробег по случайному спину

и как потом выяснялось что не для всех j,k это справедливо, например для j=0 или k=0, мы попадаем в одни и тот же участок памяти, вопрос как можно двигать по нему, но именно с индексами i,j,k. Знаю что можно обойтись одним, но генерировать надо именно i,j,k так как есть функция которая считает по ближайщим спинам.

[8Observer8]

Теперь дайте полное описание задачи. Может вы изначально неправильно решаете