# FlatHermitianOperator

Inheritance Diagram

Methods

 `FlatHermitianOperator.__init__`(npc_matvec, ...) Initialize this LinearOperator. Hermitian adjoint. Matrix-matrix or matrix-vector multiplication. `FlatHermitianOperator.eigenvectors`(*args, ...) Same as FlatLinearOperator(..., hermitian=True). Convert flat numpy vector of selected charge sector into npc Array. Convert flat vector of undetermined charge sectors into npc Array. `FlatHermitianOperator.from_NpcArray`(mat[, ...]) Create a FlatLinearOperator from a square `Array`. Create a FlatLinearOperator` from a matvec function acting on multiple legs. Matrix-matrix multiplication. Matrix-vector multiplication. Convert npc Array into a 1D ndarray, inverse of `flat_to_npc()`. Adjoint matrix-matrix multiplication. Adjoint matrix-vector multiplication. Transpose this linear operator.

Class Attributes and Properties

 `FlatHermitianOperator.H` Hermitian adjoint. `FlatHermitianOperator.T` Transpose this linear operator. `FlatHermitianOperator.charge_sector` Charge sector of the vector which is acted on. `FlatHermitianOperator.ndim`
class tenpy.linalg.sparse.FlatHermitianOperator(*args, **kwargs)[source]

Hermitian variant of `FlatLinearOperator`.

Note that we don’t check `matvec()` to return a hermitian result, we only define an adjoint to be self.

property H

Returns the Hermitian adjoint of self, aka the Hermitian conjugate or Hermitian transpose. For a complex matrix, the Hermitian adjoint is equal to the conjugate transpose.

Returns:

A_H – Hermitian adjoint of self.

Return type:

LinearOperator

property T

Transpose this linear operator.

Returns a LinearOperator that represents the transpose of this one. Can be abbreviated self.T instead of self.transpose().

Returns the Hermitian adjoint of self, aka the Hermitian conjugate or Hermitian transpose. For a complex matrix, the Hermitian adjoint is equal to the conjugate transpose.

Returns:

A_H – Hermitian adjoint of self.

Return type:

LinearOperator

property charge_sector

Charge sector of the vector which is acted on.

dot(x)

Matrix-matrix or matrix-vector multiplication.

Parameters:

x (array_like) – 1-d or 2-d array, representing a vector or matrix.

Returns:

Ax – 1-d or 2-d array (depending on the shape of x) that represents the result of applying this linear operator on x.

Return type:

array

flat_to_npc(vec)[source]

Convert flat numpy vector of selected charge sector into npc Array.

If `charge_sector` is not None, convert to a 1D npc vector with leg self.leg. Otherwise convert vec, which can be non-zero in all charge sectors, to a npc matrix with an additional `'charge'` leg to allow representing the full vector at once.

Parameters:

vec (1D ndarray) – Numpy vector to be converted. Should have the entries according to self.charge_sector.

Returns:

npc_vec – Same as vec, but converted into a npc array.

Return type:

`Array`

flat_to_npc_None_sector(vec, cutoff=1e-10)[source]

Convert flat vector of undetermined charge sectors into npc Array.

The charge sector to be used is chosen as the block with the maximal norm, not by self.charge_sector (which might be None).

Parameters:

vec (1D ndarray) – Numpy vector to be converted.

Returns:

npc_vec – Same as vec, but converted into a npc array.

Return type:

`Array`

classmethod from_NpcArray(mat, charge_sector=0, compact_flat=None)[source]

Create a FlatLinearOperator from a square `Array`.

Parameters:
• mat (`Array`) – A square matrix, with contractable legs.

• charge_sector (None | charges | `0`) – Selects the charge sector of the vector onto which the Linear operator acts. `None` stands for all sectors, `0` stands for the zero-charge sector. Defaults to `0`, i.e., assumes the dominant eigenvector is in charge sector 0.

• compact_flat (bool | None) – If True, restrict the flat array to the (only) non-zero block of given charge_sector. If False, the flat array is directly what’s represented by the npc Array’s `to_ndarray()`. Works only for fixed charge sector and if the leg of mat is blocked; None defaults to `leg.is_blocked()`.

classmethod from_guess_with_pipe(npc_matvec, v0_guess, labels_split=None, dtype=None, compact_flat=True)[source]

Create a FlatLinearOperator` from a matvec function acting on multiple legs.

