@@ -76,7 +76,7 @@ def lpu(
7676 :param diag: To return only the diagonal elements of :math:`U(Y)`.
7777 :returns: :math:`U(Y) \in \mathbb{R}^{M \times n \times n}`.
7878 """
79- return lpu_x (x .ndim - 1 , self .jac (x ), u )
79+ return lpu_x (x .ndim - 1 , self .jac (x ), u , diag )
8080
8181 @property
8282 @abstractmethod
@@ -159,7 +159,7 @@ def lpu_p(
159159 :param diag: To return only the diagonal elements of :math:`U(Y)`.
160160 :returns: :math:`U(Y) \in \mathbb{R}^{M \times n \times n}`.
161161 """
162- return lpu_p (p .ndim , self .jac_p (p , x ), u )
162+ return lpu_p (p .ndim , self .jac_p (p , x ), u , diag )
163163
164164 def lpu_x (
165165 self , p : np .ndarray , x : np .ndarray , u : np .ndarray , diag : bool = False
@@ -190,7 +190,7 @@ def lpu_x(
190190 :param diag: To return only the diagonal elements of :math:`U(Y)`.
191191 :returns: :math:`U(Y) \in \mathbb{R}^{M \times n \times n}`.
192192 """
193- return lpu_x (x .ndim - 1 , self .jac_x (p , x ), u )
193+ return lpu_x (x .ndim - 1 , self .jac_x (p , x ), u , diag )
194194
195195 @abstractmethod
196196 def estimate (
0 commit comments