biNrdZddlmZmZddlZddlmZddlm Z ddl m Z ddlmZddlmZGddeZy) a0 Copyright 2022, the CVXPY authors Licensed under the Apache License, Version 2.0 (the "License"); you may not use this file except in compliance with the License. You may obtain a copy of the License at http://www.apache.org/licenses/LICENSE-2.0 Unless required by applicable law or agreed to in writing, software distributed under the License is distributed on an "AS IS" BASIS, WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. See the License for the specific language governing permissions and limitations under the License. )ListTupleN)linalg)Atom) ConstraintceZdZdZdfd ZdZddZdeedffdZ dee e ffdZ de fd Z de fd Z de fd Zde fd Zd ZdeefdZedefdZxZS)tr_invz\ :math:`\mathrm{tr}\left(X^{-1} \right),` where :math:`X` is positive definite. returnc,tt| |y)N)superr __init__)selfX __class__s M/home/cdr/jupyterlab/.venv/lib/python3.12/site-packages/cvxpy/atoms/tr_inv.pyr ztr_inv.__init__"s fd$Q'cZtj|d|djjz dk\rtj S|d|djzdz }tj |}t|dkrtj St j|dzS)zReturns the trinv of positive definite matrix X. For positive definite matrix X, this is the trace of inverse of X. rg:0yE>) LAnormTconjnpinfeigvalshminsum)rvaluessymmeigVals rnumericztr_inv.numeric%s GGF1Iq 0 0 22 3t ;66Mq F1IKK'*T" v;! 66Mvvfbj!!rc|jdjdk(r9|jdjd|jdjdk7r*td|jdj dy)z,Verify that the argument is a square matrix.rrz The argument z& to tr_inv must be a 2-d square array.N)argsndimshape ValueErrornamers rvalidate_argumentsztr_inv.validate_arguments4sxyy|  A%1););A)>$))A,BTBTUVBW)W ! 1 1 344Z[ *Xr.ctS)z8Returns the (row, col) shape of the expression. )tupler*s rshape_from_argsztr_inv.shape_from_args;s wrcy)zCReturns sign (is positive, is negative) of the expression. )TFr*s rsign_from_argsztr_inv.sign_from_args@srcy)zIs the atom convex? Tr0r*s ris_atom_convexztr_inv.is_atom_convexEsrcy)zIs the atom concave? Fr0r*s ris_atom_concaveztr_inv.is_atom_concaveJrcy)z;Is the composition non-decreasing in argument idx? Fr0ridxs ris_incrztr_inv.is_incrOr6rcy)z;Is the composition non-increasing in argument idx? Fr0r8s ris_decrztr_inv.is_decrTr6rc*|d}tj|}tj|dkDr_tjj |j }| |z}tj|jdj gSdgS)a+Gives the (sub/super)gradient of the atom w.r.t. each argument. Matrix expressions are vectorized, so the gradient is a matrix. Args: values: A list of numeric values for the arguments. Returns: A list of SciPy CSC sparse matrices or None. rF)orderN) rrrrrinvrsp csc_arrayravel)rrr eigen_valDs r_gradz tr_inv._gradYsz 1IKKN 66) q  a ""AaALLs!34667 76Mrc(|jddz gS)z?Returns constraints describing the domain of the node. r)r%r*s r_domainztr_inv._domainos ! !""rctj|jdj|jdjjj t jt js tdt|jDcgc]}|jduc}ry|jScc}w)Nr)rtolatolz2Input array to tr_inv was not Hermitian/symmetric.) rallcloser%valuerrs ATOM_EVAL_TOLr(any parameters _value_impl)rps rrMz tr_inv.valuets{{499Q<--99Q<--//446 ! !1QR R ):;A4; <!!res/ !3a"Ta"r