@Olegg: I may have a similar problem here where I want to define constraints in JuMP involving matrix invs. Would you mind sharing how you registered a multivariate function / derivative? Or did you find an alternative approach?
inv