A Diophantine equation appearing in Diophantine approximation (Q1866440)

All Diophantine equations \(ax^2+ by^2+ cz^2= 1+dxyz\) with \(a,b,c,d\in \mathbb{N}\) and \(a\mid d\), \(b\mid d\), \(c\mid d\), with solutions \((x,y,z)\in \mathbb{N}^3\) are classified. Equations of this type occur in connection with the approximation spectrum for quaternions.