This function creates a wrapper matvec function to allow acting on a “vector” with multiple legs. The wrapper combines the legs into a `LegPipe` before calling the actual matvec function, and splits them again in the end.

Parameters:
• npc_matvec (function) – Function to calculate the action of the linear operator on an npc vector with the given split labels labels_split. Has to return an npc vector with the same legs.

• v0_guess (`Array`) – Initial guess/starting vector which can be applied to npc_matvec.

• labels_split (None | list of str) – Labels of v0_guess in the order in which they are to be combined into a `LegPipe`. `None` defaults to `v0_guess.get_leg_labels()`.

• dtype (np.dtype | None) – The data type of the arrays. `None` defaults to dtype of v0_guess (!).

• compact_flat (bool) – If True, restrict the flat array to the non-zero parts. If False, the flat array is directly what’s represented by the npc Array’s `to_ndarray()`.

Returns:

• lin_op (cls) – Instance of the class to be used as linear operator

• guess_flat (np.ndarray) – Numpy vector representing the guess v0_guess.

matmat(X)

Matrix-matrix multiplication.

Performs the operation y=A*X where A is an MxN linear operator and X dense N*K matrix or ndarray.

Parameters:

X ({matrix, ndarray}) – An array with shape (N,K).

Returns:

Y – A matrix or ndarray with shape (M,K) depending on the type of the X argument.

Return type:

{matrix, ndarray}

Notes

This matmat wraps any user-specified matmat routine or overridden _matmat method to ensure that y has the correct type.

matvec(x)

Matrix-vector multiplication.

Performs the operation y=A*x where A is an MxN linear operator and x is a column vector or 1-d array.

Parameters:

x ({matrix, ndarray}) – An array with shape (N,) or (N,1).

Returns:

y – A matrix or ndarray with shape (M,) or (M,1) depending on the type and shape of the x argument.

Return type:

{matrix, ndarray}

Notes

This matvec wraps the user-specified matvec routine or overridden _matvec method to ensure that y has the correct shape and type.

npc_to_flat(npc_vec)[source]

Convert npc Array into a 1D ndarray, inverse of `flat_to_npc()`.

Parameters:

npc_vec (`Array`) – Npc Array to be converted. If self.charge_sector is not None, this should be a 1D array with that qtotal. If self.charge_sector is not None, it should have an additional `"charge"` leg, (as returned by `flat_to_npc()`).

Returns:

vec – Same entries as npc_vec, but converted into a flat Numpy array.

Return type:

1D ndarray

rmatmat(X)

Performs the operation y = A^H * x where A is an MxN linear operator and x is a column vector or 1-d array, or 2-d array. The default implementation defers to the adjoint.

Parameters:

X ({matrix, ndarray}) – A matrix or 2D array.

Returns:

Y – A matrix or 2D array depending on the type of the input.

Return type:

{matrix, ndarray}

Notes

This rmatmat wraps the user-specified rmatmat routine.

rmatvec(x)

Performs the operation y = A^H * x where A is an MxN linear operator and x is a column vector or 1-d array.

Parameters:

x ({matrix, ndarray}) – An array with shape (M,) or (M,1).

Returns:

y – A matrix or ndarray with shape (N,) or (N,1) depending on the type and shape of the x argument.

Return type:

{matrix, ndarray}

Notes

This rmatvec wraps the user-specified rmatvec routine or overridden _rmatvec method to ensure that y has the correct shape and type.

transpose()

Transpose this linear operator.

Returns a LinearOperator that represents the transpose of this one. Can be abbreviated self.T instead of self.transpose().

eigenvectors(*args, **kwargs)[source]

Same as FlatLinearOperator(…, hermitian=True).