function [x1x2] = xprod( x1, x2, usr_par)

x1x2 = x1'*x2;