/usr/share/doc/root/test/QpRandomDriver.cxx is in root-system-doc 5.34.30-0ubuntu8.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 | #include <stdlib.h>
#include "Riostream.h"
#include "TQpDataDens.h"
#include "TQpVar.h"
#include "TQpProbDens.h"
#include "TGondzioSolver.h"
Bool_t SolutionMatches(TQpVar *vars,TQpVar *soln,TQpVar *temp,Double_t tol);
int main(int argc,char *argv[])
{
Int_t n1 = 5, m1 = 2, m2 = 2;
if (argc >= 4) {
n1 = atoi(argv[1]);
m1 = atoi(argv[2]);
m2 = atoi(argv[3]);
} else {
cout << endl
<< " Usage: QpRandomDriver n my mz " << endl
<< " where n = # primal variables, " << endl
<< " my = # equality constraints, " << endl
<< " mz = # inequality constraints " << endl << endl;
return 1;
}
Int_t nnzQ = (Int_t) .20*(n1*n1);
Int_t nnzA = (Int_t) .15*(m1*n1);
Int_t nnzC = (Int_t) .10*(m2*n1);
if (nnzQ < 3*n1) nnzQ = 3*n1;
if (nnzA < 3*m1) nnzA = 3*m1;
if (nnzC < 3*m2) nnzC = 2*m2;
TQpProbDens *qp = new TQpProbDens(n1,m1,m2);
TQpDataDens *prob;
TQpVar *soln;
qp->MakeRandomData(prob,soln,0,0,0);
// qp->MakeRandomData(prob,soln,nnzQ,nnzA,nnzC);
TQpVar *vars = qp->MakeVariables(prob);
TQpResidual *resid = qp->MakeResiduals(prob);
TGondzioSolver *s = new TGondzioSolver(qp,prob);
const Int_t status = s->Solve(prob,vars,resid);
delete s;
if (status == 0) {
cout.precision(4);
cout << endl << "Computed solution:\n\n" <<endl <<endl;
vars->fX.Print();
TQpVar *temp = qp->MakeVariables(prob);
cout << endl << "Checking the solution...";
if (SolutionMatches(vars,soln,temp,1e-4)) {
cout << "The solution appears to be correct." <<endl;
} else {
cout << endl << "The solution may be wrong "
"(or the generated problem may be ill conditioned.)" <<endl;
}
delete temp;
} else {
cout << "Could not solve this problem." <<endl;
}
delete vars;
delete soln;
delete prob;
delete qp;
return status;
}
Bool_t SolutionMatches(TQpVar *vars,TQpVar *soln,TQpVar *temp,double tol)
{
temp = vars;
// Only x is significant
temp->fX -= soln->fX;
if ((temp->fX).NormInf()/(1+(soln->fX).NormInf()) < tol)
return kTRUE;
else
return kFALSE;
}
|