本文整理汇总了C++中cppad::ADFun::Memory方法的典型用法代码示例。如果您正苦于以下问题：C++ ADFun::Memory方法的具体用法？C++ ADFun::Memory怎么用？C++ ADFun::Memory使用的例子？那么恭喜您, 这里精选的方法代码示例或许可以为您提供帮助。您也可以进一步了解该方法所在类cppad::ADFun的用法示例。
在下文中一共展示了ADFun::Memory方法的1个代码示例，这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞，您的评价将有助于系统推荐出更棒的C++代码示例。
示例1: mexFunction
//.........这里部分代码省略.........
}
//Check length of Sigma
if(mxIsEmpty(prhs[2]) || !mxIsDouble(prhs[2]) || mxIsComplex(prhs[2]) || mxGetNumberOfElements(prhs[2]) != 1)
mexErrMsgTxt("Sigma must be a real, double scalar.");
//Check length of Lambda
if(!mxIsDouble(prhs[3]) || mxIsComplex(prhs[3]) || mxIsSparse(prhs[3]) || mxGetNumberOfElements(prhs[3]) != getNoCon())
mexErrMsgTxt("Lambda must be a real, double, dense vector with ncon elements.");
//Get Sigma, Lambda
sigma = *mxGetPr(prhs[2]);
lambda = mxGetPr(prhs[3]);
imode = 5;
//Can't allocate here until we know sparsity pattern
}
else if(!strcmp(mode,"hstr")) {
imode = 6;
//Can't allocate here until we know sparsity pattern
}
else
mexErrMsgTxt("Unknown mode - options are 'obj', 'grad', 'con', 'jac', 'jacstr', 'hess' or 'hstr'");
mxFree(mode);
//Ensure we did have x for normal callbacks
if(imode != 4 && imode != 6 && nrhs < 2)
mexErrMsgTxt("You must supply the callback mode and input vector.");
//Call Req Callback
switch(imode)
{
case 0: //objective
*v = objective(xvec);
break;
case 1: //gradient
//Check if we have recorded the objective yet
if(obj.Memory()==0) { //new, tape operations
vector< CppAD::AD<double> > X(getNoVar());
memcpy(&X[0],x,getNoVar()*sizeof(double));
CppAD::Independent(X);
vector< CppAD::AD<double> > Y(1);
Y[0] = objective(X);
obj = CppAD::ADFun<double>(X, Y);
//obj.optimize();
mexAtExit(mexExit); //also register memory clear function
//mexPrintf("Evaluated Tape for Gradient\n");
}
//Evaluate "Jacobian" for gradient
memcpy(v,&(obj.Jacobian(xvec)[0]),getNoVar()*sizeof(double));
break;
case 2: //constraints
//Check if we have constraint memory yet
if(cvec.empty())
cvec.resize(getNoCon()); //allocate it
//Evaluate Constraints
constraints(xvec,cvec);
//Copy Out
memcpy(v,&cvec[0],getNoCon()*sizeof(double));
break;
case 3: //jacobian
case 4: //jacobian structure
//Check if we have recorded the constraints yet
if(con.Memory()==0){ //new, tape operations
vector< CppAD::AD<double> > X(getNoVar());
memcpy(&X[0],x,getNoVar()*sizeof(double));
CppAD::Independent(X);
vector< CppAD::AD<double> > Y(getNoCon());
constraints(X,Y);
con = CppAD::ADFun<double>(X, Y);