#pragma CODE_SECTION(iirLatticeLadderStereoProcess, ".criticalSectionInternal")
void iirLatticeLadderStereoProcess(float *restrict in, float *restrict out, iirLatticeLadder *filt, int len, int numFilters, Bool on)
{
	float *signalInL, *signalOutL, *signalInR, *signalOutR;
	float w1L, w4L, allpassL, allpoleL, sumL, inputL;
	float w1R, w4R, allpassR, allpoleR, sumR, inputR;
	int i, j;
	iirLatticeLadder *restrict filt1 = filt;

	signalInL = in;
	signalOutL = out;
	signalInR = in + len;
	signalOutR = out + len;

	// bypass filter
	if (on == FALSE)
	{
		for (j=0; j<len; j++)
		{
			*signalOutL++ = *signalInL++;
			*signalOutR++ = *signalInR++;
		}
	}
	else
	{
#pragma MUST_ITERATE(1, , )
		for (j=0; j<len; j++)
		{
			inputL = *signalInL++;
			inputR = *signalInR++;
#pragma MUST_ITERATE(1, , )
			for (i=0; i<numFilters; i++)
			{
				w1L = filt1[i].c2 * inputL - filt1[i].k2 * filt1[i].w5;	// L
				w1R = filt1[i].c2 * inputR - filt1[i].k2 * filt1[i].x5;	// R

				allpassL = filt1[i].k2 * inputL + filt1[i].c2 * filt1[i].w5;	// L
				allpassR = filt1[i].k2 * inputR + filt1[i].c2 * filt1[i].x5;	// R

				allpoleL = filt1[i].c1 * w1L - filt1[i].k1 * filt1[i].w3;	// L
				allpoleR = filt1[i].c1 * w1R - filt1[i].k1 * filt1[i].x3;	// R

				w4L = filt1[i].k1 * w1L + filt1[i].c1 * filt1[i].w3;		// L
				w4R = filt1[i].k1 * w1R + filt1[i].c1 * filt1[i].x3;		// R

				filt1[i].w5 = w4L;
				filt1[i].w3 = allpoleL;
				sumL = allpassL*filt1[i].v2+ w4L*filt1[i].v1 + allpoleL*filt1[i].v0;
				inputL = sumL;

				filt1[i].x5 = w4R;
				filt1[i].x3 = allpoleR;
				sumR = allpassR*filt1[i].v2+ w4R*filt1[i].v1 + allpoleR*filt1[i].v0;
				inputR = sumR;

			}
			*signalOutL++ = sumL;
			*signalOutR++ = sumR;
		}
	}
}