In this paper, we describe an algorithm for computing algebraic modular forms
on compact inner forms of $\mathrm{GSp}_4$ over totally real number fields. By
analogues of the Jacquet-Langlands correspondence for $\mathrm{GL}_2$, this
algorithm in fact computes Hecke eigensystems of Hilbert-Siegel modular forms
of genus 2. We give some examples of such eigensystems over $\Q(\sqrt{2})$.