code: def matrix_multiply(m0, m1):
      new_matrix = zero(len(m0), len(m1[0]))
      for i in range(len(m0)):
          for j in range(len(m1[0])):
              r = 0
              "omp parallel for reduction(+:r)"
              for k in range(len(m1)):
                  r += m0[i][k] * m1[k][j]
              new_matrix[i][j] = r
      return new_matrix
  
label_desc: Use image node_id 2 for calling a global function with example usage: matrix_multiply(m0, m1) and returns: new_matrix