//
//  pendulum_adolc.cpp
//
//  ADOLC C++ file for the vector field named: pendulum
//
//  This file was generated by the program VFGEN (Version:2.4.0)
//  Generated on 10-Jul-2008 at 11:27
//

#include <math.h>
#include "adolc/adouble.h"
#include "adolc/drivers/odedrivers.h"
#include "adolc/adalloc.h"

//
//  The vector field.
//

void pendulum_vf(short int tag, double *y_, double *f_, double *params_)
    {
    adoublev ay_(2);
    adoublev af_(2);

    const adouble Pi = M_PI;
    adouble theta, v;
    adouble g, b, L, m;

    trace_on(tag);
    ay_ <<= y_;
    theta      = ay_[0];
    v          = ay_[1];

    g          = params_[0];
    b          = params_[1];
    L          = params_[2];
    m          = params_[3];

    af_[0] = v;
    af_[1] = -g/L*sin(theta)-1.0/(L*L)*b/m*v;
    af_ >>= f_;
    trace_off(tag);
    }