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