0a7bbb0f733a33a072aa1327dc4699d317d4ce56,trunk/SUAVE/Optimization/Nexus.py,Nexus,finite_difference,#Nexus#,156

Before Change

        grad_obj = (grad_obj - obj)/(1e-8)
        
        for ii in xrange(0,conlen):
            jac_con[ii,:]  = (jac_con[ii,:]  - con[ii])/(1e-6)
            
        
        return grad_obj.astype(float), jac_con.astype(float)
    

After Change

        grad_obj = np.zeros(inplen)
        jac_con  = np.zeros((inplen,conlen))
        
        con2 = (con*np.ones_like(jac_con))
        
        for ii in xrange(0,inplen):
            newx     = x*1.0
            newx[ii] = newx[ii]+ 1e-8
            
            grad_obj[ii]  = self.objective(newx)
            jac_con[ii,:] = self.all_constraints(newx)
        
        grad_obj = (grad_obj - obj)/(1e-8)
        
        jac_con = (jac_con - con2).T/(1e-8)
        
        grad_obj = grad_obj.astype(float)
        jac_con  = jac_con.astype(float)
        

Instances