COMPDIR computes a search direction in a subspace defined by Z. % [SD,dirType] = compdir(Z,H,gf,nvars,f) returns a search direction for the % subproblem 0.5*Z'*H*Z + Z'*gf. Helper function for NLCONST. SD is Newton % direction if possible. SD is