 			}
 		}
 	}
+
+	devStressD[index] = devS3;
+	volStressD[index] = volS3;
 }
 //--------------------------------------------------------------------------------------------------------------------------------
 //without normalization

				real3 vr = R3(velMasB - velMasA);
				real3 gradW = GradW(dist3);

			}
		}
	}

		uint* cellStart,
		uint* cellEnd,
		int numBCE) {
	uint index = blockIdx.x * blockDim.x + threadIdx.x;
	if (index >= numBCE) {
		return;
	}
	//Arman take care of this
	uint BCE_Index = index + min(numObjectsD.startRigidMarkers, numObjectsD.startRigidMarkers); // updatePortionD = [start, end] index of the update portion
	uint sortedIndex = mapOriginalToSorted[BCE_Index]; //index in the sorted array

	// read particle data from sorted arrays
