FlatHermitianOperator¶
full name: tenpy.linalg.sparse.FlatHermitianOperator
parent module:
tenpy.linalg.sparse
type: class
Inheritance Diagram
Methods
|
Initialize this LinearOperator. |
Hermitian adjoint. |
|
Matrix-matrix or matrix-vector multiplication. |
|
|
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. |
|
Convert flat vector of all charge sectors into npc Array with extra “charge” leg. |
|
|
Create a FlatLinearOperator from a square |
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 |
Convert npc Array with qtotal = self.charge_sector into ndarray. |
|
Adjoint matrix-matrix multiplication. |
|
Adjoint matrix-vector multiplication. |
|
Transpose this linear operator. |
Class Attributes and Properties
Hermitian adjoint. |
|
Transpose this linear operator. |
|
Charge sector of the vector which is acted on. |
|
|
-
class
tenpy.linalg.sparse.
FlatHermitianOperator
(*args, **kwargs)[source]¶ Bases:
tenpy.linalg.sparse.FlatLinearOperator
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
¶ Hermitian adjoint.
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.
Can be abbreviated self.H instead of self.adjoint().
- 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().
-
adjoint
()¶ Hermitian adjoint.
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.
Can be abbreviated self.H instead of self.adjoint().
- 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
-
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
-
flat_to_npc_all_sectors
(vec)[source]¶ Convert flat vector of all charge sectors into npc Array with extra “charge” leg.
Deprecated since version 0.7.3: This is merged into
flat_to_npc()
withself.charge_sector = None
.- Parameters
vec (1D ndarray) – Numpy vector to be converted.
- Returns
npc_vec – Same as vec, but converted into a npc array.
- Return type
-
classmethod
from_NpcArray
(mat, charge_sector=0)[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 to0
, i.e., assumes the dominant eigenvector is in charge sector 0.
-
classmethod
from_guess_with_pipe
(npc_matvec, v0_guess, labels_split=None, dtype=None)[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 tov0_guess.get_leg_labels()
.dtype (np.dtype | None) – The data type of the arrays.
None
defaults to dtype of v0_guess (!).
- 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 byflat_to_npc()
).- Returns
vec – Same entries as npc_vec, but converted into a flat Numpy array.
- Return type
1D ndarray
-
npc_to_flat_all_sectors
(npc_vec)[source]¶ Convert npc Array with qtotal = self.charge_sector into ndarray.
Deprecated since version 0.7.3: This is merged into
npc_to_flat()
withself.charge_sector = None
.- Parameters
npc_vec (
Array
) – Npc Array to be converted. Should only have entries in self.charge_sector.- Returns
vec – Same as npc_vec, but converted into a flat Numpy array.
- Return type
1D ndarray
-
rmatmat
(X)¶ Adjoint matrix-matrix multiplication.
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)¶ Adjoint matrix-vector multiplication.
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().
-
property