Abstract In this note, we aim to review algorithms for constructing crystal representations of quiver Yangians in detail. Quiver Yangians are believed to describe an action of the BPS algebra on BPS states in systems of D-branes wrapping toric Calabi-Yau three-folds. Crystal modules of these algebras originate from molten crystal models for Donaldson-Thomas invariants of respective three-folds. Despite the fact that this subject was originally at the crossroads of algebraic geometry with effective supersymmetric field theories, equivariant toric action simplifies applied calculations drastically. So the sole pre-requisite for this algorithm’s implementation is linear algebra. It can be easily taught to a machine with the help of any symbolic calculation system. Moreover, these algorithms may be generalized to toroidal and elliptic algebras and exploited in various numerical experiments with those algebras. We illustrate the work of the algorithms in applications to simple cases of $$ \textrm{Y}\left({\mathfrak{sl}}_2\right) $$ Y sl 2 , $$ \textrm{Y}\left({\hat{\mathfrak{gl}}}_1\right) $$ Y gl ̂ 1 and $$ \textrm{Y}\left({\hat{\mathfrak{gl}}}_{\left.1\right|1}\right) $$ Y gl ̂ 1 1 .