  We have implemented a new way of computing three-point correlation functions.
It is based on a factorization of the entire correlation function into two
parts which are evaluated with open spin- (and to some extent flavor-) indices.
This allows us to estimate the two contributions simultaneously for many
different initial and final states and momenta, with little computational
overhead. We explain this factorization as well as its efficient implementation
in a new library which has been written to provide the necessary functionality
on modern parallel architectures and on CPUs, including Intel's Xeon Phi
series.
