code: def lha2matrix(values, shape):
      
      M = np.zeros(shape)
      for v in values:
          M[tuple([int(i-1) for i in v[:-1]])] = v[-1]
      return M
docs: Return a matrix given a list of values of the form
      [[1, 1, float], [1, 2, float], ...]
      referring to the (1,1)-element etc.
      `shape` is the shape of the final matrix. All elements not provided
      will be assumed to be zero. Also works for higher-rank tensors.