/**
* Copyright (c) 2008, 2009, 2010 Regents of the University of California
*
* ADModelbuilder and associated libraries and documentations are
* provided under the general terms of the "BSD" license.
*
* License:
*
* Redistribution and use in source and binary forms, with or without
* modification, are permitted provided that the following conditions are
* met:
*
* 1. Redistributions of source code must retain the above copyright
* notice, this list of conditions and the following disclaimer.
*
* 2. Redistributions in binary form must reproduce the above copyright
* notice, this list of conditions and the following disclaimer in the
* documentation and/or other materials provided with the distribution.
*
* 3. Neither the name of the University of California, Otter Research,
* nor the ADMB Foundation nor the names of its contributors may be used
* to endorse or promote products derived from this software without
* specific prior written permission.
*
* THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
* "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
* LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR
* A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT
* OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL,
* SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT
* LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE,
* DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY
* THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
* (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE
* OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
*
*/
#include
#include
double fcomp(int, dvar_vector); // Function prototype declaration
#ifdef __BCPLUSPLUS__
extern unsigned _stklen = 20000;
#endif
#ifdef __ZTC__
long _stack = 20000;
#endif
int main()
{
double f;
int i;
int nvar=20;
independent_variables x(1,nvar); // Identify the independent variables
dvector g(1,nvar); // Holds the vector of partial derivatives (the gradient)
gradient_structure gs; // must declare this structure to manage derivative
// calculations
f=fcomp(nvar,x);
gradcalc(nvar,g); // The derivatives are calculated
cout <<" The gradient vector is\n"<