It has been recently discovered by Bell, Heinle and Levandovskyy that a large class of algebras, including the ubiquitous G-algebras, are finite factorization domains (FFD for short).Utilizing this result, we contribute an algorithm to find all distinct factorizations of a given element f ∈ G, where G is any G-algebra, with minor assumptions on the underlying field.Moreover, the property of being an FFD, in combination with the factorization algorithm, enables us to propose an analogous description of the factorized Gröbner basis algorithm for G-algebras. This algorithm is useful for various applications, e.g. in analysis of solution spaces of systems of linear partial functional equations with polynomial coefficients, coming from G. Additionally, it is possible to include inequality constraints for ideals in the input.