Pinv numpy backend

tensornetwork/backends/numpy/numpy_backend.py

@@ -795,3 +795,21 @@ def eps(self, dtype: Type[np.number]) -> float:
       float: Machine epsilon.
     """
     return np.finfo(dtype).eps
+
+  def pinv(self, tensor: Tensor, rcond: float = 1E-15,
+           hermitian: bool = False) -> Tensor:
+    """
+    Compute the (Moore-Penrose) pseudo-inverse of a tensor.
+    Returns the pseudo-inverse of tensor.
+
+    Args:
+     tensor: A tensor.
+     rcond: Cutoff for small singular values.
+     hermitian(optional): If True, matrix provided is 
+                          assumed to be Hermitian (symmetric if real-valued). 
+                          Defaults to False.
+
+    Returns:
+     tensor: The pseudo inverse of tensor.
+         """
+    return np.linalg.pinv(tensor)

tensornetwork/backends/tensorflow/tensorflow_backend.py

@@ -430,3 +430,18 @@ def eps(self, dtype: Type[np.number]) -> float:
       float: Machine epsilon.
     """
     return tf.experimental.numpy.finfo(dtype).eps
+
+  def pinv(self, tensor: Tensor, rcond: float = 1E-15, hermitian: bool = False) -> Tensor:
+    """
+    Compute the (Moore-Penrose) pseudo-inverse of a tensor.
+    Returns the pseudo-inverse of tensor.
+
+    Args:
+     tensor: A tensor.
+     rcond: Cutoff for small singular values.
+     hermitian(optional): If True, matrix provided is assumed to be Hermitian (symmetric if real-valued). Defaults to False.
+
+    Returns:
+     tensor: The pseudo inverse of tensor.
+         """
+    return np.linalg.pinv(tensor)