function fg=evaluateCenterDiff(cax,centers,w,gParConst,x,pars)
    globalPar=[gParConst(3) x(1:2) 0 gParConst(1:2)];
    fv=evaluateCenters(globalPar,pars,cax);
    fg=reshape((fv-centers).*w,[],1);
end