I have an $n \times m \times p$ array, let's call it $r$.

I want to obtain an $n \times m \times p \times 5$ array (let's call it $q$) where $q$ is the same as $r$ except that every non-zero component now has a 3 in the new, fourth entry.

That is,


should return an array identical to r, and

q[[;;,;;,;;,1]],  q[[;;,;;,;;,2]], q[[;;,;;,;;,4]], q[[;;,;;,;;,5]]

should all return arrays with zeroes in all components.

This is obviously quite easy to achieve by acting directly on the array rules, but for larger arrays this becomes quite inefficient and it would be good to implement in a way which acts directly on the sparse arrays.


2 Answers 2


Another possibility is:

q = TensorProduct[r, {0, 0, 1, 0, 0}]

Is this what you're after?

  ArrayReshape[r, Append[Dimensions[r], 1]],
  Append[ConstantArray[{0, 0}, ArrayDepth[r]], {2, 2}]

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.