I’ll edit the document to have the appropriate signatures. And yes, theta and delta need only be arrays, though we may want to have them as vectors. I’ll tinker on that. I was thinking of calling the function
quadratic_optimizer – it’s a slightly different problem than solving equations.
The gitHub issue is here https://github.com/stan-dev/math/issues/888.