Jesse Grabowski

Optimize non-structured unary elementwise `Ops` on sparse matrices

3

### Description The title sucks, but the case I'm interested in is a graph like this: ``` import pytensor.sparse as pts import pytensor.tensor as pt x = pts.csr_dmatrix() out =...

`local_det_chol` doesn't work

1

### Description A simple example: ```py import pytensor.tensor as pt import pytensor x = pt.dmatrix('x') out = pt.linalg.det(pt.linalg.cholesky(x)) fn = pytensor.function([x], out) fn.dprint() Det [id A] 1 └─ Cholesky{lower=True, check_finite=False,...

### Description This came up in the context of #1674. We should consider reverting https://github.com/aesara-devs/aesara/issues/142. Sparse matrices and dense tensor are conceptually different enough that it doesn't make sense for...

Add rewrite to fuse nested `BlockDiag` Ops

2

### Description This comes up in statespace models. `BlockDiag` is variadic, so graphs like this can be canonicalized: ```py import pytensor.tensor as pt import pytensor a, b, c = pt.matrices('a',...

BLAS warnings in new github windows CI

1

### Description We're seeing BLAS warnings in the windows CI for `pymc-extras`. There are two. [One is the usual linker failure](https://github.com/pymc-devs/pymc-extras/actions/runs/18078168486/job/51437649154?pr=584): CI log snippet ``` 2025-09-28T18:44:42.2527549Z ERROR  pytensor.graph.rewriting.basic:basic.py:1757 Rewrite...

Implement pack/unpack helpers

4

## Description Adds pt.pack and pt.unpack helpers, roughly conforming to the `einops` functions of the same name. These helps are for situations where we have a ragged list of inputs...

Directly use scipy lapack functions instead of scipy wrappers

6

### Description We take small performance hits by using `scipy.linalg` functions rather than directly calling lapack functions in our `perform` methods. These wrappers include a lot of data validation checks...

### Description If aixs = None numba will try to iterate over None and error (TODO: Add MWE)

### Description Define $\text{vec}(X)$ as the ravel operator that works column-wise (that's the math convention, so in code it's `x.reshape(-1, order="F")`. Given $(A \otimes B)x = y$. If $y =...

Implement pack/unpack Ops

3

### Description This came up in the context of #1550 . There are many cases where we want to be able to ravel a sequence of inputs into a single...