Compute Deviatoric Stress at each lattice node in STL Code


In the sparse domain creation from STL file in Palabos, I want to compute the deviatoric stress at each respective lattice node. I am not yet able to figure out a way for that. The deviatoric stress has 6 components and I tried to store these components in an array using the code snippet below but it causes a deadlock while executing the code.

Array<T,3> stress_arr;

I will really appreciate if someone can help me with this.