We know that
Let $d$ be a positive square free integer and $r$ an integer satify $r^2+|r|\le d$. Suppose $x$ and $y$ are positive integers that satify $x^2-dy^2=r$. Then $\frac xy$ is a convergent to the continued fraction of $\sqrt d$.
What are the positive integer solutions of $x^2-dy^2=r$ when $r^2+|r|>d$?
All solutions to $ w^2 - 13 v^2 = 15249$ can be constructed by applying the mapping $$ (w,v) \mapsto ( \; 649 w + 2340 v \; , \; \; 180 w + 649 v \;) \; \; $$ to the first eight "SEED" solutions in the output below. For example, $$ (143,20) \mapsto ( \; 649 \cdot 143 + 2340 \cdot 20 \; , \; \; 180 \cdot 143 + 649 \cdot 20 \;) = (139607, 38720) \; \; $$
The exact same thing works for the output below that, all solutions to $ w^2 - 13 v^2 = -15249$ from those eight seed solutions. For example, $$ (26,35) \mapsto ( \; 649 \cdot 26 + 2340 \cdot 35 \; , \; \; 180 \cdot 26 + 649 \cdot 35 \;) = (98774, 27395) \; \; $$