code: def test_calculation():
      a = MutableDenseNDimArray([1] * 9, (3, 3))
      b = MutableDenseNDimArray([9] * 9, (3, 3))
  
      c = a + b
      for i in c:
          assert i == MutableDenseNDimArray([10, 10, 10])
  
      assert c == MutableDenseNDimArray([10] * 9, (3, 3))
      assert c == MutableSparseNDimArray([10] * 9, (3, 3))
  
      c = b - a
      for i in c:
          assert i == MutableSparseNDimArray([8, 8, 8])
  
      assert c == MutableDenseNDimArray([8] * 9, (3, 3))
      assert c == MutableSparseNDimArray([8] * 9, (3, 3))
  
imports: from copy import copy
  from sympy.tensor.array.dense_ndim_array import MutableDenseNDimArray
  from sympy.core.function import diff
  from sympy.core.numbers import Rational
  from sympy.core.singleton import S
  from sympy.core.symbol import Symbol
  from sympy.core.sympify import sympify
  from sympy.matrices import SparseMatrix
  from sympy.matrices import Matrix
  from sympy.tensor.array.sparse_ndim_array import MutableSparseNDimArray
  from sympy.testing.pytest import raises
label_desc: Use image node_id 7 for calling a global function with example usage: test_calculation() without return types