function [f,df] = gmmobjg(theta2)
% This function computes the GMM objective function

% Written by Aviv Nevo, May 1998.
% Modified by Bronwyn Hall, April 2005, to add gradient.

global invA theta1 theti thetj x1 IV

delta = meanval(theta2);

% the following deals with cases were the min algorithm drifts into region where the objective is not defined
if max(isnan(delta)) == 1
	f = 1e+10	   
else
    temp1 = x1'*IV;
    temp2 = delta'*IV;
    theta1 = inv(temp1*invA*temp1')*temp1*invA*temp2';
    clear temp1 temp2 
    gmmresid = delta - x1*theta1;
	temp1 = gmmresid'*IV;
	f1 = temp1*invA*temp1';
    f = f1;
    clear temp1
    if nargout>1
        load mvalold
        temp = jacob(mvalold,theta2)';
        df = 2*temp*IV*invA*IV'*gmmresid
    end
end

disp(['GMM objective:  ' num2str(f1)])