User guide

16-31
c0.bounds[lind,2] = 10*ones(rows(lind),1);
c0.bounds[tind,1] = .001*ones(rows(tind),1);
c0.bounds[tind,2] = 100*ones(rows(tind),1);
c0.output = 1;
c0.printIters = 1;
c0.trustRadius = 1;
c0.ineqProc = &ineq;
c0.covType = 1;
struct DS d0;
d0 = dsCreate;
d0.dataMatrix = loadd("maxfact");
output file = sqpfact.out reset;
struct SQPsolveMTOut out0;
out0 = SQPsolveMT(&lpr,par0,d0,c0);
lambdahat = pvUnpack(out0.par, "lambda");
phihat = pvUnpack(out0.par, "phi");
psihat = pvUnpack(out0.par, "psi");
print "estimates";
print;
print "lambda" lambdahat;
print;
print "phi" phihat;
print;
print "psi" psihat;
struct PV stderr;
stderr = out0.par;
if not scalmiss(out0.moment);
stderr = pvPutParVector(stderr,sqrt(diag(out0.moment)));
lambdase = pvUnpack(stderr, "lambda");
phise = pvUnpack(stderr, "phi");
Structures