o ְi@srddlZddlZddlmZdddZdddZeddd d Zddefd d Zdd dZ dddZ dddZ dS)N)_corecCst|tjr |j}nt|}|dkrt|}|dkr!t|}|dkrA|jt|}tj||f|j d}|| |dd<|S|dkrMt| |St d)aReturns a diagonal or a diagonal array. Args: v (array-like): Array or array-like object. k (int): Index of diagonals. Zero indicates the main diagonal, a positive value an upper diagonal, and a negative value a lower diagonal. Returns: cupy.ndarray: If ``v`` indicates a 1-D array, then it returns a 2-D array with the specified diagonal filled by ``v``. If ``v`` indicates a 2-D array, then it returns the specified diagonal of ``v``. In latter case, if ``v`` is a :class:`cupy.ndarray` object, then its view is returned. .. seealso:: :func:`numpy.diag` dtypeNzInput must be 1- or 2-d.) isinstancecupyndarrayndimnumpyarraysizeabszerosrdiagonal ValueError)vkr r retrI/home/ubuntu/.local/lib/python3.10/site-packages/cupy/_creation/matrix.pydiags    rcCs$t|r t|}t||S)aECreates a diagonal array from the flattened input. Args: v (array-like): Array or array-like object. k (int): Index of diagonals. See :func:`cupy.diag` for detail. Returns: cupy.ndarray: A 2-D diagonal array with the diagonal copied from ``v``. .. seealso:: :func:`numpy.diagflat` )r isscalarasarrayrrravel)rrrrrdiagflat/s rzint32 m, int32 kzT outzK int row = i / m; int col = i % m; out = (col <= row + k); cupy_tricCs*|dur|}tj||f|d}t|||S)aCreates an array with ones at and below the given diagonal. Args: N (int): Number of rows. M (int): Number of columns. ``M == N`` by default. k (int): The sub-diagonal at and below which the array is filled. Zero is the main diagonal, a positive value is above it, and a negative value is below. dtype: Data type specifier. Returns: cupy.ndarray: An array with ones at and below the given diagonal. .. seealso:: :func:`numpy.tri` Nr)rempty _tri_kernel)NMrroutrrrtriNs r"cCs:t|}t|jdd|td}t|||jdS)a}Returns a lower triangle of an array. Args: m (array-like): Array or array-like object. k (int): The diagonal above which to zero elements. Zero is the main diagonal, a positive value is above it, and a negative value is below. Returns: cupy.ndarray: A lower triangle of an array. .. seealso:: :func:`numpy.tril` Nrrrrrr"shapeboolwherertypemrmaskrrrtrilfs r-cCs>t|}t|jdd|dtd}t||jd|S)aReturns an upper triangle of an array. Args: m (array-like): Array or array-like object. k (int): The diagonal below which to zero elements. Zero is the main diagonal, a positive value is above it, and a negative value is below. Returns: cupy.ndarray: An upper triangle of an array. .. seealso:: :func:`numpy.triu` r#Nrr$rr%r*rrrtriu{s r.FcCst|}|jdkrtd|durt|}tjt||ft|jt d}|s3|dddddfn|}tj | ddt ||d|S)aReturns a Vandermonde matrix. Args: x (array-like): 1-D array or array-like object. N (int, optional): Number of columns in the output. ``N = len(x)`` by default. increasing (bool, optional): Order of the powers of the columns. If True, the powers increase from right to left, if False (the default) they are reversed. Returns: cupy.ndarray: A Vandermonde matrix. .. seealso:: :func:`numpy.vander` rz.x must be a one-dimensional array or sequence.Nr)r!) rrr rlenrr promote_typesrintpowerreshapearange)xr increasingrtmprrrvanders   r9)r)NF) r rrrrElementwiseKernelrfloatr"r-r.r9rrrrs   (